Spatial variability of heat-related mortality in Barcelona from 1992-2015: A case crossover study design

Abstract

Numerous studies have demonstrated the relationship between summer temperatures and increased heat-related deaths. Epidemiological analyses of the health effects of climate exposures usually rely on observations from the nearest weather station to assess exposure-response associations for geographically diverse populations. Urban climate models provide high-resolution spatial data that may potentially improve exposure estimates, but to date, they have not been extensively applied in epidemiological research. We investigated temperature-mortality relationships in the city of Barcelona, and whether estimates vary among districts. We considered georeferenced individual (natural) mortality data during the summer months (June-September) for the period 1992-2015. We extracted daily summer mean temperatures from a 100-m resolution simulation of the urban climate model (UrbClim). Summer hot days (above percentile 70) and reference (below percentile 30) temperatures were compared by using a conditional logistic regression model in a case crossover study design applied to all districts of Barcelona. Relative Risks (RR), and 95% Confidence Intervals (CI), of all-cause (natural) mortality and summer temperature were calculated for several population subgroups (age, sex and education level by districts). Hot days were associated with an increased risk of death (RR = 1.13; 95% CI = 1.10-1.16) and were significant in all population subgroups compared to the non-hot days. The risk ratio was higher among women (RR = 1.16; 95% CI= 1.12-1.21) and the elderly (RR = 1.18; 95% CI = 1.13-1.22). Individuals with primary education had similar risk (RR = 1.13; 95% CI = 1.08-1.18) than those without education (RR = 1.10; 95% CI= 1.05-1.15). Moreover, 6 out of 10 districts showed statistically significant associations, varying the risk ratio between 1.12 (95% CI = 1.03-1.21) in Sants-Montjuïc and 1.25 (95% CI = 1.14-1.38) in Sant Andreu. Findings identified vulnerable districts and suggested new insights to public health policy makers on how to develop district-specific strategies to reduce risks. ; J.B. gratefully acknowledges funding from the European Union's Horizon 2020 research and innovation programme under grant agreements No 865564 (European Research Council Consolidator Grant EARLY-ADAPT), 727852 (project Blue-Action) and 730004 (project PUCS), and from the Ministry of Science and Innovation (MCIU) under grant agreements No RYC2018-025446-I (programme Ramón y Cajal) and EUR2019-103822 (project EURO-ADAPT). V.I. acknowledges funding from the European Union's Horizon 2020 research and innovation programme under grant agreements 730004 (project PUCS). H.A. gratefully acknowledges funding from the Secretariat for Universities and Research of the Ministry of Business and Knowledge of the Government of Catalonia (grant numbers B00391 [FI-2018], B100180 [FI-2019] and B200139 [FI-2020]).