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. © 2022 The Author(s)Funding text 1: This study was supported by the Australian Research Council (DP210102076) and the Australian National Health and Medical Research Council (APP2000581). Y.W and B.W. were supported by the China Scholarship Council (nos. 202006010044 and 202006010043); S.L. was supported by an Emerging Leader Fellowship of the Australian National Health and Medical Research Council (no. APP2009866); Y.G. was supported by Career Development Fellowship (no. APP1163693) and Leader Fellowship (no. APP2008813) of the Australian National Health and Medical Research Council; J.K. and A.U. were supported by the Czech Science Foundation (project no. 20–28560S); N.S. was supported by the National Institute of Environmental Health Sciences-funded HERCULES Center (no. P30ES019776); Y.H. was supported by the Environment Research and Technology Development Fund (JPMEERF15S11412) of the Environmental Restoration and Conservation Agency; M.d.S.Z.S.C. and P.H.N.S. were supported by the São Paulo Research Foundation (FAPESP); H.O. and E.I. were supported by the Estonian Ministry of Education and Research (IUT34–17); J.M. was supported by a fellowship of Fundação para a Ciência e a Tecnlogia (SFRH/BPD/115112/2016); A.G. and F.S. were supported by the Medical Research Council UK (grant ID MR/R013349/1), the Natural Environment Research Council UK (grant ID NE/R009384/1), and the EU's Horizon 2020 project, Exhaustion (grant ID 820655); A.S. and F.d.D. were supported by the EU's Horizon 2020 project, Exhaustion (grant ID 820655); V.H. was supported by the Spanish Ministry of Economy, Industry and Competitiveness (grant ID PCIN-2017–046); and A.T. by MCIN/AEI/10.13039/501100011033 (grant CEX2018-000794-S). Statistics South Africa kindly provided the mortality data, but had no other role in the study. Y.G. A.G. M.H. and B. Armstrong set up the collaborative network. Y.G. S.L. and Y.W. designed the study. Y.G. S.L. and A.G. developed the statistical methods. Y.W. B.W. S.L. and Y.G. took the lead in drafting the manuscript and interpreting the results. Y.W. B.W. Y.G. A.G. S.T. A.O. A.U. A.S. A.E. A.M.V.-C. A. Zanobetti, A.A. A. Zeka, A.T. B. Alahmad, B. Armstrong, B.F. C.Í. C. Ameling, C.D.l.C.V. C. Åström, D.H. D.V.D. D.R. E.I. E.L. F.M. F.A. F.D. F.S. G.C.-E. H. Kan, H.O. H. Kim, I.-H.H. J.K. J.M. J.S. K.K. M.H.-D. M.S.R. M.H. M.P. M.d.S.Z.S.C. N.S. P.M. P.G. P.H.N.S. R.A. S.O. T.N.D. V.C. V.H. W.L. X.S. Y.H. M.L.B. and S.L. provided the data and contributed to the interpretation of the results and the submitted version of the manuscript. Y.G. S.L. and Y.W. accessed and verified the data. All of the authors had full access to all of the data in the study and had final responsibility for the decision to submit for publication. The authors declare no competing interests.; Funding text 2: This study was supported by the Australian Research Council ( DP210102076 ) and the Australian National Health and Medical Research Council ( APP2000581 ). Y.W and B.W. were supported by the China Scholarship Council (nos. 202006010044 and 202006010043 ); S.L. was supported by an Emerging Leader Fellowship of the Australian National Health and Medical Research Council (no. APP2009866 ); Y.G. was supported by Career Development Fellowship (no. APP1163693) and Leader Fellowship (no. APP2008813) of the Australian National Health and Medical Research Council ; J.K. and A.U. were supported by the Czech Science Foundation (project no. 20–28560S ); N.S. was supported by the National Institute of Environmental Health Sciences -funded HERCULES Center (no. P30ES019776 ); Y.H. was supported by the Environment Research and Technology Development Fund ( JPMEERF15S11412 ) of the Environmental Restoration and Conservation Agency; M.d.S.Z.S.C. and P.H.N.S. were supported by the São Paulo Research Foundation (FAPESP); H.O. and E.I. were supported by the Estonian Ministry of Education and Research ( IUT34–17 ); J.M. was supported by a fellowship of Fundação para a Ciência e a Tecnlogia ( SFRH/BPD/115112/2016 ); A.G. and F.S. were supported by the Medical Research Council UK (grant ID MR/R013349/1 ), the Natural Environment Research Council UK (grant ID NE/R009384/1 ), and the EU’s Horizon 2020 project, Exhaustion (grant ID 820655 ); A.S. and F.d.D. were supported by the EU’s Horizon 2020 project, Exhaustion (grant ID 820655 ); V.H. was supported by the Spanish Ministry of Economy, Industry and Competitiveness (grant ID PCIN-2017–046 ); and A.T. by MCIN/AEI/10.13039/501100011033 (grant CEX2018-000794-S). Statistics South Africa kindly provided the mortality data, but had no other role in the study

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

Global, regional, and national burden of mortality associated with short-term temperature variability from 2000-19: a three-stage modelling study

BACKGROUND: Increased mortality risk is associated with short-term temperature variability. However, to our knowledge, there has been no comprehensive assessment of the temperature variability-related mortality burden worldwide. In this study, using data from the MCC Collaborative Research Network, we first explored the association between temperature variability and mortality across 43 countries or regions. Then, to provide a more comprehensive picture of the global burden of mortality associated with temperature variability, global gridded temperature data with a resolution of 0.5 degrees x 0.5 degrees were used to assess the temperature variability-related mortality burden at the global, regional, and national levels. Furthermore, temporal trends in temperature variability-related mortality burden were also explored from 2000-19. METHODS: In this modelling study, we applied a three-stage meta-analytical approach to assess the global temperature variability-related mortality burden at a spatial resolution of 0.5 degrees x 0.5 degrees from 2000-19. Temperature variability was calculated as the SD of the average of the same and previous days' minimum and maximum temperatures. We first obtained location-specific temperature variability related-mortality associations based on a daily time series of 750 locations from the Multi-country Multi-city Collaborative Research Network. We subsequently constructed a multivariable meta-regression model with five predictors to estimate grid-specific temperature variability related-mortality associations across the globe. Finally, percentage excess in mortality and excess mortality rate were calculated to quantify the temperature variability-related mortality burden and to further explore its temporal trend over two decades. FINDINGS: An increasing trend in temperature variability was identified at the global level from 2000 to 2019. Globally, 1 753 392 deaths (95% CI 1 159 901-2 357 718) were associated with temperature variability per year, accounting for 3.4% (2.2-4.6) of all deaths. Most of Asia, Australia, and New Zealand were observed to have a higher percentage excess in mortality than the global mean. Globally, the percentage excess in mortality increased by about 4.6% (3.7-5.3) per decade. The largest increase occurred in Australia and New Zealand (7.3%, 95% CI 4.3-10.4), followed by Europe (4.4%, 2.2-5.6) and Africa (3.3, 1.9-4.6). INTERPRETATION: Globally, a substantial mortality burden was associated with temperature variability, showing geographical heterogeneity and a slightly increasing temporal trend. Our findings could assist in raising public awareness and improving the understanding of the health impacts of temperature variability. FUNDING: Australian Research Council, Australian National Health & Medical Research Council

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.This study was supported by the Australian Research Council ( DP210102076 ) and the Australian National Health and Medical Research Council ( APP2000581 ). Y.W and B.W. were supported by the China Scholarship Council (nos. 202006010044 and 202006010043 ); S.L. was supported by an Emerging Leader Fellowship of the Australian National Health and Medical Research Council (no. APP2009866 ); Y.G. was supported by Career Development Fellowship (no. APP1163693) and Leader Fellowship (no. APP2008813) of the Australian National Health and Medical Research Council ; J.K. and A.U. were supported by the Czech Science Foundation (project no. 20–28560S ); N.S. was supported by the National Institute of Environmental Health Sciences -funded HERCULES Center (no. P30ES019776 ); Y.H. was supported by the Environment Research and Technology Development Fund ( JPMEERF15S11412 ) of the Environmental Restoration and Conservation Agency; M.d.S.Z.S.C. and P.H.N.S. were supported by the São Paulo Research Foundation (FAPESP); H.O. and E.I. were supported by the Estonian Ministry of Education and Research ( IUT34–17 ); J.M. was supported by a fellowship of Fundação para a Ciência e a Tecnlogia ( SFRH/BPD/115112/2016 ); A.G. and F.S. were supported by the Medical Research Council UK (grant ID MR/R013349/1 ), the Natural Environment Research Council UK (grant ID NE/R009384/1 ), and the EU’s Horizon 2020 project, Exhaustion (grant ID 820655 ); A.S. and F.d.D. were supported by the EU’s Horizon 2020 project, Exhaustion (grant ID 820655 ); V.H. was supported by the Spanish Ministry of Economy, Industry and Competitiveness (grant ID PCIN-2017–046 ); and A.T. by MCIN/AEI/10.13039/501100011033 (grant CEX2018-000794-S). Statistics South Africa kindly provided the mortality data, but had no other role in the study

Seasonal variation in mortality and the role of temperature: a multi-country multi-city study

Data availability: Data have been collected within the MCC (Multi-Country Multi-City) Collaborative Research Network (https://mccstudy.lshtm.ac.uk) under a data-sharing agreement and cannot be made publicly available. The R code for the analysis is available from the first author.Copyright . Background: Although seasonal variations in mortality have been recognized for millennia, the role of temperature remains unclear. We aimed to assess seasonal variation in mortality and to examine the contribution of temperature. Methods: We compiled daily data on all-cause, cardiovascular and respiratory mortality, temperature and indicators on location-specific characteristics from 719 locations in tropical, dry, temperate and continental climate zones. We fitted time-series regression models to estimate the amplitude of seasonal variation in mortality on a daily basis, defined as the peak-to-trough ratio (PTR) of maximum mortality estimates to minimum mortality estimates at day of year. Meta-analysis was used to summarize location-specific estimates for each climate zone. We estimated the PTR with and without temperature adjustment, with the differences representing the seasonal effect attributable to temperature. We also evaluated the effect of location-specific characteristics on the PTR across locations by using meta-regression models. Results: Seasonality estimates and responses to temperature adjustment varied across locations. The unadjusted PTR for all-cause mortality was 1.05 [95% confidence interval (CI): 1.00–1.11] in the tropical zone and 1.23 (95% CI: 1.20–1.25) in the temperate zone; adjusting for temperature reduced the estimates to 1.02 (95% CI: 0.95–1.09) and 1.10 (95% CI: 1.07–1.12), respectively. Furthermore, the unadjusted PTR was positively associated with average mean temperature. Conclusions: This study suggests that seasonality of mortality is importantly driven by temperature, most evidently in temperate/continental climate zones, and that warmer locations show stronger seasonal variations in mortality, which is related to a stronger effect of temperature.This work was primarily supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI [Grant Number 19K19461]. Y.C. was supported by a Senior Research grant [2019R1A2C1086194] from the National Research Foundation of Korea (NRF), funded by the Ministry of Science, ICT (Information and Communication Technologies). V.H. received support from the Spanish Ministry of Economy, Industry and Competitiveness [Grant ID: PCIN-2017-046]. J.K. and A.U. were supported by the Czech Science Foundation [project 18-22125S]. A.S. acknowledged funding from European Union’s Horizon 2020 research and innovation programme under grant agreement No 820655 (EXHAUSTION). A.G. was supported by the Medical Research Council-UK [Grant ID: MR/R013349/1], the Natural Environment Research Council UK [Grant ID: NE/R009384/1] and the European Union’s Horizon 2020 Project Exhaustion [Grant ID: 820655]. M.H. was supported by the Japan Science and Technology Agency (JST) as part of SICORP [Grant Number JPMJSC20E4]

Seasonality of mortality under climate change: a multicountry projection study

Data sharing: All data used in our study were obtained from the MCC Collaborative Research Network under a data-sharing agreement and cannot be made publicly available. Researchers can refer to collaborators of the Network, who are listed as coauthors of this Article (primary contact: Antonio Gasparrini, [email protected]), for information on accessing the data for each country. The R code is available on request, and a reproducible example is publicly available on the personal GitHub website of the first author (https://github.com/LinaMadaniyazi).For more on the MCC see https://mccstudy.lshtm.ac.uk/Supplementary Material is available online at: https://www.sciencedirect.com/science/article/pii/S2542519623002693#sec1 .Background: Climate change can directly impact temperature-related excess deaths and might subsequently change the seasonal variation in mortality. In this study, we aimed to provide a systematic and comprehensive assessment of potential future changes in the seasonal variation, or seasonality, of mortality across different climate zones. Methods: In this modelling study, we collected daily time series of mean temperature and mortality (all causes or non-external causes only) via the Multi-Country Multi-City Collaborative (MCC) Research Network. These data were collected during overlapping periods, spanning from Jan 1, 1969 to Dec 31, 2020. We projected daily mortality from Jan 1, 2000 to Dec 31, 2099, under four climate change scenarios corresponding to increasing emissions (Shared Socioeconomic Pathways [SSP] scenarios SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5). We compared the seasonality in projected mortality between decades by its shape, timings (the day-of-year) of minimum (trough) and maximum (peak) mortality, and sizes (peak-to-trough ratio and attributable fraction). Attributable fraction was used to measure the burden of seasonality of mortality. The results were summarised by climate zones. Findings: The MCC dataset included 126 809 537 deaths from 707 locations within 43 countries or areas. After excluding the only two polar locations (both high-altitude locations in Peru) from climatic zone assessments, we analysed 126 766 164 deaths in 705 locations aggregated in four climate zones (tropical, arid, temperate, and continental). From the 2000s to the 2090s, our projections showed an increase in mortality during the warm seasons and a decrease in mortality during the cold seasons, albeit with mortality remaining high during the cold seasons, under all four SSP scenarios in the arid, temperate, and continental zones. The magnitude of this changing pattern was more pronounced under the high-emission scenarios (SSP3-7.0 and SSP5-8.5), substantially altering the shape of seasonality of mortality and, under the highest emission scenario (SSP5-8.5), shifting the mortality peak from cold seasons to warm seasons in arid, temperate, and continental zones, and increasing the size of seasonality in all zones except the arid zone by the end of the century. In the 2090s compared with the 2000s, the change in peak-to-trough ratio (relative scale) ranged from 0·96 to 1·11, and the change in attributable fraction ranged from 0·002% to 0·06% under the SSP5-8.5 (highest emission) scenario. Interpretation: A warming climate can substantially change the seasonality of mortality in the future. Our projections suggest that health-care systems should consider preparing for a potentially increased demand during warm seasons and sustained high demand during cold seasons, particularly in regions characterised by arid, temperate, and continental climates.This study was primarily supported by the Environment Research and Technology Development Fund (grant number JPMEERF20231007) of the Environmental Restoration and Conservation Agency, provided by the Ministry of the Environment of Japan. MH was supported by the Japan Science and Technology Agency as part of the Strategic International Collaborative Research Program (grant number JPMJSC20E4). AG was supported by the UK Medical Research Council (grant number MR/V034162/1) and the EU's Horizon 2020 research project Exhaustion (grant number 820655). AU and JK were supported by the Czech Science Foundation (project 22–24920S). JJKJ was supported by the Academy of Finland (grant number 310372; Global Health Risks Related to Atmospheric Composition and Weather Consortium). FS was supported by the Italian Ministry of University and Research, Department of Excellence project 2023–2027, Rethinking Data Science—Department of Statistics, Computer Science and Applications—University of Florence

Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study

© 2021 The Author(s). Background: Exposure to cold or hot temperatures is associated with premature deaths. We aimed to evaluate the global, regional, and national mortality burden associated with non-optimal ambient temperatures. Methods: In this modelling study, we collected time-series data on mortality and ambient temperatures from 750 locations in 43 countries and five meta-predictors at a grid size of 0·5° × 0·5° across the globe. A three-stage analysis strategy was used. First, the temperature–mortality association was fitted for each location by use of a time-series regression. Second, a multivariate meta-regression model was built between location-specific estimates and meta-predictors. Finally, the grid-specific temperature–mortality association between 2000 and 2019 was predicted by use of the fitted meta-regression and the grid-specific meta-predictors. Excess deaths due to non-optimal temperatures, the ratio between annual excess deaths and all deaths of a year (the excess death ratio), and the death rate per 100 000 residents were then calculated for each grid across the world. Grids were divided according to regional groupings of the UN Statistics Division. Findings: Globally, 5 083 173 deaths (95% empirical CI [eCI] 4 087 967–5 965 520) were associated with non-optimal temperatures per year, accounting for 9·43% (95% eCI 7·58–11·07) of all deaths (8·52% [6·19–10·47] were cold-related and 0·91% [0·56–1·36] were heat-related). There were 74 temperature-related excess deaths per 100 000 residents (95% eCI 60–87). The mortality burden varied geographically. Of all excess deaths, 2 617 322 (51·49%) occurred in Asia. Eastern Europe had the highest heat-related excess death rate and Sub-Saharan Africa had the highest cold-related excess death rate. From 2000–03 to 2016–19, the global cold-related excess death ratio changed by −0·51 percentage points (95% eCI −0·61 to −0·42) and the global heat-related excess death ratio increased by 0·21 percentage points (0·13–0·31), leading to a net reduction in the overall ratio. The largest decline in overall excess death ratio occurred in South-eastern Asia, whereas excess death ratio fluctuated in Southern Asia and Europe. Interpretation: Non-optimal temperatures are associated with a substantial mortality burden, which varies spatiotemporally. Our findings will benefit international, national, and local communities in developing preparedness and prevention strategies to reduce weather-related impacts immediately and under climate change scenarios. Funding: Australian Research Council and the Australian National Health and Medical Research Council.Australian Research Council; Australian National Health and Medical Research Council

Global, regional, and national burden of mortality associated with short-term temperature variability from 2000–19: a three-stage modelling study

Background: Increased mortality risk is associated with short-term temperature variability. However, to our knowledge, there has been no comprehensive assessment of the temperature variability-related mortality burden worldwide. In this study, using data from the MCC Collaborative Research Network, we first explored the association between temperature variability and mortality across 43 countries or regions. Then, to provide a more comprehensive picture of the global burden of mortality associated with temperature variability, global gridded temperature data with a resolution of 0·5° × 0·5° were used to assess the temperature variability-related mortality burden at the global, regional, and national levels. Furthermore, temporal trends in temperature variability-related mortality burden were also explored from 2000–19. Methods: In this modelling study, we applied a three-stage meta-analytical approach to assess the global temperature variability-related mortality burden at a spatial resolution of 0·5° × 0·5° from 2000–19. Temperature variability was calculated as the SD of the average of the same and previous days’ minimum and maximum temperatures. We first obtained location-specific temperature variability related-mortality associations based on a daily time series of 750 locations from the Multi-country Multi-city Collaborative Research Network. We subsequently constructed a multivariable meta-regression model with five predictors to estimate grid-specific temperature variability related-mortality associations across the globe. Finally, percentage excess in mortality and excess mortality rate were calculated to quantify the temperature variability-related mortality burden and to further explore its temporal trend over two decades. Findings: An increasing trend in temperature variability was identified at the global level from 2000 to 2019. Globally, 1 753 392 deaths (95% CI 1 159 901–2 357 718) were associated with temperature variability per year, accounting for 3·4% (2·2–4·6) of all deaths. Most of Asia, Australia, and New Zealand were observed to have a higher percentage excess in mortality than the global mean. Globally, the percentage excess in mortality increased by about 4·6% (3·7–5·3) per decade. The largest increase occurred in Australia and New Zealand (7·3%, 95% CI 4·3–10·4), followed by Europe (4·4%, 2·2–5·6) and Africa (3·3, 1·9–4·6). Interpretation: Globally, a substantial mortality burden was associated with temperature variability, showing geographical heterogeneity and a slightly increasing temporal trend. Our findings could assist in raising public awareness and improving the understanding of the health impacts of temperature variability. Funding: Australian Research Council, Australian National Health & Medical Research Council

Excess mortality attributed to heat and cold: a health impact assessment study in 854 cities in Europe

Data sharing: The exposure-response functions derived in this analysis, full results, and intermediary data are publicly available in a Zenodo repository (https://doi.org/10.5281/zenodo.7672108). The associated R code to reproduce the analysis is available in the corresponding author's GitHub page (https://github.com/pierremasselot). The mortality data have been obtained through a restricted data use agreement with each national institute and are therefore not available for public dissemination.Copyright © 2023 The Authors. Background: Heat and cold are established environmental risk factors for human health. However, mapping the related health burden is a difficult task due to the complexity of the associations and the differences in vulnerability and demographic distributions. In this study, we did a comprehensive mortality impact assessment due to heat and cold in European urban areas, considering geographical differences and age-specific risks. Methods: We included urban areas across Europe between Jan 1, 2000, and Dec 12, 2019, using the Urban Audit dataset of Eurostat and adults aged 20 years and older living in these areas. Data were extracted from Eurostat, the Multi-country Multi-city Collaborative Research Network, Moderate Resolution Imaging Spectroradiometer, and Copernicus. We applied a three-stage method to estimate risks of temperature continuously across the age and space dimensions, identifying patterns of vulnerability on the basis of city-specific characteristics and demographic structures. These risks were used to derive minimum mortality temperatures and related percentiles and raw and standardised excess mortality rates for heat and cold aggregated at various geographical levels. Findings: Across the 854 urban areas in Europe, we estimated an annual excess of 203 620 (empirical 95% CI 180 882–224 613) deaths attributed to cold and 20 173 (17 261–22 934) attributed to heat. These corresponded to age-standardised rates of 129 (empirical 95% CI 114–142) and 13 (11–14) deaths per 100 000 person-years. Results differed across Europe and age groups, with the highest effects in eastern European cities for both cold and heat. Interpretation: Maps of mortality risks and excess deaths indicate geographical differences, such as a north–south gradient and increased vulnerability in eastern Europe, as well as local variations due to urban characteristics. The modelling framework and results are crucial for the design of national and local health and climate policies and for projecting the effects of cold and heat under future climatic and socioeconomic scenarios.Medical Research Council of UK, the Natural Environment Research Council UK, the EU's Horizon 2020, and the EU's Joint Research Center. The study was funded by Medical Research Council of the UK (MR/V034162/1 and MR/R013349/1), the Natural Environment Research Council UK (NE/R009384/1), the EU's Horizon 2020 (820655), and the EU's Joint Research Center (JRC/SVQ/2020/MVP/1654). AU and JK were supported by the Czech Science Foundation (22–24920S). VH has received funding from the EU's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement (101032087