Figures
Abstract
Background
Influenza transmission is often associated with climatic factors. As the epidemic pattern varies geographically, the roles of climatic factors may not be unique. Previous in vivo studies revealed the direct effect of winter-like humidity on air-borne influenza transmission that dominates in regions with temperate climate, while influenza in the tropics is more effectively transmitted through direct contact.
Methodology/Principal Findings
Using time series model, we analyzed the role of climatic factors on the epidemiology of influenza transmission in two regions characterized by warm climate: Hong Kong (China) and Maricopa County (Arizona, USA). These two regions have comparable temperature but distinctly different rainfall. Specifically we employed Autoregressive Integrated Moving Average (ARIMA) model along with climatic parameters as measured from ground stations and NASA satellites. Our studies showed that including the climatic variables as input series result in models with better performance than the univariate model where the influenza cases depend only on its past values and error signal. The best model for Hong Kong influenza was obtained when Land Surface Temperature (LST), rainfall and relative humidity were included as input series. Meanwhile for Maricopa County we found that including either maximum atmospheric pressure or mean air temperature gave the most improvement in the model performances.
Conclusions/Significance
Our results showed that including the environmental variables generally increases the prediction capability. Therefore, for countries without advanced influenza surveillance systems, environmental variables can be used for estimating influenza transmission at present and in the near future.
Citation: Soebiyanto RP, Adimi F, Kiang RK (2010) Modeling and Predicting Seasonal Influenza Transmission in Warm Regions Using Climatological Parameters. PLoS ONE 5(3): e9450. https://doi.org/10.1371/journal.pone.0009450
Editor: Matthew Baylis, University of Liverpool, United Kingdom
Received: December 4, 2009; Accepted: February 9, 2010; Published: March 1, 2010
This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.
Funding: This study is supported by NASA Applied Sciences Public Health Program (http://nasascience.nasa.gov/earth-science/applied-sciences). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Influenza remains a global concern as estimates show that the annual epidemics may cause up to five million severe illnesses and 500,000 deaths worldwide [1]. Vaccination strategy and surveillance effort as means of prevention prevail only in temperate regions, especially in the northern hemisphere. Many countries in the subtropics and tropics, on the other hand, underutilize the prevention mechanisms [2] despite the year-round outbreaks and recent findings that (i) East and Southeast Asia act as a source of new influenza strains and transmission [3], [4], and (ii) epidemics timing in South America travels southward starting from the equator [5], [6]. The scarce surveillance data and the spatiotemporally varying transmission pattern further complicate the development of appropriate vaccination for tropical regions. The tropics consist of the equatorial regions between the Tropic of Cancer (23.4°N) and the Tropic of Capricorn (23.4°S). The regions next to the tropics with latitude less than 40°N and greater than 40°S are the subtropics. Because in this analysis the same reasoning is applicable to both the tropics and the subtropics, for brevity henceforth the term tropics will include subtropics as well.
Unlike influenza transmission in the tropics that greatly varies geographically, the consistent wintertime influenza peak in temperate regions is often associated with, arguably, the corresponding dry and cold climate [7], [8]. It is not only a condition in which aerosol-borne influenza transmission is most favorable [9], [10], but also one that promotes indoor crowding tendency which may lead to higher risk for contact transmission [11]. In addition to temperature and humidity, the El Nino Southern Oscillation (ENSO) [12] and solar radiation [13] has been implicated in influenza transmission in temperate climate. On the other hand, the role of climate on influenza in the tropics is much less understood. Several regions observe high influenza transmission that coincides with rainy season such as southern India, Vietnam and Brazil [14], [15], [16], [17]. While others such as Singapore, Thailand and Philippines, detect semi annual peaks that are not necessarily associated with rainfall [14]. In vivo study by Lowen et al. [18] emphasized the dominating effect of contact transmission in the tropics, while Alonso et al. [5], on the contrary, showed that temperature and humidity contribute more to the southward influenza spread in Brazil than contact transmission. Thus it appears that the contributing factors to influenza transmission in the tropics are region-specific due to the highly varying transmission pattern.
With the scarcely-available surveillance data for the tropics, we chose to study influenza transmission pattern in two regions that, similar to the tropics, are characterized by warm climate, but have advanced influenza surveillance system. These regions are Hong Kong, China (22° N) and Maricopa County, Arizona (33° N)–with Hong Kong in the tropics and Maricopa County in the subtropics. The objective of the study was to investigate and model the effect of climate on the transmission pattern. The resulting model, in turn, can be employed to forecast influenza epidemics in the tropics that may help to facilitate vaccination strategy development and antiviral distribution. Furthermore, the use of climatic parameters in the forecast for tropics is highly advantageous not only due to few studies focusing on the influenza prediction–as demonstrated by Viboud et al. [19])–but also due to the low surveillance efforts in these regions.
Results
Throughout this study, we utilized a time series-based model namely Autoregressive Integrated Moving Average (ARIMA), as well as SARIMA when seasonality is included. The method is briefly delineated under Methods section. In the following, results for influenza Hong Kong will be presented first, followed by Maricopa County results.
Hong Kong
As a first step in ARIMA modeling, we seek to stationarize the response series, which is the influenza weekly count in Hong Kong shown in Figure 1. Taking the log transformation of the series reduced the variances of the influenza cases, and subsequent differencing–either first or seasonal order of difference–resulted in stationary series. We further used the Auto-Correlation and Partial Auto-Correlation Function (ACF and PACF, respectively) plots to identify the order of the ARIMA model for the stationary series. For the first-order differenced series, both the ACF and PACF cut off at lag 2, whereas for the seasonally differenced series the ACF decreases very slowly with PACF cut off at lag 2.
We further fit several univariate (S)ARIMA models with different orders, and consequently excluded any models in which the residuals exhibit autocorrelation. The resulting models and the estimated coefficients are summarized in the top half of Table 1 (Please see Methods section for ARIMA model notation). As we can see in Table 1, ARIMA(2,1,2) has the best Root Mean Square Error (RMSE) for the fitted dataset, ARIMA(1,1,2) has the best predictive RMSE, while SARIMA(2,0,0)(0,1,0) has the lowest Akaike Information Criterion (AIC). Among these four univariate models, the relative differences of the worst from the best performing AIC is 6%, fit RMSE is 27% and prediction RMSE is 16%. Since the AIC differences is relatively low, and the prediction RMSE not only have smaller differences but also it is unavailable in a realistic case, we will henceforth use the model with smallest fit RMSE–that is ARIMA(2,1,2)–as a baseline univariate model for further comparison.
In order to include the environmental parameter as input series to the model, we first examined the correlations between the influenza cases and the environmental series. Our results (Table 2) show that there are significant correlations (based on the two standard error limits) with LST at lag 2 and 5, rainfall at lag 3 and relative humidity at lag 0 to 3. Subsequently (S)ARIMA models were estimated with one or more environmental variables included. The performances of these models and the estimated coefficients are shown in Table 1. For these multivariate models, the best fit RMSE is obtained from SARIMAX(0,1,2)(1,0,0) with LST, rainfall and relative humidity as covariates. SARIMAX(1,0,0)(0,1,0) with LST has the best AIC, while ARIMAX(2,1,0) with relative humidity has the best prediction RMSE. Comparing these three models with the baseline univariate model discussed previously (ARIMA(2,1,2)), we found that including the environmental input series improve the AIC by 18%, the fit RMSE by 9% and the prediction RMSE by 16% from the baseline model.
Among the three best performing multivariate models, ARIMAX(2,1,0) with rainfall as input, has the highest AIC values. Moreover, it has high p-value for estimated rainfall coefficient in this model (Table 1). Thus we can exclude this model from the selection. Between SARIMAX(1,0,0)(0,1,0) with LST only and SARIMAX(0,1,2)(1,0,0) with LST, rainfall and relative humidity, the differences in AIC is small (2%), but there is about 27% differences in the fit RMSE. Thus we choose the SARIMAX(0,1,2)(1,0,0) with LST, rainfall and relative humidity since it has the lower fit RMSE. The fitted and predicted values of this model were plotted in Figure 2, and the associated environmental variables are shown in Figure 3.
SARIMAX(0,1,2)(1,0,0) with LST lag 2 and 5, accumulated rainfall lag 3 and relative humidity lag 1.
(A) weekly average of Land Surface Temperature (LST), (B) weekly accumulated rainfall and (C) weekly average of mean relative humidity.
Maricopa County
Figure 4 showed the time series profile for positive influenza counts in Maricopa County. Similar to Hong Kong influenza, taking the log transformation of the series reduced the variation in the variances. Stationarity of the log-transformed series was achieved through first order differencing. Further plotting the ACF and PACF of the stationary series reveals significant autocorrelations at lag 2 and 21, and significant partial autocorrelation at lag 2. Based on this ACF and PACF, we fitted several univariate (S)ARIMA models and found that the two best performing models are SARIMA(1,1,0)(1,0,0) and SARIMA(0,1,0)(1,0,0). As shown in Table 3, the performances of the univariate models measured by AIC and RMSE for the fitting dataset were similar. Thus we chose SARIMA(0,1,0)(1,0,0) as the baseline model because all p-values of the estimated coefficient are relatively significant (<0.05).
We further pre-whitened the environmental variables and calculated the cross-correlation function. Variables that exhibit significant cross-correlations are shown in Table 4. As we can see, influenza cases in Maricopa County are correlated with LST at lag 3, mean air temperature at lag 1 and 7, maximum relative humidity at lag 3 and 6, minimum relative humidity at lag 6, and maximum air pressure at lag 0 and 6. The time series profiles of these variables are illustrated in Figure 5.
(A) LST weekly average and Mean Temperature from ground station, (B) weekly average maximum pressure and (C) weekly average of minimum and maximum relative humidity.
We further fitted (S)ARIMAX model with the lagged environmental variables as input series, and the results are summarized in Table 3. Overall, including one or more input series improves the model performances as compared to the baseline univariate SARIMA model previously described. Note that we kept some of the input variables even though the estimated coefficients have p-values greater than 0.05. This is mainly due to the model residuals that exhibit autocorrelation when these variables are removed. Our results indicate that incorporating the mean air temperature (Tmean) at lag 1 and 7 yield the lowest AIC value for the fitting dataset. As compared with the baseline univariate SARIMA model previously described, the addition of Tmean improves the AIC by 39% and the RMSE by 6%. Lowest RMSE for both fitting and prediction dataset is obtained by SARIMAX model with maximum air pressure (PMax) at lag 0 and 6 as covariates. Compared with the baseline univariate SARIMA, the addition of PMax improve the fit and predicted RMSE by about 9% and 0.4%, respectively. We showed in Figure 6 the fitted and predicted values produced by SARIMA model with Tmean and Pmax.
Discussion
Through the use of ARIMA models, we explored the relationship between environmental variables and influenza cases in two regions characterized by warm climate, Hong Kong (China) and Maricopa County (Arizona, USA). We first examined whether influenza cases can be modeled as a univariate (S)ARIMA where it only depends on its own past values and random errors. We found that the univariate ARIMA was capable of forecasting 1-step ahead future influenza cases relatively well. For Hong Kong influenza, the best univariate model is ARIMA(2,1,2), where influenza cases depend on cases in previous two weeks. While for Maricopa County, the best univariate model is the seasonal SARIMA(0,1,0)(1,0,0), which depends on the cases in the previous one season. In the time series plot of influenza cases in Maricopa County (Figure 4), the seasonality is very distinct compared to influenza in Hong Kong (Figure 1). Thus it is expected that influenza cases in Maricopa County are modeled best when seasonality is incorporated.
In the multivariate ARIMA models, we found that the accumulated rainfall, land surface temperature (LST) and relative humidity are significant predictors for influenza in Hong Kong. The association of rainfall with influenza is commonly observed in tropical countries such as in Brazil [17], Singapore [15] and Thailand [16]. There is yet any direct relationship connecting rainfall with either the effectiveness of influenza transmission, virus survivorship or host susceptibility. Common observation is that rainfall may cause changes in the social behavior that in turn promotes contact transmission. For example, in rainy seasons, more people prefer indoor activities that may increase the chance for social contact, and hence contact transmission. In Hong Kong, rainy season is between April to September, with heavy and persistent rain typically in May and August. Meanwhile, influenza transmission typically peaks around March–April, and June–August (Figure 1), which in general is considered to be rainy season.
Temperature and relative humidity is often associated with influenza epidemics such as in Tokyo (Japan) [8] and especially in temperate regions where influenza peak coincides with winter. The prevailing dry and cold condition during winter seems to enhance influenza transmission, though this not substantiates high influenza transmission in the tropics. Lowens et al [9] conducted in vivo study and found that at low temperature (5°C) and low relative humidity (20% to 35%), transmission through aerosol is most efficient. At 20°C, aerosol transmission efficiency varies with relative humidity [9]. On the other hand, 30°C blocks aerosol transmission but not contact transmission, which explains influenza transmission in the tropics [18]. In our Hong Kong analysis, land surface temperature (LST) is used instead of air temperature, but LST can be used as a proxy for the other.
The first influenza peak in Hong Kong typically occurs During March–April, when the temperature starts rising. The normal mean temperature at this time is 18–22°C, whereas the relative humidity is 75%–80%. According to Lowens et al. [9], 20°C temperature and a relative humidity of 65% induce relatively higher transmission, whereas no transmission occurs at 80% relative humidity. Thus Hong Kong condition in March–April, lies between intermediate to no transmission zone as concluded by Lowens et al. [9]. However note that during this time of the year, rainfall starts to increase as well. As a result, the combination of moderate aerosol transmission efficiency as determined by temperature and relative humidity, and the behavioral response that promote contact transmission due to rainfall, causes high transmission rate during this time of the year. The second peak of the influenza transmission, on the other hand, occurs during summer where mean temperature is normally between 26–28°C and relative humidity is above 80%. At this condition, aerosol transmission is theoretically blocked, but the high rainfall frequency promotes the risk of contact transmission. Thus it seems that the second influenza peak is predominantly caused by the contact transmission.
Our results for influenza in Maricopa County showed that influenza cases are correlated with air and land surface temperature, maximum air pressure, and both maximum and minimum relative humidity (Table 3). As we have previously discussed, low temperature and relative humidity induce aerosol transmission. However, the best model is obtained when only maximum pressure is included as input series (Table 4). Air pressure is an important determinant of the weather, including temperature and precipitation. In general, low pressure brings clouds and precipitation, and vice versa. Although the best model uses maximum pressure as its input variable, it should be noted that the performances of other models are not distinctly worse (Table 4). For instance, the model with mean air temperature (Tmean) as an input variable shows not only fitted RMSE that is only 2% higher than the one with maximum pressure, but also a better AIC value. In this model, the influenza incidences are positively associated with the mean air temperature. Whereas the model with minimum or maximum relative humidity shows that influenza is positively associated with humidity (both minimum and maximum). These results are in agreement with the findings that dry and cold condition enhances influenza transmission.
Both Hong Kong and Maricopa County have comparable temperatures, where it is generally warm throughout the year. Nevertheless, Maricopa County is characterized more by desert climate, where it is drier and hotter (especially in the summer) than Hong Kong's subtropical climate where rainfall is much more frequent. Therefore, the driving factors for influenza transmission may not be the same as it is reflected in our analysis. Our results show that Hong Kong influenza cases are associated with rainfall (Table 3), whereas it is not the case in Maricopa County. We had included rainfall in the model for Maricopa County but the results were unsatisfactory (not shown). The association of rainfall with influenza only in Hong Kong corroborates the finding that contact transmission is more predominant in the tropics [18], with a reasoning that rainfall induces social behavior that promotes contact transmission. Other environmental factors influencing influenza transmission seems to be common for both regions–temperature and relative humidity, although measured in different indicators for both regions–which enhance virus survivorship and aerosol transmission.
In this paper we have demonstrated the use of environmental/meteorological variables–as obtained from the satellite and ground stations–and the influenza biosurveillance data through a mathematical model, to assess the factors associated with influenza incidences. We have shown the prediction capability of the models, as measured by the RMSE of prediction dataset (Table 1 and Table 3), to forecasts the next influenza season. Presently we use the one-step ahead forecast in calculating the future influenza cases. In reality this may only be possible for cities with more advanced computer-based surveillance systems such as New York City and Hong Kong. Most of the models developed here depend on the past one to two weeks influenza cases. A more realistic approach is to predict the influenza cases using more than one-step ahead forecasts. This means that future forecasts are calculated using previously predicted number of cases instead of using the actual cases from the surveillance data (as in one-step ahead approach). However, one caveat to this approach is that more data is needed, since model selection will be based not only on the RMSE of the fitting dataset but also on the prediction dataset. When the model incorporates seasonality, the number of data that can be used for fitting process decreases significantly, especially in this study where the seasonal period is 52 weeks. In this approach one would divide the data into three: (i) for fitting process, where the coefficients are estimated, (ii) forecasting process, where future values are calculated using the predicted values, and the goodness of fit statistics will typically be used in the model selection and (iii) for validation process, to ensure that the model does not behave erratically. Thus the models in this study are a first step towards developing an early warning system for influenza.
Materials and Methods
Materials/Data
This study uses weekly count of laboratory-confirmed influenza viruses in two regions, Hong Kong and Maricopa County. We obtained Hong Kong influenza count between January 2005 to September 2008, from weekly influenza report as published by the Department of Health, Government of the Hong Kong Special Administrative [20]. Maricopa County influenza data, spanning October 2004–March 2009, were obtained from the weekly report by the Maricopa Department of Public Health [21]. Maricopa County flu reports are available during the flu season that usually begins from week 40 and ends at week 17 of the following year. We consequently assumed that flu counts outside the season to be zero.
Climatic and meteorological parameters were collected from two primary sources: ground-based and satellite-derived measurements. From the Hong Kong Observatory [22] we retrieved daily meteorological observations including temperature (maximum, mean, minimum), mean dew point temperature, mean relative humidity, global solar radiation and total evaporation. Maricopa County daily climatic observations were acquired from The Flood Control District of Maricopa County [23], where we aggregated data from 32 stations. This data includes daily mean air temperature, dew point (minimum, mean and max), relative humidity (minimum and maximum), maximum wind speed, air pressure (minimum and maximum) as well as maximum solar radiation.
Furthermore, we obtained remotely-sensed daily rainfall measurements for both Hong Kong and Maricopa County, from the instruments embarked on the Tropical Rainfall Measuring Mission (TRMM) [24], [25]. TRMM is a collaborative mission between NASA and Japan Aerospace Exploration Agency. The goal of the mission is to monitor and study tropical rainfall in order to improve the understanding of the water cycle in the climate system. Of the five instruments onboard the spacecraft, the Precipitation Radar and the TRMM Microwave Imager have the most direct relationship to measuring precipitation. As a follow-on mission for TRMM, the Global Precipitation Mission [26], [27], which is a multinational mission involving a constellation of satellites, will be launched in 2013. The TRMM data was retrieved using the GES-DISC Interactive Online Visualization ANd aNalysis Infractructure (Giovanni) as part of the NASA's Goddard Earth Sciences (GES) Data and Information and Services Center (DISC) [28]. In addition to TRMM data as rainfall measurements, we extracted Daily Land Surface Temperature (LST) from the MODerate resolution Imaging Spectroradiometer (MODIS) data set [29]. Both Terra and Aqua missions of NASA's Earth Observing System carry this instrument. The instrument has 36 bands spanning from the visible to the long wave infrared spectra, from which information related to disease transmission can be extracted. The importance of temperature and precipitation–as environmental determinants in the infectious disease transmission–have rendered TRMM and MODIS significant roles in remote sensing-based disease surveillance which has been demonstrated elsewhere [30], [31].
Analysis
The weekly influenza count data was divided into two: one set was used in the fitting process (parameter estimation), and another for prediction. We took the observations in the latest 1 year as the prediction period. Out of 194 observations in Hong Kong influenza data, we used 155 points for fitting and 39 for prediction. For Maricopa data, the fitting set consists of 159 observations and the prediction set has 75.
The influenza time series that we analyzed in this study is characterized by a strong autocorrelation, a property that commonly violates the ordinary linear regression. Thus in order to account for the autocorrelation behavior, we employed a class of time series technique namely Auto Regressive Integrated Moving Average (ARIMA) [32], [33]. We first developed a univariate ARIMA model, where the response series depends only on its past values and some random shocks, followed with multivariate ARIMA with the environmental parameters as input series/covariates. In the following we will briefly delineate the Box-Jenkins approach for ARIMA modeling that is used throughout this study.
ARIMA is based on the assumption that the response series is stationary, that is the mean and variances of the series are independent of time. Stationarity can be achieved by differencing the series, or transforming the variable so as to stabilize the variance or mean. In our analysis we took the logarithmic transformation to reduce the variances of the influenza time series, and subsequently differenced the series until it is stationary. Once the response series is stationary, we examined the Autocorrelation function (ACF) and Partial Autocorrelation Function (PACF) to determine the initial autoregressive (AR) and moving average (MA) order. An ARIMA model is notated as ARIMA(p,d,q), where p indicates the AR order, d the differencing order and q the MA order. An ARIMA model that incorporates seasonality is referred as SARIMA(p,d,q)(P,D,Q) where P,D and Q indicate the seasonal order of AR, differencing, and MA, respectively. Since the influenza time series is recorded as weekly observations, the seasonality period is 52.
Based on the ACF and PACF we fitted several ARIMA models with varying AR and MA orders. In the fitting process, the AR and MA coefficients were estimated using conditional least square method. The residuals were further inspected for autocorrelation through ACF and PACF. Models with autocorrelated residuals were discarded, else goodness of fit were examined through calculated Akaike's Information Criterion (AIC) and the Mean Square Error (RMSE). The resulting model was subsequently used to forecast (1 step ahead) the latest influenza season, and the associated RMSE were calculated.
Once we developed and selected a univariate ARIMA, we investigated the effect of the environmental variables and the corresponding lags on the influenza cases. Note that the daily environmental data were converted into weekly resolution by taking the average. The environmental series were first pre-whitened. In other words, we applied univariate ARIMA modeling, as previously discussed, such that the environmental series no longer characterized by autocorrelation. Subsequently, cross-correlations function (CCF) between the pre-whitened environmental series and the influenza cases was then calculated so as to identify the lags to be included in the model. The significance of the cross-correlations was assessed on the basis of its two standard error limits (significant at 0.05 level). Environmental variables that did not exhibit significant cross-correlations with the influenza cases were excluded from further analysis.
Similar to univariate ARIMA fitting process, we further estimated the coefficients of the AR and MA terms as well as the lagged environmental variable. In this study, (S)ARIMA model that incorporates environmental input series is referred as (S)ARIMAX. The environmental input series were first included one at a time before combining them together. Estimated coefficients with p-values greater than 0.05 are excluded when possible and the model was re-fitted.
All ARIMA modeling and the corresponding statistical tests were performed using SAS software, Version 9.1.2 of the SAS System for Windows (SAS Institute, Inc., Cary, NC).
Author Contributions
Conceived and designed the experiments: RPS RKK. Performed the experiments: RPS. Analyzed the data: RPS FA RKK. Contributed reagents/materials/analysis tools: FA RKK. Wrote the paper: RPS RKK.
Reference
- 1.
World Health Organizations (WHO) (2009) Influenza (Seasonal). Available: http://www.who.int/mediacentre/factsheets/fs211/en/. Date Accessed: January 2009.
- 2. Wong CM, Yang L, Chan KP, Leung GM, Chan KH, et al. (2006) Influenza-associated hospitalization in a subtropical city. PLoS Med 3: e121.
- 3. Rambaut A, Pybus OG, Nelson MI, Viboud C, Taubenberger JK, et al. (2008) The genomic and epidemiological dynamics of human influenza A virus. Nature 453: 615–619.
- 4. Russell CA, Jones TC, Barr IG, Cox NJ, Garten RJ, et al. (2008) Influenza vaccine strain selection and recent studies on the global migration of seasonal influenza viruses. Vaccine 26: Suppl 4D31–34.
- 5. Alonso WJ, Viboud C, Simonsen L, Hirano EW, Daufenbach LZ, et al. (2007) Seasonality of influenza in Brazil: a traveling wave from the Amazon to the subtropics. Am J Epidemiol 165: 1434–1442.
- 6. Viboud C, Alonso WJ, Simonsen L (2006) Influenza in tropical regions. PLoS Med 3: e89.
- 7. du Prel JB, Puppe W, Grondahl B, Knuf M, Weigl JA, et al. (2009) Are meteorological parameters associated with acute respiratory tract infections? Clin Infect Dis 49: 861–868.
- 8. Urashima M, Shindo N, Okabe N (2003) A seasonal model to simulate influenza oscillation in Tokyo. Jpn J Infect Dis 56: 43–47.
- 9. Lowen AC, Mubareka S, Steel J, Palese P (2007) Influenza virus transmission is dependent on relative humidity. PLoS Pathogens 3: e151.
- 10. Shaman J, Kohn M (2009) Absolute humidity modulates influenza survival, transmission, and seasonality. Proc Natl Acad Sci U S A 106: 3243–3248.
- 11. Lofgren E, Fefferman NH, Naumov YN, Gorski J, Naumova EN (2007) Influenza seasonality: underlying causes and modeling theories. J Virol 81: 5429–5436.
- 12. Viboud C, Pakdaman K, Boelle PY, Wilson ML, Myers MF, et al. (2004) Association of influenza epidemics with global climate variability. Eur J Epidemiol 19: 1055–1059.
- 13. Sagripanti JL, Lytle CD (2007) Inactivation of influenza virus by solar radiation. Photochem Photobiol 83: 1278–1282.
- 14. Hampson AW (1999) Epidemiological data on influenza in Asian countries. Vaccine 17: S19–S23.
- 15. Chew FT, Doraisingham S, Ling AE, Kumarasinghe G, Lee BW (1998) Seasonal trends of viral respiratory tract infections in the tropics. Epidemiol Infect 121: 121–128.
- 16. Chumkiew S, Srisang W, Jaroensutasinee M, Jaroensutasinee K (2007) Climatic Factosr Affecting on Influenza Cases in Nakhon Si Thammarat. Proceedings of world academy of science, engineering and technology 21: 364–367.
- 17. Moura FE, Perdigao AC, Siqueira MM (2009) Seasonality of influenza in the tropics: a distinct pattern in northeastern Brazil. Am J Trop Med Hyg 81: 180–183.
- 18. Lowen AC, Steel J, Mubareka S, Palese P (2008) High Temperature (30 C) blocks aerosol but not contact transmission of influenza virus. Journal of Virology 82: 5650–5652.
- 19. Viboud C, Boelle PY, Carrat F, Valleron AJ, Flahault A (2003) Prediction of the spread of influenza epidemics by the method of analogues. Am J Epidemiol 158: 996–1006.
- 20.
Hong Kong Department of Health Flu Express. Available: http://www.chp.gov.hk/epidemiology.asp?lang=en&id=445&pid=304&ppid=441. Date Accessed: January 2009.
- 21.
Maricopa Department of Public Health Maricopa County Weekly Influenza Summary. Available: http://www.maricopa.gov/Public_Health/epi/flu.aspx. Date Accessed: April 2009.
- 22.
Hong Kong Observatory Extract of Meteorological Observations for Hong Kong Available: http://www.hko.gov.hk/wxinfo/pastwx/extract.htm. Date Accessed: January 2009.
- 23.
The Flood Control District of Maricopa (2009) Weather Information. Available: http://www.fcd.maricopa.gov/Rainfall/Weather/weather.aspx. Date Accessed: July 2009.
- 24.
NASA Goddard Space Flight Center (2009) Tropical Rainfall Measuring Mission. Available: Date Accessed: August 24, 2009.
- 25. Kummerow C, Barnes W, Kozu T, Shiue J, Simpson J (1998) The Tropical Rainfall Measuring Mission (TRMM) Sensor Package. Journal of Atmospheric and Oceanic Technology 15: 809–817.
- 26.
NASA Goddard Space Flight Center (2009) Global Precipitation Measurement. Available: http://gpm.gsfc.nasa.gov. Date Accessed: August 24, 2009.
- 27.
Smith EA, Asrar G, Furuhama Y, A G, Adlre R, et al. (2006) International Global Precipitation Measurement (GPM) Program and Mission: An Overview. In: Levizzani V, Bauer P, Turk FJ, editors. Measuring Precipitation from Space: EURAINSAT and the Future: Kluwer Publisher.
- 28.
NASA Goddard Earth Sciences Data and Information Service Center Giovanni (GES-DISC Interactive Online Visualization and Analysis Infrastructure). Available: http://disc.sci.gsfc.nasa.gov/giovanni/overview/index.html. Date Accessed: March 2009.
- 29.
NASA Goddard Space Flight Center (2009) MODIS. Available: http://modis.gsfc.nasa.gov. Date Accessed: January 2009.
- 30. Kiang RK, Adimi F, Soika V (2006) Meteorological, environmental remote sensing and neural network analysis of the epidemiology of malaria transmission in Thailand. Geospatial Health 1: 71–84.
- 31. Xiao X, Gilbert M, Slingenbergh J, Lei F, Boles S (2007) Remote sensing, ecological variables, and wild bird migration related to outbreaks of highly pathogenic H5N1 avian influenza. Journal of Wildlife Diseases 43:
- 32. Chatfield C (2000) Time Series Forecasting: Chapman & Hall/CRC.
- 33. Box GEP, Jenkins , GM , Reinsel GC (2009) Time Series Analysis: Forecasting and Control. 4th Ed. Wiley.