GIS-aided planning of insecticide spraying to control dengue transmission
© Chu et al.; licensee BioMed Central Ltd. 2013
Received: 22 August 2013
Accepted: 24 September 2013
Published: 25 September 2013
The purpose of this paper is to integrate a multi-objective integer programming formulation and geographic information system (GIS) into dynamically planning the insecticide spraying area for preventing the transmission of dengue fever.
The optimal spraying area to combat dengue infections is calculated by the multi-objective integer programming model using the dengue epidemic in 2007 in Tainan City of southern Taiwan and is compared with the areas actually sprayed by the local health department. The dynamic epidemic indicators (i.e. frequency, intensity and duration) that identify major temporal characteristics of the dynamic process of an epidemic are all incorporated into the model.
The results indicate that the model can design the spraying area effectively when the trade-off between the coverage of dengue epidemics risk and area compactness is considered.
The model provides an alternative way to obtain a cost-effective spraying area in controlling future dengue epidemics. The proposed model in this study will be beneficial for strategically allocating dengue control resources.
KeywordsInteger programming Dengue fever Insecticide spraying planning Epidemic risk GIS Multi-objective
Dengue fever is a high-incidence, vector-borne infectious disease, and an estimated 390 million (95% credible interval 284-528 million) dengue infections occur per year around the world, chiefly distributed in tropical countries . Although a vaccine for dengue prevention is under development, there is still no good pharmaceutical intervention to prevent infection . Therefore, public health efforts focus on environmental control and surveillance of the mosquitoes’ indexes, including the House Index (HI), Container Index (CI) and Breteau Index (BI) [3, 4]. During the initial stage of an epidemic, local health departments will mobilize insecticide spraying resources within peridomestic space to prevent dengue vector breeding. The high costs associated with extensive spraying acknowledge that many residents may be inconvenienced. Based on one experimental study in southern Taiwan, the effectiveness of insecticide spraying within a radius of 100 meters can achieve a 90% reduction in the number of Aedes aegypti outdoors . In traditional practice, the spraying area focused on the infected cases near large communities. However, this strategy is feasible only when there are few, sporadic infections in the initial stage of the epidemic. When numerous dengue infections occur in close temporal proximity, spraying areas will overlap and expand while limited human resources. It will be a big challenge for the frontline public health workers to control dengue transmission without a systematic control strategy. Therefore, we develop one mathematical model which considers multiple objectives with maximal coverage of dengue fever hotspots and a minimal boundary length of spraying area at the same time. In addition, the epidemic conditions including the disease’s frequency, intensity and duration are incorporated into the strategy behind the extent of spraying . Geographic information systems (GIS) have provided useful tools for managing and surveillance of dengue fever [7, 8]. This study integrates GIS and an integer programming formulation  to dynamically plan the insecticide spraying area for prevention of the transmission of dengue fever. An integer programming approach was developed to help us to identify the insecticide spraying area that maximizes coverage of epidemic risk spots and minimizes the length of the boundaries of sprayed areas. The epidemic conditions were built into our GIS planning model. The 2007 dengue epidemic in Tainan City, located in southern Taiwan, was utilized as a case study and to compare the difference between the GIS-aided strategy and the empirical spraying area.
Material and study area
The temporal distribution of dengue infections from Jun. 2007 to Dec. 2007
Number of dengue infections
Integrating a multi-objective integer programming formulation using GIS in this study is to plan the optimal insecticide spraying area so as to maximize coverage of areas with dengue fever epidemic risk, minimize the length of the boundaries of sprayed areas, and take into account total area constraints. The epidemic risk goals can be related to (1) how often the disease occurs (frequency, α), (2) how long an epidemic persists (duration, β) and (3) how significant cases occurring in consecutive weeks are during an epidemic (intensity, γ). These three indicators, however, represent the disease component of dengue risk and enable to identify high-risk areas (hotspots) for the occurrence of dengue fever .
Where S i is epidemic risk in combination of frequency (α i ), duration (β i ) and intensity (γ i ); L ij is the length of the border between cells i and j; z is the planning area of the dengue hotspot zone; in this study, the cell size is 0.5 km2 (the average size of villages in Tainan city is 0.79 km2) and z is set as 15 or 25 km2. I i is a binary number that represents whether cell i is a hotspot or not; P ij and N ij are pseudo-binary numbers that are related to the boundary determination. P ij = 1 if I i = 1 and I i =0, P ij = 0, otherwise; N ij =1 if I j = 1 and I i =0, N ij = 0, otherwise. n is the number of cells in the study area. C i is the set of cells adjacent to cell i. w1 and w2 are the weights of epidemic risk coverage and boundary length. wα,wβ and wγ are the weights of frequency, duration and intensity indicators. For details of the optimization model, please refer to the reference . This model was programmed by LINGO (Lindo Systems, Inc., Chicago, IL, USA).
Where TW i is the total number of months during the entire epidemic period. EW i is the total number of months in which one or more cases occurred during the entire epidemic period.
where TVi is the total number of epidemic waves during the entire epidemic period.
Where IR i is the incidence rate during the defined epidemic period.
Figure 1 shows three epidemic indicators which are frequency (Figure 1b), duration (Figure 1c) and intensity (Figure 1d) in the study area, respectively. The indicators demonstrate major temporal characteristics and patterns of the epidemic dynamic process. The distribution of the frequency indicator is more widely diffused than the distribution of the other two indicators (i.e. duration and intensity). The high frequency indicators are located in the North, East and Middle-West Districts. The spatial patterns of duration and intensity indicators are consistent and highly correlated (correlation coefficient = 0.97, p-value <0.01), and were done by SPSS 10.0 (SPSS, Inc., Chicago IL). The model determines the insecticide spraying areas for reaching two objectives such as maximum coverage of environmental epidemic risk and minimizing the boundary of the areas. The compactness of the area is beneficial for environmental management and control.
The combinations of the parameters in the case studies and the epidemic risk captured by the model
The parameters in each case
Average frequency indicator in planning area
Average duration indicator in planning area
Average intensity indicator in planning area
Average Epidemic risk
Boundary length (km)
In this study, we demonstrated a useful and practical method for planning an insecticide spraying area systematically. The flexible parameters, including the user-defined size of spraying area and the weight of the boundary minimization can be adjusted according to the amounts of resources and the experience of public health workers. In addition, a dengue epidemic will evolve over time. The utilization of the three epidemic indicators (i.e. frequency, duration and intensity) to identify the temporal and spatial clusters was demonstrated to be a suitable indicator in addition to the traditional case incidence data . As a result, by including those indicators, the effectiveness at vector control and planning the appropriate area for spraying is enhanced.
The disease frequency, intensity and duration maps, however, clearly show that the main hotspot and several additional areas (e.g. in Annan District or South District) were severely affected by dengue fever in 2007. These regions are not covered by any of the optimized spraying areas (cases 1-6) in the model. If the weighting of minimal boundary in the objective function is decreased, these regions can be increasingly covered by the scattered hotspots. Thus, the model is able to aid in planning spraying when there are several scattered or clustered hotspots of dengue fever cases. Many tropical developing countries have been suffering the threat of the dengue fever and this model might help them plan the spraying area efficiently and economically. The approach can be used in any environments when the data is available. Before applying our model into practical use, there are some principal data you need to have such as dengue surveillance data (i.e. onset dates and locations of the dengue cases) and base GIS layers (i.e. administrative boundaries and road maps). In the model, the major parameters are only the weights of objective functions and the size of planning area.
The dengue control guidelines from Taiwan’s CDC  mention the standard operating procedure when a local department of health receives notification of suspected dengue infection cases from the hospitals or clinics. The public health workers need to do source reduction within 48 hours based on possible infection sites including living areas, work places and other possible infected places. Chemical space spraying will be applied when breeding sites are found in the outdoors. The strategy of space spraying is currently based on each individual case instead of a systematic preventive approach. With GIS application and our model, the compact or comprehensive spraying area would be planned based on the user’s experience during different stages of the epidemic. In addition, the temporal distribution of dengue infections has historically varied greatly in different months. Our model could be scaled down into daily or weekly units to immediately adjust to the variability of the epidemic.
Advocacy, social mobilization and legislation.
Collaboration within the health sector and with other sectors.
Evidence-based decision making.
The model proposed here mainly focuses on the third and fourth ones, especially aiding the spatial planning of the chemical vector control. Furthermore, our model also could be extended to routine vector surveillance when the epidemic indicators are replaced with entomological indicators such as the Breteau Index and House Index.
In the literature, the timing of using insecticide , the usage of the insecticide on the cover of water containers , and insecticide spraying on the breeding sites  were all shown to reduce dengue transmission. However, there is little known about how to arrange the spraying area when dengue cases occur. Hence, the core of the study used multi-objective planning to design an effective and systematic method for spraying area selection. Typically, multiple criteria simultaneously govern the hotspot selection process, such as the planning area, boundary length, and the coverage of epidemic risk. Reducing the boundary length of a hotspot area relative to its area is another critical concern in epidemic planning, owing to its influence on the economic costs of spraying.
The proposed model in this study would allow the temporal variability of dengue based on incorporating a time-varying, predefined control area. An epidemic of dengue infections occurred with a distribution over both time and space. This study incorporated three critical indicators, namely frequency, duration and intensity of epidemic cases, into the multi-objective planning of insecticide spraying. Results showed that our planned spraying area was within the real, historical spraying area of the 2007 Tainan epidemic, but the planned area was more concentrated. The model provided an alternative way to obtain a cost-effective spraying area in controlling future dengue epidemics. In the future, enhancements to the model would incorporate exposure risk and entomological indicators, which would capture the dengue risk more effectively. The quasi-experiments on comparing effectiveness between traditional and GIS-aided space spraying were needed.
This research was supported by a grant from the National Science Council, R.O.C. (NSC-102-2314-B-001-001). We would also like to express our sincere gratitude to Mr. Kent M. Suárez for his English editing.
- Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, Drake JM, Brownstein JS, Hoen AG, Sankoh O, et al: The global distribution and burden of dengue. Nature. 2013, 496 (7446): 504-507. 10.1038/nature12060.PubMedPubMed CentralView ArticleGoogle Scholar
- Konishi E: Issues related to recent dengue vaccine development. Trop Med Health. 2011, 39 (4 Suppl): 63-71.PubMedPubMed CentralView ArticleGoogle Scholar
- Seidahmed OM, Siam HA, Soghaier MA, Abubakr M, Osman HA, Abd Elrhman LS, Elmagbol B, Velayudhan R: Dengue vector control and surveillance during a major outbreak in a coastal Red Sea area in Sudan. East Mediterr Health J. 2012, 18 (12): 1217-1224.PubMedGoogle Scholar
- Sarfraz MS, Tripathi NK, Tipdecho T, Thongbu T, Kerdthong P, Souris M: Analyzing the spatio-temporal relationship between dengue vector larval density and land-use using factor analysis and spatial ring mapping. BMC Public Health. 2012, 12: 853-10.1186/1471-2458-12-853.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang C-Y, Lin C, Ge Y-F, Chien H-F, Lin C-S, Wang C-S, Jian S-W, Chen Y-H, Wang Y-Z, Teng H-J: Survey of dengue fever vector density before and after insecticide spraying. Taiwan Epidemiol Bull. 2011, 27 (19): 323-330.Google Scholar
- Wen TH, Lin NH, Lin CH, King CC, Su MD: Spatial mapping of temporal risk characteristics to improve environmental health risk identification: a case study of a dengue epidemic in Taiwan. Sci Total Environ. 2006, 367 (2–3): 631-640.PubMedView ArticleGoogle Scholar
- Aziz S, Ngui R, Lim YAL, Sholehah I, Farhana JN, Azizan AS, Yusoff WWS: Spatial pattern of 2009 dengue distribution in Kuala Lumpur using GIS application. Trop Biomed. 2012, 29 (1): 113-120.PubMedGoogle Scholar
- Chang AY, Parrales ME, Jimenez J, Sobieszczyk ME, Hammer SM, Copenhaver DJ, Kulkarni RP: Combining Google Earth and GIS mapping technologies in a dengue surveillance system for developing countries. Int J Health Geogr. 2009, 8: 49-10.1186/1476-072X-8-49. doi: 10.1186/1476-072X-8-49PubMedPubMed CentralView ArticleGoogle Scholar
- Wright J, Revelle C, Cohon J: A multiobjective integer programming-model for the land acquisition problem. Reg Sci Urban Econ. 1983, 13 (1): 31-53. 10.1016/0166-0462(83)90004-2.View ArticleGoogle Scholar
- Armstrong MP, Rushton G, Zimmerman DL: Geographically masking health data to preserve confidentiality. Stat Med. 1999, 18 (5): 497-525. 10.1002/(SICI)1097-0258(19990315)18:5<497::AID-SIM45>3.0.CO;2-#.PubMedView ArticleGoogle Scholar
- Nazri Che D, Zulkiflee Abd L, Abu Hassan A, Rodziah I, Pradhan B: Manifestation of GIS tools for spatial pattern distribution analysis of dengue fever epidemic in the City of Subang Jaya, Malaysia. EnvironmentAsia. 2012, 5 (2): 82-92.Google Scholar
- Taiwan CDC: Guidelines for Dengue Control. 2010, Taipei, Taiwan: CDC, Department of HealthGoogle Scholar
- Esu E, Lenhart A, Smith L, Horstick O: Effectiveness of peridomestic space spraying with insecticide on dengue transmission; systematic review. Trop Med Int Health. 2010, 15 (5): 619-631.PubMedGoogle Scholar
- Teng HJ, Chen TJ, Tsai SF, Lin CP, Chiou HY, Lin MC, Yang SY, Lee YW, Kang CC, Hsu HC, et al: Emergency vector control in a DENV-2 outbreak in 2002 in Pingtung City, Pingtung County, Taiwan. Jpn J Infect Dis. 2007, 60 (5): 271-279.PubMedGoogle Scholar
- World Health Organization: Handbook for integrated vector management. 2012, Geneva, Switzerland: WHO PressGoogle Scholar
- Oki M, Sunahara T, Hashizume M, Yamamoto T: Optimal timing of insecticide fogging to minimize dengue cases: modeling dengue transmission among various seasonalities and transmission intensities. PLoS Negl Trop Dis. 2011, 5 (10): e1367-10.1371/journal.pntd.0001367.PubMedPubMed CentralView ArticleGoogle Scholar
- Kroeger A, Lenhart A, Ochoa M, Villegas E, Levy M, Alexander N, McCall PJ: Effective control of dengue vectors with curtains and water container covers treated with insecticide in Mexico and Venezuela: cluster randomised trials. BMJ. 2006, 332 (7552): 1247-1252. 10.1136/bmj.332.7552.1247.PubMedPubMed CentralView ArticleGoogle Scholar
- Rizzo N, Gramajo R, Escobar MC, Arana B, Kroeger A, Manrique-Saide P, Petzold M: Dengue vector management using insecticide treated materials and targeted interventions on productive breeding-sites in Guatemala. BMC Public Health. 2012, 12: 931-10.1186/1471-2458-12-931. doi: 10.1186/1471-2458-12-931PubMedPubMed CentralView 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.