This article has Open Peer Review reports available.
A spatial simulation model for dengue virus infection in urban areas
© Karl et al.; licensee BioMed Central Ltd. 2014
Received: 5 June 2014
Accepted: 13 August 2014
Published: 20 August 2014
The World Health Organization estimates that the global number of dengue infections range between 80–100 million per year, with some studies estimating approximately three times higher numbers. Furthermore, the geographic range of dengue virus transmission is extending with the disease now occurring more frequently in areas such as southern Europe. Ae. aegypti, one of the most prominent dengue vectors, is endemic to the far north-east of Australia and the city of Cairns frequently experiences dengue outbreaks which sometimes lead to large epidemics.
A spatially-explicit, individual-based mathematical model that accounts for the spread of dengue infection as a result of human movement and mosquito dispersion is presented. The model closely couples the four key sub-models necessary for representing the overall dynamics of the physical system, namely those describing mosquito population dynamics, human movement, virus transmission and vector control. Important features are the use of high quality outbreak data and mosquito trapping data for calibration and validation and a strategy to derive local mosquito abundance based on vegetation coverage and census data.
The model has been calibrated using detailed 2003 dengue outbreak data from Cairns, together with census and mosquito trapping data, and is shown to realistically reproduce a further dengue outbreak. The simulation results replicating the 2008/2009 Cairns epidemic support several hypotheses (formulated previously) aimed at explaining the large-scale epidemic which occurred in 2008/2009; specifically, while warmer weather and increased human movement had only a small effect on the spread of the virus, a shorter virus strain-specific extrinsic incubation time can explain the observed explosive outbreak of 2008/2009.
The proof-of-concept simulation model described in this study has potential as a tool for understanding factors contributing to dengue spread as well as planning and optimizing dengue control, including reducing the Ae. aegypti vector population and for estimating the effectiveness and cost-effectiveness of future vaccination programmes. This model could also be applied to other vector borne viral diseases such as chikungunya, also spread by Ae. aegypti and, by re-parameterisation of the vector sub-model, to dengue and chikungunya viruses spread by Aedes albopictus.
Dengue is a leading cause of morbidity in tropical environments around the world  with a small proportion of infections resulting in possibly fatal dengue hemorrhagic fever (DHF) . While previous estimates of the global dengue burden were in the range of 80–100 million human infections per year, recent studies suggest a considerably higher number of infections (390 million per year) . There are also clear indications that the global range of dengue transmission is extending (e.g. in Southern Europe), with higher case numbers occurring and Aedes mosquitoes colonizing new habitats [4, 5].
Dengue epidemics can be especially severe in urban areas, where human population density is high . The vector mosquito Ae. aegypti have adapted to life in densely populated areas, using standing water as breeding sites with females feeding predominantly on humans, and are responsible for dengue transmission in urban environments .
Controlling dengue outbreaks is resource intensive and large outbreaks may overwhelm even well established and efficient control mechanisms . It is therefore important for the planning of adequate control interventions that the dynamics, frequency and scale of expected outbreaks can be predicted and simulated using suitable computational models . Furthermore, such simulation models can be used to estimate the effectiveness of alternative mitigating intervention strategies that are difficult to determine in the field .
Ae. aegypti is endemic to urban areas of northeast Queensland, Australia. Dengue viruses are often introduced to this region through viremic travellers from dengue endemic regions. This frequently causes dengue outbreaks of variable severity and the city of Cairns is often the focus of these outbreaks [7, 10].
Long-term dengue research and surveillance, including mosquito trapping studies in Cairns, has resulted in the collation of significant datasets describing dengue outbreaks in this area [7, 10–13]. These datasets were made available to the authors from the Tropical Population Health Unit (TPHU) of Queensland Health. These data have allowed for the informed development of a detailed dengue simulation model which may be used to estimate and evaluate the characteristics and determinants of dengue outbreaks and to test a range of interventions such as interior (within-house) residual spraying (IRS), the larvicidal treatment or destruction of Aedes breeding sites, potential dengue vaccination programmes and mosquito population manipulation using Wolbachia infected mosquitoes .
The aim of the present study was to develop a proof-of-concept, spatially-explicit simulation model of dengue transmission that incorporates individual humans living in the city of Cairns as well as individual mosquitoes, and uses realistic, heterogeneous human and mosquito population structures and movement patterns. A spatial modelling approach is required since many of the factors that are crucial for dengue spread are not homogeneously distributed (e.g., human population density and mosquito abundance).
The model developed in the present study builds upon previous studies which have modelled mosquito population dynamics and/or dengue transmission. The model presented here adopts many of the features of the previously developed CIMSIM/DENSIM and Skeeterbuster models, which are highly developed, complex simulation models that are focused on capturing mosquito population dynamics with great detail using a complex system of customizable mosquito breeding containers and weather data as input [15–25].
The motivation of the present study was i) to build a spatial mosquito population dynamics sub-model which is weather mediated and results in physically realistic temporal and spatial mosquito abundance patterns; ii) to build human population, dengue transmission and outbreak management sub-models with internal feedback (that is, dengue control has an effect on the mosquito population) and couple this to the mosquito population sub-model and iii) to calibrate and validate the model using detailed data available for several dengue outbreaks in the city of Cairns, Queensland.
The resulting model structure utilises new modelling techniques to capture physical system properties such as mosquito flight, rainfall dependence of the mosquito population, human movement patterns and dengue control. Here, the model is demonstrated by reproducing different outbreak scenarios recorded in Cairns. This is one of very few spatial dengue transmission model that include all of the following features i) human movement and mosquito flight, ii) geospatial estimation of mosquito breeding site abundance in an urban area, iii) internal feedback of individual-based mosquito control during a dengue outbreak on the mosquito population dynamics and iv) the use of detailed, spatial outbreak datasets to calibrate and, separately, validate the model.
The developed model is a proof-of-concept demonstration that the inherently complex phenomena observed in actual dengue outbreaks can be captured by a spatial simulation model that incorporates interlinked and interacting sub-models for each of the phenomena that determine the outcome of dengue epidemics. These are: mosquito population dynamics and movement; human population movement; transmission of dengue virus between mosquitos and humans; and vector control measures.
While the following methods section illustrates the overall approach to building, calibrating and validating the model, a detailed description of the various model components may be found in Additional file 1.
Figure 2 indicates how four types of input data were used, and these data differed between the 2003 and the 2008/2009 outbreaks: i) weather data, ii) dengue strain specific extrinsic incubation period data, iii) index case location data and iv) data describing the dengue control operations undertaken in the two outbreaks.
The overall methods applied in this study are presented in the following. A detailed desciption of all model components is given in Additional file 1.
Human population model
The human population part of the model is based on previous models designed by Milne and colleagues but with important differences e.g., in human movement behaviour as detailed below and in the Additional file 1 (pages S2-S7) [27–30].
The model of Cairns consists of approximately 52,000 cells (Additional file 1: Figures S2 and S4) and 56,000 human individuals. To assign a human population to the cells, two sets of data were used: i) the most recent population census from the Australian Bureau of Statistics (accessible at http://www.abs.gov.au, illustrated in Additional file 1: Figures S2 and S3) and ii) a geo-referenced digital cadastral dataset from Cairns provided by the Department of Natural Resources and Mines of the Queensland Government (accessible at http://www.qld.gov.au, illustrated in Additional file 1: Figures S2 - S4).
While the census data provided information on human population age distribution (illustrated in Additional file 1: Figure S1 for all of Cairns), household sizes and the number of people per area at a high spatial resolution (Additional file 1: Figure S2 Panels A and B), the cadastral dataset specified residential, commercial, industrial, educational and parkland associated properties (Additional file 1: Figure S2 Panels C and D). This made it possible to allocate households of appropriate sizes and with appropriate age distributions to each residential property. Thus, each human is assigned a home cell. Apart from their home cells, and depending on their age, humans may also frequent other cells. Properties classified as commercial and industrial in the cadastral dataset served as workplaces and shopping centres that are frequented by the human population. Data on schools (student numbers, class sizes) were sourced from the Queensland school register (accessible at http://www.education.qld.gov.au) for schools in the modelled area. Appropriate virtual classes of school children, including their teachers were built and associated with the educational establishments listed in the cadastral dataset. In addition, parkland areas are also frequented by humans for recreational purposes. A single cell in the model can contain more than one property and also properties of different types (e.g. Additional file 1: Figure S4 Panel B). We assume no further boundaries between all humans allocated to a single cell (e.g., mosquito access to humans allocated to two separate households which are located in the same cell will be the same).
In the model, time progresses in a step-wise (discrete) manner and in intervals of 6 hours, i.e. infection states and locations of humans and mosquitoes may change every 6 hours. Each 6 hour period corresponds with a specific part of the day: i) morning (3 am to 9 am), ii) daytime (9 am to 3 pm), iii) evening (3 pm to 9 pm) and iv) night (9 pm to 3 am). Humans are assumed to move between cells during the daytime and evening periods, and are assumed to be in their home cells during morning and night periods.
It has been shown that human movement is an important factor facilitating the spread of dengue . Detailed information on human movement is very sparse; we used the findings of a large survey of interpersonal human contact to derive an approximate model of daily movement of individuals, which apportions each person’s time between different locations . This model incorporates two types of human movement.
The first type is directional and always occurs between two specific cells for a given individual. These cells are i) the individual’s home cell and ii) the individual’s work or school cell. On any given weekday of the modelled period, individuals will go from their home cells (morning) to their work or school cells (daytime) and back to their home cells (evening or night), following the approach adopted in Milne et al. .
The second type of human movement used in the present study is semi-random and is incorporated in the model to account for the less predictable movement which humans may undertake within the modelled area e.g., to visit friends, go shopping, or visit recreation areas. This semi-random movement can occur on weekday evenings and during the day and evening during weekends. Based on data from a population survey by Mossong et al. it was estimated that the frequency of random human visits of this type to other cells is approximately 4–5 times per week per individual . A schematic representation of the human movement components of the present model is shown in Additional file 1: Figure S5.
Human movement in urban areas has been shown to be semi-random and to be ranked by distance i.e., the frequency of short trips is higher than the frequency of trips to destinations that are further away . In the model, this distance-dependent behaviour is modelled with a gamma distribution of the semi-random movement distances (Additional file 1: Figure S6) .
Mosquito population dynamics model
The mosquito population dynamics model used in the present study is based on previously developed models by others [17, 23, 25, 33]. The core mosquito population dynamics model (only female adult mosquitoes) is shown in Additional file 1: Figure S7. It includes egg (E), larvae (L), pupae (P) and two adult mosquito populations (A1 and A2) as determined by the length of the first (A1) versus the lengths of the remaining (A2) gonotrophic cycles, with the first gonotrophic cycle being significantly longer than the remaining ones .
The mosquito population sub-model is similar to that described in previous studies but it also differs from these in a number of ways, such as details in the implementation of larvae density dependent parameters and mosquito flight [17, 23, 24, 33]. Specific details of the mosquito population dynamics model are given in the Additional file 1 (pages S8-S23 with Figures S7 to S23).
Weather sub-model and mosquito habitat heterogeneity
The present model used a novel approach to model weather driven mosquito population dynamics that differs significantly from those used in previous studies [17, 23, 24, 33]. The weather sub-model uses temperature, rainfall and evaporation data from Cairns (available from the Australian Bureau of Meteorology http://www.bom.gov.au) as input and these weather data drive the mosquito population dynamics model component.
Temperature dependent parameters of the mosquito population dynamics model (Additional file 1: Figure S7) are determined similarly to the previous studies referenced above (described in detail in the Additional file 1 e.g., Table S1 and Figure S21). The effect of temperature on the mosquito population is assumed to be a global effect, which occurs to the same extent in all model cells. As in the previous studies, we assume that mosquito larvae develop in a density dependent manner, that is, their growth rate declines to a minimum as their density increases towards a maximum value (for further information see Additional file 1 pages S13 to S14 and Figures S12, S13, Table S2 and Equations S6 and S7).
To permit spatial heterogeneity within the mosquito population we assume that certain cells are better suited to support mosquito development than others. We define the minimum and maximum capacity of a model cell to sustain mosquito larvae (L min and L max ). The minimum capacity of a cell to sustain mosquito larvae is maintained throughout the year, irrespective of rainfall and accounts for artificially watered containers present in the cell and for those containers that are shielded from evaporation. Following rainfall (in the model this is the time when the amount of rain exceeds the amount of evaporation), a cell’s capacity to sustain mosquito larvae increases until it eventually reaches a maximum value (L max ). We assume that L max marks the point at which containers overflow and a cell is saturated with water, so that its capacity to sustain mosquito larvae does not increase further. We allow for a fraction of the immature mosquito stages to be washed out from overflowing containers (Additional file 1: Figure S11 and Equation S5). The capacity of a cell to sustain mosquito larvae thus fluctuates between the L min and L max values with rainfall and evaporation and cells can hold more or less mosquito larvae depending on rainfall, which in tropical Cairns varies between wet summers and dry winters. For detailed explanations of the weather sub-model refer to the Additional file 1 pages S9 to S15).
In contrast to previous dengue simulation models that use complicated container characteristics which need to be determined by field studies [15, 16, 18, 33], the approach adopted here is aimed at the use of geographic information system (GIS) data to estimate spatial distribution of mosquito habitats. While this may be less accurate than manual classification of each container in the modelled area, it has the benefit of allowing for transfer model to other geographic locations more readily.
It is assumed that two cell-specific geographical features are positively correlated with a cell’s capacity to sustain mosquito larvae, namely i) the degree of vegetation cover in a cell and ii) the number of dwellings per cell. This assumption is based on previous work [34, 35] and the positive association found between mosquito trapping data from Cairns and these two characteristics in the present study. We assume that the number of dwellings per cell is more important than vegetation, since it provides i) humans, who are the source for Ae. aegypti blood feeding and ii) human-made breeding habitats such as rainwater tanks, underground sumps, flower pots etc. . Previous studies have shown significant correlations between vegetation cover and the abundance of Ae. aegypti breeding sites in urban settings , in the model the presence of vegetation cover adds an additional benefit to a cell’s capacity to sustain mosquito larvae, as it may prevent increased solar exposure and enable lower local rates of evaporation. Vegetation features may also provide additional breeding sites such as fallen palm tree fronds which may collect water.
A minimum breeding site index B min is present independently from the coverage with dwellings and vegetation. Additional file 1: Figure S15 shows the distribution of the breeding site abundance indices (B) based on Equation 1 (Equation S8). B is normalized to between 0 and 1 as it is only a relative value that stands for the suitability of a cell for mosquito reproduction. We believe, that the general distribution of expected mosquito breeding sites is realistically represented by the resulting pattern (Additional file 1: Figure S15). For example, industrial areas have low breeding site indices (blue). Houses with significant amounts of surrounding vegetation, as in Parramatta Park (a suburb with a history of dengue outbreaks), have high breeding site indices (red).
As the process of assigning a mosquito breeding site abundance index to each cell does not rely on fitting-to-data, validation of this methodology was conducted using available trapping data from Cairns for the years 2006–2013, kindly provided by Scott Ritchie (James Cook University) and Peter Cook (Monash University). We relied only on the data collected using commercially available BG traps since these are considered to be reliable . Note that the available trapping data were limited to a small proportion of the modelled area (see Additional file 1: Figure S18 for trap locations). For further information on the validation process using mosquito trapping data refer to the Additional file 1 (Pages S15-S18, Figures S15-S18).
Calibration of the mosquito population dynamics model
The mosquito trapping data consisted of mean number of female Ae. aegypti captured by Biogents Sentinel traps (Figure 4; Additional file 1: Figure S20) from 09/10/2006 and 16/08/2008 . During this time, no major dengue outbreaks were recorded in the modelled area, hence no significant mosquito control activities impacted on the mosquito population. A detailed description of the calibration process is given in the Additional file 1.
The mosquito trapping data provides a measure of relative mosquito abundance, therefore, in Figure 5 (Additional file 1: Figure S22) mosquito trapping data and the calibrated mosquito abundance curve are presented without a scale, since the absolute number of mosquitoes is unknown.
Variation of L max , which is the maximum cell capacity to sustain mosquito larvae, will not change the relative shape of the calibrated curve presented in Figure 5 (Additional file 1: Figure S22). However L max can be used to scale the absolute mosquito population based on the assumption that all modelled cells should follow a similar mosquito density pattern, though some are more suitable to sustain mosquito reproduction than others. Furthermore, it is assumed that absolute mosquito density should reflect an average female mosquito number of 8–24 Ae.aegypti per property, as estimated from mosquito trapping and computer simulation studies in Cairns [39, 40]. Since the modelled area contains approximately 28,000 properties, L max was chosen to obtain an average mosquito number of around 2.24 × 105 to 6.72 × 105 mosquitoes in the modelled area (about 8–24 mosquitoes per dwelling). A resulting graph showing absolute mosquito abundance is shown in Additional file 1: Figure S23.
The calibrated version of the mosquito population sub-model allows for the calculation of temperature and rainfall dependent mosquito profiles for any year for which weather data are available. In years where there are large dengue outbreaks, the associated vector control programme will impact on the mosquito population beyond what may be calibrated with the pure mosquito population dynamics model presented here. Hence, only the coupled mosquito population dynamics, dengue transmission and dengue control model can be used to estimate mosquito population curves in years with significant dengue outbreaks.
We use a standard SEIR transmission model as done in other studies, including those modelling dengue transmission [16, 25–29]. There are two components to the transmission model i) the human one and ii) the mosquito one. The dengue transmission sub-model is described in detail in the Additional file 1 (pages S24-S26).
At any given time individual humans will be in one of four states, namely S susceptible, E exposed/infected, I infectious and R recovered/immune. The model allows for the possibility that individuals become infectious but the disease is not diagnosed or progresses asymptomatically (we do not further discriminate between the causes of non-detection; it may be due to misdiagnosis or the absence of symptoms). This results in individuals being in one of two infectious states, I S symptomatic or I A asymptomatic. A schematic representation of the human population model is presented in Additional file 1: Figure S24 (and simplified in Figure 3).
Both A1 and A2 stage mosquitoes (Additional file 1: Figure S7) bite and transmit dengue virus. The infection process in the mosquitoes does not include a recovered/immune state as it is assumed that mosquitoes remain infected until they die. A schematic representation of virus progression in the mosquito is shown in Additional file 1: Figure S25 (and simplified in Figure 3).
Note that none of the transitions between human or mosquito states described above are simple first order kinetics. Probabilities of infection from humans to mosquitoes P H→M and from mosquitoes to humans P M→H will depend on the number of humans and mosquitoes co-located in a cell, resulting in these transitions being stochastic processes based on sampling from binomial distributions. Intrinsic incubation times and human recovery rates are gamma distributed which makes early progression less likely . The extrinsic incubation periods in mosquitoes (Additional file 1: Figure S26) are sampled from temperature dependent log-normal distributions based on studies by Chan et al. .
It is assumed, that, following report of a case (i.e. that an infected, symptomatic case has occurred in a cell) IRS and ovitrap distribution would commence with an average lag time of 7 days .
IRS and larvicide treatment in the immediate vicinity of a case is known to be highly effective, killing an estimated 90% of all adult and immature mosquito stages present in the case cell and the 8 surrounding cells per day . In the model it is assumed that the effect remains active for 6 weeks and affects all mosquitoes entering the cell during that time.
Lethal ovitraps have also been shown to be effective in controlling Ae. aegypti.. Their deployment is faster than the more labour intensive IRS and larvicide treatment . It is assumed that ovitraps are deployed at a certain rate per day in the cells within a 200 m radius of a case (depending on the amount of resources committed to fight a dengue outbreak, this rate may be high or low). Ovitraps kill approximately 20-30% of adult mosquitoes per cell per day . The overall effect of lethal ovitraps is reduced due to the emergence of new mosquitoes, as we assume that lethal ovitraps only target the adult stage.
It is assumed that there are limitations on coverage and number of IRS premises treated and ovitraps deployed each day, depending on the size and number of available dengue response teams (i.e. human resources). If the required number of cells to be treated exceeds the maximum number treatable per day, the remaining cells that are left untreated are carried over to the next day.
Parameters used in the vector control part of the model
IRS + Larviciding
QH Dengue Management Plan
90% per day
QH Dengue Management Plan, 
Maximum cells treatable per day
QH Dengue Management Plan
20% per day
Maximum cells treatable per day
Additional file 2: Animation of simulated 2003 Dengue outbreak, showing mosquito density in grey and human Dengue cases in red.(MP4 700 KB)
Further information on the dengue control sub-model may be found in the Additional file 1.
Model calibration using dengue outbreak data
Sensitivity of the model to A: Main model parameters and B: Alternate settings
Mosquito to human transmission probability
Human to mosquito transmission probability
No mosquito mobility
No human mobility
No state of emergency
Index cell with low mosquito density
Following calibration of the model using the 2003 dengue outbreak data, the model was applied to a subsequent, larger epidemic which occurred in Cairns in 2008/2009 (see Figure 4). The 2008/2009 epidemic caused nearly 700 cases in the modelled area, twice as many as the 2003 outbreak. Ritchie et al.  discuss several factors that are thought to have contributed to the explosive expansion of the epidemic in late 2008 and its rapid collapse in April 2009. These contributing factors are: i) climatic factors (an unusually warm period in November 2008), ii) dengue virus related factors, specifically the shorter extrinsic incubation period of the DENV3 strain that caused the 2008/09 epidemic, compared to the 2003 DENV2 outbreak , iii) human factors, specifically the increased human movement over the Christmas period and iv) the declaration of the outbreak to be an epidemic, invoking a significant expansion of control measures in January 2009 .
Using the simulation model it was demonstrated that all of these factors may have acted together to cause the larger scale of outbreak which occurred in 2008/2009 compared to previous dengue outbreaks in Cairns, though some of these factors may have played a more significant role than others. Applying the 2008/2009 temperature profile resulted in an average predicted extrinsic incubation period that was approximately 1 day (0.75 days) shorter than that in 2003. Furthermore, the rainfall pattern in 2008/2009 did not result in a greater predicted mosquito population than in 2003, which is in agreement with mosquito trapping data . As a result, climatic factors alone did not increase the simulated case numbers significantly.
In contrast, simulation experiments demonstrated that the shorter extrinsic incubation period of the specific DENV3 strain that caused the 2008/2009 epidemic (as discussed in ) allowed a transmission cycle to be completed in approximately 10 days compared to 17 days in 2003, causing a considerable rise in case numbers and an epidemic that could not be controlled using the initial (2003-based) control measures. The modelled outbreak with significantly increased control measures activated indicated that it was the expansion of control interventions starting in January 2009, following declaration of the outbreak as an epidemic, which caused the decline in weekly case numbers and eventually a rapid collapse of the outbreak. Experiments using the model with the control measures set at the lower 2003-like level resulted in a much longer outbreak which relied on seasonal weather changes before the outbreak was contained, confirming the need for the increased response measures. Increased human movement over the Christmas period (in the model, the month of December) only slightly altered the overall case numbers (by an average of about 50 cases in total).
Differences between the 2003 and 2008/09 outbreak scenarios accounted for in the present study
Index case cell location
Onset of Control
based on 2003 weather
based on 2008 weather
Extrinsic incubation period
3-4 times per week
6-7 times per week for Dec. 2008
constant after day 43
increased after day 60 (epidemic)
Additional file 3: Animation of simulated 2008/9 Dengue outbreak, showing mosquito density in grey and human Dengue cases in red.(MP4 909 KB)
Following its development and testing, the model was used to investigate how changes to control measures affect 2003-like and 2008/2009-like outbreaks, providing guidance for the management of future outbreaks.
Effect of variation of activation of control measures
The results are in agreement with the previous study by Vasquez-Prokopec et al. , showing that variation of the start date of vector control had a major impact of the overall scale of the simulated outbreaks. In simulations of 2003-like outbreaks, initiating vector control 2 weeks earlier than that which occurred in the actual outbreak reduced the median number of cases by ~50%. Starting control 2 weeks later than the original start date approximately doubled the predicted median number of cases. In all scenarios the control measures were able to reduce weekly case numbers by a similar amount and within a similar time-frame as that observed in the original outbreaks (approximately 150 and 200 days respectively for 2003 and 2008/2009).
In the simulations of a 2008/2009-like epidemic, a 2 week earlier onset of vector control reduced the median case number to 32 (from 692), a reduction of nearly 20 fold. This large difference is due to the very early onset of control (day 13 after index case was diagnosed), resulting in onward transmission having failed to spread beyond a 200 m radius of the index case. The location of the index case was an area where control measures were actively applied, eradicating almost all mosquitoes within a couple of days (compare e.g., Figure 7 or Additional file 2 and Additional file 3). As with the 2003 simulations, lengthening the unmitigated phase of the 2008/2009-like epidemic approximately doubled case numbers but did not prolong the overall duration of the outbreak, which corresponds with previously published estimates by Vasquez-Prokopec et al. .
Effect of increased interventions
Sensitivity analyses of the model parameters such as biting rate, fraction of asymptomatic infections, mosquito to human transmission probability and human to mosquito transmission probability were conducted. The values of the relevant model parameters have been adjusted, and the simulation model rerun, to determine the sensitivity of the model to particular parameter settings. While particular parameters were adjusted, no model recalibration took place. The outcomes, for the 2003 and 2008/2009 outbreaks, are presented in Table 2.
The model has been customised for the Cairns setting, however we believe that the underlying modelling methods will be applicable to other settings. In the future, an aim will be to explore how the modelling approach presented here can be applied to other urban dengue transmission settings such as areas with endemic dengue, areas with different population structures, and to larger cities. With appropriate re-parameterisation, the model may also be applicable to other Ae. aegypti (and Ae. albopictus) vectored viruses such as chikungunya . The transfer of the model structure to other geographic locations will depend upon on the availability of data, which was extensive in the case of Cairns (a table of required data sources is given in Additional file 1: Table S4). Furthermore, a transfer of the model to a location endemic for dengue may also require changes to the model structure by introducing measures for herd immunity or properties of multiple circulating dengue strains .
This study highlights a direction of vector-borne infectious disease modelling away from compartmentalised population-based, spatially-homogenous approaches towards individual-based, spatially-explicit techniques, a move also evident in modelling studies focused on other vector-borne diseases such as malaria [48, 49].
A key goal in developing the model has been to capture the most important features of mosquito population dynamics, relevant human behaviour, dengue transmission and control without making the model overly and unnecessarily complex. It has been recognized that a challenge with developing complex models such as that presented here is to avoid the temptation to make them unnecessarily overcomplicated, as discussed by Basu et al. . The model presented here simplifies certain aspects of the physical world being modelled, usually due to limitations of data availability. These limitations include: i) lack of detailed mosquito trapping data prevented the development of a more detailed mosquito density layer, ii) simplified human movement (e.g., the use of semi-random movement is only a rough approximation of reality), iii) simplified mosquito movement (it may be that mosquitoes do not fly randomly but are directed by various factors such as availability of humans, wind direction etc.), iv) simplified effect and dynamics of vector control (lag time between a case and local control may vary on a case-by-case basis, properties may only be partially treated etc.), v) abstraction of mosquito breeding sites (we do not model every individual container that may serve as a breeding site as with the CIMSIM/DENSIM models) . The model is less detailed in the way it accounts for the presence of mosquito breeding sites compared to the CIMSIM/DENSIM approach, yet in the light of recent findings that some of the detail incorporated into the CIMSIM mosquito population dynamics model may be superfluous this approach seems justified . The validation of the model using high-quality data from a different outbreak to that used to calibrate uncertain model parameters suggests that the level of abstraction adopted is appropriate, and results in a modelling framework which will be capable of use to examine the effectiveness of new intervention strategies as well as helping to explain phenomena contributing to the conditions necessary for dengue to become endemic in Cairns-like settings.
This study describes the development of a complex, spatially-explicit and individual-based dengue transmission model, and makes use of spatio-temporal dengue outbreak data in a novel way to demonstrate that such a model is capable of simulating dengue outbreaks of the type that occur in Cairns, along with the intervention measures used to control them. The study clearly shows the need for high quality outbreak field data to inform and parameterize the model. If such data are available, as was the case in the present study, such a simulation model can be adequately calibrated.
This study presents a dengue modelling framework and highlights how key features of the physical system may be represented in a simulation model. While this modelling methodology has been validated by showing that the system as a whole reproduces the outbreak dynamics of a single dengue epidemic, the availability of new field data will permit further model validation and future model refinement, aiming at making the system more physically realistic. New field data collection will permit further validation of the system as a whole, and for the separate validation of each sub-model that makes up the overall modelling environment.
Limitations of this study include the fact that the model assumes a population centre that experiences episodic dengue outbreaks, and not populations where dengue is endemic. As such, it is assumed that the human population has no immunity to dengue, and that only one dengue strain is present during the outbreak. However the methods used to incorporate the necessary sub-models within the Cairns dengue model will also be applicable to the construction of models where dengue is endemic. Another limitation of the study was that Ae. aegypti trapping data was only available for a limited part of the Cairns urban area; availability of comprehensive mosquito trapping data would allow for further refinement of the spatially-heterogeneous vector habitat sub-models. Based on the Ae. aegypti trapping data, it was assumed rainfall was the main driver of the availability of larval habitats, and thus adult population densities. In other geographical locations (for example, where rainfall is uniform throughout the year) other factors may determine the availability of larval habitats and Ae. aegypti population dynamics. While detailed data on human movement is sparse and our use of a daily cycle of moving to work or school captures regular patterns of human movement within a community, our simple semi-random model of other types of (longer range) human movement could be refined if new reliable experimental data becomes available. Finally, we note that each sub-model which contributes to the overall model is based on the best available data and on expert-informed assumptions where experimental data is not available. Further analyses could be conducted to explore the sensitivity of the model to alternative plausible assumptions about, for example, human movement patterns, or factors determining the spatial heterogeneity of larval habitats.
This study was aimed at developing an individual-based, spatially-explicit model for dengue transmission in an urban environment. This research was facilitated by prior development of mosquito population dynamics models by others [15–25, 38, 51, 52] and research studies describing and quantifying Ae. aegypti entomology and dengue outbreaks, virus properties and the impact of vector control in Cairns [7, 10, 12, 13, 36, 39, 40, 44, 52–55]. Specific goals of the study were to produce a simulation modelling framework that could (a) make use of rich dengue outbreak data sets for calibration and validation; (b) would be capable of investigating hypotheses explaining the size of the 2008/2009 Cairns outbreak; and (c) would incorporate phenomena essential to the prediction of future dengue outbreaks and analysis of dengue control measures, viz. weather-dependent vector population dynamics, spatially heterogeneous vector habitat, spatially targeted vector control, and human host movement.
In order to achieve these goals, the resulting model contains a number of key features. These include: i) human and mosquito movement and realistic population structures, ii) weather-dependent (temperature and rainfall) spatially heterogeneous mosquito population dynamics, iii) geographically and demographically- dependent mosquito abundance, iv) spatially- explicit vector control, v) model calibration using outbreak data and vi) model validation against further outbreak data. The resulting simulation model is complex, consisting of several interacting sub-models. However, this complexity is necessary, and is required in order to represent the inherently complex physical phenomena which contribute to dengue transmission and the scale and timing of dengue epidemics.
The simulation results replicating the 2008/2009 Cairns epidemic presented in this study support several hypotheses (formulated previously) aimed at explaining the large-scale epidemic which occurred in 2008/2009 . Specifically, while warmer weather and increased human movement had only a small effect on the spread of the virus, a shorter virus strain-specific extrinsic incubation time can explain the observed explosive outbreak of 2008/2009 . In agreement with previous studies, the simulation results presented here highlight the importance of rapid diagnosis of potential index cases and prompt initiation of vector control, as presented above in Figure 9.
This study, in combination with several previous dengue modelling studies, highlights the importance of spatial and individual-based modelling approaches in order to account for local differences in mosquito and human density and differences in human movement behaviour depending on age and other factors [18, 31, 56]. This study confirms some previously known spatial phenomena, such as human movement and mosquito movement, which facilitate dengue spread (see sensitivity analyses in Table 2, where such movement has been excluded, dramatically reducing simulated case numbers) [31, 56].
The availability of this model will allow further investigation of the effect of different intervention strategies which target various stages of the transmission cycle, such as the effect of potential dengue vaccines, and may also assist with predicting the effect of the release of Wolbachia infected mosquitoes on the native mosquito population, as is currently occurring in Cairns .
Mosquito trapping data was kindly provided by Peter Cook and David Parkington (Monash University, Melbourne, Australia). Gonzalo Vasquez-Prokopec (Emory University, Atlanta, Georgia, USA) kindly provided geo-referenced and de-identified dengue case data. This study has been funded by the University of Western Australia.
- Gubler DJ, Clark GG: Dengue/dengue hemorrhagic fever: the emergence of a global health problem. Emerg Infect Dis. 1995, 1 (2): 55-57.View ArticlePubMedPubMed CentralGoogle Scholar
- Rigau-Perez JG, Clark GG, Gubler DJ, Reiter P, Sanders EJ, Vorndam AV: Dengue and dengue haemorrhagic fever. Lancet. 1998, 352 (9132): 971-977.View ArticlePubMedGoogle Scholar
- Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, Drake JM, Brownstein JS, Hoen AG, Sankoh O, Myers MF, George DB, Jaenisch T, Wint GRW, Simmons CP, Scott TW, Farrar JJ, Hay SI: The global distribution and burden of dengue. Nature. 2013, 496 (7446): 504-507.View ArticlePubMedPubMed CentralGoogle Scholar
- Giatropoulos A, Emmanouel N, Koliopoulos G, Michaelakis A: A study on distribution and seasonal abundance of Aedes albopictus (Diptera: Culicidae) population in Athens, Greece. J Med Entomol. 2012, 49 (2): 262-269.View ArticlePubMedGoogle Scholar
- Sousa CA, Clairouin M, Seixas G, Viveiros B, Novo MT, Silva AC, Escoval MT, Economopoulou A: Ongoing outbreak of dengue type 1 in the Autonomous Region of Madeira, Portugal: preliminary report. Euro Surveill. 2012, 17 (49):Google Scholar
- Pongsumpun P, Garcia Lopez D, Favier C, Torres L, Llosa J, Dubois MA: Dynamics of dengue epidemics in urban contexts. Trop Med Int Health. 2008, 13 (9): 1180-1187.View ArticlePubMedGoogle Scholar
- Ritchie SA, Pyke AT, Hall-Mendelin S, Day A, Mores CN, Cristofferson RC, Gubler DJ, Bennett SN, van den Hurk AF: An explosive epidemic of DENV-3 in Cairns, Australia. PLoS One. 2013, 8 (7): e68137-View ArticlePubMedPubMed CentralGoogle Scholar
- Anderson R, May R: Infectious diseases of humans: dynamics and control. 1991, Oxford: Oxford University Press, 768-Google Scholar
- Grassly NC, Fraser C: Mathematical models of infectious disease transmission. Nat Rev Microbiol. 2008, 6: 477-487.PubMedGoogle Scholar
- Vazquez-Prokopec GM, Kitron U, Montgomery B, Horne P, Ritchie SA: Quantifying the spatial dimension of dengue virus epidemic spread within a tropical urban environment. PLoS Negl Trop Dis. 2008, 4 (12): e920-View ArticleGoogle Scholar
- Duncombe J, Clements A, Davis J, Hu W, Weinstein P, Ritchie S: Spatiotemporal patterns of Aedes aegypti populations in Cairns, Australia: assessing drivers of dengue transmission. Trop Med Int Health. 2013, 18 (7): 839-849.View ArticlePubMedGoogle Scholar
- Ritchie SA, Long S, Smith G, Pyke A, Knox TB: Entomological investigations in a focus of dengue transmission in Cairns, Queensland, Australia, by using the sticky ovitraps. J Med Entomol. 2004, 41 (1): 1-4.View ArticlePubMedGoogle Scholar
- Russell RC, Webb CE, Williams CR, Ritchie SA: Mark-release-recapture study to measure dispersal of the mosquito Aedes aegypti in Cairns, Queensland, Australia. Med Vet Entomol. 2005, 19 (4): 451-457.View ArticlePubMedGoogle Scholar
- McMeniman CJ, O'Neill SL: A virulent Wolbachia infection decreases the viability of the dengue vector Aedes aegypti during periods of embryonic quiescence. PLoS Negl Trop Dis. 2010, 4 (7): e748-View ArticlePubMedPubMed CentralGoogle Scholar
- Ellis AM, Garcia AJ, Focks DA, Morrison AC, Scott TW: Parameterization and sensitivity analysis of a complex simulation model for mosquito population dynamics, dengue transmission, and their control. Am J Trop Med Hyg. 2011, 85 (2): 257-264.View ArticlePubMedPubMed CentralGoogle Scholar
- Focks DA, Daniels E, Haile DG, Keesling JE: A simulation model of the epidemiology of urban dengue fever: literature analysis, model development, preliminary validation, and samples of simulation results. Am J Trop Med Hyg. 1995, 53 (5): 489-506.PubMedGoogle Scholar
- Focks DA, Haile DG, Daniels E, Mount GA: Dynamic life table model for Aedes aegypti (diptera: Culicidae): simulation results and validation. J Med Entomol. 1993, 30 (6): 1018-1028.View ArticlePubMedGoogle Scholar
- Legros M, Magori K, Morrison AC, Xu C, Scott TW, Lloyd AL, Gould F: Evaluation of location-specific predictions by a detailed simulation model of Aedes aegypti populations. PLoS One. 2011, 6 (7): e22701-View ArticlePubMedPubMed CentralGoogle Scholar
- Xu C, Legros M, Gould F, Lloyd AL: Understanding uncertainties in model-based predictions of Aedes aegypti population dynamics. PLoS Negl Trop Dis. 2010, 4 (9): e830-View ArticlePubMedPubMed CentralGoogle Scholar
- Barmak DH, Dorso CO, Otero M, Solari HG: Dengue epidemics and human mobility. Phys Rev E Stat Nonlin Soft Matter Phys. 2011, 84: 011901-View ArticlePubMedGoogle Scholar
- Barmak DH, Dorso CO, Otero M, Solari HG: Modelling interventions during a dengue outbreak. Epidemiol Infect. 2014, 142 (3): 545-561.View ArticlePubMedGoogle Scholar
- Otero M, Barmak DH, Dorso CO, Solari HG, Natiello MA: Modeling dengue outbreaks. Math Biosci. 2011, 232: 87-95.View ArticlePubMedGoogle Scholar
- Otero M, Schweigmann N, Solari HG: A stochastic spatial dynamical model for Aedes aegypti. Bull Math Biol. 2008, 70 (5): 1297-1325.View ArticlePubMedGoogle Scholar
- Otero M, Solari HG: Stochastic eco-epidemiological model of dengue disease transmission by Aedes aegypti mosquito. Math Biosci. 2010, 223 (1): 32-46.View ArticlePubMedGoogle Scholar
- Otero M, Solari HG, Schweigmann N: A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. Bull Math Biol. 2006, 68 (8): 1945-1974.View ArticlePubMedGoogle Scholar
- Zeigler BP: Theory of Modeling and Simulation. 2000, London: Academic Press, 2Google Scholar
- Milne GJ, Kelso JK, Kelly HA, Huband ST, McVernon J: A small community model for the transmission of infectious diseases: comparison of school closure as an intervention in individual-based models of an influenza pandemic. PLoS One. 2008, 3 (12): e4005-View ArticlePubMedPubMed CentralGoogle Scholar
- Halder N, Kelso JK, Milne GJ: Cost-effective strategies for mitigating a future influenza pandemic with H1N1 2009 characteristics. PLoS One. 2011, 6 (7): e22087-View ArticlePubMedPubMed CentralGoogle Scholar
- Kelso JK, Halder N, Milne GJ: The impact of case diagnosis coverage and diagnosis delays on the effectiveness of antiviral strategies in mitigating pandemic influenza A/H1N1 2009. PLoS One. 2010, 5 (11): e13797-View ArticlePubMedPubMed CentralGoogle Scholar
- Milne GJ, Baskaran P, Halder N, Karl S, Kelso J: Pandemic influenza in Papua New Guinea: a modelling study comparison with pandemic spread in a developed country. BMJ Open. 2013, 3 (3): e002518-View ArticlePubMedPubMed CentralGoogle Scholar
- Stoddard ST, Forshey BM, Morrison AC, Paz-Soldan VA, Vazquez-Prokopec GM, Astete H, Reiner RC, Vilcarromero S, Elder JP, Halsey ES, Kochel DJ, Kitron U, Scott TW: House-to-house human movement drives dengue virus transmission. Proc Natl Acad Sci U S A. 2013, 110 (3): 994-999.View ArticlePubMedGoogle Scholar
- Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, Heijne J, Sadkowska-Todys M, Rosinska M, Edmunds WJ: Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008, 5: e74-View ArticlePubMedPubMed CentralGoogle Scholar
- Focks DA, Haile DG, Daniels E, Mount GA: Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. J Med Entomol. 1993, 30 (6): 1003-1017.View ArticlePubMedGoogle Scholar
- Chao DL, Halstead SB, Halloran ME, Longini IM: Controlling dengue with vaccines in Thailand. PLoS Negl Trop Dis. 2012, 6 (10): e1876-View ArticlePubMedPubMed CentralGoogle Scholar
- Vezzani D, Velazquez SM, Soto S, Schweigmann NJ: Environmental characteristics of the cemeteries of Buenos Aires City (Argentina) and infestation levels of Aedes aegypti (Diptera: Culicidae). Mem Inst Oswaldo Cruz. 2001, 96 (4): 467-471.View ArticlePubMedGoogle Scholar
- Hanna JN, Ritchie SA, Merritt AD, van den Hurk AF, Phillips DA, Serafin IL, Norton RE, McBride WJ, Gleeson FE, Poidinger M: Two contiguous outbreaks of dengue type 2 in north Queensland. Med J Aust. 1998, 168 (5): 221-225.PubMedGoogle Scholar
- Williams CR, Long SA, Russell RC, Ritchie SA: Field efficacy of the BG-Sentinel compared with CDC Backpack Aspirators and CO2-baited EVS traps for collection of adult Aedes aegypti in Cairns, Queensland, Australia. J Am Mosq Control Assoc. 2006, 22 (2): 296-300.View ArticlePubMedGoogle Scholar
- Azil AH, Long SA, Ritchie SA, Williams CR: The development of predictive tools for pre-emptive dengue vector control: a study of Aedes aegypti abundance and meteorological variables in North Queensland, Australia. Trop Med Int Health. 2010, 15 (10): 1190-1197.View ArticlePubMedGoogle Scholar
- Williams CR, Johnson PH, Ball TS, Ritchie SA: Productivity and population density estimates of the dengue vector mosquito Aedes aegypti (Stegomyia aegypti) in Australia. Med Vet Entomol. 2013, 27 (3): 313-322.View ArticlePubMedGoogle Scholar
- Ritchie SA, Montgomery BL, Hoffmann AA: Novel estimates of Aedes aegypti (Diptera: Culicidae) population size and adult survival based on Wolbachia releases. J Med Entomol. 2013, 50 (3): 624-631.View ArticlePubMedGoogle Scholar
- Wearing HJ, Rohani P, Keeling MJ: Appropriate models for the management of infectious diseases. PLoS Med. 2005, 2: e174-View ArticlePubMedPubMed CentralGoogle Scholar
- Chan M, Johansson MA: The incubation periods of Dengue viruses. PLoS One. 2012, 7 (11): e50972-View ArticlePubMedPubMed CentralGoogle 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
- Rapley LP, Johnson PH, Williams CR, Silcock RM, Larkman M, Long SA, Russell RC, Ritchie SA: A lethal ovitrap-based mass trapping scheme for dengue control in Australia: II. Impact on populations of the mosquito Aedes aegypti. Med Vet Entomol. 2009, 23 (4): 303-316.View ArticlePubMedGoogle Scholar
- : Queensland Dengue Mangement Plan 2010–2015. Fortitude Valley: Queensland HealthGoogle Scholar
- Vazquez-Prokopec GM, Chaves LF, Ritchie SA, Davis J, Kitron U: Unforeseen costs of cutting mosquito surveillance budgets. PLoS Negl Trop Dis. 2010, 4 (10): e858-View ArticlePubMedPubMed CentralGoogle Scholar
- Ruiz-Moreno D, Vargas IS, Olson KE, Harrington LC: Modeling Dynamic Introduction of Chikungunya Virus in the United States. PLoS Negl Trop Dis. 2012, 6 (11): e1918-View ArticlePubMedPubMed CentralGoogle Scholar
- Gurarie D, Karl S, Zimmerman PA, King CH, St Pierre TG, Davis TME: Mathematical modeling of malaria infection with innate and adaptive immunity in individuals and agent-based communities. PLoS One. 2012, 7 (3): e34040-View ArticlePubMedPubMed CentralGoogle Scholar
- Filipe JAN, Riley EM, Drakeley CJ, Sutherland CJ, Ghani AC: Determination of the Processes Driving, the Acquisition of Immunity to Malaria Using, a Mathematical Transmission Model. PLoS Comput Biol. 2007, 3 (12): e255-View ArticlePubMedPubMed CentralGoogle Scholar
- Basu S, Andrews J: Complexity in Mathematical Models of Public Health Policies: A Guide for Consumers of Models. PLoS Med. 2013, 10 (10): e1001540-View ArticlePubMedPubMed CentralGoogle Scholar
- Bannister-Tyrrell M, Williams C, Ritchie SA, Rau G, Lindesay J, Mercer G, Harley D: Weather-driven variation in dengue activity in Australia examined using a process-based modeling approach. Am J Trop Med Hyg. 2013, 88 (1): 65-72.View ArticlePubMedPubMed CentralGoogle Scholar
- Williams CR, Johnson PH, Long SA, Rapley LP, Ritchie SA: Rapid estimation of Aedes aegypti population size using simulation modeling, with a novel approach to calibration and field validation. J Med Entomol. 2008, 45 (6): 1173-1179.View ArticlePubMedGoogle Scholar
- Ritchie SA, Hanna JN, Selvey CE, Russell RC: To the Editor: Spatial incidence of dengue infections in Queensland, Australia. Epidemiol Infect. 2013, 141 (3): 618-View ArticlePubMedGoogle Scholar
- Ritchie SA, Rapley LP, Williams C, Johnson PH, Larkman M, Silcock RM, Long SA, Russell RC: A lethal ovitrap-based mass trapping scheme for dengue control in Australia: I. Public acceptability and performance of lethal ovitraps. Med Vet Entomol. 2009, 23 (4): 295-302.View ArticlePubMedGoogle Scholar
- Williams CR, Long SA, Webb CE, Bitzhenner M, Geier M, Russell RC, Ritchie SA: Aedes aegypti population sampling using BG-Sentinel traps in north Queensland Australia: statistical considerations for trap deployment and sampling strategy. J Med Entomol. 2007, 44 (2): 345-350.View ArticlePubMedGoogle Scholar
- Liebman KA, Stoddard ST, Morrison AC, Rocha C, Minnick S, Sihuincha M, Russell KL, Olson JG, Blair PJ, Watts DM, Kochel T, Scott TW: Spatial dimensions of dengue virus transmission across interepidemic and epidemic periods in Iquitos, Peru (1999–2003). PLoS Negl Trop Dis. 2012, 6 (2): e1472-View ArticlePubMedPubMed CentralGoogle Scholar
- Hoffmann AA, Montgomery BL, Popovici J, Iturbe-Ormaetxe I, Johnson PH, Muzzi F, Greenfield M, Durkan M, Leong YS, Dong Y, Cook H, Axford J, Callahan AG, Kenny N, Omodei C, McGraw EA, Ryan PA, Ritchie SA, Turelli M, O’Neill SL: Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission. Nature. 2011, 476 (7361): 454-U107.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/447/prepub
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.