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

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.

considering the whole range of observed values) account for more than 5 million deaths globally every year (2000-2019), with the burden attributable to cold 9 times greater than that attributable to heat 7 .Current projections indicate that heat-related deaths will steeply increase under warmer climates, while cold-related deaths will decline in most locations given the reduction in cold days 8,9 .However, these studies do not account for an aging population.Current estimates of future mortality burden mostly account only for changes in climate, while there is limited quantitative evidence on future health impacts of climate change under more complex and potentially realistic scenarios.In this respect, understanding the influence of other global changes in addition to the changing climate, on the change in temperature (heat and cold)-related mortality at different levels of warming vs. current times remains a critical issue in adequately quantifying the overall health impact of climate change.
While the climate is changing, the world is facing another public health challenge: aging.The proportion of world's population aged 65 years and above is projected to rise from 9% at present to 16% in 2050 10 .Older adults are considered among the most vulnerable populations to non-optimal temperatures 11 , due to factors such as more limited thermoregulatory responses, relatively high prevalence of chronic conditions, and a higher likelihood of social isolation 12,13 .Understanding the influence of different population demographic scenarios on the estimated temperature-related health impacts can provide important insights on plausible future health burdens.Given the increased risk present amongst older adults, it is expected that population aging will substantially amplify the future temperaturerelated impacts by increasing the vulnerability of the populations 14 .A recent review of local and regional temperature-related health impacts found that accounting for population aging could lead to a lower reduction or even an increase in the cold-related mortality burden, resulting in substantial increases in temperature-related mortality 14 .A recent nation-wide assessment in Switzerland found an increasing trend in the projected cold-related mortality under warmer climates due to an increase in the size of the population at risk when demographic changes were accounted for 1 .However, evidence on the global scale using a systematic and standardized analytical approach is limited.No study before has provided explicit estimates of the contribution of aging in future climate warming scenarios.
In this work, we assess the impact of population aging on future temperature-related excess mortality at different levels of global warming (1.5 °C: 2018-2037, 2 °C: 2032-2051, and 3 °C: 2055-2074) in 800 locations across 50 countries/areas (Fig. 1).In brief, we first estimate the location-age-specific temperature-mortality associations in a two-stage time-series analysis using quasi-Poisson regression with distributed lag nonlinear models and multivariate random metaregression using data from the Multi-Country Multi-City (MCC) Collaborative Research Network (http://mccstudy.lshtm.ac.uk/) (description provided in Supplementary Tables 1 and 2).We then combine the association estimates with bias-corrected future temperature series from 18 general circulation models (GCMs) (Coupled Model Intercomparison Project 6 -CMIP6) 15 , and future trends in age-specific baseline mortality 2 to derive excess temperature-mortality projections under three warming scenarios (1.5 °C, 2 °C, and 3 °C of global warming above preindustrial level) and for the historical period (1995-2014).Finally, we quantify the impact of population aging as the difference in the change in temperature-related mortality fractions (warming window minus historical period) between climate-only (i.e., not accounting for changes in population demographics) and climate-population scenarios.Mortality fractions refer to the percentage of all-cause deaths attributed to heat and cold over total mortality in each time window.

Projected temperature changes at different global warming levels
We used the temperature projections under the very high greenhouse gas (GHG) emission scenariothe Shared Socio-economic Pathway 5-8.5 (SSP5-8.5)and identified temperature windows consistent with such rises in global mean surface temperature 16 .To explore the influence of the timing of the temperature windows on population aging, we also conducted a sensitivity analysis using the emission scenario SSP3-7.0 (i.e., 1.5 °C: 2021-2040, 2 °C: 2037-2056, and 3 °C: 2066-2085), under which a 3 °C warming will be also crossed but under different trajectories of warming and population structure changes 16 .The 'regional rivalry' SSP3 scenario is a socio-political scenario envisioning a resurgent nationalism, concerns about competitiveness and security, and regional conflicts that result in high challenges for both mitigation and adaptation 17 .Supplementary Fig. 1A shows the multi-model ensemble-average and range of global mean surface temperature changes relative to 1995-2014 (defined as present-day in the IPCC Sixth Assessment Report) under SSP5-8.5 and SSP3-7.0scenarios from 18 GCMs (Supplementary Table 3).Trends in global mean surface temperature changes are similar under both emission scenarios before 2040 and begin to slowly diverge after 2040.By the end of the 21 st century (2081-2100), the ensemble average of 18 GCMs suggests global surface temperature changes increase by 4.6 °C and 3.7 °C under SSP5-8.5 and SSP3-7.0,respectively, relative to our historical baseline (1995-2014).We further bias-corrected the daily temperature projections at each MCC location, using the historical temperature observations and a method described previously 3,4 .Figure 1 illustrates the average temperature changes in 800 locations at 1.5 °C (2018-2037), 2 °C (2032-2051), and 3 °C (2055-2074) warming levels under the SSP5-8.5.For these time windows, temperatures in the 800 studied locations increase by 0.8 °C (range of 0.3 °C to 1.3 °C), 1.5 °C (0.7 °C to 2.2 °C), and 2.7 °C (1.3 °C to 4.2 °C) compared to the historical period (1995-2014) (Supplementary Fig. 1B).With increasing levels of global warming, larger variability can be observed in temperature changes across the locations.The largest temperature increases mostly occur in locations in Northern America (e.g., Winnipeg, Manitoba, Canada; and Fargo, North Dakota, USA) and Northern Europe (e.g., Narva linn, Estonia), with an average historical temperature below 13.8 °C.Almost identical temperature increases occur at the studied locations under SSP3-7.0,with an average of 0.8 °C (0.4 °C to 1.3 °C), 1.5 °C (0.7 °C to 2.2 °C), and 2.7 °C (1.3 °C to 4.3 °C) increase compared with the historical period at the 1.5 °C (2021-2040), 2 °C (2037-2056), and 3 °C (2066-2085) warming levels, respectively.

Population aging at different global warming levels
We estimated population aging as the change in the percentage of the population of 65 years and above under SSP5 5 .in each country during the 20-year period reaching 1.5 °C (2018-2037), 2 °C (2032-2051), and 3 °C (2055-2074) global warming relative to the historical period (Fig. 2).Across the 50 countries/areas, the average percentage of the population 65 years and above is projected to increase by 3.0%, 7.3%, and 13.8% when global warming reaches 1.5 °C, 2 °C, and 3 °C, respectively.In our study countries/areas, at 3 °C of warming (2055-2074), countries/areas in Southern Asia experience the largest increase in population aging (23.6%), followed by the countries/areas in Eastern Europe (20.1%),Latin America and the Caribbean (17.8%),Eastern Asia (17.50%),South-eastern Asia (17.2%); whereas countries/areas in North Europe (6.0%) have the lowest rate of population aging, followed by Northern America (6.7%),Western Europe (7.6%), and Australia (7.9%).In the sensitivity analysis under SSP3, a similar increasing trend in the average percentage of the population 65 years and above is expected at 1.5 °C (3.6%), 2 °C (6.4%), and 3 °C (9.1%) of warming levels, respectively (Supplementary Fig. 2).Estimates correspond to the change in temperature-related mortality fractions (heat and cold) at different levels of global warming under SSP5-8.5 and SSP3-7.0scenarios compared with the historical period (1995-2014).Heat is defined as all temperatures above the minimum mortality temperature (MMT) and cold is defined as all temperatures below the MMT (Supplementary Fig. 3).The future periods in which the 20-year running mean of global mean temperature first reaches the 1.

Changes in temperature-related mortality burden at different warming levels
We derived heat-and cold-related mortality estimates assuming: 1) constant demographic structure and size of the population consistent with the year 2015 (climate-only), and 2) time-varying population demographic structure and size according to SSP5 (fossil-fueled development with high challenge for mitigation and low challenge for adaptation) population scenario (and SSP3 population scenario in the sensitivity analysis) 18 .We derived age group-specific baseline mortality projections during each temperature window according to the selected SSP based on projected population size and mortality rates, and used the age-specific association estimates (description provided in Methods and Supplementary Text) to derive the corresponding heatand cold-related mortality.We assume no changes in population adaptation across time (i.e., constant exposure-response association).In the climate-population scenario, we calculated the heat-and cold-related mortality fractions, and total (i.e.non-optimal temperatures), and their changes in the future period relative to the historical period by applying the location and age group-specific associations (Supplementary Fig. 3) to the corresponding daily mean temperature and baseline mortality series in each 20-year period.For the "climateonly" scenario, the same methodology was used but assuming a constant baseline mortality corresponding to the average mortality in each day of the year of the historical period.Because using SSP5-8.5 and SSP3-7.0scenarios yielded similar or even identical results on the temperature-related mortality burden estimates (Table 1 and Supplementary Fig. 4), we focus on presenting the main results from the SSP5-8.5 scenario hereafter.
Under the climate-only scenario, heat-related mortality fractions increased but cold-related mortality fractions decreased overall, resulting in a small decrease in mortality burden due to non-optimal temperatures at 1.5 °C warming (−0.2% [empirical confidence interval (eCI) −0.7% to 0.1%)] and 2 °C [−0.3% (95% eCI: −1.1% to 0.5%)] (Table 1; country-specific results are shown in Figs. 3 and 4, and Supplementary Table 4), relative to the historical period.However, under 3 °C of global warming, an overall small increase in temperature-related mortality burden [0.1% (95% eCI: −2.4% to 2.3%)] can be observed due to larger increases in heat-related mortality than the decreases in cold-related mortality.
Impact of population aging on future temperature-related mortality , a constant population scenario).At 1.5 °C warming, all 50 countries/areas will experience a reduction (−1.4% in Vietnam to −0.3% in Iceland) in cold-related mortality burden driven by climate change alone (Fig. 5A).However, population aging will mostly offset this decrease and result in an increased cold-related mortality burden in 18 out of the 50 countries/areas (0.0% in Argentina to 1.9% in Kuwait).At higher levels of warming, the decline in cold-related mortality burden due to climate change will be further offset by a larger increase due to population aging, leading to 23 and 21 countries/areas having an increase in cold-related mortality burden under 2 and 3 °C warming, respectively.
For heat-related mortality (Fig. 5B), climate change alone will result in increased mortality burden in all countries/areas at 1.5 °C warming (0.02% in Iceland to 2.8% in French Guiana), 49 countries/ areas (0.03% in South Africa to 6.9% in French Guiana) except for Colombia (−0.3%) at 2 °C warming, and 48 countries/areas (0.1% in Iceland to 19.3% in French Guiana) except for Colombia [−3.4% (95% eCI: −38.6% to 10.5%)] and South Africa [−0.2% (95% eCI: −3.9% to 2.0%)] at 3 °C warming.The decreases seen in Colombia and South Africa are due to the decreased but non-significant heat-related mortality risks (Supplementary Fig. 3).Compared with climate change, population aging will have a smaller contribution to the changes in heat-related mortality, leading to an increase of 0.1-0.7% at 1.5-3 °C warming across 50 countries/areas.This translates into 1 in 5 to 1 in 4 deaths in the 0.5-2.5% increase in heat-related deaths attributable to aging (Table 1).For the non-optimal (cold and heat) temperature-related mortality, population aging alone results in an average of 0.8% (95% eCI: 0.6% to 0.9%), 1.7% (1.2% to 2.1%), and 2.6% (0.9% to 3.5%) increases in future temperature-related excess mortality when global warming reaches 1.5 °C, 2 °C, and 3 °C (Fig. 6 and Supplementary Table 7).In comparison, climate alone accounted for 3.0% (i.e., 0.08%/2.67%) of the net changes in non-optimal temperature-related mortality under the highest level (3 °C) of warming.Significantly positive estimations on mortality burden due to population aging were observed in about two-thirds (36 under all warming levels) of the countries/areas, whereas significantly negative contributions due to population aging were observed in very few countries/areas (5, 3, and 2 at 1.5 °C, 2 °C, and 3 °C, respectively).The larger role of population aging in future non-optimal temperature-related mortality may outweigh the differences in temperature rises observed between the SSP5-8.5 and SSP3-7.0scenarios under the same global warming level, culminating in similar results across these scenarios.

Discussion
In this study, we projected the future temperature-related mortality burden and quantified the contribution of population aging in 50 countries/areas at different global warming levels (1.5 °C, 2 °C, and 3 °C).Our analyses revealed an increasing temperature-related mortality burden at all warming levels studied when accounting for changes in population demographics.Population aging substantially amplifies future heat-and cold-related mortality burdens.Compared to climate change, population aging plays a dominant role in future cold-related mortality burden and would lead to a net increase in future cold-related mortality for all study locations combined.
Most of the temperature-related mortality burden can be attributable to cold rather than heat exposure, as shown in previous studies 6,7 .Thus, our findings strongly underline and support the need to account for a significant shift in the number of older individuals who will die from either cold or heat globally, regardless of whether we see large or small changes in climate.Without acknowledging the shifting population and the increasing number of people exposed to non-optimal temperatures (both heat and cold), the ability to address the health impacts of temperature extremes will be hindered.As observed in many previous studies [8][9][10] , in some locations, population exposed to warmer climates will lead to a reduction in cold-related mortality burden that could outnumber the expected increase in heat-related mortality burden.However, our study shows that accounting for population aging leads to an attenuation of this reduction or even an increase in cold-related mortality burden, resulting in a substantial increase in temperature-related mortality.Consistent with our analysis, the few previous studies considering population aging projected net increases in future temperaturerelated mortality burden under the Special Report on Emissions Scenarios (SRES) 11,12 and Representative Concentration Pathway (RCP) 13,14 , global warming scenarios that were used in the previous IPCC Fourth and Fifth Assessment Reports.
Our findings demonstrate that accounting for population aging reveals a globally consistent increase in overall temperature-related mortality burden for any increase in global warming.Due to the dominant role of population aging, the net increases in temperaturerelated mortality were evident in 45 countries/areas at the 2 °C global warming and in 40 countries/areas at the 1.5 °C global warming (Fig. 5 and Supplementary Table 6).On average, increases in heat-related mortality fractions are larger than those in cold-related mortality fractions at the 1.5 °C and 2 °C warming.The comparison between the impact of 1.5 °C and 2 °C on overall temperature-related health under climate change alone could be highly uncertain and vary geographically 15 .This poorly understood difference between the health impact of 1.5 °C and 2 °C was recognized in the IPCC Special Report on Global Warming of 1.5 °C as a key knowledge gap in climate change and health research 16 .After taking population aging into account, the projected temperature-related mortality burden in most (47/50) studied countries/areas is consistently larger at 2 °C than at 1.5 °C of global warming (Fig. 5 and Supplementary Table 6).Our study thus helps fill this knowledge gap in understanding the difference between the projected health risks at 1.5 °C and 2 °C.
To the best of our knowledge, this is the first global study to specifically account for population aging in projection of future temperature-related mortality at global warming levels below and above the goals set by the Paris Agreement.This is also one of the few studies utilizing a larger ensemble of temperature projections from the latest generation of GCMs participating in the most recent climate modeling protocol (CMIP6).Our study is also the first study to quantify the contribution of population aging through both age group-specific temperature-mortality associations and projected baseline age group-specific mortality under the SSP emission and population scenarios.Our estimates are based on an extensive dataset with more than 83 million deaths from 50 countries/areas across five continents, enabling us to assess the vulnerability of each studied population through localized temperature-mortality associations.Other strengths include a state-of-the-art epidemiological analytical method, statistical downscaling techniques, population aging projections, and an advanced approach to account for uncertainty in both temperature-mortality relationships and variability across GCMs 9 .
Our study has several limitations.First, we did not consider potential population adaptation to heat 17 given difficulties in defining and quantifying valid adaptation assumptions and mechanisms related to population-level adaptation to heat 18 .However, not accounting for adaptation allowed us to better disentangle the contribution of population aging and changes in the climate.Another limitation could be the unavailability of sub-national projections for population aging in all of our studied countries/areas.Moreover, our estimates may not necessarily represent country-specific average effects since most locations in the MCC database are urban areas and, in some cases, only one or two cities are available for a given country.Finally, due to the lack of or insufficient coverage of data, we were unable to obtain estimates of projected temperature-related mortality burden in most parts of South Asia and Africa, which will have the largest contribution of population change during the 21 st century 19,20 .Thus, our results of the population aging impact can be interpreted as conservative as we might expect larger impacts of population aging in those regions.Future international collaborations with researchers from these regions will help extend the MCC Network and address the existing data gaps.
As a major risk factor for many chronic diseases and conditions, population aging poses another threat to public health in addition to the looming climate crisis.The current global efforts to address climate change are by far insufficient to limit global warming to 1.5 °C.As of September 2021, the new and updated unconditional nationally determined contributions commitments and pledge announcements are estimated to limit global warming to 2.7 °C, inadequate to achieve the Paris Agreement's goals 21 .Even full implementation of the newly announced national net-zero emission targets could only lower global warming to 2.0-2.4 °C by 2100 22 , leaving the 1.5 °C target beyond reach.We found that population aging will play a dominant role in determining the future temperature-related mortality burden, the higher the levels (1.5 °C, 2 °C, or 3 °C) of global warming, the greater the aging-induced temperature-related mortality burden.Our findings underscore the need for ambitious and drastic climate mitigation actions to keep 1.5 °C warming within reach and targeted and efficient climate adaptation measures to prevent temperature-related health impacts under the dual threats from climate change and population aging.

Daily meteorological and mortality data
Through the MCC Collaborative Research Network, we collected observed daily meteorological and mortality data from local weather stations and health authorities in each country, respectively (Supplementary Table 1).This dataset has been applied in previous publications on the association between temperature and daily mortality and future projections of temperature-related mortality under different Representative Concentration Pathways 9,15 .Daily counts of deaths due to non-external causes (International Classification of Diseases, 9th Revision codes 0-799 or 10th Revision codes A00-R99) or all-cause deaths were used as mortality data.Although we focus here on the non-external or all-cause deaths, previous studies have shown that the majority of mortality burden attributable to heat and cold are from cardiovascular and respiratory deaths (e.g., 68.3% in Jiangsu, China 23 ).Supplementary Table 1 provides a detailed description of the data collection.

Population aging and baseline mortality projections
We obtained country-level projections of population size and mortality rates for three age groups (0-64, 65-74, and ≥75 years) under the SSP5 (Fossil-fueled development) and the SSP3 (Regional rivalry) scenarios 26 from the SSP Database -Version 2.0 (https://secure.iiasa.ac.at/web-apps/ene/SspDb/).This data was provided in 5-year bands from the year 2015 to the year 2100.For each SSP, we combined the countryage-specific mortality rate and population in each 5-year band and interpolated it into annual mortality by fitting a smoothing function (i.e., natural spline with 1 df every 10 years).We then transformed the country-age-specific mortality to the location-age-specific mortality by applying two factors: 1) an "age factor" using the projected mortality in 2015, since the MCC dataset does not include age-specific deaths, which corresponds to the share of deaths in each age category (0-64, 65-74, and ≥75 years) over the total, 2) "country-to-location" factor corresponding to the share of total deaths in each location vs. the country total.This second factor was estimated using the average annual mortality in each location using the available data between 1994-2015 and the country-specific mortality from 2015.In this way, we considered the mortality in 2015 to be held constant across the historical period between 1995-2014.In the last step, we retrieve the daily mortality from the annual series by applying a third factor "seasonal factor" corresponding to the average share of mortality by day of the year over the total annual mortality in the observed series in each location.By applying this factor, we preserve the seasonal cycle of mortality with a higher number of deaths during the winter season compared to the summer season.To disentangle the impact of population aging from climate change, we also computed the future baseline mortality without considering any population change by replicating the daily series of mortality counts in 1995-2014 along the century until 2100.

Temperature-mortality associations
We used a two-stage time-series analysis, the state-of-the-art methodology in multi-location assessments, to estimate the overall association between daily temperature and mortality in each location, i.e., the exposure-response function (ERF).In the first stage, we applied a generalized linear model with quasi-Poisson family and distributed lag non-linear models (DLNM) to estimate location-specific temperaturemortality associations using observed data included in the MCC database.Consistent with previous studies 6,17 , the time-series model included a natural cubic spline with eight degrees of freedom (df) per year to control the long-term and seasonal time trend, and indicator variables for the day of the week.We modeled the complex non-linear and delayed relationship between temperature and mortality following the DLNM methodology 9,27 .Specifically, the time-series model equation is as follows: where Y i,t is the daily counts of deaths on day t in location i; α is the intercept; cb() is the cross-basis function for temperature with a natural cubic spline with three internal knots placed at the 10 th , 75 th , and 90 th percentiles of the temperature distribution, and lag-response relationship with a natural cubic spline with three internal knots at equallyspaced log-scale values over 21 days of lag 9 ; ns() is the natural cubic spline for date t with eight df per year (a total of Nyears years) to control for both long-term trends and seasonality; and DOW i,t represents the day of the week.
In the second stage, we pooled city-specific estimates using a novel two-level hierarchical random-effects meta-regression model 28 .The model accounted for location-specific fixed-effects metapredictors including average temperature, temperature interquartile range, and country-level gross domestic product, and applied twolevel random effects with cities nested within country-specific climate zones 29 to account for heterogeneity across both.We then estimated the best linear unbiased predictions (BLUPs) for each location, which represents a trade-off between location-specific association from the first stage and the pooled association from the meta-regression 30 .Based on BLUPs, we calculated the minimum mortality temperature (MMT) for each location, which is the temperature corresponding to the lowest risk of mortality.We restricted MMT to be within the 2 nd -98 th percentile range to avoid imprecisely estimated tails of the ERF 31 .We then retrieve age-specific relationships using the output of a doseresponse multivariate meta-analytical model derived in a recent work of the consortium.The factors obtained from this meta-analytical model allowed us to predict the temperature-mortality associations for the 0-64, 65-74, and ≥75 years age groups, respectively (see Supplementary Text).

Health impact assessment
We estimated the number of deaths and the corresponding heatrelated and cold-related excess mortality fraction in each location for the historical and future 20-year periods using the standard methodology for quantification of health impacts 32 .Briefly, we applied the location-age-specific ERFs and the modeled daily series of temperature projections under SSP5-8.5 or SSP3-7.0 and baseline daily mortality projections to calculate the daily number of temperature-related excess deaths.For each 20-year period (historical and future decades corresponding to each warming level), we computed the total number of heat-related and cold-related excess deaths (EN) by summing the contributions from all the days of the series when the mean daily temperature was above or below the location-specific MMT, respectively.We then calculated the corresponding excess mortality fraction (EF -%) by dividing the estimated EN by the total number of deaths in https://doi.org/10.1038/s41467-024-45901-zNature Communications | (2024) 15:1796 each location.We also obtained mortality estimates aggregated by country-and overall across the 800 locations.
We then computed the future changes in EF as the difference between the EF estimated in each future period (corresponding to global warming targets (1.5 °C, 2 °C, and 3 °C)) and the historical period for each GCM and location.Additionally, we calculated EF in each future 20-year period assuming no changes in population growth and demographic structure (i.e.assuming the daily baseline mortality in the historical period) as done in Gasparrini et al. 9 .We finally estimated the impact of population aging by subtracting the future changes in temperature-related impacts in the constant population scenario (climate-only) from the change in EF derived under SSP5 or SSP3 mortality projections in the climate-population scenario.In other words, the climate-only scenario holds all populations constant at the baseline year 2015's size, demographic structure, and age-specific mortality rate, whereas the climate-population scenario assumes changes in population size, demographic structure, and age-specific mortality rate in the future periods under SSP5 or SSP3.We estimated the empirical confidence intervals (eCI) using a set of 1000 simulated coefficients defining the ERF derived from Monte Carlo simulations.For each GCM scenario, we estimated the corresponding distribution of impacts and then derived the lower and upper bound of the ECI as the 2.5th and 97.5th percentiles of the ensemble distribution.In this way, we account for the uncertainty in ERFs and the variability across the 18 GCMs 9 .

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Fig. 3 |
Fig. 3 | Changes in cold-related excess mortality fractions (%) at different levels of global warming by country/area under the climate-population and climateonly scenarios, compared with the historical period (1995-2014).Estimates (shown as dots) are reported as the ensemble average of 18 general circulation models under the SSP5-8.5 scenario.The whiskers represent the 95% empirical confidence intervals.The future periods in which the 20-year running mean of global mean temperature first reaches the 1.5 °C, 2 °C, and 3 °C of warming above pre-industrial level (1850-1900) are 2018-2037, 2032-2051, and 2055-2074, respectively under SSP5-8.5.

Fig. 4 |
Fig. 4 | Changes in heat-related excess mortality fractions (%) at different levels of global warming by country/area under the climate-population and climateonly scenarios, compared with the historical period (1995-2014).Estimates (shown as dots) are reported as the ensemble average of 18 general circulation models under the SSP5-8.5 scenario.The whiskers represent the 95% empirical confidence intervals.The future periods in which the 20-year running mean of global mean temperature first reaches the 1.5 °C, 2 °C, and 3 °C of warming above pre-industrial level (1850-1900) are 2018-2037, 2032-2051, and 2055-2074, respectively under SSP5-8.5.

Figure 5
Figure 5 compares the changes in cold-and heat-related excess mortality driven by climate change and population aging at different levels of global warming relative to the historical period (1995-2014) under SSP5-8.5.The climate change component corresponds to the change in temperature-(cold or heat) related mortality fractions in the climateonly scenario, while population aging is calculated from the differencebetween the change in temperature-(cold or heat) related mortality fractions estimated in climate-population scenario and the climateonly scenario (i.e., a constant population scenario).At 1.5 °C warming, all 50 countries/areas will experience a reduction (−1.4% in Vietnam to −0.3% in Iceland) in cold-related mortality burden driven by climate change alone (Fig.5A).However, population aging will mostly offset this decrease and result in an increased cold-related mortality burden in 18 out of the 50 countries/areas (0.0% in Argentina to 1.9% in Kuwait).At higher levels of warming, the decline in cold-related mortality burden due to climate change will be further offset by a larger increase due to population aging, leading to 23 and 21 countries/areas having an increase in cold-related mortality burden under 2 and 3 °C warming, respectively.For heat-related mortality (Fig.5B), climate change alone will result in increased mortality burden in all countries/areas at 1.5 °C warming (0.02% in Iceland to 2.8% in French Guiana), 49 countries/ areas (0.03% in South Africa to 6.9% in French Guiana) except for Colombia (−0.3%) at 2 °C warming, and 48 countries/areas (0.1% in Iceland to 19.3% in French Guiana) except for Colombia [−3.4% (95% eCI: −38.6% to 10.5%)] and South Africa [−0.2% (95% eCI: −3.9% to 2.0%)] at 3 °C warming.The decreases seen in Colombia and South Africa are due to the decreased but non-significant heat-related mortality risks (Supplementary Fig.3).Compared with climate change, population aging will have a smaller contribution to the changes in heat-related mortality, leading to an increase of 0.1-0.7% at 1.5-3 °C warming across 50 countries/areas.This translates into 1 in 5 to 1 in 4 deaths in the 0.5-2.5% increase in heat-related deaths attributable to aging (Table1).

Fig. 5 |
Fig. 5 | Contributions of climate change and population change to the changes in cold-and heat-related mortality at different levels of global warming.A change in cold-related excess mortality fraction (%); B change in heat-related excess mortality fraction (%).Country/area-level changes by climate change and population aging are shown at 1.5 °C, 2 °C, and 3 °C of global warming using a 20year window compared with the historical period 1995-2014 under SSP5-8.5.The future periods in which the 20-year running mean of global mean temperature first reaches the 1.5 °C, 2 °C, and 3 °C of warming above pre-industrial level (1850-1900) are 2018-2037, 2032-2051, and 2055-2074, respectively under SSP5-8.5.The impact of population aging was estimated by subtracting the future changes in temperature-related impacts in the constant population scenario ("climate-only") from the changes in temperature-related impacts under the SSP5 mortality projection in the "climate-population" scenario.Note the different scales in the x-axis used for heat and cold.

Fig. 6 |
Fig. 6 | Contributions of climate change and population change to the changes in non-optimal temperature-related (heat and cold combined) mortality at different levels of global warming.Country/area-level changes by climate change and population aging are shown at 1.5 °C, 2 °C, and 3 °C of global warming using 20year window compared with the historical period 1995-2014 under SSP5-8.5.The future periods in which the 20-year running mean of global mean temperature first reaches the 1.5 °C, 2 °C, and 3 °C of warming above pre-industrial level (1850-1900) are 2018-2037, 2032-2051, and 2055-2074, respectively under SSP5-8.5.The impact of population aging was estimated by subtracting the future changes in temperature-related impacts in the constant population scenario (climate-only) from the changes in temperature-related impacts under SSP5 mortality projection in the climate-population scenario.

Table 1 |
Changes in temperature-related excess mortality fractions (%) at different levels of global warming vs. historical period under the climate-population and climate-only scenarios