Quantifying excess deaths related to heatwaves under climate change scenarios: A multicountry time series modelling study.

BACKGROUND: Heatwaves are a critical public health problem. There will be an increase in the frequency and severity of heatwaves under changing climate. However, evidence about the impacts of climate change on heatwave-related mortality at a global scale is limited. METHODS AND FINDINGS: We collected historical daily time series of mean temperature and mortality for all causes or nonexternal causes, in periods ranging from January 1, 1984, to December 31, 2015, in 412 communities within 20 countries/regions. We estimated heatwave-mortality associations through a two-stage time series design. Current and future daily mean temperature series were projected under four scenarios of greenhouse gas emissions from 1971-2099, with five general circulation models. We projected excess mortality in relation to heatwaves in the future under each scenario of greenhouse gas emissions, with two assumptions for adaptation (no adaptation and hypothetical adaptation) and three scenarios of population change (high variant, median variant, and low variant). Results show that, if there is no adaptation, heatwave-related excess mortality is expected to increase the most in tropical and subtropical countries/regions (close to the equator), while European countries and the United States will have smaller percent increases in heatwave-related excess mortality. The higher the population variant and the greenhouse gas emissions, the higher the increase of heatwave-related excess mortality in the future. The changes in 2031-2080 compared with 1971-2020 range from approximately 2,000% in Colombia to 150% in Moldova under the highest emission scenario and high-variant population scenario, without any adaptation. If we considered hypothetical adaptation to future climate, under high-variant population scenario and all scenarios of greenhouse gas emissions, the heatwave-related excess mortality is expected to still increase across all the countries/regions except Moldova and Japan. However, the increase would be much smaller than the no adaptation scenario. The simple assumptions with respect to adaptation as follows: no adaptation and hypothetical adaptation results in some uncertainties of projections. CONCLUSIONS: This study provides a comprehensive characterisation of future heatwave-related excess mortality across various regions and under alternative scenarios of greenhouse gas emissions, different assumptions of adaptation, and different scenarios of population change. The projections can help decision makers in planning adaptation and mitigation strategies for climate change.

Health and Economic Benefits of Complying With the World Health Organization Air Quality Guidelines for Particulate Matter in Nine Major Latin American Cities.

Objectives: This study aims to estimate the short-term preventable mortality and associated economic costs of complying with the World Health Organization (WHO) air quality guidelines (AQGs) limit values for PM10 and PM2.5 in nine major Latin American cities. Methods: We estimated city-specific PM-mortality associations using time-series regression models and calculated the attributable mortality fraction. Next, we used the value of statistical life to calculate the economic benefits of complying with the WHO AQGs limit values. Results: In most cities, PM concentrations exceeded the WHO AQGs limit values more than 90% of the days. PM10 was found to be associated with an average excess mortality of 1.88% with concentrations above WHO AQGs limit values, while for PM2.5 it was 1.05%. The associated annual economic costs varied widely, between US$ 19.5 million to 3,386.9 million for PM10, and US$ 196.3 million to 2,209.6 million for PM2.5. Conclusion: Our findings suggest that there is an urgent need for policymakers to develop interventions to achieve sustainable air quality improvements in Latin America. Complying with the WHO AQGs limit values for PM10 and PM2.5 in Latin American cities would substantially benefits for urban populations.

Mortality burden and economic loss attributable to cold and heat in Central and South America.

Background: We quantify the mortality burden and economic loss attributable to nonoptimal temperatures for cold and heat in the Central and South American countries in the Multi-City Multi-Country (MCC) Collaborative Research Network. Methods: We collected data for 66 locations from 13 countries in Central and South America to estimate location-specific temperature-mortality associations using time-series regression with distributed lag nonlinear models. We calculated the attributable deaths for cold and heat as the 2.5th and 97.5th temperature percentiles, above and below the minimum mortality temperature, and used the value of a life year to estimate the economic loss of delayed deaths. Results: The mortality impact of cold varied widely by country, from 9.64% in Uruguay to 0.22% in Costa Rica. The heat-attributable fraction for mortality ranged from 1.41% in Paraguay to 0.01% in Ecuador. Locations in arid and temperate climatic zones showed higher cold-related mortality (5.10% and 5.29%, respectively) than those in tropical climates (1.71%). Arid and temperate climatic zones saw lower heat-attributable fractions (0.69% and 0.58%) than arid climatic zones (0.92%). Exposure to cold led to an annual economic loss of $0.6 million in Costa Rica to $472.2 million in Argentina. In comparison, heat resulted in economic losses of $0.05 million in Ecuador to $90.6 million in Brazil. Conclusion: Most of the mortality burden for Central and South American countries is caused by cold compared to heat, generating annual economic losses of $2.1 billion and $290.7 million, respectively. Public health policies and adaptation measures in the region should account for the health effects associated with nonoptimal temperatures.

Temperature-mortality associations by age and cause: a multi-country multi-city study.

Background: Heterogeneity in temperature-mortality relationships across locations may partly result from differences in the demographic structure of populations and their cause-specific vulnerabilities. Here we conduct the largest epidemiological study to date on the association between ambient temperature and mortality by age and cause using data from 532 cities in 33 countries. Methods: We collected daily temperature and mortality data from each country. Mortality data was provided as daily death counts within age groups from all, cardiovascular, respiratory, or noncardiorespiratory causes. We first fit quasi-Poisson regression models to estimate location-specific associations for each age-by-cause group. For each cause, we then pooled location-specific results in a dose-response multivariate meta-regression model that enabled us to estimate overall temperature-mortality curves at any age. The age analysis was limited to adults. Results: We observed high temperature effects on mortality from both cardiovascular and respiratory causes compared to noncardiorespiratory causes, with the highest cold-related risks from cardiovascular causes and the highest heat-related risks from respiratory causes. Risks generally increased with age, a pattern most consistent for cold and for nonrespiratory causes. For every cause group, risks at both temperature extremes were strongest at the oldest age (age 85 years). Excess mortality fractions were highest for cold at the oldest ages. Conclusions: There is a differential pattern of risk associated with heat and cold by cause and age; cardiorespiratory causes show stronger effects than noncardiorespiratory causes, and older adults have higher risks than younger adults.

Joint effect of heat and air pollution on mortality in 620 cities of 36 countries.

BACKGROUND: The epidemiological evidence on the interaction between heat and ambient air pollution on mortality is still inconsistent. OBJECTIVES: To investigate the interaction between heat and ambient air pollution on daily mortality in a large dataset of 620 cities from 36 countries. METHODS: We used daily data on all-cause mortality, air temperature, particulate matter ≤ 10 µm (PM10), PM ≤ 2.5 µm (PM2.5), nitrogen dioxide (NO2), and ozone (O3) from 620 cities in 36 countries in the period 1995-2020. We restricted the analysis to the six consecutive warmest months in each city. City-specific data were analysed with over-dispersed Poisson regression models, followed by a multilevel random-effects meta-analysis. The joint association between air temperature and air pollutants was modelled with product terms between non-linear functions for air temperature and linear functions for air pollutants. RESULTS: We analyzed 22,630,598 deaths. An increase in mean temperature from the 75th to the 99th percentile of city-specific distributions was associated with an average 8.9 % (95 % confidence interval: 7.1 %, 10.7 %) mortality increment, ranging between 5.3 % (3.8 %, 6.9 %) and 12.8 % (8.7 %, 17.0 %), when daily PM10 was equal to 10 or 90 µg/m3, respectively. Corresponding estimates when daily O3 concentrations were 40 or 160 µg/m3 were 2.9 % (1.1 %, 4.7 %) and 12.5 % (6.9 %, 18.5 %), respectively. Similarly, a 10 µg/m3 increment in PM10 was associated with a 0.54 % (0.10 %, 0.98 %) and 1.21 % (0.69 %, 1.72 %) increase in mortality when daily air temperature was set to the 1st and 99th city-specific percentiles, respectively. Corresponding mortality estimate for O3 across these temperature percentiles were 0.00 % (-0.44 %, 0.44 %) and 0.53 % (0.38 %, 0.68 %). Similar effect modification results, although slightly weaker, were found for PM2.5 and NO2. CONCLUSIONS: Suggestive evidence of effect modification between air temperature and air pollutants on mortality during the warm period was found in a global dataset of 620 cities.

Heat-related cardiorespiratory mortality: Effect modification by air pollution across 482 cities from 24 countries.

BACKGROUND: Evidence on the potential interactive effects of heat and ambient air pollution on cause-specific mortality is inconclusive and limited to selected locations. OBJECTIVES: We investigated the effects of heat on cardiovascular and respiratory mortality and its modification by air pollution during summer months (six consecutive hottest months) in 482 locations across 24 countries. METHODS: Location-specific daily death counts and exposure data (e.g., particulate matter with diameters ≤ 2.5 µm [PM2.5]) were obtained from 2000 to 2018. We used location-specific confounder-adjusted Quasi-Poisson regression with a tensor product between air temperature and the air pollutant. We extracted heat effects at low, medium, and high levels of pollutants, defined as the 5th, 50th, and 95th percentile of the location-specific pollutant concentrations. Country-specific and overall estimates were derived using a random-effects multilevel meta-analytical model. RESULTS: Heat was associated with increased cardiorespiratory mortality. Moreover, the heat effects were modified by elevated levels of all air pollutants in most locations, with stronger effects for respiratory than cardiovascular mortality. For example, the percent increase in respiratory mortality per increase in the 2-day average summer temperature from the 75th to the 99th percentile was 7.7% (95% Confidence Interval [CI] 7.6-7.7), 11.3% (95%CI 11.2-11.3), and 14.3% (95% CI 14.1-14.5) at low, medium, and high levels of PM2.5, respectively. Similarly, cardiovascular mortality increased by 1.6 (95%CI 1.5-1.6), 5.1 (95%CI 5.1-5.2), and 8.7 (95%CI 8.7-8.8) at low, medium, and high levels of O3, respectively. DISCUSSION: We observed considerable modification of the heat effects on cardiovascular and respiratory mortality by elevated levels of air pollutants. Therefore, mitigation measures following the new WHO Air Quality Guidelines are crucial to enhance better health and promote sustainable development.

Effect modification of greenness on the association between heat and mortality: A multi-city multi-country study.

BACKGROUND: Identifying how greenspace impacts the temperature-mortality relationship in urban environments is crucial, especially given climate change and rapid urbanization. However, the effect modification of greenspace on heat-related mortality has been typically focused on a localized area or single country. This study examined the heat-mortality relationship among different greenspace levels in a global setting. METHODS: We collected daily ambient temperature and mortality data for 452 locations in 24 countries and used Enhanced Vegetation Index (EVI) as the greenspace measurement. We used distributed lag non-linear model to estimate the heat-mortality relationship in each city and the estimates were pooled adjusting for city-specific average temperature, city-specific temperature range, city-specific population density, and gross domestic product (GDP). The effect modification of greenspace was evaluated by comparing the heat-related mortality risk for different greenspace groups (low, medium, and high), which were divided into terciles among 452 locations. FINDINGS: Cities with high greenspace value had the lowest heat-mortality relative risk of 1·19 (95% CI: 1·13, 1·25), while the heat-related relative risk was 1·46 (95% CI: 1·31, 1·62) for cities with low greenspace when comparing the 99th temperature and the minimum mortality temperature. A 20% increase of greenspace is associated with a 9·02% (95% CI: 8·88, 9·16) decrease in the heat-related attributable fraction, and if this association is causal (which is not within the scope of this study to assess), such a reduction could save approximately 933 excess deaths per year in 24 countries. INTERPRETATION: Our findings can inform communities on the potential health benefits of greenspaces in the urban environment and mitigation measures regarding the impacts of climate change. FUNDING: This publication was developed under Assistance Agreement No. RD83587101 awarded by the U.S. Environmental Protection Agency to Yale University. It has not been formally reviewed by EPA. The views expressed in this document are solely those of the authors and do not necessarily reflect those of the Agency. EPA does not endorse any products or commercial services mentioned in this publication. Research reported in this publication was also supported by the National Institute on Minority Health and Health Disparities of the National Institutes of Health under Award Number R01MD012769. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Also, this work has been supported by the National Research Foundation of Korea (2021R1A6A3A03038675), Medical Research Council-UK (MR/V034162/1 and MR/R013349/1), Natural Environment Research Council UK (Grant ID: NE/R009384/1), Academy of Finland (Grant ID: 310372), European Union's Horizon 2020 Project Exhaustion (Grant ID: 820655 and 874990), Czech Science Foundation (22-24920S), Emory University's NIEHS-funded HERCULES Center (Grant ID: P30ES019776), and Grant CEX2018-000794-S funded by MCIN/AEI/ 10.13039/501100011033 The funders had no role in the design, data collection, analysis, interpretation of results, manuscript writing, or decision to publication.

Fluctuating temperature modifies heat-mortality association around the globe.

Studies have investigated the effects of heat and temperature variability (TV) on mortality. However, few assessed whether TV modifies the heat-mortality association. Data on daily temperature and mortality in the warm season were collected from 717 locations across 36 countries. TV was calculated as the standard deviation of the average of the same and previous days' minimum and maximum temperatures. We used location-specific quasi-Poisson regression models with an interaction term between the cross-basis term for mean temperature and quartiles of TV to obtain heat-mortality associations under each quartile of TV, and then pooled estimates at the country, regional, and global levels. Results show the increased risk in heat-related mortality with increments in TV, accounting for 0.70% (95% confidence interval [CI]: -0.33 to 1.69), 1.34% (95% CI: -0.14 to 2.73), 1.99% (95% CI: 0.29-3.57), and 2.73% (95% CI: 0.76-4.50) of total deaths for Q1-Q4 (first quartile-fourth quartile) of TV. The modification effects of TV varied geographically. Central Europe had the highest attributable fractions (AFs), corresponding to 7.68% (95% CI: 5.25-9.89) of total deaths for Q4 of TV, while the lowest AFs were observed in North America, with the values for Q4 of 1.74% (95% CI: -0.09 to 3.39). TV had a significant modification effect on the heat-mortality association, causing a higher heat-related mortality burden with increments of TV. Implementing targeted strategies against heat exposure and fluctuant temperatures simultaneously would benefit public health.

Short term associations of ambient nitrogen dioxide with daily total, cardiovascular, and respiratory mortality: multilocation analysis in 398 cities.

OBJECTIVE: To evaluate the short term associations between nitrogen dioxide (NO2) and total, cardiovascular, and respiratory mortality across multiple countries/regions worldwide, using a uniform analytical protocol. DESIGN: Two stage, time series approach, with overdispersed generalised linear models and multilevel meta-analysis. SETTING: 398 cities in 22 low to high income countries/regions. MAIN OUTCOME MEASURES: Daily deaths from total (62.8 million), cardiovascular (19.7 million), and respiratory (5.5 million) causes between 1973 and 2018. RESULTS: On average, a 10 µg/m3 increase in NO2 concentration on lag 1 day (previous day) was associated with 0.46% (95% confidence interval 0.36% to 0.57%), 0.37% (0.22% to 0.51%), and 0.47% (0.21% to 0.72%) increases in total, cardiovascular, and respiratory mortality, respectively. These associations remained robust after adjusting for co-pollutants (particulate matter with aerodynamic diameter ≤10 µm or ≤2.5 µm (PM10 and PM2.5, respectively), ozone, sulfur dioxide, and carbon monoxide). The pooled concentration-response curves for all three causes were almost linear without discernible thresholds. The proportion of deaths attributable to NO2 concentration above the counterfactual zero level was 1.23% (95% confidence interval 0.96% to 1.51%) across the 398 cities. CONCLUSIONS: This multilocation study provides key evidence on the independent and linear associations between short term exposure to NO2 and increased risk of total, cardiovascular, and respiratory mortality, suggesting that health benefits would be achieved by tightening the guidelines and regulatory limits of NO2.

Geographical Variations of the Minimum Mortality Temperature at a Global Scale: A Multicountry Study.

BACKGROUND: Minimum mortality temperature (MMT) is an important indicator to assess the temperature-mortality association, indicating long-term adaptation to local climate. Limited evidence about the geographical variability of the MMT is available at a global scale. METHODS: We collected data from 658 communities in 43 countries under different climates. We estimated temperature-mortality associations to derive the MMT for each community using Poisson regression with distributed lag nonlinear models. We investigated the variation in MMT by climatic zone using a mixed-effects meta-analysis and explored the association with climatic and socioeconomic indicators. RESULTS: The geographical distribution of MMTs varied considerably by country between 14.2 and 31.1 °C decreasing by latitude. For climatic zones, the MMTs increased from alpine (13.0 °C) to continental (19.3 °C), temperate (21.7 °C), arid (24.5 °C), and tropical (26.5 °C). The MMT percentiles (MMTPs) corresponding to the MMTs decreased from temperate (79.5th) to continental (75.4th), arid (68.0th), tropical (58.5th), and alpine (41.4th). The MMTs indreased by 0.8 °C for a 1 °C rise in a community's annual mean temperature, and by 1 °C for a 1 °C rise in its SD. While the MMTP decreased by 0.3 centile points for a 1 °C rise in a community's annual mean temperature and by 1.3 for a 1 °C rise in its SD. CONCLUSIONS: The geographical distribution of the MMTs and MMTPs is driven mainly by the mean annual temperature, which seems to be a valuable indicator of overall adaptation across populations. Our results suggest that populations have adapted to the average temperature, although there is still more room for adaptation.

How urban characteristics affect vulnerability to heat and cold: a multi-country analysis.

BACKGROUND: The health burden associated with temperature is expected to increase due to a warming climate. Populations living in cities are likely to be particularly at risk, but the role of urban characteristics in modifying the direct effects of temperature on health is still unclear. In this contribution, we used a multi-country dataset to study effect modification of temperature-mortality relationships by a range of city-specific indicators. METHODS: We collected ambient temperature and mortality daily time-series data for 340 cities in 22 countries, in periods between 1985 and 2014. Standardized measures of demographic, socio-economic, infrastructural and environmental indicators were derived from the Organisation for Economic Co-operation and Development (OECD) Regional and Metropolitan Database. We used distributed lag non-linear and multivariate meta-regression models to estimate fractions of mortality attributable to heat and cold (AF%) in each city, and to evaluate the effect modification of each indicator across cities. RESULTS: Heat- and cold-related deaths amounted to 0.54% (95% confidence interval: 0.49 to 0.58%) and 6.05% (5.59 to 6.36%) of total deaths, respectively. Several city indicators modify the effect of heat, with a higher mortality impact associated with increases in population density, fine particles (PM2.5), gross domestic product (GDP) and Gini index (a measure of income inequality), whereas higher levels of green spaces were linked with a decreased effect of heat. CONCLUSIONS: This represents the largest study to date assessing the effect modification of temperature-mortality relationships. Evidence from this study can inform public-health interventions and urban planning under various climate-change and urban-development scenarios.

Mortality burden of diurnal temperature range and its temporal changes: A multi-country study.

Although diurnal temperature range (DTR) is a key index of climate change, few studies have reported the health burden of DTR and its temporal changes at a multi-country scale. Therefore, we assessed the attributable risk fraction of DTR on mortality and its temporal variations in a multi-country data set. We collected time-series data covering mortality and weather variables from 308 cities in 10 countries from 1972 to 2013. The temporal change in DTR-related mortality was estimated for each city with a time-varying distributed lag model. Estimates for each city were pooled using a multivariate meta-analysis. The results showed that the attributable fraction of total mortality to DTR was 2.5% (95% eCI: 2.3-2.7%) over the entire study period. In all countries, the attributable fraction increased from 2.4% (2.1-2.7%) to 2.7% (2.4-2.9%) between the first and last study years. This study found that DTR has significantly contributed to mortality in all the countries studied, and this attributable fraction has significantly increased over time in the USA, the UK, Spain, and South Korea. Therefore, because the health burden of DTR is not likely to reduce in the near future, countermeasures are needed to alleviate its impact on human health.

Temperature-related mortality impacts under and beyond Paris Agreement climate change scenarios.

The Paris Agreement binds all nations to undertake ambitious efforts to combat climate change, with the commitment to Bhold warming well below 2 °C in global mean temperature (GMT), relative to pre-industrial levels, and to pursue efforts to limit warming to 1.5 °C". The 1.5 °C limit constitutes an ambitious goal for which greater evidence on its benefits for health would help guide policy and potentially increase the motivation for action. Here we contribute to this gap with an assessment on the potential health benefits, in terms of reductions in temperature-related mortality, derived from the compliance to the agreed temperature targets, compared to more extreme warming scenarios. We performed a multi-region analysis in 451 locations in 23 countries with different climate zones, and evaluated changes in heat and cold-related mortality under scenarios consistent with the Paris Agreement targets (1.5 and 2 °C) and more extreme GMT increases (3 and 4 °C), and under the assumption of no changes in demographic distribution and vulnerability. Our results suggest that limiting warming below 2 °C could prevent large increases in temperature-related mortality in most regions worldwide. The comparison between 1.5 and 2 °C is more complex and characterized by higher uncertainty, with geographical differences that indicate potential benefits limited to areas located in warmer climates, where direct climate change impacts will be more discernible.

Projections of temperature-related excess mortality under climate change scenarios.

BACKGROUND: Climate change can directly affect human health by varying exposure to non-optimal outdoor temperature. However, evidence on this direct impact at a global scale is limited, mainly due to issues in modelling and projecting complex and highly heterogeneous epidemiological relationships across different populations and climates. METHODS: We collected observed daily time series of mean temperature and mortality counts for all causes or non-external causes only, in periods ranging from Jan 1, 1984, to Dec 31, 2015, from various locations across the globe through the Multi-Country Multi-City Collaborative Research Network. We estimated temperature-mortality relationships through a two-stage time series design. We generated current and future daily mean temperature series under four scenarios of climate change, determined by varying trajectories of greenhouse gas emissions, using five general circulation models. We projected excess mortality for cold and heat and their net change in 1990-2099 under each scenario of climate change, assuming no adaptation or population changes. FINDINGS: Our dataset comprised 451 locations in 23 countries across nine regions of the world, including 85 879 895 deaths. Results indicate, on average, a net increase in temperature-related excess mortality under high-emission scenarios, although with important geographical differences. In temperate areas such as northern Europe, east Asia, and Australia, the less intense warming and large decrease in cold-related excess would induce a null or marginally negative net effect, with the net change in 2090-99 compared with 2010-19 ranging from -1·2% (empirical 95% CI -3·6 to 1·4) in Australia to -0·1% (-2·1 to 1·6) in east Asia under the highest emission scenario, although the decreasing trends would reverse during the course of the century. Conversely, warmer regions, such as the central and southern parts of America or Europe, and especially southeast Asia, would experience a sharp surge in heat-related impacts and extremely large net increases, with the net change at the end of the century ranging from 3·0% (-3·0 to 9·3) in Central America to 12·7% (-4·7 to 28·1) in southeast Asia under the highest emission scenario. Most of the health effects directly due to temperature increase could be avoided under scenarios involving mitigation strategies to limit emissions and further warming of the planet. INTERPRETATION: This study shows the negative health impacts of climate change that, under high-emission scenarios, would disproportionately affect warmer and poorer regions of the world. Comparison with lower emission scenarios emphasises the importance of mitigation policies for limiting global warming and reducing the associated health risks. FUNDING: UK Medical Research Council.

Episodios de cefalea en ráfagas asociados a las fluctuaciones de temperatura en Bogotá 2009-2013 / Cluster headache attacks associated with temperature fluctuations in Bogotá 2009-2013

El calentamiento global, que conlleva un incremento en el promedio y la varianza de la temperatura a nivel mundial, se ha asociado a una gran variedad de efectos agudos y crónicos sobre la salud, los cuales han sido descritos en la ciudad de Bogotá. Este aumento de la temperatura podría estar asociado a una mayor frecuencia de eventos cardiovasculares, respiratorios y neurológicos, como la cefalea en ráfagas (CR). El objetivo de este estudio fue evaluar la asociación entre las variaciones de temperatura y la ocurrencia diaria de episodios de cefalea en ráfagas, a través del análisis ecológico de series de tiempo en el periodo comprendido entre 2009-2013 en Bogotá. Este diseño utilizó los datos diarios de consultas de urgencia por CR y su relación con los valores diarios de variables meteorológicas y contaminantes criterio durante 5 años. El análisis estadístico incluyó el resumen descriptivo de las variables, posteriormente el análisis bivariado con pruebas de correlación, autocorrelación y correlación cruzada y por último el análisis multivariado mediante una regresión dinámica de Poisson y el modelo no lineal de rezagos distribuidos. Se encontró una correlación con un comportamiento cíclico entre cefalea en ráfagas con amplitud térmica (Rho=0,054, p<0,05) y temperatura máxima (Rho=0,050, p<0,05). En el modelo multivariado, se observó una influencia alta tanto de la amplitud térmica (ß=0,11, IC95% 0,07 a 0,16), como de la temperatura máxima (ß=0,14, IC95% 0,08 a 0,20) en la incidencia de CR. La amplitud térmica y la temperatura máxima podrían ser un desencadenante de los episodios de cefalea en ráfagas, pero futuros estudios deben explorar la relación de las crisis de acuerdo a las islas de calor en Bogotá y explorar las asociaciones con estudios a nivel individual.

Global warming has been linked to the increase of both average and variance of temperature which has been associated with a variety of acute and chronic health effects and Bogotá has not been the exception. This temperature increase might be related to a greater incidence of cardiovascular, respiratory and neurological outcomes, as cluster headache. Aim of this study was to assess the association between temperature fluctuations and daily incidence of cluster headache, through a time series analysis of the cases from 2009 to 2013. An ecological design was carried out with daily emergency admissions of cluster headache and daily values of meteorological and pollutants in Bogota. Statistical analysis included a descriptive summary of variables, a bivariate analysis with correlation, crossed correlations and autocorrelations tests and multivariate analysis by Poisson regression models and distributed lag nonlinear models. A cyclic behavior was observed, with a correlation between cluster headache and thermal amplitude (Rho=0.054, p<0.05) and with maximum temperature Rho=0.050, p<0.05), which persisted in the multivariate analysis: thermal amplitude (ß=0.11, CI 95% 0.07 a 0.16), and maximum temperature (ß=0.14, IC95% 0, 08 a 0, 20). Thermal amplitude and maximum temperature might be a trigger of cluster headache, but future studies should explore +this association at individual level and in the urban heat islands.