The impact of cold and hot weather on mortality has been known for centuries. The seasonality of mortality in Sweden was studied in the 18th century using parish and national population statistics.1 Historical data for England reveal early evidence of increased mortality during severe winters, as well as after cold spells.2,3 More recently, seasonal fluctuations and heightened winter mortality, especially among the elderly, has been reported throughout Europe.4,5 Contemporary research has also focused on heat mortality especially among ageing populations.6 The frequency and duration of extreme heat events is predicted to increase.7 The 1987 heat wave in Athens resulted in about 2000 deaths,8 and the 2003 European heat wave caused some 70,000 deaths.9 Although the elderly are most susceptible to heat, the fatal effects of such events are also apparent among the middle-aged and the very youngest populations.10,11 Although not consistent, many studies suggest women face greater relative increases in mortality with heat waves.12,13
The means by which temperature extremes affect human physiology and performance are well documented.12 However, the extent and drivers of adaptation to weather and climate are not well understood. Physiologic adaptation to heat occurs within days to weeks in athletes.14 Social, societal, and behavioral adaptation can also modify a population’s vulnerability to heat and cold.15,16 Adaptation to regional ambient temperatures has been observed in a few studies of historic registers of short-term mortality caused by weather conditions.17–19 However, changes in the temperature–mortality relationship cannot be easily separated from underlying demographic changes in the population, changing prevalence of susceptibility factors (well known to have occurred during the 20th century),20 and extent of natural adaptation to the local environment.
We investigate whether sensitivity to temperature extremes of heat and cold changed in Sweden from the beginning of the 20th century to the present. We describe trends in relative risk (RR) of mortality related to temperature extremes, from all causes and stratified by sex and age.
We collected information on the daily numbers of deaths from all causes during the period 1901–2009 for present-day Stockholm County, Sweden (Table 1).21 These mortality data are complete only since 1947; from 1901 through 1946 approximately 46% of original parish records have been converted to digital format. Given that mortality distributions among parishes are similar, the available data for deaths should serve as a valid proxy for fluctuations in daily total mortality of Stockholm County.
Daily mean temperature data were provided by the Swedish Meteorological and Hydrological Institute (Table 2). Observations were made at the Stockholm Observatory, where temperature data have been recorded at almost exactly the same location since 1756. We defined heat extremes as days for which the 2-day moving average of mean temperature (“lag01”) was above the 98th percentile (20.8°C) for the entire period (1901–2009). Cold extremes were defined as days for which the 26-day moving average of mean temperature (“lag025”) was below the 2nd percentile (−6.9°C). These lags are sufficient to describe effects related to short-term variability in temperature.22 Using the same temperature criteria throughout the study period will not take into account either the increasing temperature over time or acclimatization within decades. A single criterion could lead to bias in the earlier decades due to the fact that we might be estimating mortality for a more extreme relative temperature. To address this problem, we also investigated temperature extremes occurring within each decade, using the same percentiles and lags for heat and cold extremes as defined above.
To analyze the relationship between extreme hot or cold temperatures and mortality, we used an approach that assumed the daily counts of mortality follow an overdispersed Poisson distribution. A generalized linear model was fit. Time trends in the daily number of deaths were described by a smooth function. This function was allowed 400 degrees of freedom (df) over the period, corresponding to around 4 df per year of data. The time-trend function allows baseline mortality to vary depending on demographic changes, changes in the data-collection system, and other slowly time-varying extraneous factors. However, the time-trend function cannot explain the short-term variability of daily deaths, which therefore must be explained by daily temperature fluctuations.
Day of the week was included in the model as a categorical explanatory variable. Binary variables indicated public Swedish holidays (“holiday”) and the later period when all parishes are included in register data (“from 1947 onwards”) (“complete”). Pandemic episodes of influence were included as binary variables (the Spanish flu, October 1918–November 1920; the Asian flu, October 1957–April 1958; and the Hong Kong flu, October 1969–March 1970).23
We established the model:
Model 1 was computed using cutoff temperatures for temperature extremes defined for the entire period of investigation 1901–2009, whereas for model 2 the cutoff temperatures for temperature extremes were defined within each decade.
Sensitivity analyses of the estimated heat and cold coefficients were done using 300, 350, 450, and 500 degrees of freedom to investigate the importance of the smoothing parameter. The estimates of the RRs were quite similar for the various degrees of freedom used. See eTable 1 (http://links.lww.com/EDE/A724) for the estimates of the RRs and the corresponding 95% confidence intervals (CIs) using the various degrees of freedom.
We analyzed separately to investigate differences in vulnerability to extreme temperatures between sexes.12,13 Since age has also been shown to modify susceptibility,6,12,13 we analyzed data separately for various age groups. These age groups were 0–14, 14–64, and 65+ years. These age groups are characterized by different positions in society, as well as different rates and causes of death. Due to the considerable change in age of mortality over the 110 years under study, we further separated the elderly into two categories: 65–79 and 80+.
The relative effect modification index was used to investigate whether sex and age modified the mortality effects of exposure to extreme temperatures. This index is calculated as the ratio between the RR of the group of interest and the RR from the reference group.24
To investigate time trends in the RRs, bearing in mind the substantial difference in standard errors of the RR estimates over time, we performed a weighted least squares linear regression, with the RR estimates as the dependent variable and decade as the independent variable, using the weights inversely proportional to the variance of each RR. SAS version 9.2 software (SAS Institute, Cary, NC) was used to create datasets and variables. R version 2.13.1 was used for statistical models and creation of outputs.
Number of Temperature Extremes
Figure 1 presents the number of heat and cold extremes during 1901–2009, by decade, using cutoff temperatures for temperature extremes defined for the entire period of investigation 1901–2009.
Heat extremes increased over the 20th century, with the largest upsurge from 1990 onward. The last two decades had more than half-again as many extreme heat events as any other decade of the period. The number of cold extremes per decade was much more variable, and no trend is apparent. This number during the last two decades was low, as in the 1930s when no cold extremes occurred. During the 1940s and 1980s there were large numbers of cold extremes.
RR of Mortality During Hot and Cold Weather
Figure 2 shows the RR (95% CI) of mortality associated with extreme heat (lag01) and extreme cold (lag025) compared with nonextreme days, by decade. Table 3 presents the results of the regression analyses stratified for sex and age. eTable 2 (http://links.lww.com/EDE/A724) presents the results of the regression analysis of trend over time for all investigated groups. Models 1 and 2 yielded similar risk and trend estimates; thus, if not specifically mentioned below, the results and conclusions apply to both models.
Mortality was increased with heat extremes in all decades, as well as for the entire period. The estimated RR of mortality during heat extremes over the entire period was 1.10 (95% CI = 1.08–1.13) and 1.12 (1.10–1.14) for models 1 and 2, respectively. There was a substantial decline of RR estimates for heat extremes during the period.
Although mortality with cold extremes was also increased overall (1.06 [1.03–1.08] and 1.04 [1.02–1.05] for models 1 and 2, respectively), the association was more variable across decades. During one period 1950–1959, the RR dropped below 1.00. There was little evidence of a declining trend in cold-related mortality over time.
The RRs associated with heat extremes were generally higher than those associated with cold extremes, with both sides risks stabilizing in the last three decades. Patterns and trends were similar for men and women separately (eFigures 1 and 2, http://links.lww.com/EDE/A724). Men had increased mortality during heat extremes in all but 1980–1989. Men’s RR of mortality with heat declined over time. Risk with cold also declined over time, although the trend was less marked. Results were similar for women, with a decrease in morality RR with heat over time, but no apparent change in mortality risk with extreme cold. Sex was in general not a modifier of risk, for either heat or cold extremes.
RRs for the various age groups are presented in eFigures 3–7 (http://links.lww.com/EDE/A724). An increase in RR for those above age 65 years is present with both heat and cold extremes across the entire period. Not surprisingly, this pattern is similar to overall mortality in the last three decades, when over-65 mortality accounted for eight of ten deaths.
Results for the other age groups are more difficult to interpret, given the small number of deaths in these groups. Risks were elevated during heat extremes for both the age 0–14 and 15–64 years populations in the early 20th century, and then tended to decline. Toward the end of the study period, there was an upward risk trend in the 15–64 age group, similar to another study from this region, which indicated a recent increased risk of death among persons aged 45–64 years.11
The trends over time differ somewhat in the various age groups. Although RR of mortality from heat was consistently higher than for cold, there was a significant declining trend for heat and cold extremes in the 0–14 age group but not the 15–64 age group. Time trends among the elderly followed that of total mortality, with a significant declining trend over time for heat extremes but not for cold extremes. We found no strong evidence of age being an effect modifier for changes over time.
Mortality data for the past century in Stockholm, Sweden, show a consistent elevation in RR of death after extremes of temperature, with greater risks for extremes of heat than of cold. Over time there has been a generally declining trend of mortality RR with temperature extremes, although risks remain elevated for the elderly. Age and sex are not consistent effect modifiers. Because women tend to live longer than men, the modifying effect in older age groups is hard to disentangle from the effect of sex. There are even indications that, beginning one to three decades ago, the population of Stockholm County may be becoming more susceptible to heat extremes. In particular, men and the middle-aged appear to contribute to this recent upward trend.
Recent changes are occurring at a time when the number of cold extremes is decreasing and heat extremes are becoming more frequent and intense. A decrease in cold extremes may be contributing to a shift in susceptible populations, such that mortality formerly caused by cold is more likely to occur with heat.25 A doubling of susceptibility to heat-related mortality in the last two decades has been reported in this region.25
Patterns of temperature-related mortality over a longer period may shed light on results from studies focused on shorter periods. We found reduced vulnerability to heat- and cold-related mortality over time. This decline may be attributable to better housing and health systems, behavioral and physiologic adaptation, and heightened public awareness of dangers associated with extreme temperatures. Given that Stockholm is situated in a cold climate, it is likely that many of the adaptive measures over time have aimed at reducing the adverse effects of cold temperatures.
The past 100 years has been a time of fundamental changes in both public health and mortality. Life expectancy at the national level increased from about 56 years during the first decade of the 20th century to more than 81 years in 2009.26 The epidemiologic transition20 changed mortality patterns from death primarily in childhood to death in old age. This development is connected to changes in the causes of death. The dominant cause of mortality in early industrial periods was infectious disease, often of epidemic character, whereas chronic diseases are presently responsible for most deaths.
Our results are similar to those of Carson et al17 and Ekamper et al19 who reported declining vulnerability to temperature-related mortality in England and the Netherlands over the 20th century. Davis et al18 investigated whether annual heat-related mortality rates had changed in the United States during 1964–1998. They detected a decline in heat-related vulnerability, which they attributed mainly to improvements in the health sector, population adaptation, improvements in urban design, and increased use of air conditioning. Lerchl27 found declining seasonality of mortality in Germany during 1946–1995 and concluded that this was mainly due to increased use of central heating and improvements in the public health sector.
We report results from two models—model 1 producing estimates with fixed criteria for extreme temperatures over time, and model 2 using decade-specific criteria. The results of these models should be interpreted somewhat differently. Model 2 defines temperature extremes relative to the prevailing climate, thus assuming adaptation to potentially increasing temperatures and accounting for changes in temperature extremes over time. Even though both models show increased risk of mortality during heat extremes, model 1 generally reports higher estimates than model 2 in the earlier part of the century whereas model 2 shows higher estimates in the later part of the century. This contrast is presumably due to increasing temperatures over time.
One limitation to the study is that neither intensity nor duration of the temperature extremes were taken into account when estimating the RRs of mortality. Barnett et al28found that more intense heat waves increased the risk of mortality, whereas the same was not true for cold spells. Gasparrini and Armstrong29 found heat waves lasting more than 4 days had a small effect on mortality (compared with the effect of temperature itself). Due to data availability, we report only results for all-cause mortality. Given the substantial changes in causes of mortality over the study period, and given that temperature–mortality relationships differ among vulnerable groups, it would be of interest to be able to present results for cause-specific mortality.
Previous studies have suggested that adaptation to weather has occurred over the 20th century. However, the changes observed in this and previous studies cannot be attributed to adaptation alone. None of these investigations has controlled for variables independent of adaptation, which potentially modify temperature effects and their change over time (eg, causes of death, demographic factors, and prevalence of medical conditions). Unexpectedly, we observed a recent leveling of the decreasing trend of heat susceptibility, with some groups experiencing an upward trend in temperature susceptibility in recent decades.
In conclusion, we found that susceptibility to heat extremes has decreased substantially over the last 110 years in Stockholm, Sweden. A decreasing susceptibility to cold extremes was also found, although with more variability. This apparent progress over time may indicate improvements in adaptation and general conditions, although recent increases in risk indicate possible increasing susceptibility to temperature extremes, particularly heat. To the extent that climate change brings increased episodes of extreme heat, this poses a challenge to societies and to the health sector, especially in protecting susceptible populations.
1. Wargentin P Uti hvilka månader flera människor årligen födas och dö i Sverige 1767 K. Vetenskapsac Handl Avhandling 43. Stockholm, Sweden: Kungl. Vetenskapsakademien; 1767
2. Guy WA. On temperature and its relation to mortality: an illustration of the application of the numerical method to the discovery of truth. J Statistical Soc. 1881;44:235–262
3. Farr W Vital Statistics: A Memorial Volume of Selections from the Reports and Writings of William Farr. Metuchen, NJ Scarecrow Press; 1885
4. Wilkinson P, Landon M, Armstrong B, et al. Cold Comfort: The Social and Environmental Determinants of Excess Winter Death in England, 1986-1996. 2001 York Joseph Rowntree Foundation
5. Donaldson GC, Keatinge WR. Excess winter mortality: influenza or cold stress? Observational study. BMJ. 2002;324:89–90
6. Åström DO, Forsberg B, Rocklöv J. Heat wave impact on morbidity and mortality in the elderly population: a review of recent studies. Maturitas. 2011;69:99–105
7. IPCC. Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change, Field CB, Barros V, Stocker TF, et al. 2012 Cambridge, UK, and New York, NY Cambridge University Press:582
8. Katsouyanni K, Trichopoulos D, Zavitsanos X, Touloumi G. The 1987 Athens heatwave. Lancet. 1988;2:573
9. Robine JM, Cheung SL, Le Roy S, et al. Death toll exceeded 70,000 in Europe during the summer of 2003. C R Biol. 2008;331:171–178
10. Fouillet A, Rey G, Laurent F, et al. Excess mortality related to the August 2003 heat wave in France. Int Arch Occup Environ Health. 2006;80:16–24
11. Rocklöv J, Ebi K, Forsberg B. Mortality related to temperature and persistent extreme temperatures: a study of cause-specific and age-stratified mortality. Occup Environ Med. 2011;68:531–536
12. Kovats RS, Hajat S. Heat stress and public health: a critical review. Annu Rev Public Health. 2008;29:41–55
13. Basu R. High ambient temperature and mortality: a review of epidemiologic studies from 2001 to 2008. Environ Health. 2009;8:40
14. Parsons KC. The effect of gender, acclimation state, the opportunity to adjust clothing and physical disability on requirements for thermal comfort. Energ Buildings. 2002;34:593–599
15. Moseley PL. Mechanisms of heat adaptation: thermotolerance and acclimatization. J Lab Clin Med. 1994;123:48–52
16. Stafoggia M, Forastiere F, Agostini D, et al. Vulnerability to heat-related mortality: a multicity, population-based, case-crossover analysis. Epidemiology. 2006;17:315–323
17. Carson C, Hajat S, Armstrong B, Wilkinson P. Declining vulnerability to temperature-related mortality in London over the 20th century. Am J Epidemiol. 2006;164:77–84
18. Davis RE, Knappenberger PC, Michaels PJ, Novicoff WM. Changing heat-related mortality in the United States. Environ Health Perspect. 2003;111:1712–1718
19. Ekamper P, van Poppel F, van Duin C, Garssen J. 150 years of temperature-related excess mortality in the Netherlands. Demographic Research. 2009;14:385–426
20. Omran AR. The epidemiologic transition. A theory of the epidemiology of population change. Milbank Mem Fund Q. 1971;49:509–538
21. . Sveriges Dödsbok Nr.5 1901-2009. Sweden: Sveriges Släktforskarförbund; 2010.
22. Anderson BG, Bell ML. Weather-related mortality: how heat, cold, and heat waves affect mortality in the United States. Epidemiology. 2009;20:205–213
23. Elgh F. Pandemiers påverkan på samhället. Läkartidningen. 2007;104:615–619
24. Stafoggia M, Forastiere F, Agostini D, et al. Vulnerability to heat-related mortality: a multicity, population-based, case-crossover analysis. Epidemiology. 2006;17:315–323
25. Rocklöv J, Forsberg B, Meister K. Winter mortality modifies the heat-mortality association the following summer. Eur Respir J. 2009;33:245–251
26. Historisk. Statistik för Sverige.Del 1. Befolkning. Statistics Sweden. Life Expectancy 1751–2011. 19692nd ed Stockholm SCB
27. Lerchl A. Changes in the seasonality of mortality in Germany from 1946 to 1995: the role of temperature. Int J Biometeorol. 1998;42:84–88
28. Barnett AG, Hajat S, Gasparrini A, Rocklöv J. Cold and heat waves in the United States. Environ Res. 2012;112:218–224
29. Gasparrini A, Armstrong B. The impact of heat waves on mortality. Epidemiology. 2011;22:68–73
Supplemental Digital Content
Copyright © 2013 Wolters Kluwer Health, Inc. All rights reserved.