where c = estimated effect of temperature on mortality in community c; β c = true effect of temperature on mortality in community c; JOURNAL/epide/04.02/00001648-200903000-00011/ENTITY_OV0450/v/2017-07-26T080127Z/r/image-png c = statistical variance of c; μ = true average effect across all communities; τ 2 = between-community variance of the true effect; n = number of communities.
Second-stage analysis evaluated whether community-specific variables (mean yearly and seasonal temperature and unadjusted dew point temperature) modified community-specific temperature-mortality effects.
where xj c = community-level variable j for community c: j = mean community-level variable j across communities; α 0 = average ln(relative rate) when xj c=j α1,j = change in ln(relative rate) for unit increase in (xj c − j).
Similar analysis examined how community-specific variables affect ozone mortality effects.29
Equation 4 was used to explore sensitivity of subpopulations by race, socioeconomics, urbanicity, and AC prevalence. We applied stratified models (Equations 2 and 3) to explore susceptibility by age and cause of death (respiratory, cardiovascular, and noncardiorespiratory). All analysis was performed in R-2.6.2.
Weather and mortality summary statistics are provided in eTable 1 and distributions in eFigure 1 in the supplemental material (http://links.lww.com/A720). We first calculated relative effects with each temperature metric: minimum, maximum, and mean temperature, and apparent temperature (eFig. 2, http://links.lww.com/A720). All metrics were strongly correlated (eTable 2, http://links.lww.com/A720). Mean daily temperature and apparent temperature estimates were similar. Mean daily temperature was chosen for subsequent analysis given its high correlation to other metrics as this measure provides more easily interpreted results in a policy context.
We considered lag structures of the same day and the average of the same day and up to 28 days previous. Figure 1 presents relative effect estimates and shows how this estimate changes as more days of lagged temperature are included. First the relative effects were calculated separately for each community and then combined to generate an overall effect. Although Figure 1 shows the effect of lag structure on national heat and cold effects, the lag with the strongest effect varied by community. eFigure 3 (http://links.lww.com/A720) shows these effects for the 4 largest communities. The association of heat with mortality was limited to recent days of exposure and estimates declined when longer time periods were considered. For cold-related mortality a longer exposure timeframe is relevant. Results indicate that use of identical exposure timeframes for heat and cold response is inappropriate. For subsequent analysis, we estimated heat-related mortality based on the same and previous day (Tlag0–1), and cold-related mortality based on the same day and past 25 days (Tlag0–25).
Based on the selected temperature metric and lag structures, we generated community-specific exposure–response curves. Figure 2 and eFigure 4 (http://links.lww.com/A720) provide an example, displayed as the increase in mortality risk for a given temperature (x-axis) compared with a reference temperature (60°F). eFigure 4A (http://links.lww.com/A720) shows results for a shorter exposure period (Tlag0–1), capturing the heat effect, and eFigure 4B (http://links.lww.com/A720) for a longer lag (Tlag0–25), capturing the cold effect. Figure 2 combines the heat- and cold-related portions of eFigures 4A and 4B (http://links.lww.com/A720), respectively. eFigure 5 (http://links.lww.com/A720) provides examples of other locations. This representation better captures weather-related mortality than previous research that used equal lag structures for cold and heat effects.
Slopes of exposure–response curves were summarized by comparing the risk of relative and absolute temperature changes. eFigure 6 (http://links.lww.com/A720) shows relative cold and heat effects for each community and the national effect, and eFigure 7 (http://links.lww.com/A720) provides absolute effects. The overall increase in mortality risk comparing the first and 10th percentile Tlag0–25 was 4.2% (95% posterior interval = 3.2%–5.3%). Mortality risk increased 3.0% (2.4%–3.6%) comparing the 99th and 90th percentile Tlag0–1. Results were robust to the degrees of freedom used in time splines (eTable 3, http://links.lww.com/A720). For absolute temperature changes, mortality risk at 60°F was 5.2% (3.8%–6.6%) lower than at 40°F for Tlag0–25 and 4.9% (3.8%–6.0%) lower than at 80°F for Tlag0–1.
We repeated analysis with inclusion of ozone or PM10 (Table 1 and eTable 4, http://links.lww.com/A720). A smaller dataset is available for this analysis due to the lack of pollution data for some communities and the frequency of pollution measurement. Cold effects were similar with pollution adjustment, whereas heat effects were slightly lower. eFigure 8 shows community-specific and overall results of pollution sensitivity analysis for relative heat effects (http://links.lww.com/A720).
eTable 5 (http://links.lww.com/A720) shows summary statistics for the occurrence of heat waves in this study. Table 2 presents results from the heat-wave analysis. Some communities had no heat waves of a specific definition and were thereby excluded from analysis. Mortality effects increased with intensity or duration of the heat wave. As previous research found lower heat-wave effects when longer lags were used for temperature control,26 we recalculated results with lag 0 to 2 temperature control (eTable 6, http://links.lww.com/A720). Effect estimates were very similar using this model structure.
Geographic distributions of heat and cold effects for relative (comparing risk across temperature percentiles) and absolute (comparing risk at specific temperatures) effects are shown in Figure 3 and eFigure 9 (http://links.lww.com/A720). Cold effects appear to be larger in the South than in the north (eFigs. 9A, B, http://links.lww.com/A720). Conversely, heat effects generally appear larger in the north (Fig. 3 and eFig. 9C, http://links.lww.com/A720). Heat wave and heat effects showed similar geographic patterns (eFig. 10, http://links.lww.com/A720). Spatial variation was also demonstrated by regional estimates (eFig. 11, http://links.lww.com/A720). Cold effects were more similar across regions than heat effects. Regional trends are more defined for absolute rather than relative effects (Fig. 3 and eFigs. 9, 11, http://links.lww.com/A720), which could indicate that populations acclimate to a city's weather conditions, especially heat.
We evaluated whether variation in community-specific effects could be explained by long-term average temperature and dew point temperature (Table 3). Heat effects were higher in colder communities and cold effects higher in warmer communities, consistent with the observed regional patterns. For example, a 12.9°F increase in long-term temperature was associated with a 70% decrease in absolute heat effects and a 94% increase in absolute cold effects.
Figure 4 shows relative cold and heat effects, stratified by age and cause. Respiratory mortality effects were generally higher than cardiovascular effects. Results also indicate associations for noncardiorespiratory mortality. Heat and cold effects were highest for the oldest age category (≥75 years) for all causes; however, associations were also observed for youngest age group (<65).
Heat wave effects were estimated by age and cause of death with the heat-wave definition as temperatures at the 99.5th percentile or higher and 2 days duration or longer (Table 4). Heat waves had effects on all age groups, with the largest effect for the oldest group, and for cardiovascular and noncardiorespiratory deaths, with the highest estimate for cardiovascular deaths. Results indicate an association between heat waves and respiratory mortality, although estimates are uncertain.
To investigate sensitive subpopulations, we explored the relationship between community-specific weather-mortality effects and various community-specific indicators (Table 5). Many of these community-level variables were correlated (eTable 7, http://links.lww.com/A720). Communities with higher income, unemployment, population, and urbanicity were more susceptible to heat impacts. Higher susceptibility to cold was identified for communities with a higher percentage of African Americans. A higher fraction of homes with central AC was associated with lower heat-related and higher cold-related mortality risk. Results for the 65 years and older age category were similar to results for all ages (eTables 8, 9, http://links.lww.com/A720).
A variety of modeling decisions are made when estimating the impact of weather on mortality, including the shape of the exposure–response curve, lag structure, and temperature metric. These choices affect results and comparability across studies. Earlier work investigated several exposure–response forms, such as the temperature of lowest mortality risk (minimum mortality temperature) and constant linear relationships above (heat slope) and below (cold slope) the minimum mortality temperature,6 generating a V-shaped exposure–response curve. Although useful, these methods do not fully capture the nonlinear association and are problematic for comparing across cities, for example, a comparison of a heat slope calculated for higher than 80°F versus a slope calculated for higher than 90°F. Other studies estimated constant slopes above and below city-specific threshold hot and cold temperatures.10,11,16 When comparing communities with disparate climates, this method forces a V-shaped model that may not reflect actual temperature-mortality relationships in each community.
Our spline approach allows estimation of nonlinear relationships without forcing constant slopes for specific temperature ranges or similarities among communities. Similar methods were applied in a study of 7 US cities.2 For most communities the difference in mortality risk per unit temperature decrease was fairly consistent across mild cold temperatures; however, the heat effect per unit temperature increase rose significantly at higher temperatures. Furthermore, some communities did not have a unique minimum mortality temperature. Slope approximations based on specified temperatures (eg, our absolute estimates) are useful to summarize and compare temperature-response relationships, but interpretation of results should consider that such methods reflect only a portion of the nonlinear relationship. Although the heat and cold effects generated in this study summarize only part of the temperature–mortality relationship, these methods were particularly appropriate for studying such a large range of climates. Complex nonlinear functions such as those used here and in previous studies may provide a more complete assessment of temperature and mortality risk.8,15
Previous studies of weather and mortality have used a variety of temperature measurements. Several studies have recommended apparent temperature or humidex24 because these measures incorporate humidity; others suggested minimum temperature.11 We applied minimum, maximum, and mean daily temperature, and mean daily apparent temperature, and identified heat and cold effects for all metrics. Apparent temperature effects were nearly identical to those of mean daily temperature adjusted for humidity. Heat effects were highest for mean daily temperature. A multicity European study similarly found that mean daily temperature was consistently the strongest predictor of mortality from heat and heat waves,26 compared with daily minimum, maximum, and apparent temperature. Although some differences in estimates may occur, our findings indicate that the various temperature metrics are likely to produce similar results.
We investigated lag times from same day to 28 days previous. Earlier heat-mortality studies identified risk from recent exposure (ie, same day and a few days previous).5–8,11,16 Most studies applied lags of 1 or 2 days, although some used up to 3 days.30 We found the strongest heat-related mortality association for same- and previous-day exposure. The short lag required to capture the effects of heat on mortality suggests a rapid physical response. Some of the effects observed could be the result of short-term mortality displacement, and further study is warranted.
For cold-related mortality, most US studies applied 2- to 5-day lags,1,5,6,11 whereas other researchers found cold effects after 1 or more weeks for some communities.16,31 Findings indicate that longer lags are required to capture cold's impact on mortality and that using identical lag structures for cold and heat effects is not appropriate. A limitation of longer lag structures is the introduction of more measurement error due to increased time between the exposure and event. Heat and cold effects were similar in magnitude for absolute and relative estimates, which contrasts with earlier US studies finding larger heat effects than cold effects.5,11 We hypothesize that previous studies underestimated cold-related effects through use of shorter lags. Results agree with a European study finding mortality effects occurring days to weeks after cold exposure.16 Findings suggest that cold temperatures more indirectly affect mortality than heat. Infectious diseases, which are more common in industrialized countries during colder weather (when people spend more time indoors and in proximity) could account for much of the cold-related effect. Although we found that heat effects were impacted by shorter exposures and cold effects were affected by longer exposures, the specific lag structures used here (Tlag0–25 and Tlag0–1) are intended to be representative, not to reflect the only or the exact lag measurements appropriate for temperature-mortality studies.
We took several approaches to comparing temperature effects across communities, including estimates based on each community's temperature distribution, allowing comparison despite the wide climatic range. Studies of 2 European,16 2 Mexican,30 and 50 US5 cities similarly estimated mortality risk for community-specific temperature quantiles. We also calculated the effect of a change in absolute temperature, from a relatively mild (60°F) to hot (80°F) or cold (40°F) temperatures. These 2 types of effect estimates have different interpretations with respect to acclimatization, an important consideration for climate-change studies and public health policy. Acclimatization can occur through physical adaptation, housing characteristics, or behavioral patterns (eg, staying indoors, clothing). With a high degree of acclimatization to weather, results would be similar across communities for relative temperature effects and different for absolute effects. Without a high degree of acclimatization, communities would have similar absolute effects and dissimilar relative effects. Although both absolute and relative temperature effect estimates showed variation across communities, absolute estimates exhibited larger variation, which implies some degree of acclimatization to weather conditions because a given temperature has different impact depending on location. Previous studies have also found some evidence of acclimatization.5,6
Heat effects were generally lower in communities with higher long-term temperatures. This supports the hypothesis that communities and individuals adapt, to some extent, to weather even during temperatures that are extremely warm for that area. Conversely, absolute cold effects were markedly higher in communities with higher temperatures, as observed previously in other areas.32,33 However, a similar association was not observed for relative cold effects. This indicates that those in colder cities seem to acclimatize to some degree, so they are less affected by temperatures of 40°F, but not to the extent of lessening effects at temperatures extremely cold for the community.
Most previous heat wave studies have analyzed specific extreme events (eg, Chicago 199534 or European 200312,13 heat waves). Fewer studies have considered more frequent, less severe heat waves. Our findings suggest that sustained periods of extreme heat present an elevated risk over single days of high temperatures, even if the heat wave period is as short as 2 days, and that duration and intensity of the heat wave affect mortality risk. Future studies might consider separate effects by heat wave duration, intensity, or time of occurrence during the summer.
We identified spatial heterogeneity in heat and cold effects, consistent with other US studies1,2,5–8,10,11 with larger cold effects in the southern US and smaller effects in the north.6,11 Similar to previous research,7,11 we found negligible or null effects for heat in many southeastern communities. Results emphasize the need for multicity studies because results from 1 location may not be applicable elsewhere.
Heat effects were slightly lowered when models included O3 and PM10. Cold effects were essentially unchanged. In previous studies, temperature–mortality results were also robust to particles2,16,26,30 and ozone.5,26,30,35 Earlier research found ozone-mortality associations were robust to control for temperature.20 Findings of previous studies and our work imply separate and substantial mortality effects from temperature and from air pollution; however, some studies suggest possible interaction between temperature and air pollution.3,36 Observed associations between weather and noncardiorespiratory deaths indicate that weather affects mortality beyond cardiorespiratory responses. Estimates were somewhat higher for cardiovascular and respiratory deaths, especially respiratory, compared with total deaths, consistent with earlier studies.5,37
We found higher susceptibility for older populations; however, other age categories were also subject to temperature–mortality risk. Results from stratifying effect estimates by age were consistent with earlier results based on community-level and individual-level data.11,34,37,38
We found differences in susceptibility related to socioeconomic factors and urbanicity. Previous studies also found community-level socioeconomic factors to explain some variability in communities' heat effects.6,10 This may reflect baseline health and nutrition status, access to health care, and ability to respond to extreme conditions (eg, AC). Susceptibility by urbanicity may relate to urban-heat-island effects or housing conditions. We found that heat has a lower mortality impact when communities have more central AC, as observed in smaller US studies.7,10,14 This adaptation strategy likely explains some of the regional variation in heat effects. Earlier work1,7 found that heat-related mortality decreased significantly in the southeastern US as AC prevalence increased. Over time, the number of cities without a heat effect has increased,7 especially where AC has reached almost universal prevalence. Heat wave effects, however, were not strongly associated with AC. It is possible that the protection afforded by AC is sufficient to reduce effects of high temperatures, but not to prevent more extreme heat wave effects. Heat effects were higher in communities with higher income. Although this relationship might seem surprising, a recent national study showed that median income is negatively associated with mortality in adults 65 years or younger but not with older individuals, which represents the majority of heat- and cold-related mortality.39
Socioeconomic factors were less important in explaining cold effects; however, communities with a higher percentage of black/African Americans had higher cold effects, although the relationship was uncertain. Earlier research also found insignificant impact of socioeconomic factors on cold effects,5,6 although 1 study found higher cold effects for communities with less education or a lower percentage of population identifying as black.2 A study of cold effects on the elderly found significant associations with poverty, income inequality, and deprivation rate.40 The lack of association in our study might relate to the use of community-level variables. Further investigation using individual-level data is needed. Such data could also improve exposure estimates, especially for longer lag structures.
These findings on the impact of weather on mortality have implications for policymakers and future scientific work. The identified susceptible subpopulations signify the need for targeted heat-mortality prevention efforts. The heterogeneous results across communities indicate the value of multicity research and indicate that approaches to prevent weather-related mortality might be most effective if they are community specific. Results on acclimatization and heat waves are of particular importance to research estimating weather-related mortality impacts from climate change.
1. Barnett AG. Temperature and cardiovascular deaths in the US elderly: changes over time. Epidemiology
2. O'Neill MS, Zanobetti A, Schwartz J. Modifiers of the temperature and mortality association in seven US cities. Am J Epidemiol
3. Ren C, Williams GM, Morawska L, et al. Ozone modifies associations between temperature and cardiovascular mortality analysis of the NMMAPS data. Occup Environ Med
4. Medina-Ramon M, Zanobetti A, Cavanagh DP, et al. 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
5. Medina-Ramon M, Schwartz J. Temperature, temperature extremes, and mortality: a study of acclimatization and effect modification in 50 US cities. Occup Environ Med
6. Curriero FC, Heiner KS, Samet JM, et al. Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol
7. Davis RE, Knappenberger PC, Michaels PJ, et al. Changing heat-related mortality in the United States. Environ Health Perspect
8. Braga AL, Zanobetti A, Schwartz J. The effect of weather on respiratory and cardiovascular deaths in 12 US cities. Environ Health Perspect
9. Basu R, Samet JM. Relation between elevated ambient temperature and mortality: a review of the epidemiologic evidence. Epidemiol Rev
10. Chestnut LG, Breffle WS, Smith JB, et al. Analysis of differences in hot-weather-related mortality across 44 US metropolitan areas. Environ Sci Policy
11. Kalkstein LS, Davis RE. Weather and human mortality: an evaluation of demographic and interregional responses in the United States. Ann Assoc Am Geogr
12. Filleul L, Cassadou S, Medina S, et al. The relation between temperature, ozone, and mortality in nine French cities during the heat wave of 2003. Environ Health Perspect
13. Le Tertre A, Lefranc A, Eilstein D, et al. Impact of the 2003 heatwave on all-cause mortality in 9 French cities. Epidemiology
14. O'Neill MS, Zanobetti A, Schwartz J. Disparities by race in heat-related mortality in four US Cities: the role of air conditioning prevalence. J Urban Health
15. Armstrong B. Models for the relationship between ambient temperature and daily mortality. Epidemiology
16. Pattenden S, Nikiforov B, Armstrong BG. Mortality and temperature in Sofia and London. J Epidemiol Community Health
17. Souch C, Grimmond CSB. Applied climatology: ‘heat waves.' Prog Phys Geogr
18. Seinfeld JH, Pandis SN. Atmospheric Chemistry and Physics: From Air Pollution to Climate Change.
New York: Wiley-Interscience; 1998.
19. Internet-based Health and Air Pollution Surveillance System (iHAPSS). 2007. Mortality, air pollution, and meteorological data for 108 US cities 1987–2000. Available at: http://www.ihapss.jhsph.edu/H
. Accessed November 20, 2007.
20. Bell ML, McDermott A, Zeger SL, et al. Ozone and short-term mortality in 95 US urban communities, 1987–2000. J Am Med Assoc
21. US Census Bureau. Census 1990, Summary File 1; Census 1990, Summary File 3; Census 2000, Summary File 1; Census 2000, Summary File 3
. Washington, DC: US Census Bureau.
22. US Department of Commerce, US Department of Housing and Urban Development. American Housing Survey, 1989, 1991, 1993, 1994, 1997, 1999
. Washington, DC: US Department of Commerce, US Department of Housing and Urban Development.
23. Dominici F, Peng RD, Bell ML, et al. Fine particulate air pollution and hospital admission for cardiovascular and repiratory diseases. J Am Med Assoc
24. Kalkstein LS, Valimont KM. An evaluation of summer discomfort in the United States using a relative climatological index. Bull Am Meteorol Soc
25. Samet JM, Dominici F, Zeger SL, et al. The National Morbidity, Mortality, and Air Pollution Study Part I: Methods and Methodologic Issues
. Health Effects Institute: Cambridge, MA; 2000.
26. Hajat S, Armstrong B, Baccini M, et al. Impact of high temperatures on mortality: is there an added heat wave effect? Epidemiology
27. Everson PJ, Morris CN. Inference for multivariate normal hierarchical models. J R Stat Soc Ser B
28. Kass RE, Wasserman L. The selection of prior distributions by formal rules. J Am Stat Assoc
29. Bell ML, Dominici F. Effects modification by community characteristics on the short-term effects of ozone exposure and mortality in 98 US communities. Am J Epidemiol
30. O'Neill MS, Hajat S, Zanobetti A, et al. Impact of control for air pollution and respiratory epidemics on the estimated associations of temperature and daily mortality. Int J Biometeorol
31. Braga AL, Zanobetti A, Schwartz J. The time course of weather-related deaths. Epidemiology
32. The Eurowinter Group. Cold exposure and winter mortality from ischaemic heart disease, cerebrovascular disease, respiratory disease, and all causes in warm and cold regions of Europe. Lancet
33. Barnett AG, Dobson AJ, McElduff P, et al. Cold periods and coronary events: an analysis of populations worldwide. J Epidemiol Community Health
34. Kaiser R, Le Tertre A, Schwartz J, et al. The effect of the 1995 heat wave in Chicago on all-cause and cause-specific mortality. Am J Public Health
. 2007;97(suppl 1):S158–S162.
35. Rainham DGC, Smoyer-Tomic KE. The role of air pollution in the relationship between a heat stress index and human mortality in Toronto. Environ Res
36. Nawrot TS, Torfs R, Fierens F, et al. Stronger associations between daily mortality and fine particulate air pollution in summer than in winter: evidence from a heavily polluted region in western Europe. J Epidemiol Community Health
37. Aylin P, Morris S, Wakefield J, et al. Temperature, housing deprivation and their relationship to excess winter mortality in Great Britain, 1986–1996. Int J Epidemiol
38. Stafoggia M, Forastiere F, Agostini D, et al. Vulnerability to heat-related mortality: a multicity, population-based, case-crossover analysis. Epidemiology
39. Backlund E, Rowe G, Lynch J, et al. Income inequality and mortality: a multilevel prospective study of 521 248 individuals in 50 US states. Int J Epidemiol
40. Healy JD. Excess winter mortality in Europe: a cross country analysis identifying key risk factors. J Epidemiol Community Health
Supplemental Digital Content
© 2009 Lippincott Williams & Wilkins, Inc.