Effect modification of environmental factors on influenza-associated mortality: a time-series study in two Chinese cities
© Yang et al; licensee BioMed Central Ltd. 2011
Received: 25 May 2011
Accepted: 14 December 2011
Published: 14 December 2011
Environmental factors have been associated with transmission and survival of influenza viruses but no studies have ever explored the role of environmental factors on severity of influenza infection.
We applied a Poisson regression model to the mortality data of two Chinese metropolitan cities located within the subtropical zone, to calculate the influenza associated excess mortality risks during the periods with different levels of temperature and humidity.
The results showed that high absolute humidity (measured by vapor pressure) was significantly (p < 0.05) associated with increased risks of all-cause and cardiorespiratory deaths, but not with increased risks of pneumonia and influenza deaths. The association between absolute humidity and mortality risks was found consistent among the two cities. An increasing pattern of influenza associated mortality risks was also found across the strata of low to high relative humidity, but the results were less consistent for temperature.
These findings highlight the need for people with chronic cardiovascular and respiratory diseases to take extra caution against influenza during hot and humid days in the subtropics and tropics.
Influenza used to be considered as a "cold" disease as it usually returns every cold winter in temperate countries. However, recent studies have shown that influenza can be active throughout the year in the warm tropics and subtropics and the disease burden of influenza there can be as heavy as that in temperate climates [1–4]. It has been proposed that influenza seasonality is driven by complicated interactions between antigenic drifts of virus strains, environmental factors, host susceptibility and behavior changes [5–7]. Among these potential factors, environmental factors including temperature and relative humidity have been most thoroughly explored by laboratory and observational studies [8, 9]. Recent studies raised a hypothesis that absolute humidity is one of drivers for influenza seasonality in temperate regions [10, 11], but the mechanism behind various seasonal patterns of influenza outbreaks under different climates remains unclear.
Influenza surveillance data were obtained from the Hong Kong Department of Health (DH) for the period of 1998-2006 in Hong Kong and the Guang Dong Provincial Centre for Disease Control and Prevention (GDCDC) for 2004-2006 in Guangzhou . In both cities the consultation rates for influenza-like illness (ILI) from sentinel hospitals or doctors were reported on the weekly basis, with diagnosis criteria of fever (≥38°C), cough or sore throat. Nasopharyngeal aspirate or swab was collected from ILI patients and then tested of influenza viruses by immunofluorescence and cell culture. Positive specimens were further typed (subtyped) into A/H3N2, A/H1N1 and B by haemagglutination inhibition test with strain-specific antisera provided by the World Health Organization. During 2004-2006, annual averages of 35,256 and 1,872 specimens were tested by the DH and GDCDC, respectively.
Vapor pressure will be used for absolute humidity hereafter in this paper.
Here y t denotes the number of deaths at week t. ns(t), ns(temp t ), and ns(humd t ) denote the natural cubic spline smoothing functions of time, weekly average temperature and relative humidity. The natural spline smoothing function was used to remove the small variation while maintaining the major trend of each variable, i.e. to make them smoother. The aim of smoothing was to increase the efficiency in estimating the model coefficients . Given the high correlation between temperature and vapor pressure, the smoothing function of vapor pressure was not added in order to avoid collinearity between the variables. The adequacy of this core model was evaluated by the absence of any obvious pattern in partial autocorrelation functions of its residuals (Additional file 1: Figure S1). Second, the weekly proportions of specimens positive for influenza A or B were then added into the core model as a variable for influenza virus activity to obtain a main effect model.
Here β1 was obtained from the above interaction model. To test whether our results were robust to various definitions of periods, we chose two sets of extra cutoff points: 20th and 80th, 30th and 70th percentiles of weekly average data in each city. Although we removed the autocorrelation and seasonal trends within mortality data, there are still concerns that such adjustment was inadequate and the remaining uncontrolled seasonal factors may cause interaction terms of environmental factors and virus activity to appear significant in our models. To rule out this possibility, we used accidental deaths that were expected to be unrelated to influenza infection as a control mortality group to show that our findings are unlikely the spurious results of under-adjustment of seasonal confounding factors in modeling. Some previous studies used the anomalies to assess the effects of meteorology factors on influenza .We therefore conducted a sensitivity analysis by defining the strata by anomalies, instead of absolute values of meteorology factors. The anomalies were defined as the deviations of observed metrological data from a seasonal curve with a constant and a sinusoidal pair fitted to these observed data. All the analyses were performed using the mgcv package of R software (version 2.5.1.) .
Descriptive statistics for weekly average temperature, relative humidity, vapor pressure and counts of all-cause deaths
Relative humidity (%)
All-cause deaths (< 65 years)
All-cause deaths (≥65 years)
Weekly mean of mortality numbers and proportions of specimens positive for influenza (virus %) during the periods with low, middle and high levels of environmental factors
Percentage change (%) of mortality counts associated with 1% increase in influenza virus activity during the low, middle and high periods of temperature, relative humidity and vapor pressure
Percentage change (95% CI)
0.02 (-0.50, 0.54)
0.29 (0.07, 0.50)
0.26 (-0.11, 0.63)
0.11 (-0.08, 0.30)
0.15 (0.01, 0.28)
0.24 (0.06, 0.41)
0.15 (-0.56, 0.86)
0.41 (0.11, 0.72)
0.50 (-0.03, 1.03)
0.13 (-0.14, 0.39)
0.27 (0.08, 0.47)
0.49 (0.24, 0.74)
0.14 (-1.71, 2.03)
1.37 (0.64, 2.09)
0.82 (-0.41, 2.06)
0.38 (-0.11, 0.86)
0.66 (0.33, 1.00)
0.79 (0.36, 1.23)
0.32 (-0.34, 0.98)
0.20 (-0.24, 0.63)
0.11 (-0.44, 0.67)
-0.60 (-1.50, 0.32)
0.13 (-0.13, 0.40)
0.44 (0.20, 0.68)
0.04 (-0.26, 0.33)
0.13 (0.01, 0.25)
0.29 (0.12, 0.46)
-0.79 (-2.04, 0.48)
0.32 (-0.05, 0.70)
0.57 (0.23, 0.91)
0.22 (-0.18, 0.62)
0.27 (0.09, 0.44)
0.41 (0.16, 0.66)
-1.67 (-4.95, 1.74)
0.89 (-0.02, 1.81)
1.66 (0.85, 2.48)
0.52 (-0.22, 1.26)
0.60 (0.29, 0.91)
0.76 (0.34, 1.18)
0.91 (0.00, 1.83)
0.13 (-0.27, 0.53)
0.26 (-0.29, 0.81)
-0.20 (-0.75, 0.35)
0.27 (0.04, 0.51)
0.35 (0.07, 0.62)
0.13 (-0.07, 0.32)
0.11 (-0.03, 0.25)
0.26 (0.09, 0.42)
-0.06 (-0.81, 0.70)
0.38 (0.05, 0.71)
0.54 (0.16, 0.93)
0.14 (-0.13, 0.41)
0.24 (0.04, 0.44)
0.49 (0.25, 0.72)
-0.05 (-2.03, 1.98)
1.04 (0.23, 1.86)
1.58 (0.73, 2.44)
0.45 (-0.04, 0.93)
0.71 (0.36, 1.06)
0.67 (0.26, 1.08)
0.49 (-0.16, 1.14)
0.23 (-0.22, 0.68)
-0.02 (-0.54, 0.50)
The interaction between vapor pressure levels and virus activity were found to be significant (p < 0.05) for all-cause and cardiorespiratory mortality, but not for pneumonia and influenza mortality (Table 3). Consistently higher mortality risks were found at the high levels of vapor pressure, with only exception of P&I mortality in Hong Kong. The influenza effects on mortality of cardiorespiratory, pneumonia and influenza, and all causes were found significant during the middle- and high-vapor pressure periods (with the only exception of all-cause mortality in Hong Kong), but not significant when vapor pressures remained at the relatively low levels. For the high vapor pressure periods, all-cause excess mortality counts attributable to influenza would increase by 0.35% and 0.26% for per 1% increase of virus activity, and the corresponding increases for cardiorespiratory mortality were 0.54% and 0.49% for Guangzhou and Hong Kong, respectively (Table 3).
The models with indicators for different cutoff points returned similar estimates for temperature, relative humidity and vapor pressure (Additional file 1: Figure S2). We did not observe any significant interaction between influenza and environmental factors in terms of their effects on the control mortality category of accidental mortality in Hong Kong. Influenza associated accidental mortality risks were also not statistically significant.
To ensure the standardized comparison between Guangzhou and Hong Kong, we decided to apply the same modeling approach to the data of same study period, because the core model would have been slightly different if we used the longer time series of Hong Kong. To check the robustness of our conclusions, we repeated the above analysis using a longer time series data of Hong Kong during 1998-2006 and the estimates were shown in Additional file 1: Table S1. The statistical significance of interaction terms was close to those from the study period of 2004-2006. The estimates for the low temperature periods became larger and statistically significant; those for the high relative humidity periods were smaller and comparable to the middle relative humidity; and for the low vapor pressure periods, the estimates were similar but became significant. The increasing trend across the low, middle and high levels of temperature and relative humidity was less evident.
The results of stratification analysis by anomalies are shown in Additional file 1: Table S2. The estimates were similar to those for the periods defined by absolute values of meteorological factors, in terms of magnitude and changing patterns. But the likelihood ratio tests showed more significant interaction for temperature or vapor pressure, and less for relative humidity.
In this study we quantified influenza associated mortality risks at the various ranges of temperature and humidity, and compared the results between two large subtropical cities. An increasing pattern of influenza-associated mortality risks on all-cause and cardiorespiratory along the low to high periods was found for temperature, relative and absolute humidity in both cities during the period of 2004-2006. The interaction between vapor pressure indicators and virus activity were also consistently significant for relative and absolute humidity. Although lower vapor pressure (or temperature) has been found to facilitate virus transmission and survival in the guinea pig model , our results suggested that higher vapor pressure (or temperature) was associated with a higher mortality risk attributable to influenza. Severity of seasonal influenza epidemics were not only determined by virus transmission efficiency and outdoor weather conditions, but also largely affected by host resistance, indoor living environment and social behavior . The high vapor pressure periods coincided with the summer peaks of influenza in both Guangzhou and Hong Kong, indicating that influenza could pose a higher risk when reaching its peak in seasons with high vapor pressure in subtropical cities.
Our results may help to interpret our previous findings that the effects of influenza were significantly higher in the humid and warm spring/summer period than in the dry and cold winter period in Hong Kong . The extreme low vapor pressure was usually recorded during December-January and the highest appeared during June-July, which coincided with the trough and peak periods of excess risks associated with influenza viruses. Effect modification of environmental factors was only detected in all-cause and cardiorespiratory mortality, but not in pneumonia and influenza, suggesting that the synergistic interaction between high humidity (or temperature) and virus activity may mainly lie in their similar regulation pathways in cardiovascular systems. These results are also in agreement with our previous findings that pneumonia and influenza mortality risks attributable to influenza did not exhibit a seasonal variation . Extreme heat has been documented to increase blood viscosity through evaporation of body fluid and trigger intravascular coagulation through damaging endothelial cells . Influenza infection has a similar pro-thrombotic effect by inducing inflammation around blood vessels and rupturing atherosclerotic plaques . As a result, mortality risks would be dramatically raised by the stress of both extreme weather and influenza infections. The results suggested that we need extra precautionary measures to reduce influenza infections in people with cardiovascular diseases, especially under the frequent hot and humid weather conditions experienced in the tropical and subtropical areas. Although the experiments of influenza virus transmission between guinea pig hosts found the higher transmission rates occurred under dry air (vapor pressure below 10hPA) , our results indicated that the mortality risks associated with influenza under the low vapor pressure environment were lower than the rest of study period. The reason could be the short time of exposure to the very low level of vapor pressure. During our study period there were only 2 and 13 weeks with an average vapor pressure below 10hPA in Hong Kong and Guangzhou, respectively. Most weekly average vapor pressures during the low vapor pressure period were within the range of 10-20hPA, in which the guinea pig experiments showed dramatically reduced transmission rates . In future, we may examine the seasonal variation in influenza effects in other cities to assess whether such a seasonal variation, if common in subtropical and tropical cities, is consistently determined by environmental factors, or by other factors such as host immunity and virus virulence.
The results for the different age groups suggested that the modification effects of environmental factors may mainly lie in the elderly aged over 65 years, as the consistent increasing trend over the low to high vapor pressure periods was only observed in this age group. However, since over 65% of deaths occurred in this age group for both cities (Table 1), the small numbers of weekly death counts in the younger age group (< 65 years) may not have had enough power to allow assessment of effect modification based on the data over the 3 years. A future study with a long study period or a large population may help answer whether young people maybe also expose to higher mortality risks during the hot and humid days.
In this study, we used the quartiles of weekly data in each city, to separately define the periods with normal (middle) and extreme (low and high) weather. Given the difference in weather conditions between these cities, we think that it is not appropriate to use the same cutoff points for temperature or humidity to compare their modification effects on influenza associated mortality, as people living in hot subtropical and tropical regions may adapt well to the year-round hot and humid climate and have a higher threshold for adverse effects of weather. For example, although it is widely accepted that the temperature effects on mortality exhibited a U- or V-shape curve in both temperate and tropical/subtropical areas, the turning point of this curve varied across different cities. A study conducted in 11 cities of the US found that the turning point of temperature for its effects on mortality could range from 18.4°C to 32.4°C . To our best knowledge, so far there are no studies that have ever assessed the effect modification of temperature on influenza effects. Therefore, the commonly adopted cutoff points of city-specific quartiles seem appropriate at this stage [21, 27].
There are several limitations in our study. Firstly, our study is based on 3 years of surveillance data which may not have enough power to allow assessment of exposure-response curves for the effects of environmental factors. Nevertheless, our findings did suggest an increasing trend of influenza associated mortality risks across the periods of low, middle and high vapor pressure, although such findings may be applicable only to the warm climates. Secondly, we only investigated the effect modification of environmental factors through a simple interaction model, but there were other unadjusted factors, including host susceptibility and virulence of influenza strains. These factors are unlikely to work independently with environmental factors. Other environmental factors such as ultraviolet radiation , rainfall  have been proposed to play a role in the regulation of influenza seasonality, although evidence is rather limited compared with the three factors we chose to investigate . Lastly, we did not adjust for the vaccination rate in our model. In 2003, vaccination rate was 191 doses per 1,000 total population in Hong Kong , slightly higher than the rate of 129 doses/1,000 total population in Guangzhou . However, it is not clear when people received vaccination; therefore we were unable to assess the role of vaccination in our study.
This study provides a piece of key evidence to the effect of environmental factors on severity of seasonal influenza under warm climates and helps reveal the mechanism behind global influenza seasonality. It also highlights the need for people with chronic cardiovascular and respiratory conditions to take extra caution against influenza during the hot and humid days in the subtropics.
We thank Peter Chapman for making the map and proofreading the manuscript. This study was supported by the Research Fund for the Control of Infectious Diseases/Health and Health Services Research Fund ; and the Area of Excellence Scheme of the University Grants Committee of the Hong Kong Special Administrative Region [AoE/M-12/06].
- Chow A, Ma S, Ling AE, Chew SK: Influenza-associated deaths in tropical Singapore. Emerg Infect Dis. 2006, 12: 114-121.View ArticlePubMedPubMed CentralGoogle Scholar
- Wong CM, Chan KP, Hedley AJ, Peiris JSM: Influenza-associated mortality in Hong Kong. Clin Infect Dis. 2004, 39: 1611-1617. 10.1086/425315.View ArticlePubMedGoogle Scholar
- Wong CM, Yang L, Chan KP, Leung GM, Chan KH, Guan Y, et al: Influenza-associated hospitalization in a subtropical city. PLoS Med. 2006, 3: e121-10.1371/journal.pmed.0030121.View ArticlePubMedPubMed CentralGoogle Scholar
- Viboud C, Alonso WJ, Simonsen L: Influenza in tropical regions. PLoS Med. 2006, 3: e89-10.1371/journal.pmed.0030089.View ArticlePubMedPubMed CentralGoogle Scholar
- Dowell SF: Seasonal variation in host susceptibility and cycles of certain infectious diseases. Emerg Infect Dis. 2001, 7: 369-374.View ArticlePubMedPubMed CentralGoogle Scholar
- Dowell SF, Ho MS: Seasonality of infectious diseases and severe acute respiratory syndrome-what we don't know can hurt us. Lancet Infect Dis. 2004, 4: 704-708. 10.1016/S1473-3099(04)01177-6.View ArticlePubMedGoogle Scholar
- Lofgren E, Fefferman NH, Naumov YN, Gorski J, Naumova EN: Influenza seasonality: underlying causes and modeling theories. J Virol. 2007, 81: 5429-5436. 10.1128/JVI.01680-06.View ArticlePubMedGoogle Scholar
- Schulman JL, Kilbourne ED: Airborne transmission of influenza virus infection in mice. Nature. 1962, 195: 1129-1130. 10.1038/1951129a0.View ArticlePubMedGoogle Scholar
- Lowen AC, Mubareka S, Steel J, Palese P: Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathog. 2007, 3: 1470-1476.View ArticlePubMedGoogle Scholar
- Shaman J, Kohn M: Absolute humidity modulates influenza survival, transmission, and seasonality. Proc Natl Acad Sci USA. 2009, 106: 3243-3248. 10.1073/pnas.0806852106.View ArticlePubMedPubMed CentralGoogle Scholar
- Shaman J, Pitzer VE, Viboud C, Grenfell BT, Lipsitch M: Absolute Humidity and the Seasonal Onset of Influenza in the Continental United States. PLoS Biol. 2010, 8: e1000316-10.1371/journal.pbio.1000316.View ArticlePubMedPubMed CentralGoogle Scholar
- Warren-Gash C, Bhaskaran K, Hayward A, Leung GM, Lo SV, Wong CM, et al: Circulating influenza virus, climatic factors, and acute myocardial infarction: a time series study in England and wales and Hong Kong. J Infect Dis. 2011, 203: 1710-1718. 10.1093/infdis/jir171.View ArticlePubMedPubMed CentralGoogle Scholar
- Braga AL, Zanobetti A, Schwartz J: The effect of weather on respiratory and cardiovascular deaths in 12 U.S. cities. Environ Health Perspect. 2002, 110: 859-863.View ArticlePubMedPubMed CentralGoogle Scholar
- Hajat S, Kovats RS, Lachowycz K: Heat-related and cold-related deaths in England and Wales: who is at risk?. Occup Environ Med. 2007, 64: 93-100.View ArticlePubMedGoogle Scholar
- Baccini M, Biggeri A, Accetta G, Kosatsky T, Katsouyanni K, Analitis A, et al: Heat effects on mortality in 15 European cities. Epidemiology. 2008, 19: 711-719. 10.1097/EDE.0b013e318176bfcd.View ArticlePubMedGoogle Scholar
- Yang L, Wong CM, Chan KP, Chau PY, Ou CQ, Chan KH, et al: Seasonal effects of influenza on mortality in a subtropical city. BMC Infect Dis. 2009, 9: 133-10.1186/1471-2334-9-133.View ArticlePubMedPubMed CentralGoogle Scholar
- Center for Health Protection: Sentinel surveillance. [http://www.chp.gov.hk]
- Basu R, Feng WY, Ostro BD: Characterizing temperature and mortality in nine California counties. Epidemiology. 2008, 19: 138-145. 10.1097/EDE.0b013e31815c1da7.View ArticlePubMedGoogle Scholar
- Yang L, Chiu SS, Chan KP, Chan KH, Wong HS, Peiris JSM, et al: Validation of statistical models for estimating hospitalization associated with influenza and other respiratory viruses. PLoS One. 2011, 6: e17882-10.1371/journal.pone.0017882.View ArticlePubMedPubMed CentralGoogle Scholar
- Hastie TJ, Tibshirani RJ: Generalized Additive Models. 1990, London: Chapman and HallGoogle Scholar
- Yang L, Ma S, Chen PY, He JF, Chan KP, Chow A, et al: Influenza associated mortality in the subtropics and tropics: results from three Asian cities. Vaccine. 2011Google Scholar
- Ren C, Williams GM, Mengersen K, Morawska L, Tong S: Does temperature modify short-term effects of ozone on total mortality in 60 large eastern US communities?--an assessment using the NMMAPS data. Environ Int. 2008, 34: 451-458. 10.1016/j.envint.2007.10.001.View ArticlePubMedGoogle Scholar
- R Development Core Team. R: A Language and Environment for Statistical Computing. [2.2.1]. 2005, Vienna: R Foundation for Statistical ComputingGoogle Scholar
- Bouchama A, Knochel JP: Heat stroke. N Engl J Med. 2002, 346: 1978-1988. 10.1056/NEJMra011089.View ArticlePubMedGoogle Scholar
- Madjid M, Aboshady I, Awan I, Litovsky S, Casscells SW: Influenza and cardiovascular disease: is there a causal relationship?. Tex Heart Inst J. 2004, 31: 4-13.PubMedPubMed CentralGoogle Scholar
- Curriero FC, Heiner KS, Samet JM, Zeger SL, Strug L, Patz JA: Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol. 2002, 155: 80-87. 10.1093/aje/155.1.80.View ArticlePubMedGoogle Scholar
- Medina-Ramon M, Zanobetti A, Schwartz J: The effect of ozone and PM10 on hospital admissions for pneumonia and chronic obstructive pulmonary disease: a national multicity study. Am J Epidemiol. 2006, 163: 579-588. 10.1093/aje/kwj078.View ArticlePubMedGoogle Scholar
- Sagripanti JL, Lytle CD: Inactivation of influenza virus by solar radiation. Photochem Photobiol. 2007, 83: 1278-1282. 10.1111/j.1751-1097.2007.00177.x.View ArticlePubMedGoogle Scholar
- Moura FE, Perdigao AC, Siqueira MM: Seasonality of influenza in the tropics: a distinct pattern in northeastern Brazil. Am J Trop Med Hyg. 2009, 81: 180-183.PubMedGoogle Scholar
- Tamerius J, Nelson M, Zhou S, Viboud C, Miller M, Alonso W: Global influenza seasonality: reconciling patterns across temperate and tropical regions. Environ Health Perspect. 2011, 119: 439-445.View ArticlePubMedGoogle Scholar
- The Macroepidemiology of Influenza Vaccination Study Group: The macro-epidemiology of influenza vaccination in 56 countries, 1997-2003. Vaccine. 2005, 23: 5133-5143.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/11/342/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited