An ecological time-series study of heat-related mortality in three European cities
© Ishigami et al; licensee BioMed Central Ltd. 2008
Received: 20 June 2007
Accepted: 28 January 2008
Published: 28 January 2008
Europe has experienced warmer summers in the past two decades and there is a need to describe the determinants of heat-related mortality to better inform public health activities during hot weather. We investigated the effect of high temperatures on daily mortality in three cities in Europe (Budapest, London, and Milan), using a standard approach.
An ecological time-series study of daily mortality was conducted in three cities using Poisson generalized linear models allowing for over-dispersion. Secular trends in mortality and seasonal confounding factors were controlled for using cubic smoothing splines of time. Heat exposure was modelled using average values of the temperature measure on the same day as death (lag 0) and the day before (lag 1). The heat effect was quantified assuming a linear increase in risk above a cut-point for each city. Socio-economic status indicators and census data were linked with mortality data for stratified analyses.
The risk of heat-related death increased with age, and females had a greater risk than males in age groups ≥65 years in London and Milan. The relative risks of mortality (per °C) above the heat cut-point by gender and age were: (i) Male 1.10 (95%CI: 1.07–1.12) and Female 1.07 (1.05–1.10) for 75–84 years, (ii) M 1.10 (1.06–1.14) and F 1.08 (1.06–1.11) for ≥85 years in Budapest (≥24°C); (i) M 1.03 (1.01–1.04) and F 1.07 (1.05–1.09), (ii) M 1.05 (1.03–1.07) and F 1.08 (1.07–1.10) in London (≥20°C); and (i) M 1.08 (1.03–1.14) and F 1.20 (1.15–1.26), (ii) M 1.18 (1.11–1.26) and F 1.19 (1.15–1.24) in Milan (≥26°C). Mortality from external causes increases at higher temperatures as well as that from respiratory and cardiovascular disease. There was no clear evidence of effect modification by socio-economic status in either Budapest or London, but there was a seemingly higher risk for affluent non-elderly adults in Milan.
We found broadly consistent determinants (age, gender, and cause of death) of heat related mortality in three European cities using a standard approach. Our results are consistent with previous evidence for individual determinants, and also confirm the lack of a strong socio-economic gradient in heat health effects currently in Europe.
Europe has warmed by 0.3°C per decade since the 1970s  and most countries have experienced an increase in the number of heat episodes in the past two decades . As public health measures are developed to reduce the impacts of heat-waves and hot weather, there is a need to better describe the environmental and social determinants of heat-related mortality.
In addition to physiological and clinical studies of heat stress, there is a growing literature of epidemiological studies that look at risk factors for heat-related mortality. These studies indicate that the elderly are at highest risk of heat-related mortality. However, many of these studies have used very broad age groups and do not sufficiently adjust for age when looking at gender or other subgroups, and therefore estimates are potentially subject to residual confounding. A multi-city case-crossover analysis in four Italian cities found the risk of heat related mortality was higher in women compared to men, after adjusting for age, but other studies have shown no differences by gender. Few studies have investigated the effect of socio-economic status on heat mortality in Europe. Excess mortality in Rome during the heat-wave of 2003 was lower in persons with the highest level of education , however other studies in Barcelona , Paris  and the UK  have shown little difference in impacts between high and low income groups.
We investigated the effect of heat on daily mortality in three European cities. The three cities selected were partners in the EUROHEAT project on improving public health responses to extreme weather : Budapest (Hungary), London (UK), and Milan (Italy) and represent a range of climates across Europe.
Budapest has a continental climate, with relatively high summer temperatures and recurrent heat-waves are a problem . Exposure to air pollutants is relatively high compared to other countries in Europe [10, 11]. London has a maritime climate, and a well described heat island which means that summer night time temperatures can be 1–3°C higher than the surrounding areas . Air pollution in London is relatively modest, although particulate matter and ozone levels were elevated during the 2003 heat-wave event . Recent research has shown London is more sensitive to heat-related mortality than other regions or urban areas in England and Wales . Milan and London were both severely affected by the 2003 heat-wave. Approximately 560 excess deaths (an increase of 23% in total mortality) were observed in Milan [4, 14] and 616 deaths (42%) in Greater London  during the heat event. Temperatures in Budapest were not much influenced by the heat-wave which predominantly affected western Europe. In all three cities, heat health warning systems were implemented after the 2003 event [9, 16, 17].
We investigated a range of risk factors of heat-related death using routine mortality data linked to small area indicators of environment and social status. This comparative study uses a standard time-series approach to identify the general determinants of heat-related mortality in European settings.
Mortality data were supplied by the Central Statistical Office in Budapest, the Office for National Statistics (ONS) in the UK, and the Local Health Authority in Milan. Daily counts of deaths from all-cause, cardiovascular disease (International Classification of Diseases, ICD-9 390.0–459.9 for deaths before 2001; ICD-10 I for deaths after 2001), respiratory disease (ICD-9 460.0–519.9; ICD-10 J), and external disease (ICD-9 900.0–999.9; ICD-10 S,T,V,W,X,Y,Z) were obtained for each city. We created mortality series for the following narrow age bands (0–14, 15–64, 65–74, 75–84, ≥85 years).
Linkage of mortality to deprivation indices was based on the best data available for each city. Information on income was derived from the published 2001 census in each city. In London, we computed population-weighted average score of the Multiple Deprivation Index (MDI) by Census Area Statistic ward. For Milan, we used the annual income (median) of residents by census tract after linking data from the Tax Register with the Milan Population Registry. Small area statistics were not available in Budapest so census data were linked at district-level. We used the proportion of the population who completed at least secondary school education as the indicator for socio-economic status. Other factors available from the census, such as dwelling type and proportion of elderly residents living alone, were also linked to investigate their effects on heat-related mortality.
Weather and pollution data
Characteristics of study cities
Average summer (June-Aug) temperature, °C
95th percentile daily mean temperature, °C
Particulate concentrations, μg/m3
Mean (5th to 95th percentile)
TSP: 57.0 (26.0 – 104.5)
PM10: 32.9 (18.5–60.5)
PM10: 58 (19–137)
8 hr Ozone concentration, μg/m3
Mean (5th to 95th percentile)
73.3 (26.0 – 132.5)
21.8 (5.5 – 46.0)
74 (9 – 173)
Daily number deaths
Mean (5th to 95th percentile)
75 (56 – 97)
171 (134 – 223)
30 (19 – 43)
Age at death (years), %
Cause of death, %
Each daily mortality series was examined in relation to daily temperature using Poisson generalised linear models allowing for over-dispersion, following methods used in previous analyses for England and Wales . Cubic smoothing splines of time with equally spaced knots were used to control for secular trends in the mortality series and any additional confounding by seasonally-varying factors other than temperature. The same level of seasonal control was used on each series with 7 degrees of freedom (df) per year, roughly equivalent to a two-month moving average, specified for the splines. We chose the number of df as a compromise between providing adequate control for unmeasured confounders and leaving sufficient information from which to estimate temperature effects. Sensitivity analyses were conducted to confirm that estimates were largely unchanged if other levels of seasonal control were considered.
Daily levels of PM10 and O3 (both average of the current and previous day) were incorporated into each regression model as possible confounding variables, regardless of statistical significance. Indicator variables for day-of week and public holidays were included. As cold effects can be apparent during months of even moderate temperatures, it was decided to have some degree of control for low temperatures. Cold effects are more delayed than for heat , and so cold effects were considered using a linear term of lag 0–13 (temperature averaged across values 0–13 days before the day of death) below a cut-point set at the 5th percentile of daily mean temperature.
To establish the general relationship between mortality and temperature, natural cubic splines of the temperature series (df = 3) were regressed against model residuals after controlling for the confounding factors noted above. The heat effect was modelled using lag0–1 (temperature averaged across values on the same day and the previous day). Heat effects were quantified assuming simple linear threshold models for each city. Thus, the heat effect is the log-linear increase in risk above a heat cut-point defined as the 95th percentile of daily mean temperature (lag 0–1) (Table 1). All analyses were conducted using STATA v9 .
Table 1 describes the characteristics of the mortality and environmental datasets in each of the three cities. The crude mortality rate in Budapest (4 per 1000 people in population) was high compared with London and Milan (2 per 1000), although the study periods were slightly different. Milan has a more elderly population compared to the other two cities. Summertime average of daily mean temperature was highest in Milan. Budapest and Milan have greater inter-annual temperature ranges because of their continental location. The "heat cut-points" (°C) used in this analysis were 24.4, 20.4, and 26.3, in Budapest, London and Milan, respectively. Observed PM10 and ozone levels were highest in Budapest and Milan.
Percent change in deaths (95%CI) per °C increase in temperature above cut-point
Aged &75 year
Aged ≥ 75 year
Aged &75 year
Aged ≥ 75 year
Aged &75 year
Aged ≥ 75 year
Cause of death
1.03 (1.02 to 1.05)
1.06 (1.04 to 1.07)
1.03 (1.02 to 1.04)
1.06 (1.05 to 1.07)
1.12 (1.08 to 1.16)
1.17 (1.14 to 1.20)
1.04 (0.99 to 1.08)
1.08 (1.04 to 1.12)
1.03 (1.01 to 1.04)
1.06 (1.05 to 1.07)
1.15 (1.07 to 1.24)
1.20 (1.15 to 1.25)
1.06 (0.98 to 1.14)
1.08 (1.01 to 1.15)
1.05 (1.01 to 1.08)
1.08 (1.06 to 1.10)
1.37 (1.15 to 1.62)
1.22 (1.14 to 1.31)
1.04 (0.98 to 1.11)
1.02 (0.96 to 1.09)
1.06 (1.02 to 1.10)
1.10 (1.02 to 1.18)
1.21 (1.05 to 1.40)
1.18 (1.03 to 1.34)
1 (least deprived)
1.03 (1.01 to 1.06)
1.06 (1.03 to 1.09)
1.03 (1.01 to 1.05)
1.06 (1.04 to 1.08)
1.19 (1.11 to 1.27)
1.21 (1.16 to 1.26)
1.03 (1.01 to 1.06)
1.04 (1.01 to 1.07)
1.04 (1.02 to 1.06)
1.05 (1.04 to 1.07)
1.13 (1.05 to 1.22)
1.13 (1.08 to 1.20)
1.02 (0.99 to 1.06)
1.06 (1.03 to 1.10)
1.02 (0.99 to 1.04)
1.07 (1.05 to 1.09)
1.11 (1.02 to 1.20)
1.11 (1.05 to 1.17)
1.04 (1.02 to 1.07)
1.06 (1.03 to 1.08)
1.03 (1.01 to 1.05)
1.07 (1.05 to 1.09)
1.09 (1.00 to 1.19)
1.20 (1.13 to 1.27)
5 (most deprived)
1.03 (1.01 to 1.06)
1.07 (1.04 to 1.10)
1.02 (1.00 to 1.05)
1.04 (1.03 to 1.06)
1.00 (0.90 to 1.12)
1.18 (1.10 to 1.26)
There was little evidence in Budapest and London of modification of heat effects by quintiles of the proportion of people living in flats and the proportion of elderly people (≥65 years) living alone (data not shown). These variables were not available for Milan.
We used a standard approach to investigate heat effects in the three European cities and observed broadly consistent patterns in different sub-groups. Our results confirm that elderly people are most at risk from heat-related mortality, and, within this group women seem to be at increased risk compared to men in London and Milan.
When analyzing by cause, we observed the strongest heat effect on deaths from respiratory and cardiovascular diseases, as reported by other studies [20–23]. This is important from a public health point of view since cardiovascular disease is the leading single cause of death in Europe: 49% of all deaths in Budapest, 39% in London, and 37% in Milan during the study period. Most previous studies of weather and health have focussed on cardio-respiratory causes of death. However, our results suggest that mortality from external causes is also sensitive to heat, and that this effect is apparent in both adults and the elderly. A previous study in England and Wales, also found a heat signal in deaths from external causes , and indicated a stronger heat effect in the 0–64 years age group.
The exact mechanism by which heat affects mortality from external causes needs further investigation. The risk of accidents from traffic accidents, falls, and drowning are known to be affected by summer weather [24–27]. Cause of death certification in the elderly is often not accurate and a proportion of the deaths from cardiovascular and respiratory disease are likely to be misclassified. Deaths from external causes, however, are unlikely to be misclassified. Although heat stroke deaths are underreported , deaths certified as due to classical heat illness are unlikely to contribute significantly to the number of deaths from external causes in our three cities.
We found little evidence of heat effects on mortality in children in Milan and Budapest, although there was some indication of a small effect on mortality in girls in London. A previous study in Oslo reported no effect of temperature on child mortality . Overall our study indicates that heat does not have a significant effect on mortality in children, despite their limited ability to thermoregulate. There is no clear physiological mechanism to explain a gender difference for heat effects in children, although some studies have shown gender differences in physiological heat responses in adults . More research is needed on the effects of heat on mortality and morbidity in children.
In Europe, epidemiological studies have shown little or no effect of socio-economic deprivation in the risk of heat-related mortality. Our results were broadly consistent with these reports, showing no clear gradient in heat-related mortality across the deprivation quintiles in the three cities except for the non-elderly population in Milan. It is possible that analyses of district and small area indicators are not able to detect a true but small effect on risk. In the UK, the Office for National Statistics has developed the MDI by Super Output Area (SOA) which has homogeneous characteristics (average population per SOA = 1500) from the 2001 census. There is no such published index for socio-economic deprivation in Hungary and Italy. For Budapest, we used the percentage of people who completed at least secondary school as it was more likely to have a normal distribution with a wider range than other possible indicators available in the census 2001. Area-level markers were used to represent individual circumstances and are thus subject to the ecological fallacy. Each unit area contains 32,534 households in Budapest, 4,765 households in London, and 118 households in Milan.
This study is not comparing estimated heat effects across cities, but only across subgroups within each city. There are important differences in sensitivity to heat between cities determined by weather other than temperature, air pollution, and other factors, such as population characteristics, cultural behaviours, and housing. The estimates of heat effects for subgroups in Milan appear large compared to the corresponding subgroups in the other cities. One reason is that the 95th percentile cut-point (26.3°C) in this study is rather high compared to a previous study in Milan  which suggested a cut-point of 23°C determined by model fit using data for 18 years (1985–2002). Applying a cut-point of 26.3°C results implies that only the more extreme hot days contribute to the mortality risk. Nevertheless the focus of the current study was on comparisons within cities and not between.
Our study provides consistent evidence of the environment and social determinants of heat-related mortality in three European cities. Further research should focus on identifying the spatial distribution of heat-related mortality, particularly in relation to heat island "hot spots" and housing characteristics.
We thank the Office for National Statistics, UK, NIEH Central Statistical Office, Budapest, Hungary, and Local Health Authority, Milan, Italy for supplying the data on mortality. This work was undertaken as part of "EUROHEAT: Responses to Heat-waves in Europe", funded by the European Commission DG SANCO (agreement number 2004322). AI is funded by Japanese Society for Promotion of Science. SH is funded by a Wellcome Trust Research Career Development fellowship.
- Klein Tank A, Wijngaard J, van Engelen A: Climate of Europe: Assessment of observed daily temperature and precipitation extremes. 2002, the Netherlands , KNMI, De Bilt, 36-Google Scholar
- Tank A, Konnen GP: Trends in indices of daily temperature and precipitation extremes in Europe, 1946-99. J Clim. 2003, 16 (22): 3665-3680.View ArticleGoogle Scholar
- Stafoggia M, Forastiere F, Agostini D, Biggeri A, Bisanti L, Cadum E, Caranci N, de' Donato F, De Lisio S, De Maria M, Michelozzi P, Miglio R, Pandolfi P, Picciotto S, Rognoni M, Russo A, Scarnato C, Perucci CA: Vulnerability to heat-related mortality: a multicity, population-based, case-crossover analysis. Epidemiology. 2006, 17 (3): 315-323.View ArticleGoogle Scholar
- Michelozzi P, de Donato F, Bisanti L, Russo A, Cadum E, DeMaria M, D'Ovidio M, Costa G, Perucci CA: The impact of the summer 2003 heat waves on mortality in four Italian cities. Euro Surveill. 2005, 10 (7): 161-165.Google Scholar
- Borrell C, Mari-Dell'Olmo M, Rodriguez-Sanz M, Garcia-Olalla P, Cayla JA, Benach J, Muntaner C: Socioeconomic position and excess mortality during the heat wave of 2003 in Barcelona. Eur J Epidemiol. 2006, 21 (9): 633-640.View ArticleGoogle Scholar
- Canoui-Poitrine F, Cadot E, Spira A: Excess deaths during the August 2003 heat wave in Paris, France. Rev Epidemiol Sante Publique. 2006, 54 (2): 127-135.View ArticleGoogle Scholar
- Hajat S, Kovats RS, Lachowycz K: Heat-related and cold-related deaths in England and Wales: who is at risk?. Occup Environ Med. 2007, 64 (2): 93-100.View ArticleGoogle Scholar
- WHO Europe: 1st Meeting of the project: improving public health responses to extreme weather/heat-waves - EuroHEAT: 20-22 June 2005; Rome, Italy.2005, ,Google Scholar
- Paldy A, Bobvos J, Vamos A, Kovats RS, Hajat S: The effect of temperature and heat waved on daily mortality in Budapest, Hungarym 1970-2000. Extreme weather events and public health responses. Edited by: Kirch W, Menne B, Bertollini R. 2005, NY , Springer-Verlag, 99-108.View ArticleGoogle Scholar
- Aga E, Samoli E, Touloumi G, Anderson HR, Cadum E, Forsberg B, Goodman P, Goren A, Kotesovec F, Kriz B, Macarol-Hiti M, Medina S, Paldy A, Schindler C, Sunyer J, Tittanen P, Wojtyniak B, Zmirou D, Schwartz J, Katsouyanni K: Short-term effects of ambient particles on mortality in the elderly: results from 28 cities in the APHEA2 project. Eur Respir J Suppl. 2003, 40: 28s-33s.View ArticleGoogle Scholar
- Katsouyanni K, Touloumi G, Samoli E, Gryparis A, Le Tertre A, Monopolis Y, Rossi G, Zmirou D, Ballester F, Boumghar A, Anderson HR, Wojtyniak B, Paldy A, Braunstein R, Pekkanen J, Schindler C, Schwartz J: Confounding and effect modification in the short-term effects of ambient particles on total mortality: results from 29 European cities within the APHEA2 project. Epidemiology. 2001, 12 (5): 521-531.View ArticleGoogle Scholar
- Wilby RL: Past and projected trends in London's urban heat island. Weather. 2003, 58 (7): 251-260.View ArticleGoogle Scholar
- Stedman JR: The predicted number of air pollution related deaths in the UK during the August 2003 heatwave. Atmospheric Environment. 2004, 38 (8): 1087-1090.View ArticleGoogle Scholar
- Conti S, Meli P, Minelli G, Solimini R, Toccaceli V, Vichi M, Beltrano C, Perini L: Epidemiologic study of mortality during the Summer 2003 heat wave in Italy. Environ Res. 2005, 98 (3): 390-399.View ArticleGoogle Scholar
- Johnson H, Kovats RS, McGregor G, Stedman J, Gibbs M, Walton H, Cook L, Black E: The impact of the 2003 heat wave on mortality and hospital admissions in England. Health Stat Q. 2005, 6-11.Google Scholar
- Kirchmayer U, Michelozzi P, de'Donato F, Kalkstein LS, Perucci CA: A national system for the prevention of health effects of heat in Italy. Epidemiology. 2004, 15 (4): S100-S101.View ArticleGoogle Scholar
- Department of Health: Heatwave: Plan for England - Protecting health and reducing harm from extreme heat and heatwaves (2004 edition). 2004Google Scholar
- Braga AL, Zanobetti A, Schwartz J: The time course of weather-related deaths. Epidemiology. 2001, 12 (6): 662-667.View ArticleGoogle Scholar
- Stata Corporation: Stata 9. 2006, TX , Stata CorporationGoogle Scholar
- Kovats S, Koppe C: Heatwaves: past and future impacts on health. Ingtegration of public heatlh with adaptation to climate change: Lessons learned and new directions. Edited by: Ebi KL, Smith J, Burton I. 2005, Lisse , Taylor & Fransis Group, 621-627.Google Scholar
- Nafstad P, Skrondal A, Bjertness E: Mortality and temperature in Oslo, Norway, 1990-1995. Eur J Epidemiol. 2001, 17 (7): 621-627.View ArticleGoogle Scholar
- Kunst AE, Looman CW, Mackenbach JP: Outdoor air temperature and mortality in The Netherlands: a time-series analysis. Am J Epidemiol. 1993, 137 (3): 331-341.Google Scholar
- Huynen MM, Martens P, Schram D, Weijenberg MP, Kunst AE: The impact of heat waves and cold spells on mortality rates in the Dutch population. Environ Health Perspect. 2001, 109 (5): 463-470.View ArticleGoogle Scholar
- Christoffel KK: Effect of season and weather on pediatric emergency department use. Am J Emerg Med. 1985, 3 (4): 327-330.View ArticleGoogle Scholar
- Rapport RS: Water-related injuries. The avoidable tragedies of summer. N J Med. 1999, 96 (7): 33-35.Google Scholar
- DiMaggio C, Durkin M: Child pedestrian injury in an urban setting: descriptive epidemiology. Acad Emerg Med. 2002, 9 (1): 54-62.View ArticleGoogle Scholar
- Silm S, Ahas R: Seasonality of alcohol-related phenomena in Estonia. Int J Biometeorol. 2005, 49 (4): 215-223.View ArticleGoogle Scholar
- Goldacre MJ, Duncan M, Cook-Mozaffari P, Griffith M: Mortality rates for common respiratory diseases in an English population 1979-1998: artefact and substantive trends. J Public Health (Oxf). 2004, 26 (1): 8-12.View ArticleGoogle Scholar
- Donoghue ER, Graham MA, Jentzen JM, Lifschultz BD, Luke JL, Mirchandani HG: Criteria for the diagnosis of heat-related deaths: National Association of Medical Examiners. Position paper. National Association of Medical Examiners Ad Hoc Committee on the Definition of Heat-Related Fatalities. Am J Forensic Med Pathol. 1997, 18 (1): 11-14.View ArticleGoogle Scholar
- Havenith G: Temperature, heat balance, and climatic stress. Extreme Weather Events And Public Health Responses. Edited by: Wilhelm Kirch BMRB. 2005, Berlin , Springer-Verlag, 70-80.Google Scholar
- Hajat S, Armstrong B, Baccini M, Biggeri A, Bisanti L, Russo A, Paldy A, Menne B, Kosatsky T: Impact of high temperatures on mortality: is there an added heat wave effect?. Epidemiology. 2006, 17 (6): 632-638.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.