Spatial patterns of incident malaria cases and their household contacts in a single clinic catchment area of Chongwe District, Zambia
© Pinchoff et al. 2015
Received: 8 March 2015
Accepted: 3 July 2015
Published: 7 August 2015
Reactive case detection (RACD) for malaria is a strategy that may be used to complement passive surveillance, as passive surveillance fails to identify infections that are asymptomatic or do not seek care. The spatial and seasonal patterns of incident (index) cases reported at a single clinic in Chongwe District were explored.
A RACD strategy was implemented from June 2012 to June 2013 in a single catchment area in Chongwe District. Incident (index) cases recorded at the clinic were followed up at their household, and all household contacts were tested for malaria using rapid diagnostic tests (RDTs). GPS coordinates were taken at each index household. Spatial analyses were conducted to assess characteristics related to clustering, cluster detection and spatial variation in risk of index houses. Effects of season (rainy versus dry), distance to the clinic and distance to the main road were considered as modifying factors. Lastly, logistic regression was used to identify factors associated with the proportion of household contacts testing RDT positive.
A total of 426 index households were enrolled, with 1,621 household contacts (45% RDT positive). Two space–time clusters were identified in the rainy season, with ten times and six times higher risk than expected. Significantly increased spatial clustering of index households was found in the rainy season as compared to the dry season (based on K-function methodology). However, no seasonal difference in mapped spatial intensity of index households was identified. Logistic regression analysis identified two main factors associated with a higher proportion of RDT positive household contacts. There was a 41% increased odds of RDT positive household contacts in households where the index case was under 5 years of age [OR = 1.41, 95% confidence intervals (1.15, 1.73)]. For every 500-m increase in distance from the road, there was a 5% increased odds of RDT positive household contacts [OR = 1.05 (1.02, 1.07)], controlling for season.
Areas of increased report of malaria persist after controlling for distance to the clinic and main road. Clinic-based interventions will miss asymptomatic, non-care seeking infections located farther from the road. RACD may identify additional infections missed at the clinic.
The World Health Organization estimates that there were 207 million cases of malaria globally in 2012 with 627,000 deaths . Transmission in Zambia has been declining; the 2012 malaria indicator survey found a malaria parasite prevalence of 14.9% nationally down from 16% in 2010 . Since the late 1990s, Zambia has scaled up malaria control activities, including the distribution of long-lasting, insecticide-treated nets (LLINs), indoor residual spraying, prompt diagnosis and treatment with artemisinin combination therapy (ACT) and intermittent preventive therapy in pregnant women . The surveillance system in Zambia, as in most moderate and high transmission settings, is based on passive case detection, in which symptomatic individuals seek care at the clinic or hospital. As transmission declines, it becomes critical to enhance existing surveillance methods, as many infections are asymptomatic or have low access to healthcare facilities and do not receive treatment. Reactive case detection (RACD) is one form of active case detection in which all individuals in close proximity to a passively detected ‘index’ case are tested and treated, including if they are asymptomatic [4–6].
RACD is based on the assumption that infections are clustered in space and time, within transmission ‘hotspots’ [5, 7, 8]. In Zambia, a study found that the prevalence of malaria infection was significantly higher within passively detected index case households as compared to randomly selected control households (8 vs 0.7%, respectively) . To maximize the impact and cost effectiveness of RACD, the strategy should be targeted based on geographic or demographic risk factors .
Malaria transmission is heterogeneous at varying geographical scales, and transmission can vary significantly by season. This heterogeneity is driven by a variety of ecological, biological and sociological factors. Passive surveillance may only test and treat the sub-set of infections that seek care at the facility. Describing the spatial temporal distribution of passively detected households may help to better understand care seeking and malaria transmission dynamics in this peri-urban setting.
The aims of this study were to describe characteristics of the spatial distribution of passively detected index households for malaria as identified at the clinic. Tools from the field of spatial statistics were used to assess clustering, cluster detection and spatial variation in risk for incident malaria cases. Seasonal effects were also considered. Lastly, the composition of index households was explored, to describe the rapid diagnostic test (RDT) positivity among household contacts, as related to the households’ distance from the clinic.
Study site and procedures
The study site is located in Chongwe District about 35 km from Lusaka, the capital city of Zambia. The area is considered rural to peri-urban with families living in homesteads and relying on subsistence farming. The clinic catchment area is estimated to be 14,000 individuals. A major paved road runs through the centre of the catchment area. The clinic is near the paved road about in the centre of the catchment area. A swampy area created by a natural hot spring and likely a mosquito-breeding site resides in close proximity to the clinic.
In Zambia, the single rainy season lasts from November through April, followed by a cool season from April until August, and a hot dry season through November; malaria transmission peaks during the rainy season . The primary vectors for malaria in Lusaka Province are Anopheles gambiae, Anopheles funestus and Anopheles arabiensis . Malaria transmission in the area has been declining; Lusaka Province recorded the lowest parasitaemia in the country in 2012 .
Cases were passively detected at the clinic between June 2012 and June 2013. The study procedures involved testing each febrile individual seeking diagnosis and treatment at the clinic with a RDT and/or microscopy (Giemsa thick smear) to confirm malaria infection. Those that tested positive were approached for enrolment. RDT positive patients were randomly selected for study participation. Those that were approached and gave informed consent were enrolled in the study and administered a questionnaire. If a selected patient declined to participate, then a randomly selected back-up was approached. Demographic characteristics and malaria symptoms and treatment history were collected and household GPS coordinates recorded. A household contact questionnaire and RDT was administered at the household. Data were obtained on the total number of other house members positive for malaria infection, of those that were present at the time of the survey. All participants who tested positive for malaria in the clinic or during household visits were treated according to national guidelines.
All structures in the catchment area resembling households were enumerated and digitized to provide latitude and longitude coordinates based on a high-resolution satellite image (QuickBird WorldView, DigitalGlobe Services). This methodology was piloted and applied in a study in Lusaka and found to be highly accurate . Environmental features were digitized from the satellite image and GPS coordinates taken in the field, for features such as the paved road, the swamp area and the clinic location. Data were imported into ArcGIS v10.2 and spatially integrated with data on index households.
A series of analysis were conducted to characterize the spatial distribution of the index cases of malaria infection using tools from the field of spatial statistics to assess cluster detection, spatial variation in risk and spatial clustering. The first analysis conducted was to detect spatial–temporal clusters, sub-areas within the clinic catchment area where the number of malaria index cases over an identified period of time is significantly more than what would have been expected given the distribution of all index cases. The spatial scale for this analysis was aggregated to 500-sq m grid cells overlaid onto the study area accumulating the total number of enumerated households (population) and total number of index cases within each grid cell. The temporal scale was aggregated to a week (7 days). Analysis was performed using the cluster detecting software SaTScan  specifying the Poisson model approach to accommodate the count of index cases per grid cell relative to the corresponding grid cell population. Identified space–time clusters were restricted to be no larger than 50% of the total study area population at risk (a maximum of 4-km radius) and comprise no more than 4 weeks of available data. Cluster detection analysis for grid cell index case counts was performed adjusting for each grid cell population as well as controlling for distance to the clinic and distance to the road (distances were summarized into quartiles and calculated using the grid cell centroid location).
The second analysis was conducted to jointly explore spatial clustering of malaria index cases, a different characteristic than cluster detection, and spatial variation in risk of malaria index cases . Spatial clustering of malaria index cases seeks to assess the degree to which cases tend to be near other cases, and not necessarily to the degree to which a cluster is formed, whereas spatial variation in risk explores where these cases are occurring to identify spatial trends in the risk of such an event occurring. Of particular interest was how clustering and spatial variation in risk of malaria index cases changed as a function of season (rainy or dry). The K-function, which estimates the expected number of other cases within a range of distances of an arbitrary case and plots a scaled version of these estimates as a function of distance, was used to assess spatial clustering [15, 16]. Spatial intensity, defined as the expected number of cases per unit area, was estimated in a non-parametric fashion using the kernel density approach and mapped to highlight spatial variation in the concentration of cases [15, 16]. Both K-function analysis and spatial intensity were performed using the R Statistical software spatstat .
The final analysis conducted was to identify factors associated with a higher proportion of household contacts that were RDT positive. Knowing the total number of household members (contacts) per index house and the number of other household members testing positive, logistic regression was used to model the proportion of household contacts positive for malaria. Variables explored in the regression include distance to the paved road, distance to the clinic, as well as age and gender of the index case.
This study received ethical approval from both the University of Zambia Research Ethics Committee and the University of Alabama at Birmingham Institutional Review Board.
Characteristics of households (enumerated) in the clinic catchment area and enrolled households (index) between 2012 and 2013
Index N (%)
Enumerated N (%)
Number of households
Distance to clinic (km) [median (min, max)]
9.22 (0.23, 27.54)
15 (0.04, 35.86)
Distance to swamp (km) [median (min, max)]
6.61 (0, 25.49)
12.25 (0, 32.60)
Distance to road (km) [median (min, max)]
4.94 (0.32, 24.94)
8.75 (0.01, 26.65)
Index is under 5 years of age
Index is female
Visit in rainy season
Number of total household contacts in index household [median (min, max)]
5 (0, 13)
Number of RDT+ household contacts in index household [median (min, max)]
2 (0, 10)
Logistic regression results identifying factors associated with higher proportion of RDT positive household residents
OR (95% CI)
OR (95% CI)
OR (95% CI)
Index is female
1.08 (0.89, 1.31)
1.15 (0.93, 1.41)
1.15 (0.93, 1.41)
Index under 5 years of age
1.42 (1.17, 1.73)
1.41 (1.15, 1.73)
1.41 (1.15, 1.73)
Distance to clinic (per 500 m)
1.03 (1.01, 1.05)
1.00 (0.98, 1.03)
Distance to road (per 500 m)
1.05 (1.03, 1.07)
1.05 (1.03, 1.07)
1.05 (1.02, 1.07)
1.03 (0.84, 1.25)
1.16 (0.95, 1.43)
1.17 (0.95, 1.45)
Significant spatial and seasonal variation of incident malaria infections that reported to Chinyunyu clinic was identified. The space–time clusters of index households, controlling for distance to the road and to the clinic, suggests that factors other than access to care are impacting the distribution of index malaria cases in this catchment area. Potentially, transmission near the swamp is occurring, causing spatial confounding, since the swamp is located near the clinic and along the road. The identified clusters also occurred during the rainy season, suggesting a high number of cases during this period but perhaps also highlighting that in the rainy season roads wash out, and when access is restricted cases that are reported at the clinic are only from nearby areas.
Incident malaria cases reported in the rainy season clustered significantly more than in the dry season. However, the spatial intensity maps identified no difference in the spatial intensities per season. This suggests that cases were occurring in the same general area, but as demonstrated with the K-function difference method the rainy season the pattern was more spatially compact than in the dry season. The higher level of spatial clustering of incident malaria during the rainy season is likely because houses further from the clinic are unable to travel to the clinic during this time. Future studies should include a measurement of care or treatment-seeking behaviour to apply to the area, to control for the probability of an individual in the area seeking care for malaria in a model designed to measure surveillance.
Index households located farther from the road had a significantly higher number of asymptomatic, RDT positive household contacts that were not identified during passive detection, controlling for season. Additionally, if the index case was under the age of five, the household had more RDT positive household contacts. Likely, infections in younger children experienced more severe symptoms and were brought to the clinic for care and enrolling as index cases. These findings suggest that among index households, distance to the clinic may impact reporting and treatment of household contacts. As distance from the clinic increases, asymptomatic infections are less likely to receive testing and treatment, potentially leading to an increased infectious duration . This may impact ongoing transmission in the area.
The major limitations of this study are that information was not collected at control households, and that only a random selection of clinical cases were enrolled as index cases; both limited analyses. Without control households, it was not possible to determine if the prevalence of asymptomatic infections was higher in passively detected households. However, this has been suggested in several studies [4, 18, 19]. The random selection for enrolment means that many infected persons were not enrolled, however, since cases were randomly selected it can be assumed that misclassification is unbiased. Additionally, the findings may have been spatially confounded. It is unclear if the significant clustering of cases is due to increased risk of transmission from the swamp, or if being closer to the road and clinic results in improved access. Another limitation is that potentially enumeration was incomplete, but with up to date and highly accurate satellite imagery this should not have been the case, as this methodology has previously been shown to be accurate .
Future research should examine the efficiency of RACD in treating asymptomatic infections missed by passive surveillance. In a setting such as this clinic, the incidence may be too high to recommend RACD of all incoming cases. Although our study lacked control cases and could not assess how much RACD increases detection of asymptomatic infections, previous studies have found that asymptomatic infections spatially cluster in or near index households, suggesting that RACD targeted to the index household may be useful in identifying and treating remaining infections . Our findings support the suggestion that passive surveillance alone will not be sufficient to eliminate malaria, and complementary systems are necessary. This study found that having an index case under 5 years of age increased the chances of identifying asymptomatic infections in the household. This may be because children under 5 are more likely to experience symptoms and be brought in care, identifying household locations that are at increased risk for malaria transmission. Future studies regarding RACD should account for the underlying transmission intensity and include control households in order to determine key remaining questions regarding RACD including: (1) at what transmission intensity is RACD most efficient and effective for identifying asymptomatic cases, (2) how does RACD compare to active case detection methods such as randomly visiting households, (3) how many additional infections are treated with RACD compared to just passive surveillance.
These findings suggest that there is significant spatial and seasonal variation of incident (index) cases related to proximity to the clinic and main road. Malaria cases identified from the clinic were reported from the same geographic area, regardless of the season. However, the clusters detected occurred in the rainy season, and spatial clustering was higher in the rainy season, suggesting that when smaller roads wash out, only households closer to the clinic are able to seek care for malaria. Index cases tended to be under 5 years of age and have asymptomatic household contacts, mainly older children that did not seek care. Evaluating the spatial and temporal distribution of households can determine the impact of clinic-based interventions and improve the coverage of treatment and preventive interventions.
JP assisted with development of analysis plan, carried out all statistical analyses and drafted the manuscript. GH conceived of and designed the study, participated in coordination of the study. BC assisted with data cleaning, participated in drafting of the manuscript. STR developed and implemented initial study protocol, monitoring and evaluation of the study. SH assisted with implementation of the study protocol, and oversaw data collection, entry and cleaning. BH participated in the design and implementation of the study. MH participated in the design and implementation of the study. FCC designed the statistical analysis plan, supervised the analyses, and participated in drafting of the manuscript. All authors read and approved the final manuscript.
This work was made possible by a supplemental grant award through the IeDea (International Epidemiologic Databases to Evaluate AIDS: Malaria Research Preparedness) Network. The authors would like to thank the Chongwe District Health Management Team, the staff at the clinic, and the community health workers who made data collection possible; CIDRZ staff, including Chusi Sikanyika, who assisted with study implementation and data cleaning. The authors would also like to thank Dr William Moss for his insightful comments and suggestions in developing the analysis plan. Lastly, the authors would like to thank the community and clinic for their participation and support of this research project.
Compliance with ethical guidelines
Competing interests The authors assert no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- WHO (2013) World Malaria Report. World Health Organization, GenevaGoogle Scholar
- MOH (2012) Zambia National Malaria Indicator Survey 2012. Ministry of Health Government of the Rebublic of Zambia, Lusaka, ZambiaGoogle Scholar
- Kamuliwo M, Chanda E, Haque U, Mwanza-Ingwe M, Sikaala C, Katebe-Sakala C et al (2013) The changing burden of malaria and association with vector control interventions in Zambia using district-level surveillance data, 2006–2011. Malar J 12:437PubMed CentralPubMedView ArticleGoogle Scholar
- Stresman G, Kamanga A, Moono P, Hamapumbu H, Mharakurwa S, Kobayashi T et al (2010) A method of active case detection to target reservoirs of asymptomatic malaria and gametocyte carriers in a rural area in Southern Province Zambia. Malar J 9:265–273PubMed CentralPubMedGoogle Scholar
- Sturrock HJW, Novotny JM, Kunene S, Dlamini S, Zulu Z, Cohen JM et al (2013) Reactive case detection for malaria elimination: real-life experience from an ongoing program in Swaziland. PLoS One 8:e63830PubMed CentralPubMedView ArticleGoogle Scholar
- Moonen B, Cohen J, Snow R, Slutsker L, Drakeley C, Smith D et al (2010) Tanner M. Operational strategies to achieve and maintain malaria elimination. Lancet, Targett G, p 376Google Scholar
- Bejon P, Williams T, Liljander A, Noor A, Wambua J, Ogada E et al (2010) Stable and unstable malaria hostpots in longitudinal cohort studies in Kenya. PLoS Med 7:e1000304PubMed CentralPubMedView ArticleGoogle Scholar
- Bousema T, Drakeley C, Gesase S, Hashim R, Magesa S, Mosha F et al (2010) Identification of hot spots of malaria transmission for targeted malaria control. J Infect Dis 201:1764–1774PubMedView ArticleGoogle Scholar
- Sturrock H, Hsiang M, Cohen J, Smith D, Greenhouse B, Bousema T et al (2013) Targeting asymptomatic malaria infections: active surveillance in control and elimination. PLoS Med 10:e1001467PubMed CentralPubMedView ArticleGoogle Scholar
- Kent R, Thuma P, Mharakurwa S, Norris D (2007) Seasonality, blood feeding behavior, and transmission of Plasmodium falciparum by Anopheles arabiensis after an extended drought in Southern Zambia. Am J Trop Med Hyg 76:267–274PubMed CentralPubMedGoogle Scholar
- Kasasa S, Asoala V, Gosonui L, Anto F, Adjuik M, Tindana C et al (2013) Spatio-temporal malaria transmission patterns in Navrongo demographic surveillance site, northern Ghana. Malar J 12:63PubMed CentralPubMedView ArticleGoogle Scholar
- MOH (2010) Zambia National Malaria Indicator Survey: 2010. Government of the Republic of Zambia, Ministry of Health, Central Statistics OfficeGoogle Scholar
- Lowther S, Curriero F, Shields T, Ahmed S, Monze M, Moss W (2009) Feasibility of satellite image-based sampling for a health survey among urban townships of Lusaka Zambia. Trop Med Int Health 14:7–78View ArticleGoogle Scholar
- Kulldorff M (1997) A spatial scan statistic. Commun Stat Theory Methods 26:1481–1496View ArticleGoogle Scholar
- Waller LA, Gotway CA (2004) Spatial data, in applied spatial statistics for public health data, Wiley, Hoboken, NJ, USAGoogle Scholar
- Diggle P (2003) Statistical analysis of spatial patterns, 2nd edn. Oxford University Press, New YorkGoogle Scholar
- Baddeley A, Turner R (2005) Spatstat: an R package for analyzing spatial point patterns. J Stat Softw 12:1–42Google Scholar
- Lawpoolsri S, Chavez I, Yimsamran S, Puangsa-Art S, Thanyavanich N, Maneeboonyang W et al (2010) The impact of human reservoir of malaria at a community-level on individual malaria occurrence in a low malaria transmission setting along the Thai-Myanmar border. Malar J 9:143PubMed CentralPubMedView ArticleGoogle Scholar
- Branch O, Casapia W, Gamboa D, Hernandez J, Alava F, Roncal N et al (2005) Clustered local transmission and asymptomatic Plasmodium falciparum and Plasmodium vivax malaria infections in a recently emerged, hypoendemic Peruvian Amazon community. Malar J 4:27PubMed CentralPubMedView ArticleGoogle Scholar