- Research note
- Open Access
Climatic factors influencing dengue incidence in an epidemic area of Nepal
BMC Research Notes volume 12, Article number: 131 (2019)
Geographic expansion of dengue incidence has drawn a global interest to identify the influential factors that instigate the spread of this disease. The objective of this study was to find the environmental factors linked to dengue incidence in a dengue epidemic area of Nepal by negative binomial models using climatic factors from 2010 to 2017.
Minimum temperature at lag 2 months, maximum temperature and relative humidity without lag period significantly affected dengue incidence. Rainfall was not associated with dengue incidence in Chitwan district of Nepal. The incident rate ratio (IRR) of dengue case rise by more than 1% for every unit increase in minimum temperature at lag 2 months, maximum temperature and relative humidity, but decrease by .759% for maximum temperature at lag 3 months. Considering the effect of minimum temperature of previous months on dengue incidence, the vector control and dengue management program should be implemented at least 2 months ahead of dengue outbreak season.
Dengue, caused by Dengue virus (DENV), is the most important mosquito-borne viral disease in the world, impacting half the world population [1, 2]. It is transmitted to human by Aedes aegypti and Ae. albopictus mosquitoes, and is an emerging health problem in Nepal . Although classically associated with tropical climates, occurrence of dengue was reported in Nepal for the first time in 2004 from Chitwan, a lowland district with sub-tropical climate [4, 5]. Subsequently, infection has spread towards north upland hilly districts with temperate climate [6, 7]. Dengue outbreaks in Nepal have occurred in a 3 year cyclical pattern, markedly in 2010, 2013, and 2016 . Upon its prevalence dengue has now been endemic in multiple districts of Nepal . With expansion of DENV affected territory from 9 to 30 districts, reportedly 13 million denizens residing in the lowland Terai belt are at the risk. In addition, DENV-infected patients being identified from the hilly regions in 2013 outbreak has expanded the potential risk territory beyond the lowland belt. The spread of DENV in concurrent with the evidence of stable population of primary vector Ae. aegypti at an altitude of 1310 m [9, 10] suggest the growing threat of dengue in the country.
Global spread of DENV has been attributed to socio-demographic and climatic factors [11, 12]. Vector lifecycle, viral replication and host-vector interaction are linked to temperature, rainfall and humidity making dengue a climate sensitive disease [13,14,15,16,17,18,19]. Thus, any shift in temperature to favorable condition with adequate level of precipitation can result in expansion of vector habitats and prolongation of the transmission season.
Numerous studies have been done around the world to relate dengue incidence and climatic factors [15, 20,21,22,23,24,25,26,27]. Although, a climate-based spatial risk for dengue in Nepal has been assessed  the temporal trend of dengue in relation to climatic factors has not been studied. Considered as a district with significant dengue incidence, a consistent monthly dengue data since the year 2010 was recorded only from Chitwan. Thus we deemed it best to identify climatic variables and the resultant lag periods associated with dengue in Chitwan. This research will provide information to strengthen dengue control strategy and serve as a model system to investigate other Aedes-transmitted disease such as Chikungunya which has been recently detected in Nepal.
This study was conducted in Chitwan district of Nepal. Geographically located between 27°41′ to 27°70′ North latitude and 84°26′ to 84°52′ East longitude, it is situated in the mid-western region of country at an elevation of 415 m above mean sea level. The daily temperatures range between 5 to 22 °C in winter and 24 to 38 °C in summer. The monsoon season spans from June till September with average of rainfall > 100 mm.
Monthly dengue cases from 2010 to 2017 were collected from District Public Health Office, Chitwan, Nepal. Dengue cases were reported from hospital based on detection of nonstructural protein 1 (NS1) antigen and presence of anti-dengue IgM. Daily records of maximum temperature (°C), minimum temperature (°C), rainfall (mm), and relative humidity (percentage) were obtained from the Department of Hydrology and Meteorology, Nepal. Except for the rainfall, which was calculated as cumulative monthly rainfall; monthly means were calculated from daily records of the remaining climate variables.
Maximum temperature, minimum temperature, rainfall and relative humidity were predictor variables while dengue case was response variable. Spearman correlation was assessed between dengue cases and climate factors to identify the most influencing preceding months (lag period) on the occurrence of dengue fever. The predictor variable of 1 month earlier was considered lag1 and the month corresponding to that of dengue case was lag 0. We selected lag period of one to 3 months for preceding months.
Since dengue case was count variable, the possible distribution may be Poisson or negative binomial. However, dengue cases were over dispersed with the variance greater than mean. Further, results of Kolmogorov–Smirnov (K–S) test showed that the dengue cases did not fit to the Poisson distribution (K–S Z = 7.098, p < .01). Thus, negative binomial regression for dengue cases as a response variable seemed to be appropriate for modeling. The models were attempted to include some interaction terms based on the fitting criteria. The best fit model was selected and interpreted. Similarly, incidence rate ratio (IRR) between dengue case and climatic factors was assessed to consider the relative risk of dengue incidence in relation to climatic factors. The statistical analyses were conducted using SPSS Statistics software v.21.
During the 8 years period, a total of 2176 dengue cases have been reported. Yearly distribution of dengue cases showed a cyclic pattern with major outbreaks at an interval of 3 years (Additional file 1: Figure S1). A clear seasonal pattern of dengue occurrence was observed with cases progressed from July–August, hit the highest point in September–October and declined by December. The highest dengue incidence shifted from September during 2010 outbreak to October in 2013 and 2016 outbreaks (Additional file 2: Figure S2). The dengue cases reached the peak following the months with highest temperature and rainfall (Fig. 1).
Spearman correlation showed temperature, rainfall and relative humidity have some significant correlation with dengue cases. Maximum temperature was significantly correlated with dengue cases through lag 1–3 month (P < .01) but not for the corresponding month (P = .136) (Fig. 2; Additional file 3: Table S1). The effect of minimum temperature was significant through lag 0–3 months with strong correlation at lag 2 (r = .7595). The strength of correlation with minimum temperature was higher compared to maximum temperature (Fig. 2). Rainfall was significantly correlated with dengue cases through lag 1–3 months with the strongest correlation at lag 2 month (r = .741) (Fig. 2). Relative humidity had moderate correlation at lag 0 month only (Fig. 2).
Negative binomial regression model was fitted for response variables on given predictor variables based on the fitting criteria. Both deviance test and omnibus test supported that seven models were well fitted to the data (Additional file 4: Table S2). The criteria of lowest value of AIC (or BIC) and mean deviance (− 2LL/df) close to one were used to select a better model. The first model included all 16 predictor variables but none were significant under Wald Chi square test (P < .05) (Additional file 5: Table S3). Thus, it did not meet the requirement of only significant predictors from the fitting of negative binomial regression of dengue case. Model 4 with four significant predictors (AIC = 477.15) and the model 5 with three significant predictors (AIC = 476.64) were better models.
Highly significant predictors in Model 4 under the Wald Chi square test (P < .05) (Table 1) were Minimum temperature at lag 2 month (Mini temp_2), Maximum temperature at lag 0 (Max temp), Maximum temperature at lag 3 month (Max temp_3) and Relative humidity at lag 0 month (Relative humidity). In order to make the interpretation simple and straight forward, partial regression coefficients were expressed in the exponential form, i.e., exp(bi). Exponential value greater than one indicated more effect of the predictor on the response variable and less than one indicated less effect, while one indicated no effect. Interpreting partial regression coefficients, the negative sign indicated a decrease in the response variable and the positive sign indicated an increase in the response variable per unit increase in a given predictor keeping the remaining predictors constant.
The variable Mini temp_2 has a coefficient of .363, which implies that for each one degree Celsius increase in minimum temperature at lag 2 month, the expected log of dengue case increases by .363. Similarly for Max temp with a coefficient of .176, one degree increase in maximum temperature of the corresponding month, the expected log of dengue cases increases by .176. Likewise, same implied for Relative humidity with value of .104 where one percentage increase in relative humidity, the expected log of dengue cases increases by .104.
For Max temp_3 which has a coefficient value in minus, for each one degree increase in maximum temperature at lag 3 month, the expected log of dengue cases decreases by .275. Min temp_2 [exp (.363) = 1.437], Max temp [exp (.176) = 1.193] and Relative humidity [exp (.104) = 1.109] have values greater than one, but Max temp_3 [exp (− .275) = .759] has value less than one. This means minimum temperature at lag 2 month, Maximum temperature and Relative humidity at 0 lag have shown greater impact on dengue incidence. Among them, minimum temperature at lag 2 month was the most influential.
Similarly, the percent change in the incident rate ratio (IRR) of dengue cases was a 1.437%, 1.193% and 1.109% increase for every unit increase in Mini temp_2, Max temp and relative humidity respectively, holding other predictors constant. However, for every unit increase in Maxtemp_3, the percent change in the IRR of dengue cases decreased by .759%.
Model 5 was not chosen since it included only three significant predictors and the regression coefficients values were less compared to model 4 (Additional file 4: Table S3).
In this study, important environmental determinant that influenced dengue incidence in an epidemic Chitwan District of Nepal was identified.
A persistent peak in dengue cases each year following the highest rainfall and temperature, observed in the time series graph, in conjunction with significant higher correlation of response variable with these predictor variables at lag periods (Fig. 2) indicated the influence of preceding month’s climatic factors, which was comparable to other countries [12, 28, 29]. We selected lag period of 1 to 3 months taking into consideration the time for an adult Aedes spp. to develop from egg, extrinsic and intrinsic incubation period of virus  and time for patient to visit hospital following onset of symptom.
The correlation analysis is limited to considering individual relationship and does not address the co-effect of other climate factors, thus we further performed regression analysis to elucidate the true association. Despite a significant relation between dengue cases and rainfall observed in correlation analysis, rainfall was not included as a significant predictor in the regression model. The model that included rainfall showed a low coefficient value, however this model was not selected, hence rainfall might have less influence on dengue incidence in Chitwan. Dengue transmission was not related to rainfall in numerous countries [13, 31, 32]. Influence of rainfall on dengue incidence [33, 34] has been linked to increase in vector population . However, it cannot constantly be valid since excessive rainfall has negative impact [36, 37]. Vector breeding was also possible in dry season when water stored artificial containers were available . Ae. aegypti, the major driver of dengue transmission [34, 39] usually breed in artificial containers . Their abundance in Chitwan compared to Ae. albopictus  gave an impression that water storage behavior can also be a potential contributing factor for dengue transmission which needs to be further investigated in Chitwan.
In view of Aedes spp. abundance regulated by temperature than precipitation [41, 42], temperature has profound impact on dengue prevalence. Temperature strongly influenced dengue cases in several countries [20, 24, 28, 31, 43,44,45,46]. High risk of dengue identified with temperature at lag period between 1 and 3 months [20, 28, 31] was close to the minimum temperature at lag 2 months seen in our study. The minimum temperature at lag 2 months in Chitwan was between 23 and 25 °C, which was conducive for vector survival  compared to temperature at lag 0. Besides, higher rate of dengue virus dissemination influenced by mosquito vector at 25 °C compared to 20 °C  was also in favor of Min temp_2.
Amplification of dengue virus within Ae. aegypti was augmented at 24–31 °C  and the maximum temperature of the month with peak dengue incidence in Chitwan was closer to the upper limit of this range. Thus a positive effect of maximum temperature at lag zero month might be due to the influence of vectoral capacity. However, a negative relationship of maximum temperature observed at lag 3 months can be due to temperature higher than 35 °C around this month in Chitwan. The reverse effect of lag 3 months was also observed in Vietnam . Influence of relative humidity on vector survival and dengue transmission in combination with temperature  corroborates the significant effect of maximum temperature and relative humidity in the month with lag zero.
We found dengue peak month shifted from September to October subsequently from 2013. This shift may be linked to climate change, adaptation of vector and human behavior. Though climate change is a current global concern we cannot validate its effect from 8 years study period. As the relationship between climatic factors and dengue incidence vary with countries and geographical location [15, 20, 37, 43, 51], researches to address factors owing to dengue outbreaks and geographic extension are essential for Nepal which represents a diverse climate within a short latitudinal range.
This study has identified temperature and relative humidity as potential contributors, with minimum temperature at lag 2 months being the most significant, for dengue prevalence in lowland Chitwan district of Nepal. We suggest that vector control program for dengue containment will be more effective if implemented from the months of June–July.
Climatic factors are not the sole predictor of vector borne disease like dengue. Socio-demographic components such as population growth, travel or migration rate, water storage habit should also be considered to link with incidence of dengue. This work has been limited to single district, hence in future several districts encompassing diverse geographic region should be studied.
Akaike’s information criterion
Bayesian information criterion
incident rate ratio
Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496:504–7.
Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, et al. Refining the global= spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis. 2012;6:e1760.
Pun SB. Dengue: an emerging disease in Nepal. J Nepal Med Assoc. 2011;51:203–8.
Pandey BD, Morita K, Khanal SR, Takasaki T, Miyazaki I, Ogawa T, et al. Dengue virus, Nepal. Emerg Infect Dis. 2008;14:514–5.
Pandey BD, Rai SK, Morita K, Kurane I. First case of dengue virus infection in Nepal. Nepal Med Coll J. 2004;6:157–9.
Acharya BK, Cao C, Xu M, Khanal L, Naeem S, Pandit S. Present and future of dengue fever in Nepal: mapping climatic suitability by ecological niche model. Int J Environ Res Public Health. 2018;15(2):187.
Gupta BP, Tuladhar R, Kurmi R, Manandhar KD. Dengue periodic outbreaks and epidemiological trends in Nepal. Ann Clin Microbiol Antimicrob. 2018;17:6.
DoHS. Annual Report 2071/72 (2014/2015). Kathmandu: Department of Health Services, Ministry of Health and Population, Government of Nepal; 2016.
Dhimal M, Gautam I, Kreß A, Müller R, Kuch U. Spatio-temporal distribution of dengue and lymphatic filariasis vectors along an altitudinal transect in Central Nepal. PLOS Negl Trop Dis. 2014;8:e3035.
Gautam I, Dhimal MN, Shrestha SR, Tamrakar AS. First record of Aedes aegypti (L.) vector of dengue virus from Kathmandu. Nepal. J Nat Hist Mus. 2009;24:156–64.
Gubler DJ. Dengue, urbanization and globalization: the unholy trinity of the 21(st) century. Trop Med Health. 2011;39:3–11.
Morin CW, Comrie AC, Ernst K. Climate and dengue transmission: evidence and implications. Environ Health Perspect. 2013;121:1264–72.
Bangs MJ, Larasati RP, Corwin AL, Wuryadi S. Climatic factors associated with epidemic dengue in Palembang, Indonesia: implications of short-term meteorological events on virus transmission. Southeast Asian J Trop Med Public Health. 2006;37:1103–16.
Barbazan P, Guiserix M, Boonyuan W, Tuntaprasart W, Pontier D, Gonzalez JP. Modelling the effect of temperature on transmission of dengue. Med Vet Entomol. 2010;24:66–73.
Ehelepola ND, Ariyaratne K, Buddhadasa WM, Ratnayake S, Wickramasinghe M. A study of the correlation between dengue and weather in Kandy City, Sri Lanka (2003–2012) and lessons learned. Infect Dis Poverty. 2015;4:42.
Githeko AK, Lindsay SW, Confalonieri UE, Patz JA. Climate change and vector-borne diseases: a regional analysis. Bull World Health Organ. 2000;78:1136–47.
McMichael AJ, Campbell-Lendrum DH, Corvalán CF, Ebi KL, Githeko AK, Scheraga JD, et al. Climate change and human health: risk and responses. Geneva: World Health Organization; 2003.
Reiter P. Climate change and mosquito-borne disease. Environ Health Perspect. 2001;109:141–61.
Watts DM, Burke DS, Harrison BA, Whitmire RE, Nisalak A. Effect of temperature on the vector efficiency of Aedes aegypti for dengue 2 virus. Am J Trop Med Hyg. 1987;36:143–52.
Choi Y, Tang CS, McIver L, Hashizume M, Chan V, Abeyasinghe RR, et al. Effects of weather factors on dengue fever incidence and implications for interventions in Cambodia. BMC Public Health. 2016;16:241.
Chuang TW, Chaves LF, Chen PJ. Effects of local and regional climatic fluctuations on dengue outbreaks in southern Taiwan. PLoS ONE. 2017;12:e0178698.
Karim MN, Munshi SU, Anwar N, Alam MS. Climatic factors influencing dengue cases in Dhaka city: a model for dengue prediction. Indian J Med Res. 2012;136:32–9.
Nakhapakorn K, Tripathi NK. An information value based analysis of physical and climatic factors affecting dengue fever and dengue haemorrhagic fever incidence. Int J Health Geogr. 2005;4:13.
Pham HV, Doan HTM, Phan TTT, Minh NNT. Ecological factors associated with dengue fever in a central highlands Province, Vietnam. BMC Infect Dis. 2011;11:172.
Phung D, Talukder MR, Rutherford S, Chu C. A climate-based prediction model in the high-risk clusters of the Mekong Delta region, Vietnam: towards improving dengue prevention and control. Trop Med Int Health. 2016;21:1324–33.
Ramadona AL, Lazuardi L, Hii YL, Holmner A, Kusnanto H, Rocklov J. Prediction of dengue outbreaks based on disease surveillance and meteorological data. PLoS ONE. 2016;11:e0152688.
Sharmin S, Glass K, Viennet E, Harley D. Interaction of mean temperature and daily fluctuation influences dengue incidence in Dhaka, Bangladesh. PLoS Negl Trop Dis. 2015;9:e0003901.
Hii YL, Zhu H, Ng N, Ng LC, Rocklov J. Forecast of dengue incidence using temperature and rainfall. PLoS Negl Trop Dis. 2012;6:e1908.
Minh An DT, Rocklov J. Epidemiology of dengue fever in Hanoi from 2002 to 2010 and its meteorological determinants. Glob Health Action. 2014;7:23074.
Chan M, Johansson MA. The incubation periods of dengue viruses. PLoS ONE. 2012;7:e50972.
Colon-Gonzalez FJ, Lake IR, Bentham G. Climate variability and dengue fever in warm and humid Mexico. Am J Trop Med Hyg. 2011;84:757–63.
Scott TW, Morrison AC, Lorenz LH, Clark GG, Strickman D, Kittayapong P, et al. Longitudinal studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: population dynamics. J Med Entomol. 2000;37:77–88.
Luo L, Li X, Xiao X, Xu Y, Huang M, Yang Z. Identification of Aedes albopictus larval index thresholds in the transmission of dengue in Guangzhou, China. J Vector Ecol. 2015;40:240–6.
Wijayanti SPM, Sunaryo S, Suprihatin S, McFarlane M, Rainey SM, Dietrich I, et al. Dengue in Java, Indonesia: relevance of mosquito indices as risk predictors. PLoS Negl Trop Dis. 2016;10:e0004500.
Micieli MV, Campos RE. Oviposition activity and seasonal pattern of a population of Aedes (Stegomyia) aegypti (L.) (Diptera: Culicidae) in subtropical Argentina. Mem Inst Oswaldo Cruz. 2003;98:659–63.
Baruah S, Dutta P. Seasonal prevalence of Aedes aegypti in urban and industrial areas of Dibrugarh district, Assam. Trop Biomed. 2013;30:434–43.
Campbell KM, Lin CD, Iamsirithaworn S, Scott TW. The complex relationship between weather and dengue virus transmission in Thailand. Am J Trop Med Hyg. 2013;89:1066–80.
Aziz AT, Dieng H, Ahmad AH, Mahyoub JA, Turkistani AM, Mesed H, et al. Household survey of container-breeding mosquitoes and climatic factors influencing the prevalence of Aedes aegypti (Diptera: Culicidae) in Makkah City, Saudi Arabia. Asian Pac J Trop Biomed. 2012;2:849–57.
Campbell LP, Luther C, Moo-Llanes D, Ramsey JM, Danis-Lozano R, Peterson AT. Climate change influences on global distributions of dengue and chikungunya virus vectors. Philos Trans R Soc Lond B Biol Sci. 2015;370(1665):20140135.
Tuladhar R, Singh A, Banjara MR, Gautam I, Dhimal M, Varma A, et al. Effect of meteorological factors on the seasonal prevalence of dengue vectors in upland hilly and lowland Terai regions of Nepal. Parasit Vectors. 2019;12:42.
Moore CG, Cline BL, Ruiz-Tibén E, Lee D, Romney-Joseph H, Rivera-Correa E. Aedes aegypti in Puerto Rico: environmental determinants of larval abundance and relation to dengue virus transmission. Am J Trop Med Hyg. 1978;27:1225–31.
Toma T, Sakamoto S, Miyagi I. The seasonal apperance of Aedes albopictus in Okinawajima, the Ryukyu archipelago, Japan. Mosq News. 1982;42:179–83.
Cazelles B, Chavez M, McMichael AJ, Hales S. Nonstationary influence of El Nino on the synchronous dengue epidemics in Thailand. PLoS Med. 2005;2:e106.
Gharbi M, Quenel P, Gustave J, Cassadou S, La Ruche G, Girdary L, et al. Time series analysis of dengue incidence in Guadeloupe, French West Indies: forecasting models using climate variables as predictors. BMC Infect Dis. 2011;11:166.
Johansson MA, Dominici F, Glass GE. Local and global effects of climate on dengue transmission in Puerto Rico. PLoS Negl Trop Dis. 2009;3:e382.
Lu L, Lin H, Tian L, Yang W, Sun J, Liu Q. Time series analysis of dengue fever and weather in Guangzhou. China. BMC Public Health. 2009;9:395.
Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, et al. Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasit Vectors. 2013;6:351.
Alto BW, Bettinardi D. Temperature and dengue virus infection in mosquitoes: independent effects on the immature and adult stages. Am J Trop Med Hyg. 2013;88:497–505.
Reinhold JM, Lazzari CR, Lahondere C. Effects of the environmental temperature on Aedes aegypti and Aedes albopictus mosquitoes: a review. Insects. 2018;9(4):158.
Naish S, Dale P, Mackenzie JS, McBride J, Mengersen K, Tong S. Climate change and dengue: a critical and systematic review of quantitative modelling approaches. BMC Infect Dis. 2014;14:167.
Pinto E, Coelho M, Oliver L, Massad E. The influence of climate variables on dengue in Singapore. Int J Environ Health Res. 2011;21:415–26.
RT, AS and AV conceived and designed the study. RT and DKC analyzed the data. RT, AS, AV and DKC prepared the manuscript. All authors read and approved the final manuscript.
We are thankful for this research being financed under the Climate Change Research Grants Program implemented by the Nepal Academy of Science and Technology. The programme is part of the Mainstreaming Climate Change Risk Management in Development project. This project is a component of Nepal’s Pilot Program for Climate Resilience and is executed by the Ministry of Population and Environment (Nepal), financed by the Climate Investment Funds, administered by the Asian Development Bank with technical assistance from ICEM, METCON and APTEC.
The authors declare that they have no competing interests.
Availability of data and materials
The data supporting the conclusions of this article are provided within the article and its additional files.
Consent for publication
Ethics approval and consent to participate
Ethics approval was granted by the Ethical Review Board (ERB) of the Nepal Health Research Council (reg. no. 183/2015).
This study was funded by the Nepal Academy of Science and Technology/Asian Development Bank RGP-018 provided to the first author.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Yearly distribution of dengue cases in Chitwan from 2010 to 2017.
Dengue cases by month from 2010 to 2017 in Chitwan district.
Correlation analysis between dengue cases and climate factors with lag effects of 0–3 months period.
A goodness of fit test for a response variable using negative binomial regression model.
Negative binomial regression models of dengue cases.
About this article
Cite this article
Tuladhar, R., Singh, A., Varma, A. et al. Climatic factors influencing dengue incidence in an epidemic area of Nepal. BMC Res Notes 12, 131 (2019). https://doi.org/10.1186/s13104-019-4185-4
- Negative binomial model
- Minimum temperature
- Maximum temperature
- Relative humidity