Empirical model for estimating dengue incidence using temperature, rainfall, and relative humidity: a 19-year retrospective analysis in East Delhi
Article information
Abstract
OBJECTIVES
Aedes mosquitoes are responsible for transmitting the dengue virus. The mosquito lifecycle is known to be influenced by temperature, rainfall, and relative humidity. This retrospective study was planned to investigate whether climatic factors could be used to predict the occurrence of dengue in East Delhi.
METHODS
The number of monthly dengue cases reported over 19 years was obtained from the laboratory records of our institution. Monthly data of rainfall, temperature, and humidity collected from a local weather station were correlated with the number of monthly reported dengue cases. One-way analysis of variance was used to analyse whether the climatic parameters differed significantly among seasons. Four models were developed using negative binomial generalized linear model analysis. Monthly rainfall, temperature, humidity, were used as independent variables, and the number of dengue cases reported monthly was used as the dependent variable. The first model considered data from the same month, while the other three models involved incorporating data with a lag phase of 1, 2, and 3 months, respectively.
RESULTS
The greatest number of cases was reported during the post-monsoon period each year. Temperature, rainfall, and humidity varied significantly across the pre-monsoon, monsoon, and post-monsoon periods. The best correlation between these three climatic factors and dengue occurrence was at a time lag of 2 months.
CONCLUSIONS
This study found that temperature, rainfall, and relative humidity significantly affected dengue occurrence in East Delhi. This weather-based dengue empirical model can forecast potential outbreaks 2-month in advance, providing an early warning system for intensifying dengue control measures.
INTRODUCTION
The dengue virus (DENV) is an arbovirus belonging to the Flaviviridae family, and is a cause of classical dengue fever (DF), dengue haemorrhagic fever (DHF), and dengue shock syndrome (DSS), which are major public health problems in Delhi, India [1-3]. These clinical entities are caused by four DENV serotypes (DENV-1, DENV-2, DENV-3, and DENV-4), transmitted to humans by female Aedes mosquitoes (Aedes aegypti and Aedes albopictus) [3,4]. According to the 2014 to 2015 report of the National Vector Borne Disease Control Programme, dengue is endemic in 35 states and Union territories, including Delhi [5]. Dengue is being increasingly reported in urban areas, mainly due to deficient water management, including improper water storage practices, and inconsistent attention to the elimination of vector breeding sites [5].
All four serotypes of DENV circulate in Delhi [6]. Until 2003, the predominant serotype in Delhi was DENV-2, but in 2003, Delhi became a hyperendemic state with all four DENV serotypes co-circulating. During the 2004 epidemic of DF, an abrupt shift occurred, leading to the dominance of DENV-3 at the expense of the previously circulating serotype DENV-2, followed by the complete preponderance of DENV-3 in 2005 and 2006 [6,7]. Over the next four years (2007 to 2010), DENV-1 emerged as the principal serotype [6,7]. DENV-2 was the preponderant serotype from 2011 to 2014 [8]. In 2015 all four serotypes were found to co-circulate, with DENV-2 predominant [9].
The Aedes mosquito is a climate-sensitive vector that predominantly affects tropical countries due to their climatic conditions [1,2]. Dengue cases are influenced by complex interactions of humans, vectors, the environment, and virus-related factors [10]. Studies have reported a strong and consistent relationship between the climate of a particular geographical area and the number of dengue cases [11]. Different models have been developed to predict dengue outbreaks by correlating dengue cases with climatic data [12,13]. Amongst the different climatic parameters, rainfall, temperature, and humidity have been reported to be the most important factors influencing DENV transmission. Moreover, in geographical regions where minimum thresholds of these climatic variables are adequate to sustain DENV transmission, seasonal fluctuations in these parameters act as essential determinants of the strength and period of transmission [14]. In recent years, the accuracy of predicting local weather and epidemics has improved due to advances in technology [10]. This has helped gain an understanding of the interaction between climate and the temporal-spatial distribution of infectious diseases, as well as encouraging research interest on epidemic prediction modelling [10,15]. Although dengue prediction models have been developed in many countries across the world, no such study has been reported in India [1, 12,13]. Thus, the present study was planned to develop an empirical model to predict monthly dengue cases using 19 years’ monthly data on the number of dengue cases and three major climatic factors (rainfall, temperature, and humidity) at an 1,800-bed tertiary care hospital in East Delhi.
MATERIALS AND METHODS
In the present retrospective study, the number of monthly dengue cases reported at Guru Teg Bahadur Hospital, an 1,800-bed tertiary-care hospital in East Delhi, for a period of 19 years (from January 1997 to December 2015) was obtained. Guru Teg Bahadur Hospital is the largest hospital of the government of the National Capital Territory of Delhi in the Trans-Yamuna Area (East Delhi), with a capacity of 1,800 beds. Of the population of Delhi of 12 million people, 12.07% live in East Delhi. Guru Teg Bahadur Hospital is the only Delhi Government tertiary care hospital in the Trans-Yamuna (East Delhi) area, catering to the population of East Delhi as well as patients from the adjacent districts of Noida, Meerut, Loni, Baghpat, and Bulandshahar. Hence, this hospital handles the majority of the dengue cases in East Delhi, as is reflected in the data from 2015, when Guru Teg Bahadur Hospital encountered 1,633 of the 1,737 dengue cases reported in East Delhi in 2015.
This hospital accepts all cases of suspected DF irrespective of severity. In this study, dengue cases were defined and classified according to the National Guidelines for Clinical Management for DF released by the Government of India in December 2014 [16]. The guidelines classify dengue into undifferentiated DF and severe DF based on clinical manifestations. Non-severe dengue cases include DF and DHF grades I and II, while severe dengue includes DHF grades III and IV and DSS. The clinical criteria for DF, DHF, and DSS given in the guidelines are as follows:
Clinical features of dengue fever
An acute febrile illness of two to seven days’ duration with two or more of the following manifestations: headache, retro-orbital pain, myalgia, arthralgia, rash, and haemorrhagic manifestations.
Dengue haemorrhagic fever
A case with the clinical criteria of DF plus haemorrhagic tendencies (evidenced by: positive tourniquet test or petechiae, ecchymoses, or purpura, or bleeding from mucosa, gastrointestinal tract, injection sites, or other sites) plus thrombocytopenia (<100,000 cells/mm³) plus evidence of plasma leakage due to increased vascular permeability (manifested by a rise in average haematocrit for age and sex of >20% or a drop of more than 20% in haematocrit following volume replacement treatment compared to baseline or signs of plasma leakage (pleural effusion, ascites, hypoproteinaemia [total serum protein level <6 g/dL]) ≤ 20%
Dengue shock syndrome
All the above criteria for DHF with evidence of circulatory failure manifested by rapid and weak pulse and narrow pulse pressure (mmHg) or hypotension for age, cold and clammy skin, and restlessness.
The protocol of this study was approved by the University College of Medical Sciences Institutional Review Board. Informed consent was obtained from the study subjects. A diagnosis of dengue was made on the basis of clinical findings and serology. Serological confirmation of DF during the study period was carried out using the following tests:
a) Panbio (Brisbane, Australia) Dengue Duo Cassette, a rapid diagnostic test, was used to detect anti-dengue immunoglobulin G (IgG) and immunoglobulin M (IgM) antibodies from January 1997 to December 2001.
b) Panbio (Brisbane, Australia) Dengue IgM capture enzymelinked immunosorbent assay (ELISA) was employed from January 2002 to August 2006.
c) IgM dengue capture ELISA (MAC-ELISA) kits supplied by the National Institute of Virology, Pune under the aegis of the National Vector Borne Disease Control Programme were introduced in September 2006 and have been in use since then.
d) SD Bioline Dengue Duo (non-structural protein 1 [NS1] and IgG/IgM) (Standard Diagnostics Inc., Seoul, Korea) was used from October 6, 2012 to December 31, 2013. The IgM-positive samples were confirmed by MAC-ELISA.
e) NS1 antigen capture ELISA from Panbio (Brisbane, Australia) was introduced on October 4, 2014 and used along with MAC-ELISA for the serological confirmation of DF.
Monthly climatic data (rainfall, temperature, and humidity) was collected from the Delhi Weather Station, Safdarjung, New Delhi (http://www.en.tutiempo.net/climate/ws-421820-html). The monthly climate data were correlated with the number of monthly reported dengue cases.
Average monthly rainfall, temperature, and humidity were used as independent variables and the number of dengue cases reported monthly was used as the dependent variable. One-way analysis of variance (ANOVA) was used to determine whether each of the climate variables differed significantly between seasons.
An empirical model was developed using negative binomial generalized linear model analysis. A negative binomial model with a log-link function in the generalized linear model was used to obtain the models for estimating dengue cases based on the independent variables of rainfall, temperature, and humidity. Since a significant correlation was present in the variability of dengue cases across seasons, we included season as a covariate and obtained modified models. A negative binomial model was used because our outcome variable was measured as counts and was over-dispersed; that is, the conditional variance was quite high compared to the conditional mean.
When several maximum likelihood models are available, one can compare the performance of alternative models based on several likelihood or goodness of fit measures. Two of the most regularly used measures are the Akaike information criterion (AIC) and the likelihood ratio chi-square. The model exhibiting maximum change in the value of AIC and likelihood ratio chi-square compared to the previous model was considered to be the best-fitting model.
RESULTS
Over the 19-year period (1997 to 2015), a total of 6,703 inpatient and outpatient cases of DF, including severe forms of the disease such as DHF and DSS, were reported at Guru Teg Bahadur Hospital. The number of reported dengue cases varied by year (Table 1). Over the study period, the highest number of dengue cases was reported in 2015 (n=1,633). Every year, the occurrence of dengue cases displayed a particular pattern. During the pre-monsoon season, hardly any cases of dengue were reported. Most of the cases were reported during the post-monsoon period each year, except in 2010, when the highest number of cases was reported during the monsoon. The average number of dengue cases per month (January to December) over the 19-year period (1997 to 2015) was plotted against the climatic factors (rain, temperature, and relative humidity) to assess their influence on the occurrence of DF (Figures 1 and 2). Every year, a bellwether of cases was reported in the month of July, reaching a peak in September and October, and gradually declining at the end of the year (Table 1).
Seasonal variation in the number of dengue cases each year
Average number of dengue cases and rainfall by month from 1997 to 2015.
Average number of dengue cases, temperature, and relative humidity by month from 1997 to 2015.
Temperature, rainfall, and humidity varied significantly over the pre-monsoon, monsoon, and post-monsoon periods every year (Table 2). The average monthly pre-monsoon, monsoon, and post-monsoon rainfall amounts were 28.72 mm, 150.69 mm, and 11.02 mm, respectively, over the 19-year period (1997 to 2015). The average temperature recorded during the pre-monsoon, monsoon, and post-monsoon period was 25.53˚C, 30.80˚C, and 18.20˚C, respectively. The average relative humidity recorded during the pre-monsoon, monsoon, and post-monsoon period was 49.87%, 68.68%, and 68.83%, respectively. The overall variation in all three climatic variables was statistically significant (ANOVA, p<0.001). On applying cross-correlation analysis between the three climatic factors and dengue occurrence, the best correlation was found at a time lag of 2-month for all three variables.
Distribution of climatic parameters by season (1997-2015)
Development of an empirical model
The number of dengue cases increased in the post-monsoon period, indicating a correlation between dengue infection and climatic factors (rainfall, temperature, and relative humidity), as well as providing a basis for a possible empirical model of dengue. Four models were developed using a negative binomial generalized linear model (Table 3). The first model was developed considering data from the same month. Under the influence of climatic factors, it takes 7 to 45 days for an adult mosquito to develop from an egg [8]. Hence, the influence of climatic factors was expected to manifest with a lag of 1-2 months. Therefore, three other models were generated, incorporating a lag phase of 1, 2, and 3 months, respectively. It was observed that the AIC decreased and the likelihood ratio chi-square increased from model 1 to model 2, and then from model 2 to model 3. However, the difference in both parameters was not significant between model 3 and model 4. Thus, model 3 was chosenas the final model.
Models based on a negative binomial generalized linear model relating the monthly number of dengue cases with rainfall, temperature, and relative humidity
Model 3, 2-month time lag: log Dm=rainfallm×0.001+temperaturem×0.501+humiditym×0.152-22.037; AIC=856.708, likelihood ratio chi-square=849.017.
Another set of models was prepared including season as a covariate in addition to rainfall, temperature, and relative humidity (Table 4). It was found that the AIC decreased and the likelihood ratio chi-square increased from model 1 to model 2, but the difference in the two parameters was minimal between model 2 and model 3. Hence, model 2 could be chosen as the final model if season was also considered as a covariate.
Models based on a negative binomial generalized linear model relating the monthly number of dengue cases with the rainfall, temperature, relative humidity, and season
Model 2, 1-month time lag: log Dm=Rainfallm×0.012+temperaturem×0.314+humiditym×0.058+season×4.794-23.025; AIC =1,020.722, likelihood ratio chi-square=956.592.
DISCUSSION
The growth and development of dengue vectors depends on weather conditions. As dengue is a vector-borne disease, the occurrence of dengue infections depends on the presence and density of its vector [1]. Numerous studies conducted worldwide have proposed that climatic factors, with temperature, rainfall, and humidity being the most important, are the reason for seasonal variation in the presence of both the vector Aedes aegypti and DENV [1, 3,14,17].
In the present study, we found a relationship between dengue cases and rainfall. This is consistent with some other studies [1, 3,17]. Focks & Barrera [17] asserted that vector density increases due to rainfall, causing an increase in dengue cases, which in turn because higher humidity during the rainy season provides an ideal environment for the growth and survival of mosquitoes. However, the few studies in the literature on the association between rainfall and dengue have reported contradictory findings, since the correlation depends on local characteristics [2,18].
This study demonstrated a correlation between the incidence of DF and average temperature. Similar findings were reported by Karim et al. [1] and Chandy et al. [3], while contrasting results were reported in a study by Su [13]. Temperature is a crucial limiting factor in the maturation of the dengue vector. Studies have estimated the threshold survival temperature for the dengue virus to be 11.9°C and have demonstrated that Aedes aegypti ceases to feed when the temperature falls below 17°C, along with non-amplification of the virus in the vector when the temperature falls below 18°C [2]. Therefore, at low temperatures, the mosquito does not live long enough to become infectious and transmit the virus, and the virus does not develop properly. Viral replication and the extrinsic period in the insects are shortened by high temperature. A reduced viral incubation time increases the probability of the vector living long enough to transmit the virus, thereby magnifying epidemics. Hence, temperature analysis improves our understanding of dengue epidemiology [2]. Increases in global temperatures may expand the area of involvement and number of cases of vector-borne diseases [1,13], and will increase the proportion of infective mosquitoes, giving rise to an exponential increase in dengue cases [1]. It has also been shown that a temperature increase from 26˚C-28˚C to 30˚C decreases the extrinsic incubation of the dengue virus, which may promote viral transmission [1,19,20].
The present study demonstrated an association between relative humidity and dengue cases. Relative humidity is influenced by a combination of rainfall and temperature. This important weather parameter affects the lifespan of mosquitoes and hence, viral transmission. Karim et al. [1] and Promprou et al. [21] found more dengue cases during the monsoon when the relative humidity is higher. Higher humidity during the rainy season encourages the development and propagation of mosquitoes, causing an increase in the number of infected mosquitoes [1,22]. Barbazan et al. [22] have postulated that an increase in the mosquito lifespan disproportionately augments the frequency of potential transmissions by as much as five times when the survival rate rises from 0.80 to 0.95. Hales et al. [23] have also reported that the annual average vapour pressure was the most important climatic predictor of global dengue occurrence.
The cumulative effect of temperature and humidity strongly influences the number of mosquito blood meals, survival rate of the vector, and the likelihood of a mosquito becoming infected with DENV [14]. According to studies conducted worldwide, relative humidity and temperature are the most important climatic predictors of changes in dengue transmission [14,24,25]. Hence, through their effects on the Aedes vector, rainfall, temperature, and relative humidity are essential factors determining the geographic areas within which dengue transmission can be expected to occur [14].
In the present study, significant differences were found amongst the pre-monsoon, monsoon, and post-monsoon periods in the amount of rainfall, relative humidity, and average temperature (Table 2). Dengue cases peaked during the post-monsoon period, with a time lag of 1 and 2-month in models developed including and excluding season as a covariate, respectively (Tables 3 and 4). This could be accounted for by the indirect effect of climatic factors on the incidence of dengue through their influence on the lifecycle of both the vector and virus [14]. Climatic factors affect mosquito hatching, larval and pupal development, and the emergence of adult mosquitoes, as well as virus amplification and incubation in humans, finally culminating in a dengue outbreak after a cumulative time lag [14,26]. The time lag can be explained by the influence of weather conditions on the biological development of the mosquito vector, including prolonged egg hatching periods and the propensity of Aedes eggs to survive without water for many months [10]. Depending on the time lag between the biological cycle and the clinical symptoms, a lag between climate data and dengue incidence data will emerge [14]. Our results agree with those of other studies from diverse geographical regions [11,21,27]. A season-specific pattern of dengue cases in Southeast Asia has been reported [1,28]. A similar lag phase of 2-month for explaining the occurrence of dengue cases was reported by Karim et al. [1] in Dhaka. Chen et al. [29] reported a significant positive correlation of dengue cases with temperature and relative humidity at a lag of 1-3 months in Thailand.
In this study, the greatest number of dengue cases was reported in 2015 (n=1,633). As relative humidity and temperature have been reported to be the most important climatic predictors of changes in dengue transmission, the most ideal combination of high relative humidity and high temperature was present in 2015 in comparison to all other years in the study period, thereby resulting in the greatest number of dengue cases. Several studies have reported that the climate variable indices of the Indian Ocean dipole and El Niño southern oscillation play an important role in the interannual variation in dengue transmission [30,31]. Variations in oceanic sea surface temperature have also been found to be responsible for interannual variation in dengue cases [32]. An enhanced case burden is also attributable to a general apathy to civic hygiene, unmonitored drain cleaning and silting, stagnation of rain water, highly unstructured and ill-supervised mosquito control measures, an excess focus on fogging-based vector control, and the institution of remedial measures only after the problem of an epidemic surfaces. These are supplemented by an excessive migrant population, an inadequate health-care delivery system in the districts adjacent to Delhi, and the excess congregation of ill patients accompanied by healthy caregivers in healthcare facilities, providing a milieu for vector-susceptible contact.
Studies have reported slight indications of a correlation between mosquito status and dengue spread that can be employed to predict dengue outbreaks [33]. We collected Breteau index (BI) data from 2007, 2008, and 2009 from the Municipal Corporation of Delhi (Table 5). The mean BI was very low in the pre-monsoon season. The mean BI reached its maximum levels from September to November, which is consistent with our finding of the peak number of dengue cases in September and October. The mean BI was highest in 2008 among the 3-year, which also correlates with the fact that out of those three years, the greatest number of dengue cases was reported in 2008.
Mean Breteau index (BI) before, during, and after dengue outbreaks in East Delhi (2007-2009)
Since this study suggests that temperature, rainfall, and humidity are significantly associated with DF incidence at a time lag of 2-month, mosquito control and dengue surveillance must be strengthened during the post-monsoon season. This window of opportunity provides enough time to mobilize resources for the implementation of interventional measures to minimize the impact of the epidemic [34]. It would permit vector control units to carry out their operations during the high-risk season, thus maximizing limited vector control resources, as well as offering local authorities sufficient time to alleviate a probable outbreak successfully. It also allows the avoidance of vector control measures triggered by false alarms. Disease surveillance measures enabling public health practitioners to be aware of the scale of dengue morbidity and control of mosquitoes are necessary to curtail dengue transmission. The creation of diagnostic centres is essential for facilitating surveillance and providing an early warning of the changes in dengue incidence. Dengue control requires the implementation of an integrated approach, incorporating environmental management, chemical control, and biological methods for the control of mosquitoes [13].
This study has some limitations. Only monthly data were available to us. A weekly analysis of incidence data and weather conditions would help explore the true impact of climatic variability on dengue incidence. Wind velocity is also known to have a positive effect on dengue incidence. However, such data were not available. Another limitation is the use of records from a single hospital as a surrogate for surveillance in a specific area, rather than using direct surveillance data.
This study found that climatic parameters can serve as important components for generating an uncomplicated, accurate, cost-effective, and timely dengue forecasting system. The development of a weather-based dengue forecasting model could assist local vector control, prevention, and surveillance in several ways. The model would act as an early warning system for enhancing measures of dengue control to reduce the size of an outbreak, thereby decreasing disease transmission and possibly the resulting mortality, leading to reductions in the healthcare burden and operating costs. Nonetheless, the long-term sustainability of forecast accuracy remains a potential challenge of the dengue forecasting model, since it assumes that a particular distribution pattern will be repeated in the future, while changes of dengue epidemiology in the long run are unavoidable as factors affecting dengue epidemiology evolve over time. Hence the model may have to be re-calibrated in the future to maintain long-term forecast precision by anticipating alterations in factors influencing dengue transmission and distribution patterns.
Notes
The authors have no conflicts of interest to declare for this study.