Using the SaTScan method to detect local malaria clusters for guiding malaria control programmes
© Coleman et al; licensee BioMed Central Ltd. 2009
Received: 15 December 2008
Accepted: 17 April 2009
Published: 17 April 2009
Mpumalanga Province, South Africa is a low malaria transmission area that is subject to malaria epidemics. SaTScan methodology was used by the malaria control programme to detect local malaria clusters to assist disease control planning. The third season for case cluster identification overlapped with the first season of implementing an outbreak identification and response system in the area.
SaTScan™ software using the Kulldorf method of retrospective space-time permutation and the Bernoulli purely spatial model was used to identify malaria clusters using definitively confirmed individual cases in seven towns over three malaria seasons. Following passive case reporting at health facilities during the 2002 to 2005 seasons, active case detection was carried out in the communities, this assisted with determining the probable source of infection. The distribution and statistical significance of the clusters were explored by means of Monte Carlo replication of data sets under the null hypothesis with replications greater than 999 to ensure adequate power for defining clusters.
Results and discussion
SaTScan detected five space-clusters and two space-time clusters during the study period. There was strong concordance between recognized local clustering of cases and outbreak declaration in specific towns. Both Albertsnek and Thambokulu reported malaria outbreaks in the same season as space-time clusters. This synergy may allow mutual validation of the two systems in confirming outbreaks demanding additional resources and cluster identification at local level to better target resources.
Exploring the clustering of cases assisted with the planning of public health activities, including mobilizing health workers and resources. Where appropriate additional indoor residual spraying, focal larviciding and health promotion activities, were all also carried out.
Malaria is the most important parasitic disease of humans. Over three billion people live in malarious areas and the disease causes over 500 million cases with one to three million deaths per year [1, 2]. An estimated one hundred million people in Africa are at risk of malaria epidemics . In common with most vector-borne infectious diseases, malaria is heterogeneous in its distribution in time and space [4–6], and incidence can vary greatly between districts, towns and villages. This heterogeneity is affected by patterns of malaria vector distribution, human-vector contact, human host behavioural factors, house construction, and malaria prevention methods used [5, 7–10].
Characterization of malaria heterogeneity may allow prioritization of risk areas and allow focused control interventions. The ability to identify localized malaria clusters in remote highland areas of east Africa facilitated early intervention in the absence of early warning systems, as these malaria "hotspots" remained constant in epidemic and non epidemic years .
A number of models have been developed for describing malaria spatial distribution and seasonality [12–15] transmission [8, 16–18], mosquito distribution [19, 20] and risk factors associated with malaria [21–23]. However many malaria endemic environments are resource limited and tools to assist decision support need to take account of these limitations.
Malaria remains a public health problem in north-eastern South Africa, in the low altitude provinces of KwaZulu-Natal, Limpopo and Mpumalanga [24, 25]. Malaria risk is low compared to other hyper- and holo-endemic areas of sub-Saharan Africa and naturally acquired immunity does not develop in the local population .
The burden of malaria is well described in this region as definitive diagnosis using rapid diagnostic tests (RDTs) and mandatory reporting of malaria cases is universally practised in the public health system, which manages the vast majority of malaria cases [27, 28]. All confirmed malaria cases are entered into a computerized malaria surveillance system.
This paper reports on the analysis of malaria notification data from 2002–2005 from Mpumalanga Province to determine local clustering of cases and how that was used to direct local control efforts and enhance the district outbreak identification and response system .
Mpumalanga Province, which borders Mozambique and Swaziland, is a predominantly rural area with a population of approximately four million people. Malaria occurrence in this area is seasonal during the wet and humid months of October to April. Ninety-five percent of malaria cases are due to Plasmodium falciparum infection. This area has historically experienced malaria outbreaks and epidemics with relatively high mortality [29, 30]. The main vector in this region is Anopheles arabiensis and vector control is predominantly by indoor residual spraying (IRS) with DDT [31, 32].
Original town maps acquired in 2000 were upgraded from digitized aerial photographs produced in 2002 (Cadnet, Nelspruit, South Africa) and updates were converted into MapInfo version 6.4 files (MapInfo Corporation, New York, USA) with Geomedia software (Symmetry Systems Inc., New York, USA). Unique stand (household) numbers were allocated during the digitizing process. A stand may include one or more structures belonging to a family unit on a designated piece of land.
Active malaria case investigation
Malaria case investigators followed up all malaria cases reported at local health facilities for these towns as part of malaria control program activities for the 2002/2003, 2003/2004 and 2004/2005 seasons. A malaria case is defined as a person diagnosed either by positive rapid test or blood slide. An investigator would identify the probable source of infection by confirming the travel history of the patient recorded on the notification.
Spatial and temporal clusters
SaTScan™ software, version 5.1.3, using the Kulldorf method of retrospective space-time permutation and the Bernoulli purely spatial model [34, 35], was used by malaria control programme operational staff to detect malaria clusters in individual towns in the three seasons under investigation and over the combined time period. Individual malaria cases were used and recorded against source of infection households. Households who did not seek health care for malaria or tested negative during active case detection were used as controls for the Bernoulli method. More than one case can be reported per household. A household in these areas is typically defined as a family unit with a single land owner.
The circular scan statistic is isotopic with respect to the rotation of the geographical area http://www.satscan.org. This method has previously been validated for plotting and understanding local malaria time-space-clusters [21, 36–39].
Observed cases in a cluster were compared to the distribution of expected cases if spatial and temporal locations of all cases were independent. The model adjusts for entirely spatial or entirely temporal clusters. With spatial adjustment, time remained dormant and during temporal analysis seasons were considered. The distribution and statistical significance of the clusters were explored by means of Monte Carlo replication of data sets under the null hypothesis with replications greater than 999 to ensure adequate power for defining clusters http://www.satscan.org. Clusters were prioritized for public health action according to their statistical significance. Hard copy town maps showing cluster areas were distributed to malaria field staff investigating cases allowing the coordination of intervention efforts within communities.
Cluster and outbreak profile by town, 2002/2003 – 2004/2005 season
Number of space time
Number of space
(2002/2003 to 2004/2005)
Cluster detection to guide control activities
The third season for case cluster identification, 2004/2005, overlapped with the first season of implementation of an outbreak identification and response system in the Nkomazi area . Thambokulu and Albertsnek triggered outbreak declarations during this season (Table 1). Analysis demonstrated clusters in both towns during this season. Real time cluster detection, monitoring as cases occur, triggered at the same time as the binomial outbreak identification thresholds .
The physical plotting of cases by household during this season (2004/2005) as part of the cluster identification activities at district level and on hard copy maps provided a visual distribution of risk in the towns and assisted with the timely planning of public health activities especially logistics, including active case detection, early diagnosis and treatment of positive cases in the areas of the clusters, additional indoor residual spraying, focal larviciding and health promotion activities.
Application of the SaTScan method successfully identified malaria clusters and clearly demonstrated malaria risk heterogeneity at local level. Four towns in this study experienced spatial clusters and two produced space-time clusters in the 2004/2005 malaria season, the latter resulting in targeted local control efforts. The high rates of treatment-seeking behaviour at primary health care level in Mpumalanga  combined with the use of passive notification and active case detection, strengthens disease surveillance and provides good quality data for cluster identification. SaTScan cluster identification has also proven valuable for targeting control strategies in the Kenyan highlands . Over a ten week period during a 2002 epidemic, spatial targeting of vector control interventions reduced the abundance of Anopheles mosquitoes. The investigation of causes for clusters within the areas reported here was not explored. The hypo-endemic nature of malaria in these areas does predict the unstable patterns of disease occurrence and more specifically outbreaks that can be clustered. The proximity of these areas to higher malaria risk areas such as Mozambique allows for the import of parasites and eventually local transmission.
Spatial clustering of infectious disease is enjoying renewed interest, particularly in areas of limited resources. Statistical methods have been used to investigate spatial clustering of dengue  encephalitis  and sleeping sickness but the application to malaria has been limited .
Gaudart et al  compared the oblique decision tree model, a complex statistical technique with Kulldorff's SaTScan™ cluster technique, which was used in this study. Gaudart et al  produced similar results using both methods in a village in West Africa to identify malaria risk clusters. This investigation confirmed the usefulness of the Kulldorff's scan statistic [44–46]. Due to the circular isotopic technique of Kulldorff's SaTScan™ it is a useful tool to detect clusters but has limitations on detecting irregular shaped clusters due to its fixed scan window [47, 48].
The seven Mpumalanga towns under cluster surveillance were included in a malaria outbreak identification and response system based on formal case reporting . There was strong concordance between recognized local clustering and outbreak identification in specific towns, with Albertsnek and Thambokulu both reporting malaria outbreaks in the same season as the time-space clusters. This synergy allows mutual validation of the two systems in confirming outbreaks, which demand additional resources, and cluster identification that could better target these resources within the affected town.
Special thanks and acknowledgement to the field personnel from the Malaria Control Programme in Mpumalanga province for their dedication and commitment to the cluster and outbreak activities.
- Guerra CA, Hay SI, Lucioparedes LS, Gikandi PW, Tatem AJ, Noor AM, Snow RW: Assembling a global database of malaria parasite prevalence for the Malaria Atlas Project. Malar J. 2007, 6: 17-10.1186/1475-2875-6-17.PubMed CentralView ArticlePubMedGoogle Scholar
- Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434: 214-217. 10.1038/nature03342.PubMed CentralView ArticlePubMedGoogle Scholar
- World Health Organization: Malaria epidemics: forecasting, prevention, early detection and control. From policy to practice. 2003, Leysin, Switzerland, WHOGoogle Scholar
- Carter R, Mendis KN, Roberts D: Spatial targeting of interventions against malaria. Bull World Health Organ. 2000, 78: 1401-1411.PubMed CentralPubMedGoogle Scholar
- Gamage-Mendis AC, Carter R, Mendis C, De Zoysa AP, Herath PR, Mendis KN: Clustering of malaria infections within an endemic population: risk of malaria associated with the type of housing construction. Am J Trop Med Hyg. 1991, 45: 77-85.PubMedGoogle Scholar
- Greenwood BM: The microepidemiology of malaria and its importance to malaria control. Trans R Soc Trop Med Hyg. 1989, 83: 25-29. 10.1016/0035-9203(89)90599-3.View ArticlePubMedGoogle Scholar
- Snow RW, Gouws E, Omumbo J, Rapuoda B, Craig MH, Tanser FC, le Sueur D, Ouma J: Models to predict the intensity of Plasmodium falciparum transmission: applications to the burden of disease in Kenya. Trans R Soc Trop Med Hyg. 1998, 92: 601-606. 10.1016/S0035-9203(98)90781-7.View ArticlePubMedGoogle Scholar
- Thomas CJ, Lindsay SW: Local-scale variation in malaria infection amongst rural Gambian children estimated by satellite remote sensing. Trans R Soc Trop Med Hyg. 2000, 94: 159-163. 10.1016/S0035-9203(00)90257-8.View ArticlePubMedGoogle Scholar
- Thompson R, Begtrup K, Cuamba N, Dgedge M, Mendis C, Gamage-Mendis A, Enosse SM, Barreto J, Sinden RE, Hogh B: The Matola malaria project: a temporal and spatial study of malaria transmission and disease in a suburban area of Maputo, Mozambique. Am J Trop Med Hyg. 1997, 57: 550-559.PubMedGoogle Scholar
- Trape JF, Lefebvre-Zante E, Legros F, Ndiaye G, Bouganali H, Druilhe P, Salem G: Vector density gradients and the epidemiology of urban malaria in Dakar, Senegal. Am J Trop Med Hyg. 1992, 47: 181-189.PubMedGoogle Scholar
- Ernst KC, Adoka SO, Kowuor DO, Wilson ML, John CC: Malaria hotspot areas in a highland Kenya site are consistent in epidemic and non-epidemic years and are associated with ecological factors. Malar J. 2006, 5: 78-10.1186/1475-2875-5-78.PubMed CentralView ArticlePubMedGoogle Scholar
- Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.View ArticlePubMedGoogle Scholar
- Hay SI, Snow RW, Rogers DJ: Predicting malaria seasons in Kenya using multitemporal meteorological satellite sensor data. Trans R Soc Trop Med Hyg. 1998, 92: 12-20. 10.1016/S0035-9203(98)90936-1.View ArticlePubMedGoogle Scholar
- Omumbo J, Ouma J, Rapuoda B, Craig MH, le Sueur D, Snow RW: Mapping malaria transmission intensity using geographical information systems (GIS): an example from Kenya. Ann Trop Med Parasitol. 1998, 92: 7-21.View ArticlePubMedGoogle Scholar
- Tanser FC, Sharp B, le Sueur D: Potential effect of climate change on malaria transmission in Africa. Lancet. 2003, 362: 1792-1798. 10.1016/S0140-6736(03)14898-2.View ArticlePubMedGoogle Scholar
- Kleinschmidt I, Bagayoko M, Clarke GP, Craig M, le SD: A spatial statistical approach to malaria mapping. Int J Epidemiol. 2000, 29: 355-361. 10.1093/ije/29.2.355.View ArticlePubMedGoogle Scholar
- Snow RW, Peshu N, Forster D, Bomu G, Mitsanze E, Ngumbao E, Chisengwa R, Schellenberg JR, Hayes RJ, Newbold CI, Marsh K: Environmental and entomological risk factors for the development of clinical malaria among children on the Kenyan coast. Trans R Soc Trop Med Hyg. 1998, 92: 381-385. 10.1016/S0035-9203(98)91056-2.View ArticlePubMedGoogle Scholar
- Yang G, Zhou X, Malone JB, McCarroll JC, Wang T, Liu J, Gao Q, Zhang X, Hong Q, Sun L: GIS prediction model of malaria transmission in Jiangsu province. Zhonghua Yu Fang Yi Xue Za Zhi. 2002, 36: 103-105.PubMedGoogle Scholar
- Kitron U, Pener H, Costin C, Orshan L, Greenberg Z, Shalom U: Geographic information system in malaria surveillance: mosquito breeding and imported cases in Israel, 1992. Am J Trop Med Hyg. 1994, 50: 550-556.PubMedGoogle Scholar
- Sithiprasasna R, Linthicum KJ, Liu GJ, Jones JW, Singhasivanon P: Use of GIS-based spatial modeling approach to characterize the spatial patterns of malaria mosquito vector breeding habitats in northwestern Thailand. Southeast Asian J Trop Med Public Health. 2003, 34: 517-528.PubMedGoogle Scholar
- Bostoen K, Chalabi Z: Optimization of household survey sampling without sample frames. Int J Epidemiol. 2006, 35: 751-755. 10.1093/ije/dyl019.View ArticlePubMedGoogle Scholar
- Gemperli A, Vounatsou P, Kleinschmidt I, Bagayoko M, Lengeler C, Smith T: Spatial patterns of infant mortality in Mali: the effect of malaria endemicity. Am J Epidemiol. 2004, 159: 64-72. 10.1093/aje/kwh001.View ArticlePubMedGoogle Scholar
- Hu H, Singhasivanon P, Salazar NP, Thimasarn K, Li X, Wu Y, Yang H, Zhu D, Supavej S, Looarecsuwan S: Factors influencing malaria endemicity in Yunnan Province, PR China (analysis of spatial pattern by GIS). Geographical Information System. Southeast Asian J Trop Med Public Health. 1998, 29: 191-200.PubMedGoogle Scholar
- leSueur D, Sharp BL, Gouws E, Ngxongo S: Malaria in South Africa. S Afr Med J. 1996, 86: 936-939.Google Scholar
- Sharp BL, le SD: Malaria in South Africa – the past, the present and selected implications for the future. S Afr Med J. 1996, 86: 83-89.PubMedGoogle Scholar
- Richie TL, Saul A: Progress and challenges for malaria vaccines. Nature. 2002, 415: 694-701. 10.1038/415694a.View ArticlePubMedGoogle Scholar
- Durrheim DN, La Grange JJ, Govere J, Mngomezulu NM: Accuracy of a rapid immunochromatographic card test for Plasmodium falciparum in a malaria control programme in South Africa. Trans R Soc Trop Med Hyg. 1998, 92: 32-33. 10.1016/S0035-9203(98)90944-0.View ArticlePubMedGoogle Scholar
- Durrheim DN, Harris BN, Speare R, Billinghurst K: The use of hospital-based nurses for the surveillance of potential disease outbreaks. Bull World Health Organ. 2001, 79: 22-27.PubMed CentralPubMedGoogle Scholar
- Coleman M, Coleman M, Mabuza A, Kok G, Coetzee M, Durrheim D: Evaluation of an operational malaria outbreak identification and response system in Mpumalanga Province, South Africa. Malar J. 2008, 7: 69-10.1186/1475-2875-7-69.PubMed CentralView ArticlePubMedGoogle Scholar
- Govere J, Durrheim DN, Baker L, Hunt R, Coetzee M: Efficacy of three insect repellents against the malaria vector Anopheles arabiensis. Med Vet Entomol. 2000, 14: 441-444. 10.1046/j.1365-2915.2000.00261.x.View ArticlePubMedGoogle Scholar
- Kager PA: Malaria control: constraints and opportunities. Trop Med Int Health. 2002, 7: 1042-1046. 10.1046/j.1365-3156.2002.00981.x.View ArticlePubMedGoogle Scholar
- Kapp C: Hazard or help?. Lancet. 2004, 364: 1113-1114. 10.1016/S0140-6736(04)17122-5.View ArticlePubMedGoogle Scholar
- Booman M, Durrheim DN, La Grange K, Martin C, Mabuza AM, Zitha A, Mbokazi FM, Fraser C, Sharp BL: Using a geographical information system to plan a malaria control programme in South Africa. Bull World Health Organ. 2000, 78: 1438-1444.PubMed CentralPubMedGoogle Scholar
- Kulldorff M, Athas WF, Feurer EJ, Miller BA, Key CR: Evaluating cluster alarms: a space-time scan statistic and brain cancer in Los Alamos, New Mexico. Am J Public Health. 1998, 88: 1377-1380. 10.2105/AJPH.88.9.1377.PubMed CentralView ArticlePubMedGoogle Scholar
- Kulldorff M, Heffernan R, Hartman J, Assuncao R, Mostashari F: A space-time permutation scan statistic for disease outbreak detection. PLoS Med. 2005, 2: e59-10.1371/journal.pmed.0020059.PubMed CentralView ArticlePubMedGoogle Scholar
- Boscoe FP, McLaughlin C, Schymura MJ, Kielb CL: Visualization of the spatial scan statistic using nested circles. Health Place. 2003, 9: 273-277. 10.1016/S1353-8292(02)00060-6.View ArticlePubMedGoogle Scholar
- kulldorf m, Feuer EJ, Miller BA, Freeman LS: Breast Cancer in northeastern United States: a geographical analysis. Am J Epidemiol. 1997, 146: 161-170.View ArticleGoogle Scholar
- Nkhoma ET, Hsu CE, Hunt VI, Harris AM: Detecting spatiotemporal clusters of accidental poisoning mortality among texas counties, United States, 1980–2001. Int J Health Geogr. 2004, 3: 25-10.1186/1476-072X-3-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Odoi A, Martin SW, Michel P, Middleton D, Holt J, Wilson J: Investigation of clusters of giardiasis using GIS and a spatial scan statistic. Int J Health Geogr. 2004, 3: 11-10.1186/1476-072X-3-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Govere J, Durrheim D, La GK, Mabuza A, Booman M: Community knowledge and perceptions about malaria and practices influencing malaria control in Mpumalanga Province, South Africa. S Afr Med J. 2000, 90: 611-616.PubMedGoogle Scholar
- Morrison AC, Getis A, Santiago M, Rigau-Perez JG, Reiter P: Exploratory space-time analysis of reported dengue cases during an outbreak in Florida, Puerto Rico, 1991–1992. Am J Trop Med Hyg. 1998, 58: 287-298.PubMedGoogle Scholar
- Chadee DD, Kitron U: Spatial and temporal patterns of imported malaria cases and local transmission in Trinidad. Am J Trop Med Hyg. 1999, 61: 513-517.PubMedGoogle Scholar
- Gaudart J, Poudiougou B, Ranque S, Doumbo O: Oblique decision trees for spatial pattern detection: optimal algorithm and application to malaria risk. BMC Med Res Methodol. 2005, 5: 22-10.1186/1471-2288-5-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Hjalmars U, Kulldorff M, Gustafsson G, Nagarwalla N: Childhood leukaemia in Sweden: using GIS and a spatial scan statistic for cluster detection. Stat Med. 1996, 15: 707-715. 10.1002/(SICI)1097-0258(19960415)15:7/9<707::AID-SIM242>3.0.CO;2-4.View ArticlePubMedGoogle Scholar
- Kulldorf M: A spatial scan statistic. Communications in Statistics: Theory and Methods. 1997, 26: 1481-1496. 10.1080/03610929708831995.View ArticleGoogle Scholar
- Shaheen P, Clart JD, Williams D: Small area analysis: a review and analysis of North American literature. J Health Polit Policy Law. 1987, 12: 741-809. 10.1215/03616878-12-4-741.View ArticleGoogle Scholar
- Aamodt G, Samuelsen SO, Skrondal A: A simulation study of three methods for detecting disease clusters. Int J Health Geogr. 2006, 5: 15-10.1186/1476-072X-5-15.PubMed CentralView ArticlePubMedGoogle Scholar
- Patil GP, Taille C: Geographic and network surveillance via scan statistics for critical area detection. Stat Sci. 2003, 18: 457-465. 10.1214/ss/1081443229.View ArticleGoogle Scholar
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.