The impact of weather on human health is a matter of increasing concern, especially in light of climate change and the 2003 heat wave in Western Europe.1,2 Climatologists forecast that temperature across Europe will rise over coming decades and the frequency of periods characterized by extremely high temperatures will double.3 Optimal health protection requires an understanding of the nature of the effect of weather conditions on health.
High temperatures have been examined in relation to natural (nonaccidental) mortality, to certain specific causes of death, and to other health effects such as hospitalizations.1,4 Time-series studies of the effects of heat on mortality have been conducted for several cities in the United States and in various European areas.5 The association between temperature and mortality has often been shown as a J- or V-shaped function, with the lowest mortality rates at moderate temperatures and rising progressively as temperatures increase or decrease.6 Variations in the temperature at which minimum mortality occurs may reflect population differences in acclimatization.6,7
In epidemiologic studies of the effects of heat, temperature has been generally used as the exposure of interest, sometimes with adjustment for humidity. More recently, several studies have considered the effect on mortality of apparent temperature, which is a thermal discomfort index that combines air temperature and humidity.7–9
The present work was developed within the PHEWE project (Assessment and Prevention of acute Health Effects of Weather conditions in Europe).10 In this project, separate analyses have been performed on the health effects of meteorologic conditions during warm and cold seasons. In the present paper, we report findings concerning the association between maximum apparent temperature and mortality during the warm season, defined as the period from 1 April to 30 September. The choice of a 6-month season was made with the aim of being able to detect the effect of early and late warm-season exposures. Our primary purpose was to investigate the relationship between heat and mortality across Europe, with particular attention to a comparison among cities. We also explored the delayed effect of heat and the assumption of a constant heat effect over the summer period.
The effect of weather on mortality during the warm season was explored in 15 cities: Athens, Barcelona, Budapest, Dublin, Helsinki, Ljubljana, London, Milan, Paris, Prague, Rome, Stockholm, Turin, Valencia, and Zurich. Data were collected for the period 1990–2000. Population size, percentage of elderly (≤75 years) and study period are reported in Table 1. The total population studied was about 30 million. The number of study years ranged from 5 in Athens to 11 in Dublin, Helsinki, Milan, and Stockholm.
All cities provided daily counts of deaths for all causes excluding external causes (natural mortality, ICD-9: 1–799), cardiovascular diseases (ICD-9: 390–459) and respiratory diseases (ICD-9: 460–519). Mortality counts were subdivided by age groups (15–64 yrs, 65–74 yrs, 75+ yrs). The daily average number of deaths for all natural causes during the warm season (April-September) ranges from 6.3 in Ljubljana (263,290 inhabitants) to 149 in London (6,796,900 inhabitants). Daily averages range from 2.6 to 61.0 per day for deaths from cardiovascular causes and from 0.4 to 23.7 for mortality due to respiratory diseases (Table 1).
Meteorologic and Air Pollution Data
All enrolled cities provided 3-hour meteorologic data (air temperature and dew point temperature, wind speed, and barometric pressure at sea level). We retrieved this information from the nearest airport weather station.
We focused on the effect of apparent temperature on mortality. Apparent temperature is a measure of relative discomfort due to combined heat and high humidity. It was developed by Steadman11 on the basis of physiologic studies on evaporative skin cooling. It can be calculated as a combination of air temperature (temp) and dew point (dew), according to the following formula12:
Daily values of apparent temperature, barometric pressure, and wind speed were computed for each city from 3-hour measurements. For all cities except Barcelona, we used daily maximum apparent temperature defined as the maximum of 3-hour values. On Barcelona, we calculated the daily average apparent temperature from the daily mean values of temperature and dew point, given the lack of 3-hour data. Sensitivity analyses using daily minimum apparent temperature as exposure variable were performed.
A large variability in climate characteristics was observed among cities. Daily average of maximum apparent temperature during April-September ranged from 14.3°C in Helsinki to 27.9°C and 29.5°C in Athens and Valencia, respectively (Table 1).
Air pollution data coming from the urban monitoring network were collected for each city. We used the maximum hourly value of nitrogen dioxide (NO2) as an indicator of the overall daily air pollution level for all the cities, except in Dublin, where we considered the daily average level of black smoke.
A negligible percentage of isolated missing values (less than 1% for most cities, with a maximum of 4% for Turin) characterized meteorologic and air pollution data. We imputed these through the average of values observed in the same month. For a few cities (Barcelona, Turin, and London), meteorologic information was lacking for series of consecutive days (1 or even 2 months). In this case, the missing periods were not included in the analysis (Table 1). Sensitivity analyses excluding years with higher number of missing measurements produced similar results (results not reported).
Rationale for the City-Specific Modeling Approach
In the larger project, we have investigated the effect of cold and heat separately. Restriction to 6-month seasons permitted us to reduce the complexity of controlling for temporal confounding and enabled us to use different model specifications by season. A common a priori definition of warm season, from 1 April to 30 September, was adopted for all the cities.
As data from each city were represented by several disjoint 6-month daily time series (clusters), we used a generalized estimating equations (GEE) approach, which extends generalized linear models to the analysis of longitudinal data, when the observations on different clusters can be assumed to be independent and the observations within clusters are correlated.13 For each city, the data consisted of an outcome variable (number of deaths) and several covariates (confounders and apparent temperature), observed day by day for each available summer. We assumed that observations during a single summer were correlated, while observations from different summer periods were independent. A similar approach has been suggested by Schwartz and Dockery14 in the context of the analysis of short-term effects of air pollution on health and recently by Hajat et al15 in studying the health effect of heat waves in 3 major European cities. Other examples are in the papers by Fouillet et al16 and by Gorjanc et al.17
Since the number of clusters is mostly lower than the number of observations within clusters, specifying the appropriate autocorrelation structure is a relevant point to assure consistency of the standard error estimates.18 We established the correct specification of the error structure taking advantage of an exploratory analysis based on dynamic models.19 The dynamic models are linear regression models with lagged regressors and correlated errors. To highlight a common autocorrelation structure, we defined dynamic models with different correlations between errors and, for each season and for each city, we employed an automatic method for the best model selection. The results suggested that we should use a first order autocorrelation structure. We then used the model-based variance estimator, as recommended in the presence of few large clusters.13
City-Specific Model Specification
For each city, we specified a GEE model, assuming the dependent variable was Poisson distributed. Confounders included in the models were dummy variables for holidays, day of the week and calendar month; linear and quadratic terms to pick up long-term time trend; linear terms for barometric pressure (lag 0–1) and wind speed. Air pollution was found to confound the effect of high temperature.20,21 We adjusted for maximum hourly NO2 concentrations (lag 0–1). The choice of NO2 as the air pollution indicator was based on comparability, availability, and completeness of daily measurements among cities.
We studied the effect of lag 0–3 exposure (the average of the current and the previous 3 days maximum apparent temperature).
First, we used a flexible parametric approach to describe the exposure-response curve, including in the model a cubic regression spline for apparent temperature with 1 knot every 8°C.22 Second, to summarize and to address heterogeneity, we further described the relationship between apparent temperature and mortality during the warm season by 2 linear terms constrained to join at a common point, which we call the threshold. The threshold is the value of apparent temperature, which corresponds to a change in the effect estimate. For a V-shaped curve, this is the value of apparent temperature associated with the minimum mortality rate. The city-specific thresholds were obtained by the maximum likelihood approach proposed by Muggeo.23 The slopes above the breakpoint were used as effect estimates.
A drawback of this method is that threshold estimation can be strongly dependent on the algorithm starting point.23 We used the same starting point for all the cities to avoid post hoc data dredging, but a sensitivity analysis showed that in our data the threshold estimate is substantially independent of the starting point.
Delayed Effect of Exposure
The lagged effect of daily maximum apparent temperatures above the city-specific threshold was investigated up to 40 days using a constrained distributed-lag model.24 We assumed that all coefficients from lag 0 to lag 40 fell on a curve of a 5th-order polynomial, to obtain a reasonably flexible distribution of lags and parsimony in the number of parameters estimated. We considered alternative constraint definitions and performed analyses leaving the current day effect unconstrained to check the robustness of estimates to outlying immediate effects. The various approaches produced negligible differences (results not shown).
Time-Varying Effect of Maximum Apparent Temperature
We relaxed the season definition extending our GEE model to allow a time-varying effect of lag 0–3 maximum apparent temperature during the warm season. We assumed that the effect of exposure above the threshold was smooth over time. We used sine and cosine terms to describe this pattern25:
where t indexes time, from 1 April to 30 September (t = 1, 2, …, 183), βt is the effect of maximum apparent temperature and (α0, α10, α11, α20, α21) are unknown coefficients to be estimated. This analysis provided a way for checking whether the slope above the threshold is constant over the warm season. At the same time, it gave information about the effective period when the effect is significant.
All city-specific analyses were conducted using STATA/SE 9.0 (StataCorp, College Station, TX) and R Software 1.2.1 (Team R Development Core).
City-specific results were pooled into 2 groups defined on the basis of meteorologic and geographic criteria10:
- “Mediterranean” cities: Athens, Rome, Barcelona, Valencia, Turin, Milan, and Ljubljana.
- “north-continental” cities: Prague, Budapest, Zurich, Paris, Helsinki, Stockholm, London, and Dublin.
The aim was to control part of the heterogeneity among cities and to compare the 2 regions. The definition of the geographic subgroups was done a priori.
We obtained overall exposure-response curves fitting GEE regression models with a cubic regression spline for maximum apparent temperature (with one knot every 8°C) on the pooled data sets.22 We did the analysis separately for Mediterranean and north-continental cities. Fixed city-specific intercepts and interaction terms between confounders and city indicators were introduced in the model. We assumed independence among years and first-order autocorrelation structure within season.
To produce regional averages of city-specific estimates of the effect parameters, we specified Bayesian hierarchical models.26 Let 1, …, c, …, C be the vectors of the city-specific estimates of interest and V1, …, Vc, …, VC the corresponding variance-covariance matrix estimates. In combining results from the distributed lag models and from the time-varying models, c is a vector of model parameters. In pooling estimates of threshold or slope, c is a scalar. A multivariate Gaussian distribution with unknown mean vector λc and variance-covariance matrix Vc was assumed for each c:
and a multivariate normal distribution with common mean vector and unknown variance-covariance matrix Σ was assumed on λc:
Σ models the intercity variability. Vague proper priors were elicited on hyper-parameters. Inference was based on an algorithm proposed by Everson and Morris, implemented in the TLNise library of R software or on MCMC methods with WinBUGS 14.27
The relationship between maximum apparent temperature (lag 0–3) and log mortality rates were V or J shaped for most cities. This indicated an excess of risk for exposures to apparent temperature above a threshold that varies among cities (Fig. 1).
We combined the city-specific curves except Barcelona, where the exposure was defined in terms of daily mean apparent temperature. The curves in Figure 2 represent the overall estimate of the exposure-response relationship, separately for Mediterranean and north-continental cities. The pooled curves had a minimum around higher apparent temperatures and the effect of heat appeared slightly stronger in Mediterranean cities than in north-continental cities. When we used minimum apparent temperature as the exposure indicator, the combined curve was similar in shape, but, as expected, with a minimum around lower apparent temperature values (Fig. 2).
We report the heat effect as percent change in mortality associated with a 1°C increase in maximum apparent temperature above the city-specific threshold. City-specific and overall meta-analytic estimates of thresholds and percent change are reported in Table 2. For Ljubljana, Stockholm, and Zurich estimated thresholds were just below 22°C whereas for Athens, Milan, and Rome estimates were over 30°C. The overall meta-analytic value of the threshold was 29.4°C (95% credibility interval [CrI] = 25.7 to 32.4) for Mediterranean cities (excluding Barcelona) and about 6 degrees lower for North-continental cities (23.3°C; 22.5 to 24.0). The posterior distributions in Figure 3A indicate larger heterogeneity in thresholds among Mediterranean cities than among north-continental cities. Heterogeneity among Mediterranean cities remained high even if Ljubljana was removed, with the city-specific threshold estimates ranging in this case between 27.0°C (Turin) and 32.7°C (Athens). Overall meta-analytic percent change per degree of above-threshold apparent temperature equaled 3.1 (0.6 to 5.7) and 1.8 (0.1 to 3.6) for Mediterranean and north-continental cities, respectively. Heterogeneity was larger among Mediterranean cities (Fig. 3B).
Meta-analytic estimates of percent change in mortality per degree of above-threshold apparent temperature by cause and age class are reported in Table 3. We estimated different city-specific thresholds for each cause of death, which were close to those obtained for all natural mortality (not reported). The overall meta-analytic estimates of percent change in cardiovascular mortality per degree of above-threshold temperature were 3.7 (0.4 to 7.0) for Mediterranean cities and 2.4 (−0.1 to 5.3) for north-continental cities. Higher associations were found between heat and mortality due to respiratory diseases, with estimated percent changes equal to 6.7 (2.4 to 11.3) and 6.1 (2.6 to 11.1) for Mediterranean and north-continental cities, respectively.
The effect of heat was particularly large in the elderly. For people aged 75 and older, we estimated that a 1°C increase in maximum apparent temperature above the threshold was associated with an increase in mortality for all natural causes of 4.2% for the Mediterranean region and of 2.1% for the north-continental region. The same effect estimates were 8.1% and 6.6%, respectively, when only deaths for respiratory causes were considered.
Figure 4 shows the meta-analytic distributed lag curves by geographic region. A strong effect of high apparent temperatures was evident within the first week. The excess mortality declined in subsequent days, became negative, then returned to the level of baseline mortality-more slowly in the Mediterranean cities, faster in the north-continental cities. Table 4 summarizes the posterior distributions of the cumulative effects calculated at various lags from unrestricted models. Cumulative effects were larger for the Mediterranean cities than for the north-continental cities. They declined after the first 5 days from the day of above-threshold apparent temperature, becoming immediately not significant for north-continental cities and after 15 days for Mediterranean cities.
The results of the time-varying coefficients analysis indicated an effect of apparent temperature exposure above the threshold between June and August. Evidence of a stronger association between early above-threshold exposure and total mortality was found, although there were a small number of days during which maximum apparent temperature exceeded the threshold in April and May (Fig. 5).
Study Design and Modeling Issues
The exposure-response curves using daily maximum versus minimum apparent temperature were substantially similar in shape although (of course) shifted. The minimum apparent temperature curve was less precise and had a reduced slope above the threshold. This result could indicate a misclassification bias with the minimum apparent temperature: perhaps minimum (usually overnight) temperature is a less sensitive indicator of risk than the maximum apparent temperature.
We based our analysis on a GEE approach, assuming independence of the observations from different summer periods. This assumption implies the absence of long-term harvesting-like phenomena; ie, we are excluding the possibility that high mortality levels during one summer induce low mortality levels in the next summer, and vice versa, thus focusing on the effect of heat only within season.
For each city, the effect of heat was summarized by a minimum risk threshold and a linear term describing the effect of exposures exceeding the threshold. In principle, the algorithm for segmented regression would allow more than one threshold,23 but the introduction of 2 cut-points would have made comparison among cities more difficult (despite a possible gain in terms of fitting).
The Muggeo algorithm for threshold estimation can be unstable, especially in small sample or when cut-points are not well defined.23 In some studies, threshold estimation has been avoided and the effect of heat calculated by comparing the mortality rate at 2 different points over a flexible curve that describes the mortality-temperature association.9 This approach provides effect estimates independent from the threshold, but dependent on the 2 points selected for comparison. This arbitrariness can be a drawback in comparing results from cities characterized by different ranges of exposure.
We did not make the strong assumption of a common threshold for all the cities, but instead obtained city-specific thresholds. In the same way, we estimated different thresholds for each cause of death, because in principle the minimum of the exposure-response curves can be different for different specific causes of death. In either situation, an incorrect assumption of common threshold could produce biased estimates of the slope. The same reasoning is in principle valid also for the age-specific analyses, but in this case estimating different thresholds for each age group is problematic due to the small number of daily events among younger people. This would worsen the performance of the estimation algorithm.
We investigated the effect of heat on a large scale across a variety of climate conditions and socioeconomic and demographic characteristics. We applied a standardized methodology to data from 15 European cities, with a gain of power and assessment of variability across city-specific effect estimates.
We did not focus on heat wave episodes. Our results provide evidence of an effect of heat on daily mortality. Many investigators have reported V- or J-shaped associations between temperature and mortality.6,7,9,15 We confirmed these findings.
Previous studies have shown that the temperature level corresponding to the minimum mortality rate (threshold) varies from city to city, and across different latitudes according to the local climate.6,7,28 We found that the apparent temperature threshold in the Mediterranean cities was higher than in the north-continental cities, indicating that residents of north-continental cities are susceptible at lower values of apparent temperature.
As far as the effect of exposures that exceed the threshold is concerned, stronger percent changes in all natural mortality and in respiratory and cardiovascular mortality were observed in the Mediterranean region. Steeper slopes were found for mortality for respiratory causes, consistent with results of previous studies.29
The associations with mortality were stronger after 74 years of age than earlier in life, suggesting that the elderly are particularly vulnerable.9,29,30
Knowledge of the lag time between exposure to extreme weather conditions and negative health outcomes is important for public health authorities and health care providers in developing prevention plans. The effect of high apparent temperature is immediate both in Mediterranean and north-continental cities. This result is consistent with findings of previous studies indicating that the impact of heat on mortality reaches its maximum in less than a week.31,32
There is contradictory evidence as to whether the increase of mortality is followed by a deficit that partly compensates for the negative effect (harvesting). In some cases, a harvesting effect is found to fully balance the observed excess, whereas in others only part of the excess is compensated.32,33 We found evidence of harvesting in both the Mediterranean and north-continental regions. The mortality displacement partially compensates for the effect of heat observed during the first week after exposure. In the Mediterranean cities, the harvesting phenomenon is more prolonged, but in both regions the cumulative effect at lag 25 is around 30% of the cumulative effect at lag 3. The evidence of harvesting suggests the presence of subgroups of susceptible individuals for which heat precipitates deaths by a few days to weeks.
An overall heat effect was clearly present from June to August. However, the time-varying coefficients models suggested a greater effect of earlier exposures. This evidence was based on few observations (there were few days during the first months of the warm season with a maximum apparent temperature above the threshold), but was largely consistent among cities. The stronger effect of early exposure could have various explanations. From a physiologic point of view, the human organism may react better to later heat exposure due to acclimatization,34 although this finding could also be explained by a changing composition of the population at risk over time, due to a harvesting-like phenomenon.
We adjusted for the confounding effect of air pollution. The discrepancy between adjusted and unadjusted effect estimates varied among cities. For example, in Athens, adjustment changed the percent increase in total mortality with 1°C rise in maximum apparent temperature from 6.2 to 5.5, while in Stockholm, adjustment had a negligible effect.
Future studies should investigate possible effect modifiers that could explain heterogeneity in threshold and slope estimates, both between and within regions. Varying acclimatization to weather changes could explain the observed heterogeneity.35 However the observed differences among cities and between regions could also reflect demographic, cultural, socioeconomic, and technological circumstances that produce different proportions of susceptible persons in the enrolled cities.21 In this sense, the estimated exposure-response curves can be considered a mixture of curves, each of which describes the relationship between apparent temperature and mortality in a subgroup of the population. Also meteorologic and geographic characteristics should be studied as possible contributors to heterogeneity among cities, insofar as they can modify mechanisms leading to acclimatization or interact with temperature, modifying the effect of heat.
1. McMichael AJ, Woodruff RE, Hales S. Climate change and human health: present and future risks. Lancet
2. Kosatsky T. The 2003 European heat waves. Eurosurveillance
3. WHO/WMO/UNEP. Climate change and human health: risk and responses.
Geneva, Swizterland: World Health Organization; 2003.
4. Schwartz J, Samet JM, Patz JA. Hospital admissions for heart disease: the effects of temperature and humidity. Epidemiology
5. Basu R, Samet JM. Relation between elevated ambient temperature and mortality: a review of the epidemiologic evidence. Epidemiol Rev
6. Curriero FC, Heiner KS, Samet JM, et al. Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol
7. Michelozzi P, De Sario M, Accetta G, et al. Temperature and summer mortality: geographical and temporal variations in four Italian cities. J Epidemiol Community Health
8. Choi GY, Choi JN, Kwon HJ. The impact of high apparent temperature on the increase of summertime disease-related mortality in Seoul: 1991–2000. J Prev Med Pub Health
9. Stafoggia M, Forastiere F, Agostini D, et al. Vulnerability to heat-related mortality: a multicity, population-based, case–crossover analysis. Epidemiology
10. Michelozzi P, Kirchmayer U, Katsouyanni K, et al. Assessment and prevention of acute health effects of weather conditions in Europe, the PHEWE project: background, objectives, design. Environ Health
11. Steadman RG. The assessment of sultriness. Part II: Effects of wind, extra radiation and barometric pressure on apparent temperature. J Appl Meteorol.
12. Kalkstein LS, Valimont KM. An evaluation of summer discomfort in the United States using a relative climatological index. Bull Am Meteorol Soc
13. Diggle P, Liang KY, Zeger SL. Analysis of Longitudinal Data.
Oxford-New York-Tokyo: Oxford University Press; 1994.
14. Schwartz J, Dockery DW. Increased mortality in Philadelphia associated with daily air pollution concentrations. Am Rev Respir Dis
15. Hajat S, Armstrong B, Baccini M, et al. Impact of high temperatures on mortality: is there an added heat wave effect? Epidemiology
16. Fouillet A, Rey G, Jougla E, et al. A predictive model relating daily fluctuations in summer temperatures and mortality rates. BMC Public Health
17. Gorjanc ML, Flanders WD, VanDerslice J, et al. Effects of temperature and snowfall on mortality in Pennsylvania. Am J Epidemiol
18. Xie M, Yang Y. Asymptotics for generalized estimating equations with large cluster sizes. Ann Statist
19. Chiogna M, Gaetan C. Mining epidemiological time series: an approach based on dynamic regression. Statist Model
20. 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
21. Medina-Ramón M, Schwartz J. Temperature, temperature extremes, and mortality: a study of acclimatization and effect modification in 50 United States cities. Occup Environ Med
22. Green PJ, Silverman BW. Nonparametric Regression and Generalized Linear Models.
London: Chapman and Hall; 1994.
23. Muggeo VMR. Estimating regression models with unknown breakpoint. Statist Med
24. Almon S. The distributed lag between capital appropriations and expenditures. Econometrica
25. Peng RD, Dominici F, Pastor-Barriuso R, et al. Seasonal analyses of air pollution and mortality in 100 US cities. Am J Epidemiol
26. Sutton AJ, Abrams KR. Bayesian methods in meta-analysis and evidence synthesis. Stat Methods Med Res
27. Everson PJ, Morris CN. Inference for multivariate normal hierarchical models. J R Stat Soc Ser A Stat Soc
28. Keatinge WR, Donaldson GC, Cordioli E, et al. Heat related mortality in warm and cold regions of Europe: observational study. BMJ
29. Hajat S, Kovats RS, Lachowycz K. Heat-related and cold-related deaths in England and Wales: who is at risk? Occup Environ Med
30. 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
31. Braga A, Zanobetti A, Schwartz J. The effect of weather on respiratory and cardiovascular deaths in 12 US cities. Environ Health Perspect
32. Hajat S, Armstrong BG, Gouveia N, et al. Mortality displacement of heat-related deaths: a comparison of Delhi, Sao Paulo, and London. Epidemiology
33. Braga A, Zanobetti A, Schwartz J. The time course of weather-related deaths. Epidemiology
34. Bouchama A, Knochel JP. Heat stroke. N Engl J Med
35. Ellis FP. Heat illness. III. Acclimatization. Trans R Soc Trop Med Hyg