Identifying Individual Risk Factors and Documenting the Pattern of Heat-Related Illness through Analyses of Hospitalization and Patterns of Household Cooling

Abstract

Background

As climate change increases the frequency and intensity of extreme heat events researchers and public health officials must work towards understanding the causes and outcomes of heat-related morbidity and mortality. While there have been many studies on both heat-related illness (HRI), there are fewer on heat-related morbidity than on heat-related mortality.

Objective

To identify individual and environmental risk factors for hospitalizations and document patterns of household cooling.

Methods

We performed a pooled cross-sectional analysis of secondary U.S. data, the Nationwide Inpatient Sample. Risk ratios were calculated from multivariable models to identify risk factors for hospitalizations. Hierarchical modeling was also employed to identify relationships between individual and hospital level predictors of hospitalizations. Patterns of air conditioning use were analyzed among the vulnerable populations identified.

Results

Hospitalizations due to HRI increased over the study period compared to all other hospitalizations. Populations at elevated risk for HRI hospitalization were blacks, males and all age groups above the age of 40. Those living in zip-codes in the lowest income quartile and the uninsured were also at an increased risk. Hospitalizations for HRI in rural and small urban clusters were elevated, compared to urban areas.

Conclusions

Risk factors for HRI include age greater than 40, male gender and hospitalization in rural areas or small urban clusters. Our analysis also revealed an increasing pattern of HRI hospitalizations over time and decreased association between common comorbidities and heat illnesses which may be indicative of underreporting.

About the Authors:

Michael T. Schmeltz

* E-mail: [email protected]

Affiliation: School of Public Health, City University of New York (CUNY), New York, United States of America

Grace Sembajwe

Affiliations School of Public Health, City University of New York (CUNY), New York, United States of America, Hunter College, City University of New York (CUNY), New York, United States of America

Peter J. Marcotullio

Affiliations Hunter College, City University of New York (CUNY), New York, United States of America, CUNY Institute for Sustainable Cities, New York, United States of America

Jean A. Grassman

Affiliations School of Public Health, City University of New York (CUNY), New York, United States of America, Health and Nutrition Sciences, Brooklyn College, City University of New York (CUNY), Brooklyn, NY, United States of America

David U. Himmelstein

Affiliations School of Public Health, City University of New York (CUNY), New York, United States of America, Hunter College, City University of New York (CUNY), New York, United States of America

Stephanie Woolhandler

Affiliations School of Public Health, City University of New York (CUNY), New York, United States of America, Hunter College, City University of New York (CUNY), New York, United States of America

Introduction

Heat from high environmental temperatures is a natural hazard that can adversely affect human health [1–3]. Extreme heat events account for more fatalities in the United States than any other weather hazard, with approximately 650 deaths per year [4]. As climate change increases the frequency and intensity of extreme heat events researchers and public health officials must work aggressively towards understanding the causes and outcomes of heat-related morbidity and mortality [5–7]. While there have been many studies on both heat-related morbidity and mortality [8–10], there are fewer studies on heat-related morbidity than those which focus on heat-related mortality [11–15]. The effect of temperature on morbidity is a significant public health concern that will escalate as global mean temperature continues to increase and more episodes of extreme heat events affect large geographic regions [16].

Human populations have acclimatized to a variety of local climates physiologically and behaviorally, though there are limits to the extremes in temperature the human body can tolerate. An individual’s ability to thermoregulate protects the body from ambient temperatures by maintaining a core body temperature around 36.0°C to 37.5°C [17]. Body temperatures outside of this range are classified as either hypothermia (body temperature below 35.0°C) or hyperthermia (body temperature above 37.5°C) [18,19]. Heat-related illnesses (HRIs) are a spectrum of conditions ranging from minor to life threatening. They include heat stress, heat exhaustion and heat stroke. Heat stress is a perceived discomfort and physiological strain; heat exhaustion, a mild to moderate illness; and heat stroke, a severe illness characterized by a core body temperature above 40.0°C [19]. Clinically, severe heat illness can be seen as a combination of systematic inflammatory response and cytotoxicity, which left untreated results in multiple organ failure [17,19,20]. Infants are at increased risk of heat-related illness as are the elderly, individuals taking medications that interfere with thermoregulation, and those with underlying medical conditions such as cardiac, respiratory, renal disease or diabetes [3,21,22].

Studies of heat mortality during specific events, such as the 1995 Chicago heat wave and the 2003 heat wave that affected much of Europe, highlight characteristics of populations that are at risk of dying during extreme heat events [9,23]. These studies described elevated heat mortality among the elderly, persons of color or of lower socioeconomic status, individuals confined to bed or unable to care for themselves and individuals with social risk factors such as not leaving the home daily, living alone, and limited social networks [23,24]. Studies have also identified underlying medical conditions as risk factors for heat mortality, including cardiovascular, cerebrovascular, respiratory, renal and neurologic diseases, and diabetes [25,26]. However, the diagnosis of a specific heat-related illness is not always documented during heat events and diagnoses for cause of death are usually listed as the underlying medical condition [27–29]. Some studies of heat-related morbidity have focused on emergency dispatch data and emergency department visits during individual heat events, with specific HRI diagnoses [30,31], while others have analyzed excess hospitalizations during specific extreme heat events [13,26,32]. Few investigations have investigated HRI hospitalizations across different geographic regions, perhaps because classification of heat events often varies across regions [15,33]. Additionally the underreporting of HRI diagnoses may limit our understanding of health-related illness [34].

One of the strongest protectors against HRIs is ownership and use of an air conditioner [35–37]. As of 2009, 87% of homes in the United States have air conditioning, although the prevalence of air conditioning is lower (and less efficient cooling systems such as fans, are higher) among apartment dwellers and low income households [38]. A recent Morbidity and Mortality Weekly Report highlighting heat illness and deaths in New York City indicated that a majority of heat illness decedents did not have a working air conditioner, although information on cooling practices of decedents was limited [39]. Research is needed regarding a number of economic and behavioral issues that might improve access to air conditioning and reduce heat mortality and morbidity [40,41].

In addition to individual level determinants and behaviors around cooling during extreme heat events, the environment plays a key role in exacerbating or protecting us from extreme heat events. While the urban heat island (UHI) is not a new concept, it has been used to explain higher rates of HRI in urban areas [42,43]. The average ambient urban temperature is increased where vegetation has been replaced by heat-retaining materials such as asphalt and concrete. Recent studies using remote sensing and spatial analysis attempt to create heat vulnerability maps and indexes based on environmental factors [44–49]. An integrated understanding of environmental factors, individual-level factors and socio-behavioral determinants is needed to accurately identify populations most vulnerable to extreme events due to climate change. Health impacts associated with heat waves disproportionately harm vulnerable populations. However, with the increase in frequency and magnitude of climate related disasters, extreme heat events may ultimately disrupt and affect all populations. This research examines individual and environmental risk factors for HRI hospitalization using a large nationally-representative administrative database. In addition we document patterns of household cooling, with particular focus on populations at risk for HRI hospitalization.

Materials and Methods

Our primary source of data is the 2001 to 2010 Healthcare Cost and Utilization Project (HCUP) Nationwide Inpatient Sample (NIS), which was developed by the Agency for Healthcare Research and Quality (AHRQ). The NIS is an all-payer database (i.e. contains patients without regard to their health insurance), with data from approximately 8 million hospitalizations per year (about 20% of all hospital discharges in the United States) [50]. The AHRQ provides a weighting variable for descriptive and regression analyses to provide a better estimate of error and statistical significance for national estimates. The NIS includes data on patient demographics and hospital characteristics, primary and secondary diagnoses coded according to the International Classification of Diseases, Ninth Revision (ICD-9) (with up to 15 diagnoses available for the years 2001–2008 and up to 25 diagnoses for 2009–2010), primary and secondary procedures (ICD-9 coding up to 15 diagnoses available for the years 2001–2008 and up to 25 diagnoses for 2009–2010), type and source of admission, discharge disposition, primary payer type, total hospital charges, and length of stay. Stratification and weighting variables allow calculation of national estimates that account for the complex sampling design [50,51]. These data include information on urban vs. rural hospital location. Additional spatial data for geocoded hospital location was obtained from ESRI based on American Hospital Association (AHA) identifiers, and merged to the NIS dataset to create a detailed and updated hospital location indicator. Using the 2010 Census Urban and Rural classification, geocoded hospitals were mapped and designated as being in urbanized areas of 50,000 or more people, in small urban clusters of 2,500 to 50,000 people or in rural areas which encompassed all populations outside urbanized areas and small urban clusters [52]. All spatial joinings were done using ArcGIS (version 10; ESRI, Redlands, CA).

Data on air conditioner use was derived from the Residential Energy Consumption Survey (RECS), a national household survey conducted by the United States Energy Information Administration (EIA) starting in 1978. In 2005 and 2009 surveys were collected from approximately 4,000 and 12,000 households, respectively, which were nationally representative of the 110 million residential U.S. households. We analyzed RECS data on air conditioner ownership and use examining households below 100% of the federal poverty line, and below 150% of the federal poverty line, employment status, and age of householder (head of household) [38].

The study population for HRI consists of patients in the 2001–2010 NIS with at least one diagnosis of a heat-related illness (ICD-9 codes 992.0–992.9). The control group consisted of all other NIS hospitalizations during this time period. The HRI variable was coded as a binary indicator variable. The initial cohort consisted of 15,885 patient records for heat-related illness (HRI) and 79,140,110 patient records for controls. Since approximately 95% of all HRI hospitalizations occurred in the summer months (May-September), hospitalizations in other months were excluded yielding a final cohort of 14,949 HRI patient records and 37,019,792 controls.

The outcome variable was heat-related illness. First, a descriptive bivariate analysis of baseline characteristics was performed for HRI and control populations. Chi-square tests and t-tests were used to examine differences between these populations. Second, multivariable analysis was performed to determine risk factors associated with heat morbidity. A log-binomial model, which directly models risk ratios (RRs), was used. Analyses were weighted for national estimates and to enhance likelihood ratios [53]. Risk factors for heat-related illnesses were explored, including patient characteristics (gender, age, race, zip-code income quartile, insurance, comorbidities), and hospital characteristics (hospital size—based on number of beds, hospital location, hospital region). Patient comorbidities were analyzed using the AHRQ algorithms based on methods developed by Elixhauser et al [54]. Third, a multilevel (i.e. hierarchical) model was used to better account for the correlation between patient and hospital-level characteristics. Weighting was not used for multilevel analysis as recommended by the HCUP hierarchical modeling report [55]. Lastly air conditioner data was used to analyze the ownership and use of air conditioning among vulnerable populations as identified in the hospitalization data.

All statistical analyses were performed using SAS 9.3. All data used in this study are publically available and had been anonymized and de-identified prior to our analyses. The study was a secondary data analysis and data were de-identified prior to being obtained by the authors and did not require patient consent. This study was also approved by the Institutional Review Board of Hunter College, CUNY and conforms to the HCUP data use agreement for data protection.

Results

The characteristics of patients with HRI and of all hospitalized patients during the five “summer” months from 2001 to 2010 in the NIS data are shown in Tables 1 and 2. After weighting there were 73,185 patient discharges for HRI and 181,094,795 for all other discharges during this time period. Compared to all non-HRI hospitalizations, more HRI admissions were emergencies. Emergency department (ED) admission (as source of admission) accounted for 59% for HRI hospitalizations and 34% of non-HRI hospitalizations, while emergency admission (as type of admission) accounted for 73% for HRI and 41% of non-HRI hospitalizations. The mean age was also higher for HRI than non-HRI hospitalizations, 55 years compared to 48 years. These differences were statistically significant (p < 0.05). The number of HRI hospitalizations also increased faster than non-HRI hospitalizations, albeit with large year-to-year variation (Fig. 1).

[Figure omitted. See PDF.]

Table 1. Demographic Characteristics of Heat-Related Illness (HRI) patients, 2001–2010 (Summer).

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

[Figure omitted. See PDF.]

Table 2. Hospital Characteristics of Heat-Related Illness (HRI) patients, 2001–2010 (Summer).

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

[Figure omitted. See PDF.]

Fig 1. Heat-Related Illness versus All Hospitalizations Rates, Summer 2001–2010.

The dotted line represents the rate of all hospitalizations during the study time period which remained approximately level over time. The solid line represents heat-related illness hospitalizations which although varied by year have steadily increased over the study period.

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

Multivariable analyses identifying risk factors for patients with an HRI diagnosis and for any other hospitalization are shown in Table 3. Blacks had an elevated risk ratio of 1.21 (95% CI: 1.19–1.23) compared to Whites; males had a markedly higher risk ratio of 3.54 (95% CI: 3.49–3.58) compared to females; as did each age group above the age of 40 (ages 40–64) RR 1.34 (95% CI: 1.32–1.37), (ages 65–74) RR 1.28 (95% CI: 1.25–1.30), (ages 75+) RR 1.52 (95% CI: 1.47–1.55) compared to ages 18–39. Patients who lived in zip codes with the lowest median income quartile were more likely to be hospitalized for a HRI, RR 1.19 (95% CI: 1.17–1.21) as were those in the second lowest income quartile RR 1.08 (95% CI: 1.07–1.10) compared to the highest zip code income quartile. Lack of health insurance was a particularly strong social predictor of being hospitalized due to HRI, RR 2.33 (95% CI: 2.29–2.37).

[Figure omitted. See PDF.]

Table 3. Multivariable model of risk factors of hospitalization for heat-related illnesses, United State 2001–2010 (Summer).

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

Patients with neurological disorders and psychoses were more likely to be hospitalized with HRI, with a relative risk of 1.79 (95% CI: 1.75–1.82) and 1.83 (95% CI: 1.79–1.87) respectively. Other comorbidities such as cardiovascular diseases, diabetes, renal failure and pulmonary diseases did not predict HRI relative to non-HRI hospitalizations with risk ratios of 0.80 (95% CI: 0.78–0.82) for congestive heart failure, 0.72 (95% CI: 0.71–0.82) for chronic pulmonary diseases, 0.90 (95% CI: 0.89–0.92) for diabetes, 1.00 (95% CI: 0.99–1.01) for hypertension and 0.61 (95% CI: 0.60–0.63) for renal failure. Examination of hospital characteristics indicate that patients in the Western U.S. region had higher risk of HRI hospitalization, RR 1.23 (95% CI: 1.20–1.25), compared to those in the Northeastern US, as did patients at small and medium sized hospitals (RR 1.39 95% CI: 1.37–1.42 and RR 1.25 95% CI: 1.23–1.27 respectively), compared to large hospitals. Hospitals located in rural areas, (RR 1.79 95% CI: 1.74–1.83) and in small urban clusters (RR 2.09 95% CI: 2.06–2.12) had increased risk of HRI hospitalizations as compared to those in urban areas.

The multilevel analysis (Table 4), found risk ratios similar to the single-level multivariable analysis: Blacks (RR 1.17, 95% CI: 1.08–1.27), males (RR 3.55, 95% CI: 3.35–3.77) and patients above the age of 40 (ages 40–64) (RR 1.29, 95% CI: 1.19–1.40), (ages 65–74) (RR 1.15, 95% CI: 1.04–1.27), (ages 75+) (RR 1.40, 95% CI: 1.28–1.53). Being in the lowest zip code income quartile significantly predicted of HRI at RR 1.13 (95% CI: 1.03–1.25) as did being uninsured (RR 2.51, 95% CI: 2.31–2.72). Hospital level characteristics for the multilevel analysis also had similar results with small hospital size having a relative risk of 1.35 (95% CI: 1.19–1.53) and medium hospital size a relative risk of 1.17(95% CI: 1.05–1.30). Risk of HRI hospitalization in rural areas and small urban clusters also remained higher than in urban areas. Analyses of the energy use data revealed that overall air conditioner use for all U.S. households increased from 82% in 2005 to 87% in 2009. However, households below 100% to 150% of the federal poverty line had 3% to 6% decrease in using an air conditioner from 2005 to 2009. Air conditioner use among householders who were unemployed, retired or employed part-time was 25% less likely in both 2005 and 2009, than among householders employed full-time. Additionally, air conditioner use increased among householders of most age groups from 2005 to 2009, but decreased by 13% among householders aged 75 and older.

[Figure omitted. See PDF.]

Table 4. Multi-level multivariable model of risk factors of hospitalization for heat-related illnesses, United State 2001–2010 (Summer).

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

Discussion

Most of the risk factors for heat related morbidity that we identified were similar to those found in previous heat-related morbidity and mortality studies [27,56,57]. For instance persons identified as Black, elderly and of lower socioeconomic status were at a higher risk for HRI hospitalizations. A recent study examining United States emergency department visits for acute heat illnesses indicated that males were at a higher risk, as did we [56].

However, our nationally representative data yielded some surprising findings as well. First, all older age groups were at an increased risk of heat-related illness compared to 18–39 year olds. Prior studies have frequently highlighted only infants and elderly as being most vulnerable, especially to heat mortality [58,59]. This increased risk of hospitalization among middle-aged adults, not just the elderly may signal a lack of risk perception among the population in regards to heat-related illness. Studies have shown that adults and many elderly populations do not perceive themselves as at risk for a HRI [60]. A study of older people in London and Norwich, United Kingdom found that while people are able to recognize warnings and provide examples of what to do to reduce the effects of heat, many did not apply these warnings or behaviors as they failed to consider themselves at risk [61].

Second, much research has been done on the urban heat island effect and its implications on heat-related morbidity and mortality [62]; our study indicates that heat morbidity is not restricted to urban environments. Risk of hospitalization due to HRI was actually higher in small urban clusters (populations 2,500 to 50,000) and rural areas as compared to urban areas (populations 50,000+). This is in contrast to current research indicating that the urban heat island is exacerbating warming in urban areas. A recent study from China indicated that summer mortality rates have increased in large urban areas like Shanghai compared to suburban locales [63]. Additional research is needed to understand why risk of HRI hospitalization is higher in small urban clusters and rural areas, despite the fact that urban areas have higher temperatures, and hold heat for longer periods of time. Hospital travel distance, absent heat-health warning systems and social factors such as higher rates of manual labor and lower health literacy outside of urban areas may explain these differences.

Third, surprisingly, we found that common serious comorbidities, aside from neurological disorders and psychoses, were not a significant risk factor for being hospitalized for a heat-related illness compared to all other hospitalizations. This is contrary to previous research indicating that these common comorbidities are exacerbated during extreme heat events [12,24]. We interpret this cautiously as it may indicate only that other hospitalizations (i.e. our comparator) are very frequent for persons with “chronic” medical conditions. Additionally artifacts may have been introduced if the severity of the underlying condition and comorbidities led providers to under-diagnose or under-code HRI. Such underreporting obstructs our understanding of risk factors for heat-related illnesses.

Lastly, our results also indicate that individuals from poorer neighborhoods and especially persons without health insurance are at an increased risk of hospitalizations due to heat-related illnesses. Although low socioeconomic status has been shown to be a predictor for heat morbidity our insurance status indicator measures a specific aspect of socioeconomic status among hospitalizations of heat-related illnesses.

As pervasive as air conditioning use has become in the United States, many among these vulnerable populations still lack access to this protective factor [64]. Our air conditioning data showed that poor and elderly householders, as well as those who are unemployed or underemployed have lower rates of air conditioning use. These same populations are at the highest risk of being hospitalized for heat-related illnesses. In some communities air conditioning is still considered a luxury, but the increase in extreme heat events may make air conditioning a necessity. Steps are needed to improve the availability and affordability of air conditioning and to improve access to cooling centers. One study from California (where temperature variation is less extreme than in other parts of the U.S) indicated that while central air conditioning had a protective effect on heat morality, central air conditioning use was not associated with income [36]. However, the study did not analyze the use of room air conditioning or use of cooling centers which may impact lower socioeconomic populations. While air conditioning may prevent heat-related illness hospitalizations, equitable policies surrounding energy consumption during summer months and during heat events will become an important topic as the demand for energy increases and further stresses energy distribution systems [65].

To our knowledge this is one of the few studies to assess HRI nationally [46,56]. However, it has limitations. Administrative data is less accurate than clinical data and HRI is almost certainly under-coded. Because the NIS does not capture information regarding hospital readmission, an individual can be counted multiple times if hospitalized repeatedly for HRI. We were unable to account for clustered patient data and thus may slightly underestimate the width of our confidence intervals, although given the large data set we believe it should make little difference in the overall results. The NIS is a 20% sample survey and may be subject to sampling error. A sub-population may seriously be over-represented or under-represented. To overcome this limitation in the sampling design we included the use of weights that allow for national estimates. Reporting for some hospital level variables also varied from state to state. Additionally, the RECS data does not link with NIS records; thus comparisons could only be made at an ecological level concerning air conditioning use.

Heat is the number one cause of death from natural disasters in the United States [4], yet is often overlooked due to the limited visibility of heat as a hazard. Tornados, hurricanes, floods and blizzards are more tangible and dramatic than heat waves. The latest Intergovernmental Panel on Climate Change (IPCC) states that global mean temperature is projected to increase between 1.5 and 2.3°C by mid-century and over 4°C by the end of the century [66]. As temperatures increase the intensity and frequency of extreme heat events will also increase. In turn this will likely increase the number of hospitalizations of heat-related illnesses. Our research is pertinent in understanding the predictors for HRI hospitalizations and identifies three key risk factors for heat morbidity, including ages greater than 40, hospitalization in rural areas and small urban clusters and lack of health insurance as a predictor. Our analysis also revealed an increasing pattern of HRI hospitalizations over time. Our findings reinforce previous research on heat morbidity and explore the prevalence and use of air conditioning among vulnerable populations. Future research should employ spatial analytic techniques to identify environmental risk factors with focus on the differences between rural and urban heat-related morbidity. Additionally research should address common comorbidities associated with HRI and HRI's impact on hospital resource use and costs.

Author Contributions

Conceived and designed the experiments: MS GS JG PM DH SW. Performed the experiments: MS GS. Analyzed the data: MS GS. Contributed reagents/materials/analysis tools: PM JG DH SW. Wrote the paper: MS GS PM JG DH SW.

Citation: Schmeltz MT, Sembajwe G, Marcotullio PJ, Grassman JA, Himmelstein DU, Woolhandler S (2015) Identifying Individual Risk Factors and Documenting the Pattern of Heat-Related Illness through Analyses of Hospitalization and Patterns of Household Cooling. PLoS ONE 10(3): e0118958. https://doi.org/10.1371/journal.pone.0118958

References

1. Li T, Horton RM, Kinney PL. Projections of seasonal patterns in temperature-related deaths for Manhattan, New York. Nature Climate Change. 2013;3(8): 717–721

2. Hansen AL, Bi P, Ryan P, Nitschke M, Pisaniello D, Tucker G. The effect of heat waves on hospital admissions for renal disease in a temperate city of Australia. Int J Epidemiol. 2008;37: 1359–1365. pmid:18710886

3. Curriero FC, Heiner KS, Samet JM, Zeger SL, Strug L, Patz J. Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol. 2002;155: 80–87. pmid:11772788

4. Centers for Disease Control and Prevention (CDC). Heat-related deaths after an extreme heat event—four states, 2012, and united states, 1999–2009. MMWR Morbidity and Mortality Weekly Reports. 2013;62: 433–436. pmid:23739336

5. Kovats R, Menne B, Ahern M, Patz J. National assessments of health impacts of climate change: A review. Climate Change and Human Health: Risks and Responses (McMichael AJ, Campbell-Lendrum D, Corvalan CF, Ebi KL, Githeko A, Scheraga JD, Woodward A, eds).Geneva: WHO/World Meteorological Organization/United Nations Environment Programme. 2003;181–203.

6. Bell ML, O’Neill MS, Ranjit N, Borja-Aburto VH, Cifuentes LA, Gouveia NC. 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. 2008;37: 796–804. pmid:18511489

7. Medina-Ramón M, Zanobetti A, Cavanagh DP, Schwartz J. Extreme temperatures and mortality: Assessing effect modification by personal characteristics and specific cause of death in a multi-city case-only analysis. Environ Health Perspect. 2006;114: 1331–1336. pmid:16966084

8. Knowlton K, Lynn B, Goldberg RA, Rosenzweig C, Hogrefe C, Rosenthal JK, et al. Projecting heat-related mortality impacts under a changing climate in the New York City region. American Journal of Public Health. 2007;97(11): 2028–2034. pmid:17901433

9. Semenza JC, Rubin CH, Falter KH, Selanikio JD, Flanders WD, Howe HL, et al. Heat-related deaths during the July 1995 heat wave in Chicago. N Engl J Med. 1996;335: 84–90. pmid:8649494

10. Vaneckova P, Beggs PJ, de Dear RJ, McCracken KW. Effect of temperature on mortality during the six warmer months in Sydney, Australia, between 1993 and 2004. Environ Res. 2008;108: 361–369. pmid:18774130

11. Kovats RS, Hajat S. Heat stress and public health: A critical review. Annu Rev Public Health. 2008;29: 41–55. pmid:18031221

12. Kenny GP, Yardley J, Brown C, Sigal RJ, Jay O. Heat stress in older individuals and patients with common chronic diseases. Can Med Assoc J. 2010;182: 1053–1060.

13. Mastrangelo G, Fedeli U, Visentin C, Milan G, Fadda E, Spolaore P. Pattern and determinants of hospitalization during heat waves: An ecologic study. BMC Public Health. 2007;7(1): 200.

14. Li B, Sain S, Mearns LO, Anderson HA, Kovats S, Ebi KL, et al. The impact of extreme heat on morbidity in Milwaukee, Wisconsin. Clim Change. 2012;110: 959–976.

15. Ye X, Wolff R, Yu W, Vaneckova P, Pan X, Tong S. Ambient temperature and morbidity: A review of epidemiological evidence. Environ Health Perspect. 2012;120(1): 19–28. pmid:21824855

16. Meehl GA, Tebaldi C. More intense, more frequent, and longer lasting heat waves in the 21st century. Science. 2004;305: 994–997. pmid:15310900

17. Lugo-Amador NM, Rothenhaus T, Moyer P. Heat-related illness. Emerg Med Clin North Am. 2004;22(2): 315–328. pmid:15163570

18. Conlon KC, Rajkovich NB, White-Newsome JL, Larsen L, O’Neill MS. Preventing cold-related morbidity and mortality in a changing climate. Maturitas. 2011;69(3): 197–202. pmid:21592693

19. Bouchama A, Knochel JP. Heat stroke. N Engl J Med. 2002;346(25): 1978–1988. pmid:12075060

20. Yeo TP. Heat stroke: A comprehensive review. AACN Advanced Critical Care. 2004;15(2): 280–293. pmid:15461044

21. Basu R. High ambient temperature and mortality: A review of epidemiologic studies from 2001 to 2008. Environ Health. 2009;8(1): 40.

22. Kravchenko J, Abernethy AP, Fawzy M, Lyerly HK. Minimization of heatwave morbidity and mortality. Am J Prev Med. 2013;44(3): 274–282. pmid:23415125

23. Fouillet A, Rey G, Laurent F, Pavillon G, Bellec S, Guihenneuc-Jouyaux C, et al. Excess mortality related to the august 2003 heat wave in France. Int Arch Occup Environ Health. 2006; 80(1): 16–24. pmid:16523319

24. Semenza JC, McCullough JE, Flanders WD, McGeehin MA, Lumpkin JR. Excess hospital admissions during the July 1995 heat wave in Chicago. Am J Prev Med. 1999;16(4): 269–277. pmid:10493281

25. Khalaj B, Lloyd G, Sheppeard V, Dear K. The health impacts of heat waves in five regions of New South Wales, Australia: A case-only analysis. Int Arch Occup Environ Health. 2010;83(7): 833–842. pmid:20464412

26. Knowlton K, Rotkin-Ellman M, King G, Margolis HG, Smith D, Solomon G, et al. The 2006 California heat wave: Impacts on hospitalizations and emergency department visits. Environ Health Perspect. 2009;117(1): 61–67. pmid:19165388

27. Stafoggia M, Forastiere F, Agostini D, Biggeri A, Bisanti L, Cadum E, et al. Vulnerability to heat-related mortality: A multicity, population-based, case-crossover analysis. Epidemiology. 2006;17(3): 315–323. pmid:16570026

28. D'Ippoliti D, Michelozzi P, Marino C, de'Donato F, Menne B, Katsouyanni K, et al. Research the impact of heat waves on mortality in 9 European cities: Results from the EuroHEAT project. Environ Health. 2010;9.

29. Schwartz J. Who is sensitive to extremes of temperature?: A case-only analysis. Epidemiology. 2005;16(1): 67–72. pmid:15613947

30. Bassil KL, Cole D, Moineddin R, Gournis E, Schwartz B, Craig AM, et al. Development of a surveillance case definition for heat-related illness using 911 medical dispatch data. Can J Public Health. 2008;99(4): 339–343. pmid:18767283

31. Josseran L, Caillère N, Brun-Ney D, Rottner J, Filleul L, Brucker G, et al. Syndromic surveillance and heat wave morbidity: A pilot study based on emergency departments in France. BMC Med Inform and Decision Making. 2009;9(1): 14.

32. Zhang Y, Nitschke M, Bi P. Risk factors for direct heat-related hospitalization during the 2009 Adelaide heatwave: A case crossover study. Sci Total Environ. 2013; 442: 1–5. pmid:23168533

33. Vaneckova P, Bambrick H. Cause-specific hospital admissions on hot days in Sydney, Australia. PloS ONE. 2013;8(2): e55459. pmid:23408986

34. Kilbourne EM. Heat-related illness: Current status of prevention efforts. Am J Prev Med. 2002;22(4): 328–329. pmid:11988389

35. Kondo M, Ono M, Nakazawa K, Kayaba M, Minakuchi E, Sugiomoto K, et al. Population at high-risk of indoor heatstroke: The usage of cooling appliances among urban elderlies in Japan. Environ Health and Preven Med. 2013;18(3): 251–257. pmid:23160849

36. Ostro B, Rauch S, Green R, Malig B, Basu R. The effects of temperature and use of air conditioning on hospitalizations. Am J Epidemiol. 2010;172(9): 1053–1061. pmid:20829270

37. Nunes B, Paixão E, Dias CM, Nogueira P, Falcão JM. Air conditioning and intrahospital mortality during the 2003 heatwave in Portugal: Evidence of a protective effect. Occup Environ Med. 2011;68(3): 218–223. pmid:20921273

38. United States Energy Information Administration (EIA). Residential energy consumption survey (RECS). About the RECS. Residential Energy Consumption Survey 2013. Available: http://www.eia.gov/consumption/residential/about.cfm.

39. Centers for Disease Control and Prevention (CDC). Heat illness and deaths—New York City, 2000–2011. MMWR Morb Mortal Wkly Rep. 2013;62: 617–621.

40. Yun GY, Steemers K. Behavioural, physical and socio-economic factors in household cooling energy consumption. Appl Energy. 2011;88(6): 2191–2200.

41. Richard L, Kosatsky T, Renouf A. Correlates of hot day air-conditioning use among middle-aged and older adults with chronic heart and lung diseases: The role of health beliefs and cues to action. Health Educ Res. 2011; 26(1): 77–88. pmid:21068164

42. Imhoff ML, Zhang P, Wolfe RE, Bounoua L. Remote sensing of the urban heat island effect across biomes in the continental USA. Remote Sens Environ. 2010;114(3): 504–513.

43. Stone B, Hess JJ, Frumkin H. Urban form and extreme heat events: Are sprawling cities more vulnerable to climate change than compact cities? Environ Health Perspect. 2010;118(10): 1425–1428. pmid:21114000

44. Reid CE, Mann JK, Alfasso R, English PB, King GC, Lincoln RA, et al. Evaluation of a heat vulnerability index on abnormally hot days: An environmental public health tracking study. Environ Health Perspect. 2012;120(5): 715–720. pmid:22538066

45. Kershaw SE, Millward AA. A spatio-temporal index for heat vulnerability assessment. Environ Monit Assess. 2012;184(12): 7329–7342. pmid:22270590

46. Reid CE, O’Neill MS, Gronlund CJ, Brines SJ, Brown DG, Diez-Roux AV, et al. Mapping community determinants of heat vulnerability. Environ Health Perspect. 2009; 117(11): 1730–1736. pmid:20049125

47. Kestens Y, Brand A, Fournier M, Goudreau S, Kosatsky T, Maloley M, et al. Modelling the variation of land surface temperature as determinant of risk of heat-related health events. Intrnal J of Health Geographics. 2011;10: 7. pmid:21251286

48. Chow WT, Chuang W, Gober P. Vulnerability to extreme heat in metropolitan phoenix: Spatial, temporal, and demographic dimensions. The Professional Geographer. 2012;64(2): 286–302.

49. Johnson D, Webber J, Beerval Ravichandra K, Lulla V, Stanforth A. Spatio-temporal variations in heat-related health risk in three midwestern US cities between 1990 and 2010. Geocarto International. 2014;29: 65–84.

50. Healthcare Cost and Utilization Project (HCUP). Overview of the nationwide inpatient sample. 2013. Avilable: http://www.hcup-us.ahrq.gov/nisoverview.jsp.

51. Steiner C, Elixhauser A, Schnaier J. The healthcare cost and utilization project: An overview. Effective clinical practice: ECP. 2002;5(3): 143–151. pmid:12088294

52. United States Census Bureau. 2010 census urban and rural classification and urban area criteria. 2013. Available: https://www.census.gov/geo/reference/ua/urban-rural-2010.html.

53. Deddens J, Petersen M. Approaches for estimating prevalence ratios. Occup Environ Med. 2008;65(7): 501–506.

54. Elixhauser A, Steiner C, Harris DR, Coffey RM. Comorbidity measures for use with administrative data. Med Care. 1998;36(1): 8–27. pmid:9431328

55. Houchens R, Chu B, Steiner C. Hierarchical modeling using HCUP data. HCUP Methods Series: Report# 2007–01. 2007. Available: http://www.hcup-us.ahrq.gov/reports/methods/2007_01.pdf.

56. Hess JJ, Saha S, Luber G. Summertime acute heat illness in U.S. emergency departments from 2006 through 2010: Analysis of a nationally representative sample. Environ Health Perspect. 2014;122(11): 1209–1215. pmid:24937159

57. Kovats RS, Hajat S, Wilkinson P. Contrasting patterns of mortality and hospital admissions during hot weather and heat waves in greater London, UK. Occup Environ Med. 2004;61(11): 893–898. pmid:15477282

58. Xu Z, Etzel RA, Su H, Huang C, Guo Y, Tong S. Impact of ambient temperature on children's health: A systematic review. Environ Res. 2012;117: 120–131 pmid:22831555

59. Oudin Åström D, Bertil F, Joacim R. Heat wave impact on morbidity and mortality in the elderly population: A review of recent studies. Maturitas. 2011; 69(2): 99–105. pmid:21477954

60. Kalkstein AJ, Sheridan SC. The social impacts of the heat–health watch/warning system in Phoenix, Arizona: Assessing the perceived risk and response of the public. Int J Biometeorol. 2007; 52(1): 43–55. pmid:17262221

61. Abrahamson V, Wolf J, Lorenzoni I, Fenn B, Kovats S, Wilkinson P, et al. Perceptions of heatwave risks to health: Interview-based study of older people in London and Norwich, UK. Journal of Public Health. 2009;31(1): 119–126. pmid:19052099

62. Wilby RL. Past and projected trends in London's urban heat island. Weather. 2003;58(7): 251–260.

63. Tan J, Zheng Y, Tang X, Guo C, Li L, Song G, et al. The urban heat island and its impact on heat waves and human health in Shanghai. Int J Biometeorol. 2010;54(1): 75–84. pmid:19727842

64. O’Neill MS, Zanobetti A, Schwartz J. Disparities by race in heat-related mortality in four US cities: The role of air conditioning prevalence. Journal of Urban Health. 2005;82(2): 191–197. pmid:15888640

65. Isaac M, Van Vuuren DP. Modeling global residential sector energy demand for heating and air conditioning in the context of climate change. Energy Policy. 2009;37(2): 507–521.

66. IPCC. (2014) Climate change 2014: Impacts, adaptation, and vulnerability. part A: Global and sectoral aspects. Contribution of working group II to the fifth assessment report of the intergovernmental panel on climate change [Field, C.B., V.R. Barros, D.J. Dokken, K.J. Mach, M.D. Mastrandrea, T.E. Bilir, M. Chatterjee, K.L. Ebi, Y.O. Estrada, R.C. Genova, B. Girma, E.S. Kissel, A.N. Levy, S. MacCracken, P.R. Mastrandrea, and L.L. White (eds.)]. Cambridge university press, Cambridge, United Kingdom and New York, NY, USA, 1132 pp.