The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project

Background The present study aimed at developing a standardized heat wave definition to estimate and compare the impact on mortality by gender, age and death causes in Europe during summers 1990-2004 and 2003, separately, accounting for heat wave duration and intensity. Methods Heat waves were defined considering both maximum apparent temperature and minimum temperature and classified by intensity, duration and timing during summer. The effect was estimated as percent increase in daily mortality during heat wave days compared to non heat wave days in people over 65 years. City specific and pooled estimates by gender, age and cause of death were calculated. Results The effect of heat waves showed great geographical heterogeneity among cities. Considering all years, except 2003, the increase in mortality during heat wave days ranged from + 7.6% in Munich to + 33.6% in Milan. The increase was up to 3-times greater during episodes of long duration and high intensity. Pooled results showed a greater impact in Mediterranean (+ 21.8% for total mortality) than in North Continental (+ 12.4%) cities. The highest effect was observed for respiratory diseases and among women aged 75-84 years. In 2003 the highest impact was observed in cities where heat wave episode was characterized by unusual meteorological conditions. Conclusions Climate change scenarios indicate that extreme events are expected to increase in the future even in regions where heat waves are not frequent. Considering our results prevention programs should specifically target the elderly, women and those suffering from chronic respiratory disorders, thus reducing the impact on mortality.


Introduction
Several studies in the United States and Europe have shown marked short-term increases in the number of deaths during heat waves, however the role of these extreme events on mortality needs to be better clarified [1]. Climate change projections for Europe show that over the next century, heat waves will become more frequent, intense and will last longer, not only in Mediterranean regions, but also in Northern areas currently not characterized by heat wave events [2]. These changes could contribute to the burden of disease and premature deaths, particularly in vulnerable populations with limited adaptation resources [3].
In meteorological terms, a heat wave is defined as a prolonged period of unusually hot weather. To date, a standard definition of heat wave has not been agreed upon and different definitions have been used to evaluate its impact on health [4][5][6][7].
During summer 2003, Europe experienced one of the worst heat wave events in recent history [8], with an estimated excess mortality varying between 25.000 and 70.000 deaths in Western Europe [9,10]. The comparison of the impact of the 2003 heat wave between countries is hampered by the substantial differences in the methodologies employed to define heat wave events and to estimate the excess mortality [4,6,[11][12][13].
The present paper presents the main results of the EuroHEAT project (Improving Public Health Responses to extreme weather/heat-waves), which aimed to develop a standardized definition of a heat wave event and to compare the impact of heat waves in European cities. The study took advantage of the database and the findings of a previous European collaborative study (Assessment and Prevention of acute Health Effects of Weather conditions in Europe -PHEWE), which estimated the impact of high temperatures on mortality using a time series approach [14]. The methodology developed within EuroHEAT meets the need to investigate the specific contribution of single heat wave events on mortality and allows for a valid comparison between cities.
Given the evidence that heat waves do not affect the entire population but susceptibility is greater among the elderly [3,4,7], EuroHEAT was focused on people aged 65 years and over.
The main objectives of this paper are: • to provide a unique definition of "heat wave" to be applied across cities in the subsequent analyses; • to estimate the health impacts of heat wave episodes on the elderly population, by gender, specific causes and age groups; • to evaluate the role of specific heat wave characteristics, namely duration, intensity and timing within the season on mortality; • to compare the impact of the 2003 heat wave across nine European cities using the same methodological approach.

Study population and setting
The EuroHEAT project involves nine European cities (Athens, Barcelona, Budapest, London, Milan, Munich, Paris, Rome, Valencia), with a total population of around 25 million citizens, which represent a variety of climatic, socio-economic, and air pollution characteristics. Available data for daily mortality, meteorological and air pollution were provided for each city between 1990 and 2004. Only summer months (June-August) were included in the study. Summer 2003 was analyzed separately to assess the impact of this exceptional heat wave episode that affected most of the European cities, and results were compared with heat waves from other years included in the study period.

Exposure definition
Exposure to heat waves considered both extreme day time values in terms of maximum apparent temperature (Tappmax) and high night time temperatures through minimum temperature (Tmin).
Tappmax is a discomfort index based on air and dew point temperature, calculated using the following formula [15,16]: Heat waves were thus defined as 1) periods of at least two days with Tappmax exceeding the 90 th percentile of the monthly distribution or 2) periods of at least two days in which Tmin exceeds the 90 th percentile and Tappmax exceeds the median monthly value. Daily counts of total and cause specific mortality (cardiovascular, cerebrovascular and respiratory) in the heat wave periods were considered as outcome variable, stratifying by gender and age groups (65-74 yrs, 75-84 yrs, 85 + yrs).

City specific analysis
As the first stage, a city specific analysis was conducted using Generalized Estimating Equations (GEE models) [17] to analyze longitudinal data. A Poisson distribution was assumed for the outcome variable (mortality) and days were characterized as "heat wave" or "non heat wave" days as the exposure variable to estimate the effect on mortality.
Observations from different years were assumed to be independent, while observations within the same summer were correlated. A similar approach has already been suggested in other studies [18][19][20][21]. Since the number of clusters (summers) was small, and equal to the number of years in the study period, we used the model-based estimator for the coefficients' standard errors, as recommended in the presence of few large clusters [22]. A first order autoregressive structure within each year was chosen, based on an exploratory analysis similar to the one described by Chiogna & Gaetan [23,24]. . .

T T air dewpt
A common model was specified for each city taking into account the following as potential confounders: holidays, day of the week and calendar month, linear terms for barometric pressure (lag 0-3) and wind speed, linear and quadratic terms for time trend and the maximum 1hour daily value of NO 2 (μg/m 3 ) at (lag 0-1). NO 2 was chosen to adjust for air pollution as it has shown to be a good indicator of traffic in a large European collaborative project [25] which assessed the impact of air pollution on mortality, using meteorological variables as potential effect modifiers.
The effect was estimated as the percent increase in daily mortality during heat wave compared to non heat wave days. The effect of each heat episode was also investigated with respect to specific heat wave characteristics such as duration, intensity and timing within the season.
Duration was categorized using the city specific median values of heat wave duration (number of consecutive heat wave days) as the cut-point: short heat wave if duration was shorter than the median, long heat wave if duration was equal to or longer than the median. Intensity was also categorized into two levels according to the extreme values of Tappmax reached during the heat wave: low intensity heat wave if Tappmax was below the monthly 95 th percentile, high intensity heat wave if Tappmax was equal to or above the monthly 95 th percentile. Timing within the season was defined according to the time interval between different heat waves: the first heat wave of each summer was considered separately, heat waves that occurred between 1 and 3 days after the previous one, heat waves that occurred three or more days after the previous one.

Pooled analysis
To summarize city-specific results, cities were grouped into two regions, according to geographical and climatological criteria in order to control for heterogeneity, as defined in the PHEWE study [18]: "Mediterranean" (Athens, Barcelona, Milan, Rome, and Valencia) and "North-Continental" (Budapest, London, Munich and Paris) Cityspecific estimates were combined through a random effect meta-analysis using the method described by Der-Simonian and Laird [26]. To estimate the impact in each region, a GEE model was used, similar to the city-specific model, but adding a city indicator variable and interaction terms of the exposure variable with the confounders.

Heat wave definition
The appropriateness of the heat wave definition was tested on the 2003 time series data, as this was the year when heat waves were observed in almost all cities. For all cities, the definition was able to capture consecutive days with high temperatures, excluding only isolated peaks. Table 1 presents datasets and the city characteristics. The number of natural deaths per day in the age group 65 + ranged from 12 in Valencia to 117 in London. Cardiovascular causes accounted for about 40-50% of the total, except in Paris (34.1%), Budapest (60.9%) and Athens (51.6%). Cerebrovascular mortality was about 10% of total mortality, with the lowest proportion observed in Munich and Valencia (8.3%) and the highest in Athens (17.7%). Great heterogeneity was observed in respiratory deaths, generally accounting for less than 10% of the total mortality, with the highest proportion in London (16.2%) and the lowest in Budapest and Munich (4.3%). Mortality data for specific causes of death were not available for Paris in 2003.
Temperature indicators showed a large variability of the meteorological conditions. Considering all years in study, the range for minimum temperature was from 15.9 °C in Munich to 26.4 °C in Athens, and for maximum apparent temperature from 27.1 °C in London to 39.9 °C in Valencia. In every city the number of heat waves days was heterogeneous among years. In 2003, an increase in both minimum and maximum apparent temperature was observed in all cities, except in Athens. The number of heat wave days was generally higher in 2003 in all cities except in Athens, Budapest and London. Table 2 showed the heat wave effect in the study period (1990-2002 and 2004), estimated as the percentage increase in mortality per day with respect to non heat waves days for all natural deaths, respiratory, cardiovascular and cerebrovascular mortality. During heat waves days a significant increase in total daily mortality was observed in all cities, with the greatest increase in Milan (+ 33.6%, 95%CI + 28.5%-39.0%) and the lowest in Munich (+ 7.6, 90% C.I. + 3.8-11.5). A significant increase in daily mortality was observed for all the causes of death considered, except for respiratory mortality in Munich. In most cities the highest increment was observed for respiratory causes. Table 3 summarize the pooled results for North continental and Mediterranean cities, by gender, age groups, and cause of death (all causes, respiratory, cardiovascular, and cerebrovascolar causes). The effect was higher in Mediterranean cities than in North-Continental cities for all groups of causes considered. In both two groups of cit- ies the highest increase in daily mortality was observed for respiratory causes. The effect of heat waves increased with age and a statistically significant higher impact among females than males was observed in Mediterranean cities in the 75-84 age group.

Role of specific heat wave characteristics
The effect of specific characteristics of heat waves in terms of duration and a combination of duration with intensity was investigated and results are shown in Figure  1. A higher effect of heat waves of longer duration and high intensity is shown in most cities. Results give evidence for duration to play a bigger role than intensity. In the pooled analysis, a two-fold impact of longer and high intense heat waves is observed in both European regions.
The analysis of the role of timing of heat waves did not produce homogeneous results across cities. The first heat wave of the summer had a greater impact only in some cities. In case of subsequent heat waves, those occurring after less than 3 days after the end of the preceding one had a weaker effect (data not shown).

Discussion
This study is the first in Europe to compare the impact of heat waves on mortality in different cities using a common heat wave definition and a standardized methodological approach. Heat wave episodes were defined taking into account both temperature and humidity levels during the day as well as high night-time temperatures, thus accounting for the physiological impact of heat on health [27]. Moreover, the impact on mortality was also evaluated considering specific heat wave characteristics; results showed that the intensity, duration and timing of heat waves influenced the risk of mortality during extreme events. Heat waves of long duration had the greatest impact on mortality, and resulted in 1.5 to 3 times higher daily mortality than for other heat waves.
These results together with findings from the previous study (PHEWE project) conducted in the same cities [18], using a time series approach, provide a complete overview of the effect of high temperature in European cities.
Results showed great heterogeneity in the effect of heat waves on daily mortality. In the study period, the strongest impact was observed in the Mediterranean cities where heat waves are more frequent and characterised by higher temperatures, confirming results from the PHEWE project time series analysis [18].
In contrast, in 2003 the greatest impact on mortality was observed in cities where heat wave episodes are rare events or were characterized by temperatures largely outside the usual meteorological conditions. In several cities, the 2003 episodes were characterized by more extreme temperatures, with both maximum and minimum temperature several degrees above the average values. However, differences in exposure alone do not seem to explain the greater impact observed in North-Continental than in Mediterranean cities. A poor ability to adapt to high  temperatures among residents in North-Continental cities may be the explanation for the highest impact on these populations, as observed in France in 2003 [28]. We found that the impact on daily mortality increases with age. Moreover, gender was among the factors that increases individual vulnerability to heat waves; we observed a higher susceptibility of females even after stratifying by age groups. Other authors observed similar results [9,10,12] and this finding may be attributable to social conditions of elderly women living alone and to physiological differences, such as a reduced sweating capacity that affects the ability to respond to heat stress [9]. It may also be due to residual confounding of age, as within the older age groups there will be a greater number of women than men.
In most cities, a greater effect of heat waves on respiratory than on cardiovascular mortality was observed. This is in agreement with other studies [7,29] and probably reflects the fact that health status of people suffering from chronic respiratory diseases rapidly deteriorates during hot periods or that the capacity to cope with heat stress is generally diminished in particular population subgroups, such as the elderly [30,31]. However, the underlying mechanisms still remain unclear [32].
Several population characteristics like social isolation and income level have been identified as factors that affect the susceptibility of urban populations to heat, as have the adaptive measures and actions in place [33]. However, information on these characteristics were not available in the present study.
A limit of the analysis is that neither the lag effect, nor the role of harvesting could be investigated. A relevant question that still remains unclear is whether a reduction in mortality is observed shortly after a severe heat episode. Results from the PHEWE study [18] showed some evidence of a harvesting effect both in the Mediterranean and in the North-Continental cities, although the observed excess in mortality associated with high temperatures was only partially compensated for by mortality displacement. Previous studies on heat wave episodes have also reported a compensatory decrease in overall mortality during the weeks following a heat wave [9,34-

36]
; however, the study on the impact of the 2003 heat wave in 9 French cities [28] did not report any short-term mortality displacement. Differences in exposure intensity and the differential impact of mortality on susceptible subgroups may explain why short-term mortality displacement is reported in some studies [9,[34][35][36] and not in others [28].
Considering the latest IPCC Assessment Report, climate change predictions for Europe show an increase in the frequency and the intensity of heat waves, especially in central, southern and eastern areas [3] and as consequence heat-related mortality will become a relevant threat even in cities usually not exposed to extreme hot temperatures. Considering our results prevention programs should specifically target the elderly, especially women, and those suffering from chronic respiratory disorders, in order to reduce in the future the burden of heat-related mortality.