Skip to main content

Meteorological factors affecting dengue incidence in Davao, Philippines



Dengue fever is a major public health concern in the Philippines, and has been a significant cause of hospitalizations and deaths among young children. Previous literature links climate change to dengue, and with increasingly unpredictable changing climate patterns, there is a need to understand how these meteorological variables affect dengue incidence in a highly endemic area.


Weekly dengue incidences (2011–2015) in Davao Region, Philippines were obtained from the Department of Health. Same period of weekly local meteorological variables were obtained from the National Climatic Data Center (NCDC) and the National Oceanic and Atmospheric Administration (NOAA). Wavelet coherence analysis was used to determine the presence of non-stationary relationships, while a quasi-Poisson regression combined with distributed lag nonlinear model (DLNM) was used to analyze the association between meteorological variables and dengue incidences.


Significant periodicity was detected in the 7 to 14-week band between the year 2011–2012 and a 26-week periodicity from the year 2013–2014. Overall cumulative risks were particularly high for rainfall at 32 mm (RR: 1.67, 95% CI: 1.07–2.62), while risks were observed to increase with increasing dew point. On the other hand, lower average temperature of 26 °C has resulted to an increased RR of dengue (RR: 1.96, 95% CI: 0.47–8.15) while higher temperature from 27 °C to 31 °C has lower RR.


The observed possible threshold levels of these meteorological variables can be integrated into an early warning system to enhance dengue prediction for better vector control and management in the future.

Peer Review reports


Dengue fever is a common mosquito-borne viral disease of humans transmitted by the bite of an infected female adult mosquito namely the Aedes aegypti as the primary vector and Aedes albopictus as the secondary vector. Aedes aegypti is a holometabolous insect with a life cycle consisting of four stages, namely: egg, instars, pupa and adult [1]. It is characterized with a round or globular head structure, a white flat scales in the middle vertex and a slender, black, long and cylindrical proboscis [2]. Dengue can cause fatal complications such as dengue haemorrhagic fever (DHF) and dengue shock syndrome [3] . Dengue virus belongs to the single-stranded RNA virus of the Flaviviridae family that has four viral serotypes namely; DEN-1, DEN-2, DEN-3, DEN-4 [4]. The adult female Aedes aegypti, the primary vector is a small, black mosquito with white markings around its body. Frequently, mosquito vectors lay its eggs in places within or near human dwellings; with only female adult mosquitoes transmitting the dengue virus [3, 5].

Dengue remains a major public health concern in tropical and subtropical areas [3]. Over the last 50 years, dengue incidence has increased by 30-fold and around 2.5 billion people live in areas where dengue is endemic. Moreover, an estimated 50–200 million cases of dengue infections occur annually in the world [6]. The spread of dengue may be partly due to the increase of international travel, unplanned urbanization, rapid increase in population growth, lack of effective vector management, climate change and extreme weather events, and poor socio-economic status [7,8,9].

Dengue disease has risen in an alarming state in the Philippines in recent years. From January 1st to August 6th of 2016, the Philippines’ Department of Health (DOH) reported an estimated suspected dengue cases of 84,085 in the country, which is 15.8% higher compared to the same period of last year in 2015 with only 72,627 reported cases; out of this, 372 resulted to death [10]. Out of the 10 Association of South East Asian Nations (ASEAN), the Philippines ranked fourth for having the highest number of dengue cases as of 2012 [11]. This alarming rate is partly due to several factors such as environmental degradation, climatic condition, lack of clean water supply, inappropriate waste disposal and management, rapid urbanization, increasing population, and poor mosquito surveillance and control system all contributed to the increasing number of dengue cases in the country [12].

The increasing dengue incidence worldwide is caused by several factors and one of them, which is our primary focus, are the meteorological factors. Change in these factors is believed to influence people’s health through the spread of vector-borne diseases [8]. For example, meteorological factors such as temperature, rainfall, and humidity influence the life stages of female adult Aedes mosquitoes. A warm temperature is important to adult mosquitos’ behavior and maturation, especially the larval development rate is shortened [13, 14]. In addition, rainfall provides plenty of breeding sites for mosquito vectors such as puddles, while humidity affects the adult mosquitoes’ survival and biting frequency [15].

Many countries have conducted studies on the relationship between meteorological factors and dengue cases. For instance, in the temperature and dengue relationship studies, different lagged effect was observed. An increase in RR of dengue was reported to be related with an increase in minimum and maximum temperature by up to 1 to 2-month lag period in Brazil [16] and a lag of 1 month for maximum temperature in Mexico [17]. Meanwhile, a longer lag of up to 3 to 4 months was observed in Australia [18] and Barbados [19]. Furthermore, cumulative rainfall and dengue have also been found to have a varying lag effect; such as a 2-week lag in Mexico [20], a 4-week lag in Thailand [21], and a 10-week lag in Taiwan [22]. The differences on the effects of weather on dengue incidences could be due to the different variations of the amount of rainfall or the range of temperatures in the different regions with respect to their geographical locations [23].

Various approaches have been utilized to estimate the risks and the associated delayed risks of various local meteorological variables on dengue incidence, with a variety of linear [24] and non-linear models [25]. Recent methodological advancements have resulted to the utilization of a distributed lag nonlinear model (DLNM) [26], taking into account the bi-dimensionality of the risks, exposure and lag components, evident in previous studies [27,28,29,30]. Further methodological specifications of DLNM are extensively discussed elsewhere [26, 31, 32]. Here, we aim to elucidate the effects of the local meteorological variables, one of the significant driving forces on dengue transmission [33]. This will enable us to determine the period of high risk of dengue infection since we hypothesized that, local meteorological factors may have an impact on the dengue incidences in Davao Region, Philippines. Clarifying the effects of these meteorological factors on dengue could provide an insight into the seasonal mechanisms of the disease thereby rendering information to understand the complex relationship between meteorological factors and health. Moreover, this study could provide valuable information to health officers for a more effective management of the disease for its control and prevention.


Study area

The study site is located in Davao Region, which is situated at the southeastern portion of the third largest island in the Philippines called Mindanao at 7°05′ latitude North and 125°35′ longitude East (Fig. 1). It is comprised of five provinces namely Compostela Valley, Davao del Norte, Davao del Sur, Davao Oriental, Davao Occidental [34]. It has a total land area of 20,357 km2 and with a total population of 4,893,318 in 2015 [35], making it the country’s fifth fastest growing in terms of the population and 5% in the whole country [36]. The climate is tropical with two distinct seasons: a dry and wet season. According to the Philippine Atmospheric Geophysical and Astronomical Services Administration (PAGASA), wet season occurs from the months of April to October, while dry season occurs from the months of November to March [37]. There are four existing climate types in the Philippines based on the country’s modified Coronas classification, namely Type I (dry from November to April, wet from May to October), Type II (seasonal rainfall from November to December), Type III (same as Type I, but with maximum rainfall from May to October) and Type IV (even distribution of rainfall year-round) [38]. Davao region is classified under Type IV with an annual average temperature of 28 °C; a temperature average shared with other major cities across the country [39].

Fig. 1
figure 1

Geographical location of Davao Region, Philippines. Red-colored area indicates the location of Davao Region, which situated in the Southernmost Island of Mindanao, Philippines

Data collection

We conducted a retrospective time-series study of weekly dengue incidences and meteorological variables in Davao Region from 2011 to 2015. Weekly secondary data of dengue incidences were obtained from the Department of Health-Davao Region (DOH-Region XI) [40] and was freely available through formal data request. The case classification for the dengue incidences used in this study were based on the standard definitions of DOH and were not laboratory confirmed. The meteorological variables of interest, daily average temperature (°C), cumulative rainfall per week (mm), and dew point (°C) reported by Francisco Bangoy International Airport in Davao City, Philippines, were extracted from the climate prediction center of the US National Weather Service website NOAA and NCDC (, which were then aggregated to weekly measures. Daily average temperature was computed as the mean of daily maximum and minimum temperatures. Weekly average temperature, dew point, and weekly cumulative rainfall were generated from the average of the seven-day daily observations for each meteorological variable.

Statistical analyses

Wavelet coherence analysis

This was performed to detect changes in the periodicity in the dengue incidence time series. According to Cazelles et al. [41] wavelet approach extracts the time and frequency of a time series and it is the most efficient method in studying non-stationary data. It transforms a time series into a wave [42]. Furthermore, wavelet was used in this study to determine whether the presence of a particular periodic frequency at a given time in dengue incidences corresponded to the presence of the same periodic frequency at the same time in the given meteorological variables [43].

Distributed lag nonlinear model

A quasi-Poisson time series model coupled with DLNM, which was introduced by Armstrong in 2006 [44] and expanded by Gasparrini in 2010 [31, 32], was used to assess the impact of meteorological variables on dengue incidences. The relationship between meteorological factors and diseases are seldom linear, and the former’s effect are commonly delayed [45]. DLNM uses a cross-basis function, which enables users to explore the association between the exposure variable and its lag, which is in weeks in this instance [46]. We assumed that the reported dengue incidences, y, at week t, follows an over-dispersed quasi-Poisson distribution and can be written as shown below:


$$ \log \left(\mathrm{E}\left({y}_t\right)\right]=\upalpha +{\beta}_1{Rainfall}_{t,l}+ ns\left( Week,4\times 5\right) $$

Where y t refers to the weekly dengue incidences on week t, while t denotes the week of the observation, α is the model intercept, Rainfall t,l is the cross-basis matrix for cumulative rainfall as a covariate on time (t) with 4 degrees of freedom (df), l is the lag in weeks, ns is the smoothing parameter specified on a natural cubic spline (NCS), Week is assigned to control for the seasonal variation.

Moreover, we conducted the DLNM analysis for a maximum lag of 12 weeks. For other meteorological variables, we replaced the cross-basis of the rainfall in Eq. 1 with the respective cross-basis of the other meteorological variables. Model selection was facilitated using three performance indicators namely, Quasi-Akaike Information Criterion (QAIC), Root Mean Squared Error (RMSE) and R-squared as indicated in Table 1. For comparison, we also assumed a simple linear regression with the same specifications.

Table 1 The QAIC, RMSE and pseudo R-squared values for the selected models

All the statistical analyses were performed using R software program, version 3.3.1 [47], while the map figure was created using QGIS version 3.0 [48].

Sensitivity analyses

The passive surveillance system, utilized by the country, is prone to under-reporting and may have considerable impact on the effects estimate [49]. To establish the robustness of the effects estimates, we applied an expansion factor (EF) of 7.0 for the suspected number of cases [11, 50]. EFs have been utilized to adjust for the underreporting of the total number of cases from surveillance systems [50]. For model performance, we further analyzed for a multivariate model with DLNM and simple linear models (in Additional file 1: Table S1).


Dengue incidence distribution

Table 2 shows the descriptive statistics of both weekly dengue incidences and local meteorological variables in Davao Region from 2011 to 2015. From 2011 to 2015, there were 38,169 dengue cases reported at the DOH-Davao Region. In the study period, dengue increased from 4115 cases in 2011 to 9507 cases in 2012, and reached a peak in the year 2013 with 10,762 cases, the highest reported number of dengue cases within the study period. Then, it gradually decreased to 8643 in 2014 and at 5142 cases in 2015 (Fig. 2). The highest reported number of dengue was at the 6th week of the year 2013 around the month of February (Additional files 2 and 3: Figures S2 and S3A) with 442 cases while the lowest was in the 13th week around the month of March 2015 with 30 cases.

Table 2 Descriptive statistics of the weekly dengue incidences and local meteorological variables in Davao Region, 2011–2015
Fig. 2
figure 2

Annual distribution of dengue incidences in Davao Region, 2011–2015. Dengue incidences increased from 4115 cases in 2011 to 9507 cases in 2012. The highest number of dengue incidences among the study period was recorded around the year 2013 with 10,762 cases. Dengue incidence gradually decreased in 2014 with 8643 cases towards 2015 with 5142 cases

Figure 3 represents the time trend of the weekly dengue incidences and local meteorological variables in Davao Region. The data showed that the number of dengue incidences were high from the year 2012, 2013 and 2014 (Fig. 3a). A higher number of dengue generally occurred from the rainy months of June throughout October (Additional file 3: Figure S3A), which consequently falls within the same period of the annual monsoon rains.

Fig. 3
figure 3

Time series plot of weekly dengue incidences and local meteorological variables in Davao Region from 2011 to 2015. (a) Dengue incidences, (b) Average temperature (°C), (c) Cumulative rainfall (mm) and (d) Dew point (°C)

Temperature and rainfall distribution

The highest weekly average temperature recorded within the study period was 30.9 °C in the 18th week of 2014 around the month of May (Fig. 3b and Additional file 2: Figure S2), which is the second hottest month of the year after April, followed by 30.3 °C in the 18th week of 2013 on the same month. Meanwhile, the lowest weekly average temperature recorded was 25.9 °C in the 2nd week of 2014 around the month of January (Additional file 3: Figure S3B).

Between the study period from 2011 to 2015, 2012 has the largest volume of cumulative rainfall recorded at around 823.0 mm, while the smallest volume of rainfall was in 2015 with 604.0 mm (Table 3 and Additional file 4: Figure S1). The highest weekly cumulative rainfall reported was around the 51st week of 2013 around the month of December with 102.3 mm of rainfall (Fig. 3c and Additional file 3: Figure S3C).

Table 3 Annual distribution of dengue incidences and local meteorological variables in Davao Region, 2011–2015

Coherence between local meteorological variables and dengue

Wavelet coherence analysis was used to detect periodicity between dengue incidences and local meteorological variables. In Fig. 4, the coherence between dengue incidences and local meteorological variables varied at different periods and the periodicity of the signals were different through time. For the average temperature and dengue incidences, there were significant periodicities detected in the 7- to 14-week band between the year 2011–2012 and a 26-week periodicity was also observed from the year 2013–2014. There were also several mild periodicities at 2- to 6-week periods that appear in an occasional pattern (Fig. 4a). Likewise, for the cumulative rainfall and dengue incidences, wavelet analysis showed a significant 32-week periodicity band from 2011 to 2012 and an 18- to 30-week periodicity from 2013 to 2014. Furthermore, several high periodicities were also observed in an infrequent pattern at around 2- to 10-week bands (Fig. 4b). For dew point and dengue incidences, wavelet results revealed a significant periodicity of a 20- to 26-week from 2012 to 2013 as well as a 50- to 60-week band from 2011 to 2015. There were also several high frequency periodicities at 2- to 6-week (Fig. 4c).

Fig. 4
figure 4

Wavelet coherence analysis between weekly dengue incidences and local meteorological variables in Davao Region from 2011 to 2015. Blue color represents low coherence, while red color shows high coherence. The black lines represent the areas where coherence is significantly high (95% confidence level). The cone of influence shown as a white curve indicates the region not influenced by edge effects. (a) Average temperature (°C), (b) Cumulative rainfall (mm), and (c) Dew point (°C)

Relationship of climate variables and dengue

To analyze the relationship between local meteorological variables and dengue incidences, a quasi-Poisson time series combined with DLNM was used. The effect of rainfall on dengue were different between low and high cumulative rainfall. As shown in the overall cumulative plot, risk was gradually increasing from 20 to 40 mm of rainfall, with a peak at 32 mm (RR: 1.67, 95% CI: 1.07–2.62). On the other hand, lower average temperature of 26 °C has resulted to an increased RR of dengue (RR: 1.96, 95% CI: 0.47–8.15) while higher temperature from 27 °C to 31 °C has lower RR (Fig. 5b). Meanwhile, for dew point, risks were increasing 25.3 °C. The highest RR was observed at dew point value of 26 °C (RR: 3.10, 95% CI: 1.20–8.06) (Fig. 5c).

Fig. 5
figure 5

Estimated overall effects of the RR of dengue incidences relative to the meteorological variables. (a) Rainfall (mm), (b) Average temperature (°C) and (c) Dew point (°C) with a maximum lag of 12 weeks, using a “natural cubic B-spline-natural cubic spline” DLNM. The center red line represents the estimated spline curve, while the shadow shows 95% CI


In this study, we have observed significant but varying non-stationary periodicities between the local meteorological variables and that of dengue incidence. Further analysis indicated the effects of the varying levels of meteorological variables on dengue incidence. Findings from our study can be utilized for an integrated dengue early warning system, relevant for disease control and management.

In Fig. 4, we have observed significant, but varying inter-annual periodicities of dengue incidence with respect to the specific meteorological variables of interest. In particular, there are consistent longer significant bands from the 16th to 32nd week, for all the three local variables, which coincides with the summer month of April and gradually transitions towards the rainy month of August. There are also small pockets of significant mild periodicities in the 4th to 16th week across the years amongst the three meteorological variables, however, we have observed a significant inter-annual periodicity (2010–2012) from the 6th week to the 15th week, dry season, apparent only with average temperature. Though some studies have noted that mosquito activity would be high in the rainy periods, there are indications that even in dry season mosquito activity may be heightened. Wai et al. [51] observed that vector breeding was enhanced in the dry season in the Philippines due to the presence of water storage, conducive for mosquito growth and development. Tsuzuki et al. [52] also noted the potential of dengue transmission even in hot-dry periods in Nha Trang, Viet Nam, with similarly possible linkage to unchecked and left out indoor water containers.

Beyond the non-stationary relationship observed in the wavelet coherence patterns, we further investigated the impact of local meteorological variables on the dengue morbidity, with the primary focus on identifying the lag (weeks) in which dengue cases has occurred. Based on DLNM, there was a positive association between rainfall, average temperature and dew point with dengue cases. We found that, a moderate amount of rainfall has resulted with increasing RR and gradually decreased as the amount of rainfall increased.

Our findings are consistent with previous studies, whereby a high occurrence of dengue in the few weeks were also observed after moderate rainfall. Ehelepola et al. [53] have found that regular rain favors an increase in dengue but not heavy rain. Similarly, Sarfraz [54] noted that heavy rainfalls may flush away the eggs, larvae and pupae of dengue mosquitoes, which could have consequently affected the mosquito abundance [55]. On the other hand, we have observed that moderate amount of rainfall (20–30 mm) was related to higher dengue incidences in Davao Region. Eastin et al. [56] and Hii et al. [5] noted that light to moderate rainfall can increase the usage of water containers, which are conducive breeding sites for the mosquito.

According to DOH, Davao City which is located near Davao del Sur and a highly urbanized area in Davao Region, has always had the highest number of dengue, comprising around 70% of the cases [57] while the rest are from the neighboring provinces. Dengue is usually higher in highly urbanized areas like cities, where there is overcrowding and poor environmental waste management [58]. In 2013, there was a spike in dengue incidence in Davao City, which was assumed to be linked to the lack of cleanliness drive, unpredictable weather conditions and floods in the area [59] despite the increasing exertion of the DOH and the local health units in Davao Region on the implementation of dengue vector control programs throughout the years. One of the initiated clean-up activities by the DOH is the 4 o’clock habit, which entails the search and destruction of possible mosquito breeding sites, usually done in four in the afternoon [59].

Average temperature at 26 °C has resulted to an increased RR, while higher temperature from 27 °C to 31 °C has lower RR. This was suggestive of the established fact that mosquito development has an optimum range of 25 to 27 °C [60]. This optimum temperature for mosquito strongly enhanced the development from larva to adult, the biting frequency in humans, and the extrinsic incubation period of dengue virus in the mosquito. The decrease in the RR of temperatures from 27 °C to 31 °C is indicative that higher temperature above the optimum range for the Aedes mosquito development brings about a protective effect on dengue transmission [61], as observed in Fig. 5. Higher temperatures may have a negative effect on adult life span of mosquitoes, thereby affecting consequent transmission [62]. In particular, reduced vector competence and activity may result from an increased temperature [63].

Beyond a dew point temperature of 25.3 °C, risks were apparently increasing. Similar observations of significant relationship by dew point on dengue incidence was observed in Brazil [64]. Mechanisms of how dew point affects dengue incidence maybe related to that of the mechanisms posed by humidity. Mathematically and theoretically, dew point and humidity have a nearly linear relationship [65]. Taking this into account, we can observe an identical J-shaped pattern of absolute humidity and RR in a study done in Singapore [29] compared to Fig. 4. High humidity favors an increased longevity of adult mosquitoes as well as the shortening of viral incubation period, thereby allowing an increased transmission intensity [66].

We acknowledge that under-estimation of the burden poses one of the few challenges in establishing the robustness of the effects estimates in this study. However, even after applying an EF of 7, effects estimates as well as the risk curve remained the same, making the estimates robust (in Additional file 5: Figure S4). Furthermore, a recent study in the Philippines noted a high proportion of laboratory confirmed cases (86.1%) from the suspected cases, which thereby indicate that accuracy of clinical diagnosis at admission [67].

This is the very first study in the Philippines which extensively described the association between local meteorological variables and dengue incidences. This study could help improve the dengue surveillance system in the country by taking into account the underlying mechanisms which could be framed in the context of non-stationary relationship as well as the candidate threshold levels, for better dengue prediction. Furthermore, we also acknowledge some limitations in this study. First, the weekly dengue data used in this study were notified suspected dengue cases from clinics and hospitals and are not laboratory confirmed. Second, we did not take into account the mosquito density, population immunity, age classification, social behavior, and socioeconomic conditions for these data are unavailable.


The meteorological variables, though may have varying effects on dengue incidence, exhibit associations which coincide with the plausible biological pathways. Risks were particularly higher in meteorological events with moderate rain, low temperature and high dew point. Furthermore, the observed possible threshold levels of these meteorological variables can be integrated into an early warning system to enhance dengue prediction for better vector control and management in the future.



Association of South East Asian Nations


Dengue haemorrhagic fever


Distributed lag nonlinear model


Department of Health


Expansion Factor


National Climatic Data Center


Natural cubic spline


National Oceanic and Atmospheric Administration


Philippine Atmospheric Geophysical and Astronomical Services Administration


Quasi-Akaike Information Criterion


Root Mean Squared Error; RR – Relative risk


  1. Farnesi LC, Brito JM, Linss JG, Pelajo-Machado M, Valle D, Physiological RGL. Morphological aspects of Aedes aegypti developing larvae: effects of the chitin synthesis inhibitor novaluron. PLoS One. 2012;7:e30363.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Andrew J, Bar A. Morphology and morphometry of Aedes aegypti adult mosquito. Annu Rev Res Biol. 2013;3:52–69.

    Google Scholar 

  3. Gubler DJ. Dengue and dengue hemorrhagic fever. Clin Microbiol Rev. 1998;11:480–96.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Johnson BW, Russell BJ, Lanciotti RS. Serotype-specific detection of dengue viruses in a Fourplex real-time reverse transcriptase PCR assay. J Clin Microbiol. 2005;43:4977–83.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Hii YL, Rocklov J, Ng N, Tang CS, Pang FY, Sauerborn R. Climate variability and increase in intensity and magnitude of dengue incidence in Singapore. Glob Health Action. 2009;2

  6. Murray NEA, Quam MB, Wilder-Smith A. Epidemiology of dengue: past, present and future prospects. Clinical Epidemiology. 2013;5:299–309.

    PubMed  PubMed Central  Google Scholar 

  7. Morin CW, Comrie AC, Climate EK. Dengue transmission: evidence and implications. Environ Health Perspect. 2013;121:1264–72.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Ebi KL, Nealon J. Dengue in a changing climate. Environ Res. 2016;151:115–23.

    Article  PubMed  CAS  Google Scholar 

  9. 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.

  10. DOH: 372 deaths due to dengue since January. Accessed 11 Apr 2017.

  11. Edillo FE, Halasa YA, Largo FM, Erasmo JNV, Amoin NB, Alera MTP, et al. Economic cost and burden of dengue in the Philippines. Am J Trop Med Hyg. 2015;92:360–6.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Bravo L, Roque VG, Brett J, Dizon R, L’Azou M. Epidemiology of dengue disease in the Philippines (2000-2011): a systematic literature review. PLoS Negl Trop Dis. 2014;8:e3027.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 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:1003–17.

    Article  PubMed  CAS  Google Scholar 

  14. Patz JA, Martens WJ, Focks DA, Jetten TH. Dengue fever epidemic potential as projected by general circulation models of global climate change. Env Heal Perspect. 1998;106:147–53.

    Article  CAS  Google Scholar 

  15. Banu S, Guo Y, Hu W, Dale P, Mackenzie JS, Mengersen K, et al. Impacts of El Niño southern oscillation and Indian Ocean dipole on dengue incidence in Bangladesh. Sci Rep. 2015;5 NOVEMBER:16105. doi:

  16. Gomes AF, Nobre AA, Cruz OG. Temporal analysis of the relationship between dengue and meteorological variables in the city of Rio de Janeiro, Brazil, 2001-2009. Cad Saude Publica. 2012;28:2189–97.

    Article  PubMed  Google Scholar 

  17. Chowell G, Sanchez F. Climate-based descriptive models of dengue fever: the 2002 epidemic in Colima, Mexico. J Environ Health. 2006;68:40–4. 55

    PubMed  Google Scholar 

  18. Bi P, Tong S, Donald K, Parton KA, Hobbs J. Climate variability and the dengue outbreak in Townsville, Queensland, 1992-93. Environ Health. 2001;1:54.;dn=229296751516574;res=IELHEA

    Google Scholar 

  19. Depradine C, Lovell E. Climatological variables and the incidence of dengue fever in Barbados. Int J Environ Health Res. 2004;14:429–41.

    Article  PubMed  Google Scholar 

  20. Brunkard JM, Cifuentes E, Rothenberg SJ. Assessing the roles of temperature, precipitation, and ENSO in dengue re-emergence on the Texas-Mexico border region. Salud Publica Mex. 2008;50:227–34.

    Article  PubMed  Google Scholar 

  21. Jeefoo P, Tripathi NK, Souris M. Spatio-temporal diffusion pattern and hotspot detection of dengue in Chachoengsao province, Thailand. Int J Environ Res Public Health. 2011;8:51–74.

    Article  PubMed  Google Scholar 

  22. Chen MJ, Lin CY, Wu YT, Wu PC, Lung SC, Su HJ. Effects of extreme precipitation to the distribution of infectious diseases in Taiwan, 1994-2008. PLoS One. 2012;7:e34651.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Sun W, Xue L, Xie X. Spatial-temporal distribution of dengue and climate characteristics for two clusters in Sri Lanka from 2012 to 2016. Sci Rep. 2017;7:12884.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Johansson MA, Reich NG, Hota A, Brownstein JS, Santillana M. Evaluating the performance of infectious disease forecasts: a comparison of climate-driven and seasonal dengue forecasts for Mexico. Sci Rep. 2016;6:33707.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Xiang J, Hansen A, Liu Q, Liu X, Tong MX, Sun Y, et al. Association between dengue fever incidence and meteorological factors in Guangzhou, China, 2005–2014. Environ Res. 2017;153:17–26.

    Article  PubMed  CAS  Google Scholar 

  26. Gasparrini A. Modeling exposure-lag-response associations with distributed lag non-linear models. Stat Med. 2014;33:881–99.

    Article  PubMed  Google Scholar 

  27. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Horta MA, Bruniera R, Ker F, Catita C, Ferreira AP. Temporal relationship between environmental factors and the occurrence of dengue fever. Int J Environ Health Res. 2014;24:471–81.

    Article  PubMed  CAS  Google Scholar 

  29. Xu HY, Fu X, Lee LKH, Ma S, Goh KT, Wong J, et al. Statistical modeling reveals the effect of absolute humidity on dengue in Singapore. PLoS Negl Trop Dis. 2014;8:e2805.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Chuang TW, Chaves LF, Chen PJ. Effects of local and regional climatic fluctuations on dengue outbreaks in southern Taiwan. PLoS One. 2017;12:e0178698.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Gasparrini A. Distributed lag linear and non-linear models in R: the package dlnm. J Stat Softw. 2011;43:1–20.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Gasparrini A, Armstrong B, Kenward MG. Distributed lag non-linear models. Stat Med. 2010;29:2224–34.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Hales S, de Wet N, Maindonald J, Woodward A. Potential effect of population and climate changes on global distribution of dengue fever: an empirical model. Lancet. 2002;360:830–4.

    Article  PubMed  Google Scholar 

  34. PSA. Philippine Standard Geographic Code. Accessed 9 Apr 2017.

  35. PSA. Census of Population (2015): Highlights of the Philippine Population 2015 Census of Population. Accessed 9 Apr 2017.

  36. PSA. Population and Annual Growth Rates for The Philippines and Its Regions, Provinces, and Highly Urbanized Cities. 2010 Census and Housing Population. and Annual Growth Rates for The Philippines and Its Regions%2C Provinces%2C and Highly Urbanized Cities Based on 1990%2C 2000%2C and 2010 Censuses.pdf. Accessed 9 Apr 2017.

  37. PAGASA. Accessed 12 Apr 2017.

  38. Seposo XT, Dang TN, Honda Y. Exploring the effects of high temperature on mortality in four cities in the Philippines using various heat wave definitions in different mortality subgroups. Glob Health Action. 2017;10:1368969.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Seposo XT, Dang TN, Honda Y. Effect modification in the temperature extremes by mortality subgroups among the tropical cities of the Philippines. Glob Health Action. 2016;9:31500.

    Article  PubMed  Google Scholar 

  40. DOH. DOH - Region XI. Accessed 12 Apr 2017.

  41. Cazelles B, Chavez M, McMichael AJ, Hales S. Nonstationary influence of El Ni??O on the synchronous dengue epidemics in Thailand. PLoS Med. 2005;2:0313–8.

    Article  Google Scholar 

  42. Johansson MA, Cummings DAT, Glass GE. Multiyear climate variability and dengue - El Nino southern oscillation, weather, and dengue incidence in Puerto Rico, Mexico, and Thailand: a longitudinal data analysis. PLoS Med. 2009;6:e1000168.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Torrence C, Compo GPA. Practical guide to wavelet analysis. Bull Am Meteorol Soc. 1998;79:61–78.

    Article  Google Scholar 

  44. Armstrong B. Models for the relationship between ambient temperature and daily mortality. Epidemiology. 2006;17:624–31.

    Article  PubMed  Google Scholar 

  45. Imai C, Armstrong B, Chalabi Z, Mangtani P, Hashizume M. Time series regression model for infectious disease and weather. Environ Res. 2015;142:319–27.

    Article  PubMed  CAS  Google Scholar 

  46. Talagala T. CLIMATIC FACTORS AND DENGUE INCIDENCE IN COLOMBO DISTRICT distributed lag nonlinear modelling approach to identify relationship between climatic factors and dengue incidence in Colombo District. Sri Lanka Epidemiol Biostat Public Heal. 2015;12

  47. R Core Team. R development Core team. R: A Language and Environment for Statistical Computing. 2017;55:275–86:

    Google Scholar 

  48. QGIS Development Team. QGIS geographic information system. Open Source Geospatial Foundation Project. 2015;

  49. L’Azou M, Moureau A, Sarti E, Nealon J, Zambrano B, Wartel TA, et al. Symptomatic dengue in children in 10 Asian and Latin American countries. N Engl J Med. 2016;374:1155–66.

    Article  PubMed  CAS  Google Scholar 

  50. Undurraga EA, Halasa YA, Shepard DS. Use of expansion factors to estimate the burden of dengue in Southeast Asia: a systematic analysis. PLoS Negl Trop Dis. 2013;7:e2056.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Wai KT, Arunachalam N, Tana S, Espino F, Kittayapong P, Abeyewickreme W, et al. Estimating dengue vector abundance in the wet and dry season: implications for targeted vector control in urban and peri-urban Asia. Pathog Glob Health. 2012;106:436–45.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Tsuzuki A, Duoc VT, Higa Y, Yen NT, Takagi M. High potential risk of dengue transmission during the hot-dry season in Nha Trang City, Vietnam. Acta Trop. 2009;111:325–9.

    Article  PubMed  Google Scholar 

  53. Ehelepola NDB, Ariyaratne K, Buddhadasa WMNP, Ratnayake S, Wickramasinghe M, Promprou S, et al. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  54. Sarfraz MS, Tripathi NK, Tipdecho T, Thongbu T, Kerdthong P, Souris M. Analyzing the spatio-temporal relationship between dengue vector larval density and land-use using factor analysis and spatial ring mapping. BMC Public Health. 2012;12:853.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Seidahmed OME, Eltahir EABA. Sequence of flushing and drying of breeding habitats of Aedes aegypti (L.) prior to the low dengue season in Singapore. PLoS Negl Trop Dis. 2016;10:e0004842.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Eastin MD, Delmelle E, Casas I, Wexler J, Self C. Intra- interseasonal autoregressive prediction of dengue outbreaks using local weather and regional climate for a tropical environment in Colombia. Am J Trop Med Hyg. 2014;91:598–610.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Philippines: “The situation of dengue fever in Davao Region is alarming” - Outbreak News Today. Accessed 22 May 2017.

  58. Mulligan K, Dixon J, Joanna Sinn C-L, Elliott SJ. Is dengue a disease of poverty? A systematic review. Pathog Glob Health. 2015;109:10–8.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Dengue cases reach 5,757 in Davao City; floods, poor sanitation blamed | Inquirer News. Accessed 23 May 2017.

  60. Gama ZP, Nakagoshi N. Climatic variability and dengue Haemaorrhagic fever incidence in Nganjuk District, east java. Indonesia Acta Biol Malaysiana. 2013;2:31–9.

    Google Scholar 

  61. Focks DA, Barrera R. Dengue Transmission Dynamics: Assessment And Implications For Control. Working Paper for the Scientific Working Group on Dengue Research. 2006;2006:1–5.

  62. Alto BW, Temperature BD. Dengue virus infection in mosquitoes: independent effects on the immature and adult stages. Am J Trop Med Hyg. 2013;88:497–505.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Hardy JL, Rosen L, Kramer LD, Presser SB, Shroyer DA, Turell MJ. Effect of rearing temperature on transovarial transmission of St. Louis encephalitis virus in mosquitoes. Am J Trop Med Hyg. 1980;29:963–8.

    Article  PubMed  CAS  Google Scholar 

  64. Correia Filho WLF. Influence of meteorological variables on dengue incidence in the municipality of Arapiraca, Alagoas, Brazil. Rev Soc Bras Med Trop. 2017;50:309–14.

    Article  Google Scholar 

  65. Lawrence MG. The relationship between relative humidity and the dewpoint temperature in moist air: a simple conversion and applications. Bull Am Meteorol Soc. 2005;86:225–33.

    Article  Google Scholar 

  66. Mutheneni SR, Morse AP, Caminade C, Upadhyayula SM. Dengue burden in India: recent trends and importance of climatic parameters. Emerg Microbes Infect. 2017;6:e70.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Capeding MRZ, L’Azou M, Manalaysay M, Vince-Woo CR, Rivera RG, Sy AK, et al. Laboratory-confirmed dengue in children in three regional hospitals in the Philippines in 2009-2010. Pediatr Infect Dis J. 2015;34:1145–51.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We would also like to thank the National Oceanic and Atmospheric Administration for the meteorological data, as well as the Department of Health-Davao Region (Region XI) regional director Dr. Abdullah B. Dumama Jr., for generously providing the weekly dengue data report. Without their generosity, this study wouldn’t be possible.


This research is supported by the Environment Research and Technology Development Fund (S-10 and S-14) of the Ministry of the Environment, Japan, and by the Global Research Laboratory grant (Grant K21004000001-10A0500–00710) through the National Research Foundation of Korea, which is funded by the Ministry of Education, Science and Technology, Korea.

Availability of data and materials

The availability of data is subject to both request from the first author, and the approval from the respective data handler; in this case the DOH – Davao Region.

Author information

Authors and Affiliations



JAI collected and analyzed the data. JAI wrote the manuscript. XTS analyzed the data. XTS and YH provided technical support, drafted and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xerxes T. Seposo.

Ethics declarations

Ethics approval and consent to participate

Since secondary data were utilized in this study, no ethics approval was necessary.

Consent for publication

No data or information in this study included an individual person’s data.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Sensitivity analysis with the combination of DLNM and simple linear models. The table depicts that no matter the combination of the linear and DLNM models, there is not much improvement in the model performance, with similar observations in the previous univariate linear and DLNM models. (DOCX 14 kb)

Additional file 2:

Figure S2. Weekly distribution of dengue incidence, average temperature, and cumulative rainfall from 2011 to 2015 in Davao Region. Red line is the cumulative rainfall, dotted blue line are the dengue incidences, and the green dot-and-line is the average temperature. The right-hand side y-axis is in degrees Celsius (for the average temperature), while the left-hand side y-axis is for the dengue incidence and rainfall levels. (JPEG 38 kb)

Additional file 3:

Figure S3. Distribution of monthly dengue incidences and local meteorological variables from 2011 to 2015. (A) Box plot of the monthly dengue incidences, (B) Average temperature (°C), (C) Cumulative rainfall (mm), and (D) Dew point (°C). The horizontal line in the middle of each box is the mean, while the top and bottom borders of the box represent the 25th and 75th percentiles, respectively and the whiskers indicates the 10th and 90th percentiles. (TIF 234 kb)

Additional file 4:

Figure S1. Annual distribution of cumulative rainfall (mm) in Davao Region, 2011–2015. Annual variations in the rainfall, with 2012 recording the highest level (of rainfall). (TIFF 6784 kb)

Additional file 5:

Figure S4. EF-applied dose-response risk curve. Relatively the same risk patterns can be observed, with heightened risks for moderate rain, low temperature and high dew point. (TIFF 1551 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Iguchi, J.A., Seposo, X.T. & Honda, Y. Meteorological factors affecting dengue incidence in Davao, Philippines. BMC Public Health 18, 629 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: