Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Assessing the Short-Term Effects of Heatwaves on Mortality and Morbidity in Brisbane, Australia: Comparison of Case-Crossover and Time Series Analyses

  • Shilu Tong ,

    s.tong@qut.edu.au

    Affiliation School of Public Health and Social Work, Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia

  • Xiao Yu Wang,
  • Yuming Guo

Abstract

Background

Heat-related impacts may have greater public health implications as climate change continues. It is important to appropriately characterize the relationship between heatwave and health outcomes. However, it is unclear whether a case-crossover design can be effectively used to assess the event- or episode-related health effects. This study examined the association between exposure to heatwaves and mortality and emergency hospital admissions (EHAs) from non-external causes in Brisbane, Australia, using both case-crossover and time series analyses approaches.

Methods

Poisson generalised additive model (GAM) and time-stratified case-crossover analyses were used to assess the short-term impact of heatwaves on mortality and EHAs. Heatwaves exhibited a significant impact on mortality and EHAs after adjusting for air pollution, day of the week, and season.

Results

For time-stratified case-crossover analysis, odds ratios of mortality and EHAs during heatwaves were 1.62 (95% confidence interval (CI): 1.36–1.94) and 1.22 (95% CI: 1.14–1.30) at lag 1, respectively. Time series GAM models gave similar results. Relative risks of mortality and EHAs ranged from 1.72 (95% CI: 1.40–2.11) to 1.81 (95% CI: 1.56–2.10) and from 1.14 (95% CI: 1.06–1.23) to 1.28 (95% CI: 1.21–1.36) at lag 1, respectively. The risk estimates gradually attenuated after the lag of one day for both case-crossover and time series analyses.

Conclusions

The risk estimates from both case-crossover and time series models were consistent and comparable. This finding may have implications for future research on the assessment of event- or episode-related (e.g., heatwave) health effects.

Introduction

Heatwaves or excessive ambient heat exposures have significant impacts on mortality and morbidity [1][6]. For example, during the 1995 Chicago heatwave, there were over 700 excess deaths in a single day [7]. The well-known 2003 heatwaves led to 15,000 excess deaths in France alone [8], [9], and over 70,000 deaths across Europe [10], [11]. The 2006 California heatwave resulted in an increase in morbidity which included 16,166 excess emergency department visits and 1,182 excess hospitalizations state-wide [12]. Heat-related impacts may have greater public health implications as climate change continues. It is important to appropriately characterize the relationship between heatwaves and health outcomes.

Two common epidemiologic methods have been frequently used to assess the heat-related health effects. Time series analysis has been used to investigate the health impact of time varying environmental exposures (eg, air pollution and temperature) for many years [13], [14]. Recently, a case-crossover design (introduced by MaClure in 1991) has been increasingly used to examine an association between a transient exposure (eg, temperature or air pollution) and acute health outcomes [15], [16]. This design controls for time-invariant confounders by study design itself [17]. Therefore, it has some advantages compared with commonly-used time series analysis. However, some methodological issues in the use of case-crossover analysis have attracted much research attention. For example, unidirectional case-crossover design was initially applied and the referent period was designated by specific time period(s) before the case period [18]. Recently, ambidirectional and time-stratified case-crossover analyses have been assumed as ideal approaches because unidirectional design has often produced biased results [18][20]. The previous research mainly focused on the risk assessment of time-varying exposures (eg, air pollution and temperature) using relatively long time series datasets. However, little information is available on whether these findings are applicable to the assessment of event- or episode-related (eg, heatwave) health effects. Since time series and case-crossover methods are often viewed as two competing analytical approaches, this study examined whether these methods produced equivalent risk estimates in the assessment of the health effects of heatwaves in Brisbane, Australia.

Materials and Methods

Data collection

Brisbane, Australia's third largest city, is located in the south-east corner of the Queensland state (27°29′S, 153°8′E) and has a sub-tropical climate. The population increased from 896,649 on 30 June 2001 to 991 260 on 30 June 2006. 18% of the residents in Brisbane were aged 0–14, while 11% of them were aged 65+[21].

We obtained emergency hospital admissions (EHAs) data during 1st January 1996 to 31st December 2005, and mortality data during 1st January 1996 to 30th November 2004. Daily data on mortality and EHAs were provided by the Office of Economic and Statistical Research of the Queensland Treasury and the Health Information Centre of Queensland Health, respectively. Non-external causes (NEC) mortality and EHAs were categorised according to the International Classification of Diseases (revisions 9 and 10) (ICD 9, <800; and all ICD 10 codes excluding S00–U99 for external causes).

Daily data on maximum temperature and relative humidity data were obtained from the Australian Bureau of Meteorology during January 1996 to December 2005. The daily average values of climatic variables were calculated from five monitoring stations. We retrieved daily air pollution data from the Queensland Department of Environment and Resource Management (formerly Queensland Environmental Protection Agency), including ambient 24-hour average concentrations of particulate matter with diameter less than 10 µm (PM10), daily maximum 1-hour average nitrogen dioxide (NO2) and ozone (O3). Daily air pollution concentrations were averaged from seventeen monitoring stations in Brisbane. When data were missing for a particular monitoring station on a given day, the data recorded from other monitoring stations were used to calculate the daily average values.

Data analysis

According to the local heatwave definition (ie, daily maximum temperature higher than 37°C for two or more consecutive days) developed in the previous research [22], [23], three heatwaves occurred (20 and 21 January 2000; 24 to 26 December 2001; 21 and 22 February 2004) during the whole study period. We examined the short-term effects of heatwaves on mortality and EHAs within three different periods, 84 days (28 days as a strata length for time-stratified case-crossover analysis), summer season (from December to February) and whole study period (1996–2005). Spearman's correlation coefficients were used to evaluate the interrelations between air pollutants and climate variables in these three periods.

In this study, both time-series and case–crossover analyses were used to examine the heatwave effects on NEC mortality and EHAs. Poisson generalised additive model (GAM) was used to perform time series analyses on three different periods. We also used time-stratified case–crossover with a stratum length of 28 days, and matched case-control days using day of the week (one case day was matched with three control days). Because only three heatwaves were defined during the whole study period, three strata (84 days) were used for the case–crossover analysis. Heatwave days were categorised as 1, while non-heatwave days were categorised as 0. We adjusted for humidity and air pollutants (PM10, NO2 and O3) in these models. Additionally, lagged effects (ie, lag 0 to lag 5 and moving average of lags 0–5) were assessed using the same methods. When lagged effects of heatwaves on NEC mortality and EHAs were assessed, the same lagged effects of humidity and air pollution were controlled for. Relative risks (RRs) for time series analysis, odds ratios (ORs) for time-stratified case–crossover analysis, and 95% confidence intervals (CIs) were calculated in each model. R software (version 2.12.2) and the “mgcv” package were applied to fit the time series GAM and case–crossover analyses. The R codes are provided in Information S1 using an example dataset (Information S2), and results of estimation using different degree of freedom were presented in Figure S1.

Results

There were a total of 51,233 deaths and 488,005 EHAs recorded in Brisbane during the whole study period. Table 1 shows summary statistics of climatic variables, air pollutants, mortality and EHAs for three different periods, as well as the total number of different health outcomes used in each model. There was some variation in these variables during three periods. The mean concentrations of PM10, NO2 and O3 differed slightly across the different periods. The daily average deaths (17) and EHAs (142) during 84 days were higher than those during summer days (15 deaths and 128 EHAs) or whole study period (16 deaths and 134 EHAs).

thumbnail
Table 1. Summaries of daily weather, air pollutants and health outcomes in Brisbane, Australia.

https://doi.org/10.1371/journal.pone.0037500.t001

Table 2 shows that maximum temperatures positively correlated with most air pollutants in all periods except for NO2 during the whole study period. Moderate to high correlations were observed for temperature and PM10, NO2 and ozone. However, humidity inversely correlated with most air pollutants in different periods, although most of these correlations were week. There were significant associations between air pollutants.

thumbnail
Table 2. Spearman correlation between climatic variables and air pollutants in Brisbane, Australia.

https://doi.org/10.1371/journal.pone.0037500.t002

In time-stratified case-crossover analysis we used 28 days as a stratum length. Since three heatwaves occurred during the study period, therefore 84 days (ie, three strata: 3×28 days) were used in data analysis. Table 3 shows the estimated odds ratios (ORs) and relative risks (RRs) for NEC mortality and EHAs during heatwave days compared with non-heatwave days during these 84 days. There was a broadly consistent and statistically significant increase in ORs or RRs for both mortality and EHAs during heatwaves.

thumbnail
Table 3. Odds ratios (ORs) and Relative risks (RRs) of mortality and EHAs during heatwaves in Brisbane (84 days).

https://doi.org/10.1371/journal.pone.0037500.t003

Figure 1 shows the similar results as Table 3 when both methods was used to adjust for humidity, PM10, NO2 and O3. ORs of mortality and morbidity during heatwaves were 1.62 (95% CI: 1.36–1.94) (Figure 1A) and 1.22 (95% CI: 1.14–1.30) (Figure 1B) at lag 1, respectively. RRs of mortality and morbidity during heatwaves were 1.72 (95% CI: 1.40–2.11) (Figure 1A) and 1.14 (95% CI: 1.06–1.23) (Figure 1B) at lag 1, respectively. The risk estimates also generally alleviated after lag 1 for both mortality and EHAs.

thumbnail
Figure 1. Odds ratios and Relative risks of daily mortality (1A) and emergency hospital admissions (1B) during heatwaves (84 days), using both case crossover and time series analyses.

https://doi.org/10.1371/journal.pone.0037500.g001

Table 4 shows the delayed effects of heatwaves on NEC mortality and EHAs for summer season and whole study period using time series models to adjust for different sets of confounders. There were also broadly consistent increased risks for both mortality and EHAs during heatwaves across different periods. The highest risk estimates were observed at lag 1 for both mortality and EHAs after adjusting for confounders.

thumbnail
Table 4. Relative risks (RRs) of daily mortality and EHAs during heatwaves in Brisbane (summer season and whole study period), using time series analysis.

https://doi.org/10.1371/journal.pone.0037500.t004

Figure 2 reveals the similar results as Table 4 when time series GAM models were used with different periods to adjust for humidity, PM10, NO2 and O3. The results show that RRs of mortality and EHAs ranged from 1.77 (95% CI: 1.53–2.04) to 1.81 (95% CI: 1.56–2.10) in Figure 2A and 1.28 (95% CI: 1.21–1.36) (Figure 2B) at lag 1, respectively. The risk estimates generally attenuated after the lag of one day for both mortality and EHAs.

thumbnail
Figure 2. Relative risks of daily mortality (2A) and emergency hospital admissions (2B) during heatwaves (summer season and whole study period), using time series analyses.

https://doi.org/10.1371/journal.pone.0037500.g002

Discussion

An increase in the frequency, duration and intensity of heatwaves is one of the most certain impacts of global climate change [24], and therefore, it is important to characterise the heat-related health risks. In this study, we found consistent and significant risks of NEC deaths and EHAs during heatwaves using both time series and case-crossover methods. For time series analyses, a similar pattern of heatwave-related risks was observed even though different periods of data were used. In general, the effect estimates using time series analyses were quite comparable to those estimated by the case-crossover method. We also found that heatwaves had acute effects on mortality and EHAs, and the highest impact usually occurred at lag 1 in Brisbane – a subtropical city in Australia.

In a recent study of the relationship between temperature and mortality among the elderly, the results from a case-crossover study design using three different approaches for the selection of referent periods were compared with those from the time-series analysis [19]. Consistent results were found by using either the time-series or time-stratified case-crossover analysis. Our findings support and extend the previous study, demonstrating that time-stratified case-crossover design can also be used to estimate the event or episode-related (e.g., heatwave) health effects, and both time series and time-stratified case-crossover methods produced robust and comparable results.

Previous studies reported that increases in mortality and morbidity usually occurred within a short time frame after a heatwave [5], [25][28]. For example, in the 2003 French heatwave, relative risk for mortality increased rapidly in most cities during a 4-day period from 10 to 13 August [9]. Similarly, in the Chicago heatwave, there was the highest dose–response relationship for lags of 1 and 2 days of high minimum or maximum temperatures before all-cause mortality occurred, with an RR of 1.20 (95% CI: 1.15–1.27) for high temperature on lag day 1 and an indication of a cumulative effect (RR = 1.40; 95% CI:1.30–1.51) when temperatures stayed high both on lag days 1 and 2 [15]. In our study, we also found that exposure to heatwaves at the same day, lag 1 and lags 0–5 was strongly associated with risks of mortality and EHAs. The strongest heatwave effects appeared at lag 1 and lagged effects were gradually decreased after lag 1 for both EHAs and mortality (Figures 1 and 2).

Previous studies including our own research show an immediate impact of high temperature on the patients with cardiovascular diseases [29], suggesting that cardiovascular deaths tend to happen quickly during hot days because their ability to cope with heat is already compromised. Several pathophysiological mechanisms may explain this. Firstly, heat stress can reduce cerebral blood velocity and markedly impair orthostatic tolerance in humans [23], [30]. Secondly, water loss and reduced plasma volume during hot days may facilitate the release of platelets into circulation and increase red and white cell counts, blood viscosity, and plasma cholesterol levels, which may be connected to the increased mortality from arterial thrombosis in hot weather [31]. Finally, results from a mice study indicate that heat stress will stimulate cells of living organisms to generate heat shock proteins which may cause systematic damages in the body [32].

This study has three main strengths: (1) it is the first study to compare time series and case crossover analyses in the examination of heatwave effects on mortality and EHAs. Time-stratified case-crossover design is comparable to time series analysis in estimating event or episode-related health risks; (2) sophisticated statistical methods were used to assess lagged and cumulative average effects of heatwaves on both mortality and EHAs after adjustment for confounding factors; and (3) the datasets used is this study are quite comprehensive, with no missing values.

Several limitations of this study must also be acknowledged. We only focused on one city, so the results might not be generalisable to other areas. However, the approaches applied in this study can be used in further research in other areas. We only compared time series and case crossover analyses using the data on NEC mortality and EHAs. We did not use cause-specific mortality and morbidity because the number for each category was much smaller. Additionally, there might be exposure misclassification, as we used exposure data from fixed monitors rather than individual exposure data.

In conclusion, the results of this study demonstrate that both time-stratified case-crossover and time series analyses produced a similar pattern of the relationship between heatwave and health outcomes. This finding may have implications for future studies of event or episode-related health effects.

Supporting Information

Figure S1.

The coefficients of temperature effects on mortality using different degrees of freedom for day of the year.

https://doi.org/10.1371/journal.pone.0037500.s001

(TIFF)

Acknowledgments

We thank the Queensland Health, Department of Environment and Resource Management, Office of Economic and Statistical Research of the Queensland Treasury, and Australia Bureau of Meteorology for providing the relevant data.

Author Contributions

Conceived and designed the experiments: ST. Performed the experiments: XW YG. Analyzed the data: XW YG. Contributed reagents/materials/analysis tools: XW YG. Wrote the paper: ST XW YG. Had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis: ST.

References

  1. 1. Anderson GB, Bell ML (2011) Heat Waves in the United States: Mortality Risk during Heat Waves and Effect Modification by Heat Wave Characteristics in 43 U.S. Communities. Environ Health Perspect 119: 210–218.
  2. 2. Gasparrini A, Armstrong B (2011) The impact of heat waves on mortality. Epidemiology 22: 68–73.
  3. 3. Fouillet A, Rey G, Wagner V, Laaidi K, Empereur-Bissonnet P, et al. (2008) Has the impact of heat waves on mortality changed in France since the European heat wave of summer 2003? A study of the 2006 heat wave. Int J Epidemiol 37: 309–317.
  4. 4. Kovats RS, Hajat S, Wilkinson P (2004) Contrasting patterns of mortality and hospital admissions during hot weather and heat waves in Greater London, UK. Occupational and Environmental Medicine 61: 893–898.
  5. 5. Michelozzi P, Accetta G, De Sario M, D'Ippoliti D, Marino C, et al. (2009) High temperature and hospitalizations for cardiovascular and respiratory causes in 12 European cities. Am J Respir Crit Care Med 179: 383–389.
  6. 6. Argaud L, Ferry T, Le QH, Marfisi A, Ciorba D, et al. (2007) Short- and long-term outcomes of heatstroke following the 2003 heat wave in Lyon, France. Arch Intern Med 167: 2177–2183.
  7. 7. Semenza JC, Rubin CH, Falter KH, Selanikio JD, Flanders WD, et al. (1996) Heat-related deaths during the July 1995 heat wave in Chicago. The New England Journal Of Medicine 335: 84–90.
  8. 8. Fouillet A, Rey G, Jougla E, Frayssinet P, Bessemoulin P, et al. (2007) A predictive model relating daily fluctuations in summer temperatures and mortality rates. BMC Public Health 7: 114.
  9. 9. Le Tertre A, Lefranc A, Eilstein D, Declercq C, Medina S, et al. (2006) Impact of the 2003 heatwave on all-cause mortality in 9 French cities. Epidemiology 17: 75–79.
  10. 10. Conti S, Meli P, Minelli G, Solimini R, Toccaceli V, et al. (2005) Epidemiologic study of mortality during the Summer 2003 heat wave in Italy. Environ Res 98: 390–399.
  11. 11. Johnson H, Kovats RS, McGregor G, Stedman J, Gibbs M, et al. (2005) The impact of the 2003 heat wave on mortality and hospital admissions in England. Health Stat Q 6–11.
  12. 12. Knowlton K, Rotkin-Ellman M, King G, Margolis HG, Smith D, et al. (2009) The 2006 California heat wave: impacts on hospitalizations and emergency department visits. Environ Health Perspect 117: 61–67.
  13. 13. Anderson BG, Bell ML (2009) Weather-related mortality: how heat, cold, and heat waves affect mortality in the United States. Epidemiology 20: 205–213.
  14. 14. Hajat S, Armstrong B, Baccini M, Biggeri A, Bisanti L, et al. (2006) Impact of high temperatures on mortality: is there an added heat wave effect? Epidemiology 17: 632–638.
  15. 15. Kaiser R, Le Tertre A, Schwartz J, Gotway CA, Daley WR, et al. (2007) The effect of the 1995 heat wave in Chicago on all-cause and cause-specific mortality. Am J Public Health 97: Suppl 1S158–162.
  16. 16. Stafoggia M, Forastiere F, Agostini D, Biggeri A, Bisanti L, et al. (2006) Vulnerability to Heat-Related Mortality: A Multicity, Population-Based, Case-Crossover Analysis. Epidemiology 17: 315–323.
  17. 17. Maclure M (1991) The case-crossover design: a method for studying transient effects on the risk of acute events. Am J Epidemiol 133: 144–153.
  18. 18. Lumley T, Levy D (2000) Bias in the case – crossover design: implications for studies of air pollution. Environmetrics 11: 689–704.
  19. 19. Basu R, Dominici F, Samet JM (2005) Temperature and mortality among the elderly in the United States - A comparison of epidemiologic methods. Epidemiology 16: 58–66.
  20. 20. Levy D, Lumley T, Sheppard L, Kaufman J, Checkoway H (2001) Referent selection in case-crossover analyses of acute health effects of air pollution. Epidemiology 12: 186–192.
  21. 21. QLD Queensland Government Population and Housing Fact sheet - Brisbane City Council (2009) Available: http://www.oesr.qld.gov.au/queensland-by-theme/demography/population-characteristics/profiles/pop-housing-fact-sheets-lga/pop-housing-fact-sheets-brisbane-200908.pdf. Accessed 2009 Nov 26.
  22. 22. Tong S, Wang XY, Barnett AG (2010) Assessment of heat-related health impacts in Brisbane, Australia: comparison of different heatwave definitions. PLoS One 5: e12155.
  23. 23. Wang XY, Barnett AG, Yu W, FitzGerald G, Tippett V, et al. (2012) The impact of heatwaves on mortality and emergency hospital admissions from non-external causes in Brisbane, Australia. Occup Environ Med 69: 163–169.
  24. 24. IPCC (2007) Intergovernmental Panel on Climate Change. Climate Change 2007: Synthesis Report. Geneva, Switzerland. Available: http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr.pdf. Accessed 2009 Nov 26.
  25. 25. Basu R, Ostro BD (2008) A multicounty analysis identifying the populations vulnerable to mortality associated with high ambient temperature in California. Am J Epidemiol 168: 632–637.
  26. 26. Bell ML, O'Neill MS, Ranjit N, Borja-Aburto VH, Cifuentes LA, et al. (2008) Vulnerability to heat-related mortality in Latin America: a case-crossover study in Sao Paulo, Brazil, Santiago, Chile and Mexico City, Mexico. Int J Epidemiol 37: 796–804.
  27. 27. Ishigami A, Hajat S, Kovats RS, Bisanti L, Rognoni M, et al. (2008) An ecological time-series study of heat-related mortality in three European cities. Environmental Health 7: 5.
  28. 28. Baccini M, Biggeri A, Accetta G, Kosatsky T, Katsouyanni K, et al. (2008) Heat effects on mortality in 15 European cities. Epidemiology 19: 711–719.
  29. 29. Huynen MM, Martens P, Schram D, Weijenberg MP, Kunst AE (2001) The impact of heat waves and cold spells on mortality rates in the Dutch population. Environ Health Perspect 109: 463–470.
  30. 30. Wilson TE, Cui J, Zhang R, Crandall CG (2006) Heat stress reduces cerebral blood velocity and markedly impairs orthostatic tolerance in humans. Am J Physiol Regul Integr Comp Physiol 291: R1443–1448.
  31. 31. Keatinge WR, Coleshaw SR, Easton JC, Cotter F, Mattock MB, et al. (1986) Increased platelet and red cell counts, blood viscosity, and plasma cholesterol levels during heat stress, and mortality from coronary and cerebral thrombosis. American Journal of Medicine 81: 795–800.
  32. 32. Li SQ, Li RF, Xi SM, Hu S, Jia ZQ, et al. (2012) Systematical analysis of impacts of heat stress on the proliferation, apoptosis and metabolism of mouse hepatocyte. J Physiol Sci 62: 29–43.