Interactive effects of ambient fine particulate matter and ozone on daily mortality in 372 cities: two stage time series analysis

Objective To investigate potential interactive effects of fine particulate matter (PM2.5) and ozone (O3) on daily mortality at global level. Design Two stage time series analysis. Setting 372 cities across 19 countries and regions. Population Daily counts of deaths from all causes, cardiovascular disease, and respiratory disease. Main outcome measure Daily mortality data during 1994-2020. Stratified analyses by co-pollutant exposures and synergy index (>1 denotes the combined effect of pollutants is greater than individual effects) were applied to explore the interaction between PM2.5 and O3 in association with mortality. Results During the study period across the 372 cities, 19.3 million deaths were attributable to all causes, 5.3 million to cardiovascular disease, and 1.9 million to respiratory disease. The risk of total mortality for a 10 μg/m3 increment in PM2.5 (lag 0-1 days) ranged from 0.47% (95% confidence interval 0.26% to 0.67%) to 1.25% (1.02% to 1.48%) from the lowest to highest fourths of O3 concentration; and for a 10 μg/m3 increase in O3 ranged from 0.04% (−0.09% to 0.16%) to 0.29% (0.18% to 0.39%) from the lowest to highest fourths of PM2.5 concentration, with significant differences between strata (P for interaction <0.001). A significant synergistic interaction was also identified between PM2.5 and O3 for total mortality, with a synergy index of 1.93 (95% confidence interval 1.47 to 3.34). Subgroup analyses showed that interactions between PM2.5 and O3 on all three mortality endpoints were more prominent in high latitude regions and during cold seasons. Conclusion The findings of this study suggest a synergistic effect of PM2.5 and O3 on total, cardiovascular, and respiratory mortality, indicating the benefit of coordinated control strategies for both pollutants

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.Massimo Stafoggia, Francesca K. de’ Donato, Masna Rai and Alexandra Schneider were partially supported by the European Union’s Horizon 2020 Project Exhaustion (Grant ID: 820655). Jan Kyselý and Aleš Urban were supported by the Czech Science Foundation project (22-24920S). Joana Madureira was supported by the Fundação para a Ciência e a Tecnologia (FCT) (grant SFRH/BPD/115112/2016). Masahiro Hashizume was supported by the Japan Science and Technology Agency (JST) as part of SICORP, Grant Number JPMJSC20E4. Noah Scovronick was supported by the NIEHS-funded HERCULES Center (P30ES019776). South African Data were provided by Statistics South Africa, which did not have any role in conducting the study. Antonio Gasparrini was supported by the Medical Research Council-UK (Grants ID: MR/V034162/1 and 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)

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

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

The burden of heat-related mortality attributable to recent human-induced climate change

Medical Research Council-UK (Grant ID: MR/M022625/1); Natural Environment Research Council UK (Grant ID: NE/R009384/1); European Union’s Horizon 2020 Project Exhaustion (Grant ID: 820655); N. Scovronick was supported by the NIEHS-funded HERCULES Center (P30ES019776); Y. Honda was supported by the Environment Research and Technology Development Fund of the Environmental Restoration and Conservation Agency, Japan (JPMEERF15S11412); J. Jaakkola was supported by Academy of Finland (Grant No. 310372); V. Huber was supported by the Spanish Ministry of Economy, Industry and Competitiveness (Grant ID: PCIN-2017-046) and the German Federal Ministry of Education and Research (Grant ID: 01LS1201A2); J Kysely and A. Urban were supported by the Czech Science Foundation (Grant ID: 20-28560S); J. Madureira was supported by the Fundação para a Ciência e a Tecnologia (FCT) (SFRH/BPD/115112/2016); S. Rao and F. di Ruscio were supported by European Union’s Horizon 2020 Project EXHAUSTION (Grant ID: 820655); M. Hashizume was supported by the Japan Science and Technology Agency (JST) as part of SICORP, Grant Number JPMJSC20E4; Y. Guo was supported by the Career Development Fellowship of the Australian National Health and Medical Research Council (#APP1163693); S. Lee was support by the Early Career Fellowship of the Australian National Health and Medical Research Council (#APP1109193)

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

Impact of population aging on future temperature-related mortality at different global warming levels

Data availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Data were collected within the MCC Collaborative Research Network under a data sharing agreement and cannot be made publicly available.Code availability: A sample of the analysis code is available from https://github.com/CHENlab-Yale/MCC_ProjAging_Temp .Supplementary information is available online at: https://link-springer-com.ezproxytest.brunel.ac.uk/article/10.1038/s41467-024-45901-z#Sec15 .Older adults are generally amongst the most vulnerable to heat and cold. While temperature-related health impacts are projected to increase with global warming, the influence of population aging on these trends remains unclear. Here we show that at 1.5 °C, 2 °C, and 3 °C of global warming, heat-related mortality in 800 locations across 50 countries/areas will increase by 0.5%, 1.0%, and 2.5%, respectively; among which 1 in 5 to 1 in 4 heat-related deaths can be attributed to population aging. Despite a projected decrease in cold-related mortality due to progressive warming alone, population aging will mostly counteract this trend, leading to a net increase in cold-related mortality by 0.1%–0.4% at 1.5–3 °C global warming. Our findings indicate that population aging constitutes a crucial driver for future heat- and cold-related deaths, with increasing mortality burden for both heat and cold due to the aging population.We acknowledge the World Climate Research Programme, which, through its Working Group on Coupled Modeling, coordinated and promoted CMIP6. We thank the climate modeling groups for producing and making available their model output, the Earth System Grid Federation (ESGF) for archiving the data and providing access, and the multiple funding agencies who support CMIP6 and ESGF. K.C. was supported by the Yale Planetary Solutions Project seed grant. A.G., A.S., and S.R. were supported by the European Union’s Horizon 2020 Project Exhaustion grant (820655). A.G. was also supported by the Medical Research Council UK grant (MR/V034162/1). J.M. received funding from the Fundação para a Ciência e a Tecnlogia Grant (SFRH/BPD/115112/2016). A.T. was supported by the MCIN/AEI/10.13039/501100011033 grant (CEX2018-000794-S). A.U. and J.K. were supported by the Czech Science Foundation (22-24920S). F.S. was supported by the Italian Ministry of University and Research (MUR), Department of Excellence project 2023-2027 ReDS ‘Rethinking Data Science’ - Department of Statistics, Computer Science and Applications - University of Florence. MNM. was supported by the European Commission (H2020-MSCA-IF-2020) under REA grant agreement no. 101022870. A.V.C. acknowledges the support of the Swiss National Foundation (TMSGI3_211626). V.H. received funding from the European Union’s Horizon 2020 research and innovation program (Marie Skłodowska-Curie Grant Agreement No.: 101032087)

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

Geographical Variations of the Minimum Mortality Temperature at a Global Scale

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

Mortality risk attributable to wildfire-related PM<inf>2·5</inf> pollution: a global time series study in 749 locations

Data sharing: Data used in this study were collected by collaborators within the MCC Network under a data sharing agreement and cannot be made available publiclyCopyright © 2021 The Author(s). Background: Many regions of the world are now facing more frequent and unprecedentedly large wildfires. However, the association between wildfire-related PM2·5 and mortality has not been well characterised. We aimed to comprehensively assess the association between short-term exposure to wildfire-related PM2·5 and mortality across various regions of the world. Methods: For this time series study, data on daily counts of deaths for all causes, cardiovascular causes, and respiratory causes were collected from 749 cities in 43 countries and regions during 2000–16. Daily concentrations of wildfire-related PM2·5 were estimated using the three-dimensional chemical transport model GEOS-Chem at a 0·25° × 0·25° resolution. The association between wildfire-related PM2·5 exposure and mortality was examined using a quasi-Poisson time series model in each city considering both the current-day and lag effects, and the effect estimates were then pooled using a random-effects meta-analysis. Based on these pooled effect estimates, the population attributable fraction and relative risk (RR) of annual mortality due to acute wildfire-related PM2·5 exposure was calculated. Findings: 65·6 million all-cause deaths, 15·1 million cardiovascular deaths, and 6·8 million respiratory deaths were included in our analyses. The pooled RRs of mortality associated with each 10 μg/m3 increase in the 3-day moving average (lag 0–2 days) of wildfire-related PM2·5 exposure were 1·019 (95% CI 1·016–1·022) for all-cause mortality, 1·017 (1·012–1·021) for cardiovascular mortality, and 1·019 (1·013–1·025) for respiratory mortality. Overall, 0·62% (95% CI 0·48–0·75) of all-cause deaths, 0·55% (0·43–0·67) of cardiovascular deaths, and 0·64% (0·50–0·78) of respiratory deaths were annually attributable to the acute impacts of wildfire-related PM2·5 exposure during the study period. Interpretation: Short-term exposure to wildfire-related PM2·5 was associated with increased risk of mortality. Urgent action is needed to reduce health risks from the increasing wildfires.Australian Research Council, Australian National Health & Medical Research Council