- Research article
- Open Access
- Open Peer Review
Mortality on extreme heat days using official thresholds in Spain: a multi-city time series analysis
BMC Public Health volume 12, Article number: 133 (2012)
The 2003 heat wave had a high impact on mortality in Europe, which made necessary to develop heat health watch warning systems. In Spain this was carried-out by the Ministry of Health in 2004, being based on exceeding of city-specific simultaneous thresholds of minimum and maximum daily temperatures. The aim of this study is to assess effectiveness of the official thresholds established by the Ministry of Health for each provincial capital city, by quantifying and comparing the short-term effects of above-threshold days on total daily mortality.
Total daily mortality and minimum and maximum temperatures for the 52 capitals of province in Spain were collected during summer months (June to September) for the study period 1995-2004. Data was analysed using GEE for Poisson regression. Relative Risk (RR) of total daily mortality was quantified for the current day of official thresholds exceeded.
The number of days in which the thresholds were exceeded show great inconsistency, with provinces with great number of exceeded days adjacent to provinces that did not exceed or rarely exceeded. The average overall excess risk of dying during an extreme heat day was about 25% (RR = 1.24; 95% confidence interval (CI) = [1.19-1.30]). Relative risks showed a significant heterogeneity between cities (I2 = 54.9%). Western situation and low mean summer temperatures were associated with higher relative risks, suggesting thresholds may have been set too high in these areas.
This study confirmed that extreme heat days have a considerable impact on total daily mortality in Spain. Official thresholds gave consistent relative risk in the large capital cities. However, in some other cities thresholds
According to the most recent World Health Organisation executive committee report, the 2003 summer heat wave caused a total of 44,000 excess deaths in Europe as a consequence of the extremely high temperatures recorded that summer . Since that, many European countries have drawn up heat-surveillance systems to alert and prevent the negative effects of extreme heat temperatures on their health . The Spanish Government, through their Ministry of Health, started implementing in 2004 the "National plan for preventive actions against the effects of excess temperatures on health" . This Plan sets city-specific thresholds for extreme heat days during summer months (from June to September) for each of the 52 provincial capital cities. Although, extreme heat temperatures have been associated with an increase of mortality [4–7] and morbidity [8, 9] in some Spanish cities, the health impact of extreme heat days established under the Ministry's official definition has not yet been assessed. Therefore, this could be considered the first national study conducted in Spain.
The objective of this study is to assess effectiveness of the official temperature thresholds established by the Ministry of Health for each provincial capital city in identifying high-risk days by quantifying short-term effects of above-threshold days on total daily mortality in each city over the ten years prior to the operation of the Plan. In particular we explored whether the implicit above-threshold mortality risk level has been similar across cities. We also explore determinants of variation in the city-specific risks.
Spain is the second largest country in Western Europe (504,030 km2) with a total population of 46 million people. Its mainland is bordered to the South and East by the Mediterranean Sea; to the North by France, and the Bay of Biscay; and to the North-west and West by the Atlantic Ocean and Portugal. Spanish territory also includes the Balearic Islands in the Mediterranean, the Canary Islands in the Atlantic Ocean off the African coast, and two autonomous cities in Northern Africa, Ceuta and Melilla. It is divided in 17 Autonomous Regions, each composed from one up to nine provinces, and each province has a capital city (Additional file 1 shows their location and Additional file 2 reports major geographic and socio-demographic characteristics for each provincial capital city).
Due to Spain's geographical situation and conditions, the climate is extremely diverse. It can be roughly divided into five main zones: Mediterranean climate extends along the southern and eastern coasts up to the Pyrenees; Semiarid Mediterranean climate in the South-east; Continental Mediterranean climate in the inland areas of the Peninsula; Oceanic climate in the North-west and the coastal strip near the Bay of Biscay; and a Subtropical climate in the Canary Islands .
Data on total daily all-cause mortality, excluding accidents (International Classification of Diseases-9th revision/ICD-9: 1-799), for the 52 capital cities were provide by the Spain National Institute of Statistics for summer months (from 1st June to 30th September) from 1995 to 2004.
Temperature data and extreme heat days
Daily minimum and maximum temperatures (in C) for the 52 capital cities were collected from the Spain National Meteorology Agency for the same study period. In two capital cities, Palencia and Ceuta, daily temperatures were not available.
Using these data, city-specific extreme heat days were defined under the official thresholds established by the "National Plan for Preventive Actions against the Effects of Excess Temperatures on Health" of the Ministry of Health. This Plan sets city-specific thresholds as follows: "The criteria have been established based on studies sponsored by the Ministry of Health, the observations made by the Autonomous Regions and the technical details of National Meteorology Agency. Extreme heat days are based on the 95 th percentile of the 20-year historical series of daily summer (from June to September) of minimum and maximum temperatures for each of the 52 cities, which are capitals of provinces. As exceptions, for the mild climate stations with low daily temperature variations (mainly sea areas), Northern and North-western Spain, the threshold for the maximum temperature correspond to the 95 th percentile of the series of absolute maximum temperatures of each summer. Similarly, in continental climate stations, the threshold considered for the minimum temperature corresponds to the 95 th percentile of the series of higher minimum temperatures of summer. The values obtained through these calculations have been rounded to nearest integer, because the error in predicting temperature maximum and minimum is about one Celsius degree. In capital cities where the threshold for minimum temperature were below 20 C and below 33 C for maximum temperature, these values have been assigned because minimum temperatures lower than 20 C and maximum lower than 33 C do not have consequences for the purposes of this Plan"  (Author's translation). This means that original thresholds given by 95th percentiles were revised upwards in cooler cities. The activation of this Plan occurs when both thresholds, for minimum and maximum daily temperatures, are simultaneously exceeded.
Following methodology used in studies previously conducted in other European cities [2, 11, 12] the risk of mortality on extreme heat days was quantified using Generalised Estimating Equation models (GEE) with the number of daily deaths being assumed to follow a Poisson distribution. We modelled the marginal relationship between total daily mortality and a binary variable for extreme heat days in each capital city, assuming independence among summers and treating serial dependence of daily number of deaths within each summer as a nuisance parameter by specifying a first-order autoregressive correlation structure in the Poisson GEE model. Long-term trends were modelled using natural cubic splines of the variable time allowing one degree of freedom (df) for every 5-year of data. To allow for within summer seasonal variation not explained by extreme heat days, we fit natural cubic splines of day-in-year (4 df) constrained to be the same over all years.
A binary variable (0,1) was defined to identify extreme heat days for each capital city using the official thresholds, as previously described. To avoid inclusion in the baseline of days with temperatures just below the thresholds, which could still include some affected by heat, we used as a comparison category those days on which the 50th percentile of both, minimum and maximum temperatures in summer months, was not simultaneously exceeded (i.e. days of below-median minimum and maximum temperature).
Overall estimates of the extreme heat relative risks were obtained by combining the city-specific estimates with a random-effects meta-analysis, when estimates between cities were heterogeneous . To explore heterogeneity of effects we used random-effects meta-regression  to identify whether size of extreme heat relative risk was associated with city-specific geographic, socio-demographic and climatic characteristics.
All analyses were done using Stata 11.2 (StataCorp, College Station, TX, 2010).
Extreme heat days
Table 1 shows the distribution of daily temperature and the number of days when official thresholds for extreme heat days were exceeded for the study period and Figure 1 shows its geographical distribution identifying if exceptions to the use of 95th percentile had been used. These range from 0 (in Teruel, Santander, Guadalajara, Soria and Lugo) to 112 days (in Tarragona). In half of these cities the threshold was exceeded less than 10 days, while a fourth part of the cities was exceeded between 10 and 30 days, otherwise in more than 30 days. The year 2003 was when the limit was exceeded in a greater number of capital cities (35), followed by the year 1995 (29 cities), whilst the years where the threshold was exceeded by a smaller number of capital cities were in 1997 and 1996 (6 and 14, respectively).
Risk of mortality
Figure 2 shows the relative risk (RR) of dying on days on which the official thresholds, for both minimum and maximum temperatures, were exceeded with respect to a baseline comprising those days on which the maximum and minimum temperatures were simultaneously below the 50th percentile. The overall excess risk of dying during an extreme heat day was about 25% (RR = 1.24; 95% confidence interval (CI) = [1.19-1.30]).
Estimates between cites showed a significant heterogeneity (I2 = 54.9%, p < 0.001). The highest city-specific risks were found for Zamora (RR = 2.70; 95%CI = [1.50-4.86]), Oviedo (RR = 2.08; 95%CI = [1.30-3.31]) and Cuenca (RR = 1.87; 95%CI = [1.26-2.78]). The largest cities showed similar risk estimates, close to 30% (Madrid: RR = 1.31; 95%CI = [1.24-1.39], Barcelona RR = 1.32; 95%CI = [1.22-1.43], and Seville RR = 1.32; 95%CI = [1.25-1.40]). The geographical pattern (Figure 3) seems to be a trend from East, with lower risks, to West, with higher risks. Comparing mean risk estimates between Autonomous Regions showed higher values for central regions and lower values for those in the Mediterranean area.
To explain the variation in risk of mortality due to above-threshold days across cities we considered geographic (latitude, longitude, and altitude), socio-demographic (population density, proportion of people older than 65 years, and per capita personal income) and climatic (mean temperature, mean humidity, and total hours of sunshine in summer) characteristics for each capital city. Although many of these were significant determinants of extreme heat risk when considered singly, in a multiple meta-regression we found that a model including only latitude and mean temperature in summer explained most of the heterogeneity between city-specific risk estimates (residual heterogeneity, I2 = 28.3%), leaving no other factors significantly associated with heat risk. Western situation and low mean summer temperatures were associated with higher relative risks. We split latitude and mean temperature in two levels, below and above their respective median values of -3.4° (running slightly East of Madrid) and 23 C. The cross-tabulation of both defined four groups that roughly identify the climatic zones in the Iberian Peninsula (Figure 4): Oceanic climate (RR = 1.57: 95%CI = [1.39-1.79]), Semiarid Mediterranean climate (RR = 1.27: 95%CI = [1.19-1.34]), Continental Mediterranean (RR = 1.29: 95%CI = [1.23-1.36]), and Mediterranean climate (RR = 1.12: 95%CI = [1.04-1.20]).
Our results confirm that extreme heat days have a clear impact on total mortality in Spain. This has also been found in other multi-city European Studies [12, 15, 16], even in temperate climates . To our knowledge, this is the first study using official thresholds defined by a heat health warning system.
These have been mainly based on the 95th percentile (though with modifications in some scenarios, as previously described) of the 20-years historic series of minimum and maximum daily summer temperatures in each capital city. However, this criterion was only previously associated with health effects in the largest cities of Madrid , Barcelona  and Seville , whereas in others was not validated. Although the number of excesses was relatively homogeneous between the provincial capital cities in each Autonomous Region, there were substantial differences between regions. There were also a few with a surprisingly lower or higher than expected in comparison to other cities of the same Region.
There are several reasons that cities might show atypical relative risks during the supra-threshold days. Firstly, the official thresholds may not be a good indicator, for that city, of days at high risk. Specifically, the 95th percentile may not be a good indicator of risk or the exceptions made to the 95th percentile in the official thresholds criteria were inadvisable. Some of the cities with unusually high risks do appear to have been such exceptions. For example, Oviedo (RR = 2.08, 95%CI = [1.30-3.31]) had both the minimum and maximum threshold revised upwards to their 'floor' levels of 20 C and 33 C, and Zamora (RR = 2.70, 95%CI = [1.50-4.86]) had a minimum temperature threshold (22 C) appreciably higher that the 95th percentile (18.6 C), indicating that the special rule for continental climates had operated. More generally, the meta-regression analysis showed that low mean summer temperature, which was where the exception rules applied, was a strong predictor for high relative risks. Together, these suggest that the upwards revision of thresholds in cooler cities lower than would be appropriate if the intention was to identify days which the same excess risk as other cities on supra-threshold days.
To shed further light on performance of the official thresholds in identifying days of greatest heat-related events, we compared the percentage increase in risk of death on extreme heat days in 2003, calculated in our study (results not shown) with a previous study, by Martinez-Navarro et al. , where observed mortality for the summer of 2003 (from 1st June to 31st August) was compared with the mortality that would have been expected on the basis of historic series since 1990. We found a very weak correlation (r = 0.09) between estimates of both studies, perchance due to imprecision in the two estimates. However, if we only consider the largest cites (over 400,000 inhabitants) this correlation becomes considerably larger (r = 0.59). Indicating that in all these the original thresholds given by 95th percentiles were not revised upwards.
Identifying heat-health thresholds is a practical decision that should respond to credibility, precision and cost criteria . In this study we have not sought to explore all aspects of this nor review all the ways this can be done, so that we retain focus on the simple question of whether the official thresholds have been successful in identifying, in each provincial capital city, days with elevated risks of mortality and whether extent of elevation was about the same level across cities. For the same reason we did not explore factors that might modify the extreme heat impact, such period of time in the summer and duration of the extreme heat events, nor did we compare extremely hot to all other days, as done in some studies. Therefore our estimates are not entirely comparable with other published studies that quantified health effects of extreme hot temperatures.
The geographical variability in the distribution of the extreme heat-mortality effects roughly identifies the climatic zones in the Iberian Peninsula. This could partially be explained by the fact that city-specific risk estimates are calculated using official thresholds, which are at different percentiles of the temperature distribution, which may have added 'artificial' heterogeneity between climatic zones. However, when we used a fixed 95% percentile thresholds heterogeneity in relative risks did not differ substantially (I2 = 51.4% vs. 54.9%). Thus there may still be a case for different thresholds for provincial capital cities based on their climatic region. Effects of high daily temperatures on mortality in English regions have been predicted from the region's climate at the 93 rd percentile of summer mean temperature . Similarly, regional thresholds have been defined for heat-health warning systems in France  and Australia . Age, urban environment and low socio-economic level have been identified as risk factors associated with heat-related mortality [22, 23].
A major strength of the current study was the availability of long time series data sets from all provincial capital cities of Spain, providing enough power to identify and estimate extreme heat effects in all climatic regions of Spain. This study also includes the 2003 summer, which was an unusually hot summer in Europe. For this reason, we tested the sensitivity of our results by removing 2003 from the analysis. However, the overall risk estimate did not differ substantially (RR = 1.21; 95%CI = [1.16-1.27]). This may be explained by the fact that 1996 summer also recorded similar high temperatures in Spain, especially in Central Regions. Overall risk estimate for 1996 (RR = 1.29; 95%CI = [1.08-1.53]) was close to 2003 (RR = 1.31; 95%CI = [1.22-1.40]). Therefore, in a 10 years period we found another extremely hot summer, besides 2003, with a high impact on mortality. This agrees with the latest report issued by the Intergovernmental Panel for Climate Change indicates that climate change will lead to an increase in the frequency and intensity of heat waves . Predictions for the Iberian Peninsula, using general circulation models, indicate a uniform increase in temperature over the course of the 21st century, with an average upward trend every 10 years between 0.6 C to 0.7 C in summer, a greater number and higher frequency of days with extreme heat temperatures in summer .
However, ecological time-series studies can be affected by various biases, in particular potential confounding due temporal patterns. In this study well-established methods were used to control for trend and seasonality. Additionally, some air pollutant levels, which are higher in summer in many European areas , are associated with an increase in mortality and may be associated with hot temperatures as well. In our case, it is difficult to ascertain whether the triggering factor for death in extreme heat days is exclusively due to an increase in temperature or to an increase of air pollution levels that tend to accompany such high temperatures. Unfortunately, due to lack of data our study was unable to assess the influence of air pollutants and other weather variables, such as humidity (although there is little evidence that it is associated with mortality) .
Although official thresholds gave consistent relative risk in the large capital cities, mainly because from previous studies [4–7] the 95th percentile seems to be a 'natural' threshold of triggering mortality for heat days, in some other cities heat thresholds should be updated to better represent impact on health. In this sense, a very recent study conducted in Castilla La Mancha region showed how these 'natural' heat thresholds for maximum temperature range from the 92nd percentile in Cuenca up to 97th percentile in Albacete and Toledo . Therefore, criterion for the choice of thresholds must be based on epidemiological instead only on meteorological basis. In this sense, further research needs to be done to describe and understand how is the relationship of heat and mortality in these cities.
This study confirmed that extreme heat days, as defined by the official thresholds of the "National plan for preventive actions against the effects of excess temperatures on health", have a considerable impact on total daily mortality in Spain. However, relative risks on above-threshold days varied across the climatic zones of the Iberian Peninsula being higher in Northern cities with cooler climates, probably because thresholds were reset to values above the original 95th percentiles in these cities. These should be updated to identify days of equally high risk.
World Health Organisation Executive Board: Climate change and health. 2008, Geneva
Michelozzi P, Kirchmayer U, Katsouyanni K, Biggeri A, McGregor G, Menne B: Assessment and prevention of acute health effects of weather conditions in Europe, the PHEWE project: background, objectives, design. Environ Health. 2007, 6: 12-10.1186/1476-069X-6-12.
Ministerio de Sanidad y Consumo: Plan Nacional de acciones preventivas de los efectos del exceso de temperaturas sobre la salud. 2009, Madrid, [http://www.msc.es/ciudadanos/saludAmbLaboral/planAltasTemp/2009/home.htm, last accessed 1st July 2011]
Díaz J, García R, Velázquez de Castro F, Hernández E, López C, Otero A: Effects of extremely hot days on people older than 65 years in Seville (Spain) from 1986 to 1997. Int J Biometeorol. 2002, 46: 145-149. 10.1007/s00484-002-0129-z.
Díaz J, Linares C, Tobías A: Impact of extreme temperatures on daily mortality in Madrid (Spain) among the 45-64 age-group. Int J Biometeorol. 2006, 50: 342-348. 10.1007/s00484-006-0033-z.
Tobías A, de Olalla PG, Linares C, Bleda MJ, Caylà JA, Díaz J: Short-term effects of extreme hot summer temperatures on total daily mortality in Barcelona, Spain. Int J Biometeorol. 2010, 54: 115-117. 10.1007/s00484-009-0266-8.
Mirón IJ, Montero JC, Criado-Alvarez JJ, Díaz J, Linares C: Efectos de los extremos térmicos sobre la mortalidad diaria en Castilla-La Mancha: evolución temporal 1975-2003. Gac Sanit. 2010, 24: 117-122. 10.1016/j.gaceta.2009.10.016.
Linares C, Díaz J: Impact of high temperatures on hospital admissions: comparative analysis with previous studies about mortality in Madrid. Eur J Public Health. 2008, 18: 317-322. 10.1093/eurpub/ckm108.
García-Pina R, Tobías Garcés A, Sanz Navarro J, Navarro Sánchez C, García-Fulgueiras A: Effect of weather temperature on hospital emergencies in the Region of Murcia, Spain, throughout the 2000-2005 and its use in epidemiological surveillance. Rev Esp Salud Publica. 2008, 82: 153-166. 10.1590/S1135-57272008000200002.
Prieto L, García-Herrera R, Díaz J, Hernández Tesoa T: Minimum extreme temperatures over Peninsular Spain. Glob Planet Chang. 2004, 44: 59-71. 10.1016/j.gloplacha.2004.06.005.
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: 632-638. 10.1097/01.ede.0000239688.70829.63.
D'Ippoliti D, Michelozzi P, Marino C, De'donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA: The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health. 2010, 9: 37-10.1186/1476-069X-9-37.
DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7: 177-188. 10.1016/0197-2456(86)90046-2.
Berkey CS, Hoaglin DC, Mosteller F, Colditz GA: Random effects regression model for meta-analysis. Stat Med. 1995, 14: 396-411.
Stafoggia M, Forastiere F, Agostini D, Caranci N, de' Donato F, Demaria M, Michelozzi P, Miglio R, Rognoni M, Russo A, Perucci CA: Factors affecting in-hospital heat-related mortality: a multi-city case-crossover analysis. J Epidemiol Community Health. 2008, 62: 209-215. 10.1136/jech.2007.060715.
Almeida SP, Casimiro E, Calheiros J: Effects of apparent temperature on daily mortality in Lisbon and Oporto. Portugal Environ Health. 2010, 9: 12-
Armstrong BG, Chalabi Z, Fenn B, Hajat S, Kovats S, Milojevic A, Wilkinson P: Association of mortality with high temperatures in a temperate climate: England and Wales. J Epidemiol Community Health. 2011, 65: 340-345. 10.1136/jech.2009.093161.
Martinez-Navarro F, Simon-Soria F, Lopez Abente G: Valoracion del impacto de la ola de calor del verano de 2003 sobre la mortalidad. Gac Sanit. 2004, 18 (Supl 1): 250-258.
Koppe C, Kovats RS, Jendritzky G, Menne B: Heat-waves: risks and responses. 2004, Health and Global Environmental Change, series No. 2. Geneva, WHO
Fouillet A, Rey G, Wagner V, Laaidi K, Empereur-Bissonnet P, Le Tertre A, Frayssinet P, Bessemoulin P, Laurent F, De Crouy-Chanel P, Jougla E, Hémon D: Has the impact of heat waves on mortality changed in France since the European heat wave of summer 2003? A study of the 2006 heat wave. Int J Epidemiol. 2008, 37: 309-317. 10.1093/ije/dym253.
Nicholls N, Skinner C, Loughnan M, Tapper N: A simple heat alert system for Melbourne, Australia. Int J Biometeorol. 2008, 52: 375-84. 10.1007/s00484-007-0132-5.
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: 93-100.
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: 315-323. 10.1097/01.ede.0000208477.36665.34.
Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report: Climate Change 2007. 2007, Geneva
Castro M, Marín-Vide J, Alonso S: El Clima en España: Pasado, presente y escenarios de clima para el siglo XXI. Evaluación Preliminar de los Impactos en España por Efecto del Cambio Climático. Edited by: Moreno JME. 2005, Madrid: Ministerio de Medio Ambiente
Montero JC, Mirón IJ, Criado-Álvarez JJ, Linares C, Díaz J: Influence of local factors in the relationship between mortality and heat waves: Castile-La Mancha (1975-2003). Sci Total Environ. 2012, 414: 73-80.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/12/133/prepub
This study was funded by grant SEPY 1106/07 of the Instituto de Salud Carlos III (ISCIII) and by grant 200930I008 of the Consejo Superior de Investigaciones Científicas (CSIC). To Fernando Simon, from the Alert and Emergency Unit of the National Centre of Epidemiology, to provide data.
The authors declare that they have no competing interests.
AT, IZ, CL, and JD conceived the study. AT analyzed the data in consultation with BA, and AG. AT, and IZ wrote the draft version and revisions of the manuscript according to the contribution of BA, AG, CL, and JD. All authors read and approved the final version of the manuscript.