- Open Access
Mapping hotspots of malaria transmission from pre-existing hydrology, geology and geomorphology data in the pre-elimination context of Zanzibar, United Republic of Tanzania
Parasites & Vectorsvolume 8, Article number: 41 (2015)
Larval source management strategies can play an important role in malaria elimination programmes, especially for tackling outdoor biting species and for eliminating parasite and vector populations when they are most vulnerable during the dry season. Effective larval source management requires tools for identifying geographic foci of vector proliferation and malaria transmission where these efforts may be concentrated. Previous studies have relied on surface topographic wetness to indicate hydrological potential for vector breeding sites, but this is unsuitable for karst (limestone) landscapes such as Zanzibar where water flow, especially in the dry season, is subterranean and not controlled by surface topography.
We examine the relationship between dry and wet season spatial patterns of diagnostic positivity rates of malaria infection amongst patients reporting to health facilities on Unguja, Zanzibar, with the physical geography of the island, including land cover, elevation, slope angle, hydrology, geology and geomorphology in order to identify transmission hot spots using Boosted Regression Trees (BRT) analysis.
The distribution of both wet and dry season malaria infection rates can be predicted using freely available static data, such as elevation and geology. Specifically, high infection rates in the central and southeast regions of the island coincide with outcrops of hard dense limestone which cause locally elevated water tables and the location of dolines (shallow depressions plugged with fine-grained material promoting the persistence of shallow water bodies).
This analysis provides a tractable tool for the identification of malaria hotspots which incorporates subterranean hydrology, which can be used to target larval source management strategies.
It is becoming more apparent that if we are to achieve the ambitious goal of eliminating malaria, it will be necessary to complement current priority interventions for controlling adult mosquitoes that feed and rest indoors, such as long-lasting insecticide treated nets (LLINs) and indoor residual house spraying (IRS) by targeting vector mosquitoes at aquatic habitats where they are most vulnerable. This can be achieved through larval source management techniques such as environmental management or larviciding [1-8].
Zanzibar has one of the most intensive and successful programs for malaria control and elimination in Africa with LLINs and IRS applied to 90% of all dwellings [9,10]. While dramatically reducing the numbers of anthropophagic vectors, specifically Anopheles gambiae s.s. and An. funestus, species that are not suited to control by indoor interventions, such as An. arabiensis, persist . Larval source management interventions that tackle such behaviourally evasive vectors at their aquatic habitats would be most effective if the habitats to be targeted could be identified beforehand, with simple replicable tools .
The distribution and characteristics of aquatic habitats used by malarial mosquitoes for oviposition and larval development are controlled by hydrological and geomorphological processes [5,12-20]. The explicit control of hydrology and geomorphology on the formation of malaria vector habitats means that their distribution will be heterogeneous, enabling the identification of hotspots. This has been exploited in studies that associate malaria transmission rates to proximity to marshland or swamps which tend to be found in areas of topographic convergence, such as valley bottoms, for example in the Western Kenyan Highlands and the Usambara Mountains, Tanzania [21-28].
The main assumption behind the hydrological processes that govern this relationship is that the local water table lies parallel to the ground surface  and that subsurface flow pathways follow surface flow directions. However, malaria is also prevalent in limestone dominated (karst) landscapes, including southeast Asia, the south Pacific and coastal east Africa [30-36], where the hydrology is markedly different from previous examples given for the Western Kenyan Highlands and the Usambara Mountains, Tanzania. This is due to the high permeability of limestone rock enabling water to infiltrate into the bedrock to form cave systems, rather than converging in topographic depressions or river channels, resulting in a variable relationship between terrain and water table depth. These environments need to be studied if we are to understand the physical controls behind the formation and persistence of vector habitats in malarial regions where limestone is dominant. This is particularly poignant in east Africa where a large proportion of the population live in cities built in coastal areas dominated by coral limestones.
In karst landscapes, water reacts with calcium carbonate in limestone, to create a highly corrosive substance (carbonic acid) that exploits natural cracks and crevasses in the rock. Through dissolution widening of fractures, preferential flowpaths and conduits develop in a positive feedback loop eventually leading to the development of sink holes and cave systems (Figure 1A). Rivers and streams intercepting this landscape will be diverted underground (Figure 1A) forming groundwater aquifers, which can re-emerge downstream to feed surface pools of water as natural springs. Corrosion in karst landscapes can lead to the formation of closed bowl-shaped depressions up to 1 km in diameter known as dolines providing foci for local surface drainage [37,38] (Figure 1B). Fine-grained soil and weathered material known as terra rossa, is transported fluvially into dolines  and is deposited to form localised areas of low infiltration (Figure 1B) and, in some cases, the development of shallow water bodies [38-40].
Although the hydrology of karst landscapes is well understood, no reported attempt has been made to link this information to aquatic malaria vector habitats in sub-Saharan Africa. In this paper we examine wet and dry season spatial patterns of malaria infection rates reported by health facilities as part of routine surveillance mechanisms on Unguja, Zanzibar which is in the pre-elimination phase and thus making our findings immediately applicable to malaria elimination efforts. Spatial patterns of malaria infection rates are analysed in relation to the physical geography of the island, including its karst geomorphology, hydrology, geology and land cover in order to identify hot spots of persisting malaria transmission.
The aim of this study is to explore the influence of geology, geomorphology and hydrology on patterns of malaria transmission in a karst landscape. The focus of this study is on the limestone dominated island of Unguja, Zanzibar, United Republic of Tanzania, which is in the pre-elimination phase of its Malaria Control Programme.
Unguja is the largest island (1600 km2) of the Zanzibar Archipelago, located 40 km off the east coast of Tanzania (Figure 2). The island is underlain by Miocene sandy clay marl (Figure 3A). Alluvial deposits and laterites are found on the northwest part of the island up to 130 m above sea-level (Figure 3A). This area supports a small number of perennial rivers and numerous seasonally active streams (Figure 3A), which tend to divert into the ground once they intercept the porous Quaternary coralline limestone reef terraces (Figure 1A) that dominate the rest of the island, particularly in the east and southeast  (Figure 3). Over these limestones the landscape is typical of karst environments with the development of sink holes, caves, and doline features. A description of these features and other geographical terms, including rock types, can be found in Table 1.
Unguja receives between 1000 and 2250 mm of rainfall per year. Rainfall is strongly seasonal, typically with dry and hot weather during January and February, heavy rains from March to May, a dry season during June to September and light rains during October to December . The vegetation of the island comprises secondary forest, mangrove swamps, and degraded fallow bush. Agriculture is mainly dominated by root crops, vegetables and rain-fed rice plantations . In some areas, low infiltrating soils and the availability of groundwater supply supports irrigated rice and sugar cane plantations. Intense deforestation in advance of agricultural development has left little original vegetation with the exception of the Jozani Forest National Park  which is located within a shallow basin where the water table is relatively high.
Malaria infection rates
Data on malaria cases for the 2011 wet season (April-June) and dry season (July-September) are based on weekly summary reports of numbers of suspected, tested and treated malaria cases amongst patients from 144 health facilities (49 located on Unguja) located across the Zanzibar islands . Transmission data since 2011 was not available. Data prior to this year was not considered in the analysis because patterns of malaria transmission will be strongly affected by a campaign distributing treated bed nets and indoor residual spraying. At the health facilities, patients presenting signs and symptoms of malaria, including history of fever, were tested for malaria using either quality assured microscopic examination of Giemsa-stained thick blood films or a Rapid Diagnostic Test (Paracheck® or SD Bioline®) . To enable binomial analysis of the data, malaria positivity rates (number of patients tested positive for malaria/the number of patients tested) were classified as hotspot or non-hotspot based on whether they are above or below the mean rate, following the example of Bousema et al. . Information regarding the condition of the facility is also recorded and used as a control variable in subsequent analyses. Facility condition was categorised as very bad, bad, good and very good by the Zanzibar Ministry of Health according to the infrastructure (electricity, safe water supply, telephone), services and care offered (basic equipment, laboratory capacity, infection control), and treatments administered (HIV/AIDS, mental health, diabetes, tuberculosis) at the facility.
Physical geography data
Datasets describing the physical geography of Unguja were obtained following a study into the water resources of the island , including geology (Figure 3) , the location of streams and dolines  (Figure 4A), and soil infiltration rate  (Figure 4A). Additional geological and water resource information, including land cover type (Figure 4B), was obtained from the Zanzibar Water Authority. Elevation data from the Shuttle Radar Topography Mission (SRTM) was downloaded from United States Geological Survey online data archive Earth Explorer  (Figure 4C). This was used to calculate slope angle (Figure 4D) using ArcMap 10.1 . To facilitate spatial analyses all data were converted to raster datasets with a spatial resolution of 90 m, the rivers and dolines data were summarised using the Distance to Feature tool in ArcMap. Streams and rivers can provide aquatic habitats for malaria vector larval development [14,15,17] meaning that proximity of households to river features may increase the risk of malaria transmission. There have been no documented attempts to link the presence of doline landforms to malaria incidence, though their characteristics, as shallow inwardly draining basins [37,38], make them conducive to supporting vector larval development. The distance to feature rasters for both dolines and rivers will test whether proximity to these features, and their potential for providing malaria vector habitats, affects the intensity of malaria transmission. Values for each of the physical geography raster datasets were extracted at the location of each health facility using the Extract Multi Values to Points tool in ArcMap using bilinear interpolation of the raster cells adjacent to the central cell at the location of each facility. The resulting dataset is summarised in Table 2.
Some of the variables listed in Table 2 do not necessarily have a direct biological explanation for the presence or absence of malaria transmission hotspots. For instance, the maximum values for the variables distance to dolines and distance to perennial streams (36 km and 26.7 km respectively) far exceed the maximum flight distance of an adult anopheline mosquito . Rather, these distances represent the highest single values for any distance separating these landscape features from the nearest health facility. The statistical tests employed in this study rely on all the data encompassing the full range of distances between pairs of geographic locations where these dependent and independent variables were observed. We cannot therefore infer that any single geographic feature has a statistically significant impact on malaria risk in its own right.
Despite this, the human malaria infection indicator used this study was diagnostic positivity as recorded at health facilities where patients self-reported at their own discretion, actively travelling to their own preferred choices among these units of observation over a range of distances. The distribution of observed hot spots is therefore not equivalent to the distribution of transmission exposure determined by mosquito proliferation and dispersal alone, but also to human mobility. As such, the landscape features examined in this study could affect malaria infection risk observed at health facilities beyond the normal flight range of mosquitoes.
The relationship between the malaria positivity rates and the physical geography descriptors was explored using boosted regression tree (BRT) analysis fitted in R  using the ‘dismo’ package . This type of modelling generates a robust estimate of response variables by combining a large number of simple models, or regression trees [52-54] which has advantages over Generalised Additive Models or multiple linear regression models that relate the response variable to a number of predicting variables via a single model with well documented drawbacks, including bias in parameter estimation and inconsistencies in model selection [53,55-57].
The binary wet and dry season malaria infection hotspot data were linked to the physical geography variables using a binomial BRT model fitted with a Bernoulli distribution. The BRT models were fitted using a tree complexity of 5 and a learning rate of 0.002 which fitted a minimum recommended number of trees (>1000) while providing optimal training and cross-validation AUC scores [53,56]. A bag fraction of 0.5 was used following Elith et al. .
The BRT procedure returns a list of the predicting variables ranked by their relative contribution (%) to the model which is used to evaluate the influence of each predicting variable on the model. The relative contribution represents the number of times the variable was selected for splitting, weighed by the improvement to the model as a result of each split averaged over all trees [56,58]. Model predictive performance is evaluated using the area under the receiver operating curve (AUC) reported for training of the model and subsequent cross-validation [53,56].
The BRT modelling procedure is used in this instance to form a correlative model between hotspots of malaria infection and metrics summarising the physical geography of the island. As such, we may only infer but not directly imply causal relationships between these metrics and the abundance of malaria vectors, and of course any such statistical associations may be spurious, particularly those with modest levels of significance. Nevertheless, this robust statistical approach can be used to illustrate inherent spatial heterogeneity in the distribution of malaria transmission through a process-based understanding of the physical hydrological and biological mechanisms taking place in the landscape and how this can be exploited in future integrated malaria control initiatives.
The malaria infection data used in this study was collected as part of a routine surveillance initiative by the Zanzibar Malaria Control Program and was not collected for research purposes and is fully anonymised with no personal identifiers. These records are centred on health facilities so residence is also anonymised.
The total number of patients tested for malaria infection in the wet season (April-June) of 2011 was 40,173 with 3.2% being diagnosed with malaria. This rate decreased in the dry season (July-September) of 2011 with 1.5% of 25,858 patients being tested positive for malaria. The mean malaria infection rate across the health facilities was 4.5% for the wet season and 2.6% for the dry season. Malaria positivity rates varied in health facilities across Unguja (Figure 5) with a tendency for higher rates to be found in the south of the island and lower rates found in the north. The spatial distribution of malaria infection rates was similar in both the wet and dry season with a Pearson’s correlation, corrected for spatial autocorrelation [59,60], of 0.71 (P < 0.001).
The BRT model for wet season hotspots of malaria transmission was fitted with 1600 trees and the dry season model was fitted with 1350 trees. AUC scores were significantly better than random with a cross-validation AUC of 0.89 and 0.8 for the wet and dry season models respectively (Table 3) indicating that hotspots of malaria infection were successfully modelled using the physical geography variables outlined in Table 2.
The relative influence of each variable on the wet and dry season BRT models were similar with distance to dolines and slope angle having a high influence and soil infiltration rate having little or no influence (Table 4). Malaria transmission hotspots were negatively associated with distance to dolines which is supported by Pearson correlation coefficients of −0.33 and −0.34 (P value < 0.05) for malaria positivity rates in the dry and wet season respectively. Although, malaria transmission hotspots were not significantly univariately correlated with distance to streams, a positive relationship (more malaria further from streams) had a high influence in the BRT model. A spatially corrected Pearson’s correlation indicated that the metrics distance to dolines and distance to streams were not significantly autocorrelated with a correlation coefficient of −0.42 (P > 0.05).
The fitted functions for the predicting variables were similar for both wet and dry season models (Figure 6) although some small differences existed for the landcover variable. Specifically, areas of bushland were less associated with malaria infection hotspots in the dry season compared to the wet season. Conversely, cultivated land was more associated to malaria infection hotspots during the dry season, relative to the wet season. The scrubland landcover type had the greatest association with hotspots of malaria transmission. Geology had a similar level of influence in both the wet and dry season models. Of the different geological types, coralline and reef limestone (Q2) and the mixture of crystalline, reef and detrital limestone, with marine and fluvial sands and sandstone (Q2Q3M1) was shown to have the clearest relationship with hotspots of malaria transmission. The condition of the health facility had no influence on either the wet or dry season models meaning that subsequent inferences are independent from the state of facilities at which malaria incidence is reported.
Predicted maps of malaria hotspot probability are shown in Figure 7. Predictions for both the wet and dry seasons indicate similar patterns, with an increase in predicted hotspots in the south of the island and the northeast also demonstrating potential for hotspot occurrence. Primarily, hotspots are concentrated in areas with steep slopes close to dolines. The model follows observed hotspots for most locations across the island, although some apparent false negatives occur. For instance, hotspots occur at Fumba on northwest coast and Donge Vijibweni in the far north of the island (see Figure 2 for locations) but are not predicted by the model. Some false positives also occur, particularly in areas with steep coastal cliffs, such as Pwani, where slope angle has a high influence on both the wet and dry season BRT models.
Relatively high malaria infection rates were reported in health facilities towards the south and central part of the island with relatively low rates in the north (Figure 5). The south and central part of the island is dominated by limestones with high infiltration rates where surface water tends to be diverted underground, illustrated by the influence of coralline and reef limestone (Q2) in both the wet and dry season BRT models (Figure 6) which forms the island’s main aquifer complex . The dominance of groundwater in this area provides a temporally stable supply of water to malaria vector habitats independently from rainfall. Additionally, this area includes outcrops of dense limestone which maintain relatively high water tables feeding surface water bodies including cave wells and hand dug wells [45,61] (Figure 3B). These outcrops are represented by the rock complex Q2Q3M1 (a mixture of crystalline, reef and detrital limestone, with marine and fluvial sands and sandstone) which was shown to influence the BRT models (Figure 6). These outcrops are characterised by steep slopes, resulting in the variable slope angle being one of the most influential variables for predicting dry and wet season malaria infection hotspots (Table 4).
Distance to dolines was also an important variable for predicting hotspots of malaria infection (Table 4) as high concentrations of these landform features occur in the central and southern part of the island, providing localised depressions for the formation of surface water bodies . The locations of these landforms often coincide with areas of scrub vegetation, which was shown to be positively associated with hotspots of malaria infection (Figure 6). This landcover type typically occurs at the fringes of forested areas, such as the groundwater-fed Jozani National Forest [42,43,45], and is characterised by a fine-grained red-brown sandy top soil  that often occupies crevices in the limestone rock  providing a focus for local drainage and the development of shallow pools [38-40].
Elevation had a marginal influence on both the wet and dry season BRT models of malaria infection hotspots (Table 4). This reflects the influence of subterranean processes governing the movement of water on Unguja through the karst landscape via an aquifer network or the concentration of water bodies in localised areas of low infiltration, such as doline features infilled with fine grained material. This is in contrast to other, non-karst landscapes, such as the Western Kenyan Highlands and the Usambara Mountains, Tanzania, are dominated by the overland, or near-surface movement of water, resulting in an association between elevation and/or terrain and malaria transmission [21-28].
Distance to streams was a relatively important predicting variable (Table 4) with hotspots of transmission being associated with large distances from perennial streams and rivers (Figure 6). The positive relationship suggests that while streams on Unguja may support malaria vector larval development, their suitability is low and is probably outweighed by their ability to provide habitats in the surrounding floodplains due to flooding. This is likely to be particularly important in the wet season when increased flows not only increase the suitability of the floodplains, but also make habitats located within river channels less suitable for anopheline larvae due to habitat flushing [62,63] and an intolerance of fast flowing water [2,63,64]. Additionally, pollution, due to the absence of sewerage infrastructure , may also negatively influence the ability of rivers to support vector larval development . Additional work focussed on this habitat type is needed to establish their ability to support malaria vector larvae particularly where channels of similar dimensions (approximately 10 m in width and 5 m deep) in the Kilombero Valley southern Tanzania have been found to support vector larval development within chains of pools that form in the river bed once the river ceases to flow in the dry season . The contrast between these two distinct ecosystems merits consideration: In contrast to Unguja, in the Kilombero Valley floodplain, extensive aquatic habitats are formed during the wet season, when the floodplains of these channels are inundated with water during peak river flows, while in the dry season, river channels provide suitable habitats in chains of pools that form in the river bed once the river ceases to flow . Specifically, this function (the formation of habitats in river channels during the dry season) has the potential to support dry season transmission, making it a crucial process for maintaining malaria endemicity throughout the year. As such, the potential for streams and rivers on Unguja to provide productive malaria vector habitats cannot be dismissed and requires further detailed investigation.
We found that the spatial pattern of malaria infection remained relatively static during the 2011 wet and dry seasons. Given the pronounced seasonality to rainfall patterns in the Zanzibar islands this suggests that aquatic malaria vector habitats are chiefly fed by groundwater sources which persist throughout the year. The influence of groundwater processes dominating the formation and persistence of water bodies means that geology and hydrological processes need to be understood if the dynamics of aquatic malaria vector habitats are to be mapped and targeted for larval source interventions.
Malaria infection rates on Unguja demonstrated spatial heterogeneity which can reduce the effectiveness of control strategies where resources are wasted on areas with low or negligible rates of malaria . It is therefore crucial to identify geographic foci of disease transmission for targeting interventions [44,68]. Bejon et al.  noted the uncertainty surrounding the temporal stability of disease hotspots at a vector dispersal level (0.5-1 km) in Kenya. For Unguja, patterns of malaria infection were similar in both the 2011 wet and dry seasons, increasing the potential for efficient targeting of interventions. Targeting malarial habitats in this way may help to eliminate vector species, such as An. arabiensis, that are not vulnerable to conventional indoor-based interventions  and help to eliminate malaria on Unguja.
Relatively low malaria infection rates were found in the central northern part of the island in 2011. The availability of surface water in this region of the island, owing to the dominance of low infiltrating sandy clay soils and recent alluvial deposits (Q1, Q3 and M3) [45,46], is exploited by widespread cultivation of crops including sugar cane and rice. Saturated rice fields often support high vector abundances, particularly for species such as An. arabiensis which are the dominant vector on the island . Despite this, wards within this region demonstrated low infection rates in both the wet and dry seasons (Figure 5). Reasons behind the reduction in malaria positivity rate in this region of island remain unclear. Some authors have purported a ‘paddies paradox’ whereby reduced malaria transmission has been associated with irrigated crop production in areas with stable transmission due to a focus of indoor-based vector control interventions [69,70]. However, we do not believe this is the case on Unguja where there are no distinct geographical differences in the use of LLINs and IRS [9,10]. From our study we hypothesise that these habitats are in fact relatively unproductive for vector mosquitoes, but this needs to be tested by entomological survey. Although irrigated areas showed no association with malaria, hotspots were more likely to be found in areas of cultivated land during the dry season. These areas are likely to be occupied at night during certain points of the year to protect seedlings from birds, or farmers working late for ploughing or harvesting, increasing their risk of exposure to malarial mosquitoes . The apparent seasonal independence of the spatial distribution of malaria transmission on Unguja is testament to the notion that dry season transmission is more important for sustaining endemicity than wet season transmission and should, therefore, provide encouraging opportunities for focussing seasonal malaria control activities at exactly the time of year when they are easiest to implement [72,73].
Despite the general trend of relatively low malaria transmission rates in the northern part of Unguja, two hotspots of malaria infection occur in this region: at Fumba on the northwest coast (Figure 8) and Donge Vijibweni in the far north of the island (Figure 9). Both health facilities are located close to relatively small doline (400 m wide) features which have not been identified in the environmental layers obtained from the Zanzibar Water Authority. In addition, both sites are located in close proximity to streams, the latter having a negative association with malaria infection rates for most hotspots on the island. Such areas could be mapped using high spatial resolution (<2 m) aerial/satellite imagery, or through using remote sensing systems such as L-band Radar , with contextual information to differentiate from anthropogenically induced wetted areas, such as irrigated rice paddies. Interestingly, both these locations are known hotspots of malaria transmission highlighted at the onset of the Zanzibar Malaria Control Programme’s Malaria Epidemic Early Detection System  and remain despite extensive distribution of indoor-based interventions.
This study identified two key processes that are enabling the development and persistence of malaria vector habitats, particularly in the south and central part of the island where malaria infection rates remain high. Firstly, outcrops of dense crystalline, reef and detrital limestone (M1), intercepting the main groundwater aquifer complex, maintain relatively high water tables, increasing the potential for surface water bodies to persist [43,45]. Secondly, dolines, which are naturally occurring depressions in limestone dominated karst landscapes, can become plugged with fine-grained material which is washed in from the surrounding area leading to pockets of low infiltrating soils providing a focus for local drainage and the development of shallow pools [38-40] that are likely to be suitable for vector oviposition and larval development . Using available static geological maps we were able to identify areas where these two processes exist representing a significant tool for informing larval source intervention strategies.
The accumulation of fine-grained material and subsequent development of shallow water bodies is also prominent at the fringes of the Jozani National Forest [42,43,45]. Deforestation, which has been prevalent on the island , is likely to increase the availability of fine-grained material , potentially increasing the coverage of low infiltrating soils with a propensity to support surface water bodies and potential habitats. This process been linked to a rise in historical malaria cases in ancient Rome [77,78], as well as more recent examples [79,80] and should provide a lesson for present and future land managers in malarial regions [79,81].
The malaria infection rates used in this study are based on summary reports from health facilities. This makes the assumption that the landscape, and therefore physical geography variables, is the same where the patient lives. Given the heterogeneous nature of the landscape across Unguja this is likely to introduce uncertainty. Using information regarding the specific location of patient dwellings would provide a more precise indication of the role of physical geography on the risk of being infected with malaria. Furthermore, this may help to identify specific breeding habitats prolific in malaria vector production which can be used to target larval source management interventions. Data from a greater number of sites will also improve the spread across the different physical geography components across the island. For instance, few health facilities were located close to mangroves and therefore this landcover type is not identified in our study as being as related to malaria infection. In this specific case, the analysis may not take into account malaria transmission due to salt-tolerant vectors found in mangrove forests, such An. merus . This is not assumed to have a significant impact on the study findings because An. merus contributes a relatively small proportion (5%) of adult mosquitos sampled on Unguja in 2010 . Additionally, the data used in this study does not take into account patients that contracted malaria away from their dwelling. Information regarding the travel history of the patient should be gathered to take this aspect into account, although research has indicated a relatively low risk of importing malaria from mainland Tanzania .
Our study indicates that proximity to streams is negatively correlated with hotspots of malaria infection. Despite this correlation we do not suggest that habitats located in river channels do not represent important breeding sites. In particular, whereas steep sides, strong flows, an abundance of predators and other contributing factors have the potential to make streams less likely to support malaria vector habitats during the wet season [2,63,64], streams during periods of low flow have been shown to support vector larval habitats, particularly in during the dry season in ephemeral channels when the river stops flowing, forming small, shallow water bodies in the river bed . Additionally, a number of sources indicate the importance of shallow pools of water forming along the periphery of rivers [5,13-18,84], particularly following periods of overbank flooding, with the potential to support vector larval development. To this end, the function of streams on Unguja on supporting malaria vector populations needs to be explored further.
This study has demonstrated the effects of hydrology and geology on the distribution of malaria transmission hotspots. It shows that patterns of malaria infection gathered using a routine programmatic surveillance platform demonstrate disease heterogeneity on Unguja, Zanzibar, United Republic of Tanzania. Furthermore, spatial patterns of malaria infection can be predicted using static data, such as geology and slope angle, which underpin the hydrology of the island, helping to develop a framework for targeting malaria vector habitats for intervention strategies. These methods can be applied at the national scale where similar surveillance programmes exist, such as Zambia  and the Solomon Islands , to provide fine-scale information for malaria elimination campaigns. All of the data used in this work was freely available, providing an example of the great potential for replicating the methods outlined in this paper for similar landscapes burdened by malaria.
Ferguson H, Domhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla M, et al. Ecology: a prerequisite for malaria elimination and eradication. PLoS Med. 2010;7(8):1–7.
Ndenga BA, Simbauni JA, Mbugi JP, Githeko AK, Fillinger U. Productivity of malaria vectors from different habitat types in the Western Kenya highlands. PLoS One. 2011;6(4):e19473.
Gouagna LC, Rakotondranary M, Boyer S, Lemperiere G, Dehecq JS, Fontenille D. Abiotic and biotic factors associated with the presence of Anopheles arabiensis immatures and their abundance in naturally occurring and man-made aquatic habitats. Parasites Vectors. 2012;5(1):96.
Imbahale SS, Githeko A, Mukabana WR, Takken W. Integrated mosquito larval source management reduces larval numbers in two highland villages in western Kenya. BMC Public Health. 2012;12(1):362.
Zhou G, Munga S, Minakawa N, Githeko AK, Yan G. Spatial relationship between adult malaria vector abundance and environmental factors in western Kenya highlands. Am J Trop Med Hyg. 2007;77(1):29–35.
Killeen GF, Seyoum AK, Knols BGJ. Rationalizing Historical successes of malaria control in Africa in terms of mosquito resource availability management. Am J Trop Med Hyg. 2004;71(2 suppl):87–93.
Fillinger U, Lindsay SW. Suppression of exposure to malaria vectors by an order of magnitude using microbial larvicides in rural Kenya. Trop Med Int Health. 2006;11(11):1629–42.
Mbare O, Lindsay S, Fillinger U. Aquatain(R) Mosquito Formulation (AMF) for the control of immature Anopheles gambiae sensu stricto and Anopheles arabiensis: dose-responses, persistence and sub-lethal effects. Parasites Vectors. 2014;7(1):438.
Haji KA, Khatib BO, Smith S, Ali AS, Devine GJ, Coetzee M, et al. Challenges for malaria elimination in Zanzibar: pyrethroid resistance in malaria vectors and poor performance of long-lasting insecticide nets. Parasites Vectors. 2013;6(1):82.
Zanzibar Malaria Control Programme. Zanzibar malaria epidemic early detection system biannual report. 2011.
Smith DL, Perkins TA, Tusting LS, Scott TW, Lindsay SW. Mosquito population regulation and larval source management in heterogeneous environments. PLoS One. 2013;8(8):e71247.
Hardy AJ, Gamarra JGP, Cross DE, Macklin MG, Smith MW, Kihonda J, et al. Habitat hydrology and geomorphology control the distribution of malaria vector larvae in Rural Africa. PLoS One. 2013;8(12):1–13.
Smith MW, Macklin MG, Thomas CJ. Hydrological and geomorphological controls of malaria transmission. Earth Sci Rev. 2013;116:109–27.
Amerasinghe PH, Amerasinghe FP, Konradsen F, Fonseka KT, Wirtz RA. Malaria vectors in a traditional dry zone village in Sri Lanka. Am J Trop Med Hyg. 1999;60(3):421–9.
Mutuku FM, Alaii JA, Bayoh MN, Gimnig JE, Vulule JM, Walker ED, et al. Distribution, description, and local knowledge of larval habitats of Anopheles Gambiae s.l. in a village in Western Kenya. Am J Trop Med Hyg. 2006;74(1):44–53.
Oesterholt M, Bousema J, Mwerinde O, Harris C, Lushino P, Masokoto A, et al. Spatial and temporal variation in malaria transmission in a low endemicity area in northern Tanzania. Malar J. 2006;5(1):98.
Van Der Hoek W, Konradsen F, Amerasinghe PH, Perera D, Piyaratne M, Amerasinghe FP. Towards a risk map of malaria for Sri Lanka: the importance of house location relative to vector breeding sites. Int J Epidemiol. 2003;32(2):280–5.
Bøgh C, Lindsay SW, Clarke SE, Dean A, Jawara M, Pinder M, et al. High spatial resolution mapping of malaria transmission risk in the Gambia, West Africa, using Landsat TM satellite imagery. Am J Trop Med Hyg. 2007;76(5):875–81.
Thomas CJ, Cross DE, Bøgh C. Landscape movements of anopheles gambiae malaria vector mosquitoes in Rural Gambia. PLoS One. 2013;8(7):e68679.
Bomblies A. Modeling the role of rainfall patterns in seasonal malaria transmission. Climatic Change. 2012;112(3–4):673–85.
Dobson M. “Marsh fever”—the geography of malaria in England. J Hist Geogr. 1980;6(4):357–89.
Staedke SG, Nottingham EW, Cox J, Kamya MR, Rosenthal PJ, Dorsey G. Short report: proximity to mosquito breeding sites as a risk factor for clinical malaria episodes in an urban cohort of Ugandan children. Am J Trop Med Hyg. 2003;69(3):244–6.
Ernst KC, Lindblade KA, Koech D, Sumba PO, Kuwuor DO, John CC, et al. Environmental, socio‐demographic and behavioural determinants of malaria risk in the western Kenyan highlands: a case–control study. Trop Med Int Health. 2009;14(10):1258–65.
Protopopoff N, Van Bortel W, Marcotty T, Van Herp M, Maes P, Baza D, et al. Spatial targeted vector control in the highlands of Burundi and its impact on malaria transmission. Malar J. 2007;6(1):158.
Balls MJ, Bødker R, Thomas CJ, Kisinza W, Msangeni HA, Lindsay SW. Effect of topography on the risk of malaria infection in the Usambara Mountains, Tanzania. Trans R Soc Trop Med Hyg. 2004;98(7):400–8.
Cohen J, Ernst C, Lindblade K, Vulule J, John C, Wilson M. Local topographic wetness indices predict household malaria risk better than land-use and land-cover in the western Kenya highlands. Malar J. 2010;9(328):1–10.
Minakawa N, Seda P, Yan G. Influence of host and larval habitat distribution on the abundance of African malaria vectors in western Kenya. Am J Trop Med Hyg. 2002;67(1):32–8.
Minakawa N, Munga S, Atieli F, Mushinzimana E, Zhou G, Githeko AK, et al. Spatial distribution of anopheline larval habitats in Western Kenyan highlands: effects of land cover types and topography. Am J Trop Med Hyg. 2005;73(1):157–65.
Hewlett JD, Hibbert AR. Factors affecting the response of small watersheds to precipitation in humid areas. In: Sopper WE, Lull HW, editors. International symposium on forest hydrology. Oxford: Pergammon; 1967. p. 275–90.
Kobayashi J, Nambanya S, Miyagi I, Vanachone B, Manivong K, Koubouchan T, et al. Collection of anopheline mosquitos in three villages endemic for malaria in Khammouane, Lao PDR. Southeast Asian J Trop Med Public Health. 1997;28:615–20.
Obsomer V, Defourny P, Coosemans M. The Anopheles dirus complex: spatial distribution and environmental drivers. Malar J. 2007;6(1):26.
Toma T, Miyagi I, Okazawa T, Kobayashi J, Saita S, Tuzuki A, et al. Entomological surveys on malaria in Khammouane Province, Lao PDR, in 1999 and 2000. Southeast Asian J Trop Med Public Health. 2002;33(3):532–46.
Manguin S, Garros C, Dusfour I, Harbach R, Coosemans M. Bionomics, taxonomy, and distribution of the major malaria vector taxa of Anopheles subgenus Cellia in Southeast Asia: an updated review. Infect Genet Evol. 2008;8(4):489–503.
O'Loughlin S, Somboon P, Walton C. High levels of population structure caused by habitat islands in the malarial vector Anopheles scanloni. Heredity. 2007;99(1):31–40.
Bugoro H, Cooper RD, Butafa C, Iro'ofa C, Mackenzie DO, Chen C-C, et al. Bionomics of the malaria vector Anopheles farauti in Temotu Province, Solomon Islands: issues for malaria elimination. Malar J. 2011;10:133.
Mwangangi JM, Mbogo CM, Orindi BO, Muturi EJ, Midega JT, Nzovu J, et al. Shifts in malaria vector species composition and transmission dynamics along the Kenyan coast over the past 20 years. Malar J. 2013;12:13.
Summerfield MA. Global geomorphology: an introduction to the study of landforms: Longman Scientific & Technical. Wiley; 1991.
Sweeting M. Karst landforms. London: The Macmillan Press Ltd; 1972.
Ford DC, Williams P. Karst hydrogeology and geomorphology. Chichester: John Wiley & Sons; 2007.
Piñol J, Beven K, Freer J. Modelling the hydrological response of Mediterranean catchments, Prades. Catalonia. The use of distributed models as aids to hypothesis formulation. Hydrol Processes. 1997;11(9):1287–306.
Johnson J. A review of the hydrogeology of Zanzibar (Tanzania). New York: United Nations; 1984.
Siex KS, Struhsaker TT. Ecology of the Zanzibar red colobus monkey: demographic variability and habitat stability. Int J Primatol. 1999;20(2):163–92.
Sikat L. Assessing the spatial and temporal characteristics of groundwater recharge in Zanzibar: towards the optimal management of grounwater resources. Twente, The Netherlands: University of Twente; 2011.
Bousema T, Drakeley C, Gesase S, Hashim R, Magesa S, Mosha F, et al. Identification of hot spots of malaria transmission for targeted malaria control. J Infect Dis. 2010;201(11):1764–74.
Colbert G, Wagner BH, Pinther M. Hydrogeological map of Zanzibar. New York, USA: United Nations; 1987.
Hettige ML. Land evaluation and land suitability classification - Unguja and Pemba islands. Rome: UN Food and Agricultural Organisation; 1990.
Earth Explorer [http://earthexplorer.usgs.gov]
ESRI. ArcGIS desktop: release 10.2.1. Redlands, CA: Environmental Systems Research Institute; 2013.
Costantini C, LI SG, Torre AD, Sagnon NF, Coluzzi M, Taylor CE. Density, survival and dispersal of Anopheles gambiae complex mosquitoes in a West African Sudan savanna village. Med Vet Entomol. 1996;10(3):203–19.
R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2012.
Hijmans RJ, Phillips S, Leathwick J, Elith J. Species distribution modeling. R package version 0.8-11. 2013.
Buston PM, Elith J. Determinants of reproductive success in dominant pairs of clownfish: a boosted regression tree analysis. J Anim Ecol. 2011;80(3):528–38.
Leathwick J, Elith J, Francis M, Hastie T, Taylor P. Variation in demersal fish species richness in the oceans surrounding New Zealand: an analysis using boosted regression trees. Mar Ecol Prog Ser. 2006;321:267–81.
Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, et al. The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis. Parasites Vectors. 2010;3:117.
Whittingham MJ, Stephens PA, Bradbury RB, Freckleton RP. Why do we still use stepwise modelling in ecology and behaviour? J Anim Ecol. 2006;75(5):1182–9.
Elith J, Leathwick JR, Hastie T. A working guide to boosted regression trees. J Anim Ecol. 2008;77(4):802–13.
Burnham KP, Anderson DR. Model selection and multi-model inference: a practical information-theoretic approach. Verlag: Springer; 2002.
Friedman JH, Meulman JJ. Multiple additive regression trees with application in epidemiology. Stat Med. 2003;22(9):1365–81.
Clifford P, Richardson S, Hemon D. Corrected Pearson’s correlation for spatial autocorrelation. Biometrics. 1989;49:305–14.
Osorio F, Vallejos R. Package for analysis of spatial data. R package version 0.2. 2012.
Vuai S. Geochemical characteristics of spaleotherm formation in caves from Zanzibar Island, Tanzania. vol. 1: OMICS; 2012: 1-5.
Paaijmans KP, Wandago MO, Githeko AK, Takken W. Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007;2(11):e1146.
Gillies MT, De Meillon B. The Anophelinae of Africa South of the Sahara. Johannesburg: South African Institute for Medical Research; 1968.
Thomas C, Lindsay S. Local-scale variation in malaria infection amongst rural Gambian children estimated by satellite remote sensing. Trans R Soc Trop Med Hyg. 2000;94(2):159–63.
Moynihan MA, Baker DM, Mmochi AJ. Isotopic and microbial indicators of sewage pollution from Stone Town, Zanzibar, Tanzania. Mar Pollut Bull. 2012;64(7):1348–55.
Sattler MA, Mtasiwa D, Kiama M, Premji Z, Tanner M, Killeen GF, et al. Habitat characterization and spatial distribution of Anopheles sp. mosquito larvae in Dar es Salaam (Tanzania) during an extended dry period. Malar J. 2005;4(1):4.
Bejon P, Williams T, Liljander A, Noor A, Wambua J, Ogada E, et al. Stable and unstable malaria hotspots in longitudinal cohort studies in Kenya. PLoS Med. 2010;7(7):1–14.
Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, et al. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012;9(1):e1001165.
Ijumba J, Shenton F, Clarke S, Mosha F, Lindsay S. Irrigated crop production is associated with less malaria than traditional agricultural practices in Tanzania. Trans R Soc Trop Med Hyg. 2002;96(5):476–80.
Ijumba J, Lindsay S. Impact of irrigation on malaria in Africa: paddies paradox. Med Vet Entomol. 2001;15(1):1–11.
Dunn CE, Le Mare A, Makungu C. Malaria risk behaviours, socio-cultural practices and rural livelihoods in southern Tanzania: implications for bednet usage. Soc Sci Med. 2011;72(3):408–17.
Killeen GF. A second chance to tackle African malaria vector mosquitoes that avoid houses and don’t take drugs. Am J Trop Med Hyg. 2013;88(5):809–16.
Alonso PL, Bassat Q, Binka F, Brewer T, Chandra R, Culpepper J, et al. A research agenda for malaria eradication: drugs. PLoS Med. 2010;8(1):1–9.
Washino RK, Wood BL. Application of remote sensing to vector arthropod surveillance and control. Am J Trop Med Hyg. 1993;50:134–44.
Zanzibar Malaria Control Programme. Zanzibar malaria epidemic early detection system biannual report. 2009.
Sidle RC, Ziegler AD, Negishi JN, Nik AR, Siew R, Turkelboom F. Erosion processes in steep terrain - truths, myths, and uncertainties related to forest management in Southeast Asia. Forest Ecol Manag. 2006;224(1–2):199–225.
O'Sullivan L, Jardine A, Cook A, Weinstein P. Deforestation, mosquitoes, and ancient Rome: lessons for today. Bioscience. 2008;58(8):756–60.
Angel JL. Ecology and population in the eastern Mediterranean. World Archaeol. 1972;4(1):88–105.
Vittor AY, Gilman RH, Tielsch J, Glass G, Shields T, Lozano WS, et al. The effect of deforestation on the human-biting rate of Anopheles darlingi, the primary vector of falciparum malaria in the Peruvian Amazon. Am J Trop Med Hyg. 2006;74(1):3–11.
Yasuoka J, Levins R. Impact of deforestation and agricultural development on anopheline ecology and malaria epidemiology. Am J Trop Med Hyg. 2007;76(3):450–60.
Pattanayak S, Dickinson K, Corey C, Murray B, Sills E, Kramer R. Deforestation, malaria, and poverty: a call for transdisciplinary research to support the design of cross-sectoral policies. Sustainability. 2006;2(2):45–56.
Kipyab P, Khaemba B, Mwangangi J, Mbogo C. The bionomics of Anopheles merus (Diptera: Culicidae) along the Kenyan coast. Parasites Vectors. 2013;6(1):37.
Tatem AJ, Qiu Y, Smith DL, Sabot O, Ali AS, Moonen B. The use of mobile phone data for the estimation of the travel patterns and imported Plasmodium falciparum rates among Zanzibar residents. Malar J. 2009;8:287.
Bøgh C, Clarke S, Jawara M, Thomas C, Lindsay S. Localized breeding of the Anopheles gambiae complex (Diptera: Culicidae) along the River Gambia, West Africa. Bull Entomol Res. 2003;93(4):279–87.
Steketee RW, Sipilanyambe N, Chimumbwa J, Banda JJ, Mohamed A, Miller J, et al. National malaria control and scaling up for impact: the Zambia experience through 2006. Am J Trop Med Hyg. 2008;79(1):45–52.
Kelly GC, Hii J, Batarii W, Donald W, Hale E, Nausien J, et al. Modern geographical reconnaissance of target populations in malaria elimination zones. Malar J. 2010;9:289.
This work was funded by the UK Natural Environmental Research Council (NERC) grant number NE/H022740/1 and the Bill & Melinda Gates Foundation. We would like to thank all our colleagues in Zanzibar, notably the Zanzibar Malaria Control Programme and the Zanzibar Water Authority. In the UK we would like to thank the support and guidance of Javier G P Gamarra, Dónall Cross and Richard Lucas. At USAID we thank David Ngilangwa and Peter McElroy. A special thanks to Leonila Bron Sikat for providing data concerning the water resources of the island.
The authors declare that they have no competing interests.
AH, ZM, SD, GK, MGM, SM, MS, and CT designed the study, prepared the manuscript and interpreted the results. AA, AA-M and MM designed and implemented the disease transmission surveillance. AH and ZM obtained the geographical data and carried out statistical analyses. All authors read and approved the final version of this manuscript.