Advertisement
The Lancet Planetary Health
Articles| Volume 6, ISSUE 5, e410-e421, May 2022

Download started.

Ok

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

Open AccessPublished:May, 2022DOI:https://doi.org/10.1016/S2542-5196(22)00073-0

      Summary

      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.

      Introduction

      Climate change is a major public health concern in the 21st century. Climate change affects both the global mean surface temperature and its variability, resulting in more frequent extreme weather events and unstable weather.
      • Denman KL
      • Brasseur G
      • Chidthaisong A
      • et al.
      Couplings between changes in the climate system and biogeochemistry.
      • Meehl GA
      • Tebaldi C
      More intense, more frequent, and longer lasting heat waves in the 21st century.
      • McBean G
      Climate change and extreme weather: a basis for action.
      Globally, non-optimum temperatures have been identified as an important indicator of climate change given the largely recognised warming trend, and as one of the leading causes of the global burden of diseases.
      • Stocker TF
      • Qin D
      • Plattner G-K
      • et al.
      Climate change 2013: the physical science basis. Working group I contribution to the fifth assessment report of the intergovernmental panel on climate change.
      • Murray CJ
      • Aravkin AY
      • Zheng P
      • et al.
      Global burden of 87 risk factors in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019.
      • Watts N
      • Amann M
      • Arnell N
      • et al.
      The 2020 report of The Lancet Countdown on health and climate change: responding to converging crises.
      However, temperature variability, another challenging aspect of climate change that reflects weather stability, has less public awareness and has been less investigated compared with non-optimum temperatures.
      • Shi L
      • Kloog I
      • Zanobetti A
      • Liu P
      • Schwartz JD
      Impacts of temperature and its variability on mortality in New England.
      • Cheng J
      • Xu Z
      • Bambrick H
      • Su H
      • Tong S
      • Hu W
      The mortality burden of hourly temperature variability in five capital cities, Australia: time-series and meta-regression analysis.
      Research in context
      Evidence before this study
      In the past few years, attention has been increasingly paid to the negative health effects of short-term temperature variability, an indicator of unstable weather that causes challenges to human adaption. We searched MEDLINE (via PubMed), Web of Science, and Google Scholar from database inception until Nov 6, 2021, for articles in English. We used a combination of search terms, including exposure terms (“temperature variability”, “temperature change”, “temperature fluctuation”, “temperature range”) and health outcome terms (“mortality”, “mortality burden”, “death”, “excess death”). Although many studies have reported the association between temperature variability and premature deaths, most evidence has been obtained from studies in single cities or countries; the systematic evaluation of results is challenged by differences in modelling, parameterisation, and publication bias. One study evaluated the mortality risk associated with temperature variability based on data from 12 countries but did not provide the attributable burden of temperature variability-related mortality. The global burden of mortality associated with temperature variability remains unknown.
      Added value of this study
      To our knowledge, this is the largest and first study, using global gridded observation data at a spatial resolution of 0·5° × 0·5°, to systematically estimate the global burden of mortality associated with temperature variability and explore its temporal trend over 20 years. This study identified an increasing trend of temperature variability at the global level from 2000–19. Globally, 1 753 392 deaths were associated with temperature variability per year, accounting for 3·4% of all deaths. The percentage excess in mortality increased about 4·6% per decade from 2000–19. Disparate geographical variations were also found.
      Implications of all the available evidence
      This study provides robust epidemiological evidence of the impact of temperature variability on mortality globally, and across and within countries or regions. Our findings suggest that regions with a higher percentage excess in mortality due to temperature variability (eg, the whole of Asia, Australia and New Zealand, and northern Africa) are of great importance to contribute to coordinated actions for health, and more targeted policies should be implemented in regions with a higher increasing rate of temperature variability-related mortality burden (eg, Australia and New Zealand, Europe, and the Americas).
      Temperature variability can be measured in many ways, such as diurnal temperature range (DTR), reflecting intra-day temperature variability, and temperature change between two adjacent days (TCN), reflecting inter-day temperature variability.
      • Lee W
      • Kim Y
      • Sera F
      • et al.
      Projections of excess mortality related to diurnal temperature range under climate change scenarios: a multi-country modelling study.
      • Zhan Z
      • Zhao Y
      • Pang S
      • Zhong X
      • Wu C
      • Ding Z
      Temperature change between neighboring days and mortality in United States: a nationwide study.
      Both DTR and TCN are associated with various health outcomes. Given that unstable weather is a continuous process that takes time for human beings to adapt to, a composite index of temperature variability has been developed using the SD of daily minimum and maximum temperatures (Tmin and Tmax) during several exposure days, to assess the combined effect of both intra-day and inter-day temperature variability and their lagged effects.
      • Guo Y
      • Gasparrini A
      • Armstrong BG
      • et al.
      Temperature variability and mortality: a multi-country study.
      • Zhao Q
      • Coelho MSZS
      • Li S
      • et al.
      Spatiotemporal and demographic variation in the association between temperature variability and hospitalizations in Brazil during 2000–2015: a nationwide time-series study.
      Some studies have reported adverse health impacts of short-term exposure to temperature variability, showing a significant association between temperature variability and mortality risk.
      • Guo Y
      • Gasparrini A
      • Armstrong BG
      • et al.
      Temperature variability and mortality: a multi-country study.
      • Zhao Q
      • Coelho MSZS
      • Li S
      • et al.
      Spatiotemporal and demographic variation in the association between temperature variability and hospitalizations in Brazil during 2000–2015: a nationwide time-series study.
      • Yang J
      • Zhou M
      • Li M
      • et al.
      Vulnerability to the impact of temperature variability on mortality in 31 major Chinese cities.
      • Ma C
      • Yang J
      • Nakayama SF
      • Honda Y
      The association between temperature variability and cause-specific mortality: evidence from 47 Japanese prefectures during 1972–2015.
      Our previous study based on Multi-Country Multi-City (MCC) Collaborative Research Network data showed significant but varied associations between temperature variability and mortality risk across 12 countries with various climate patterns, indicating that temperate variability can affect the entire population, but poses a higher risk to particular population groups (eg, children, older people, and people with pre-existing illness).
      • Guo Y
      • Gasparrini A
      • Armstrong BG
      • et al.
      Temperature variability and mortality: a multi-country study.
      • Zhao Q
      • Coelho MSZS
      • Li S
      • et al.
      Spatiotemporal and demographic variation in the association between temperature variability and hospitalizations in Brazil during 2000–2015: a nationwide time-series study.
      • Tian Y
      • Liu H
      • Si Y
      • et al.
      Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: a national time-series study.
      Although several studies have investigated the association between temperature variability and mortality,
      • Guo Y
      • Gasparrini A
      • Armstrong BG
      • et al.
      Temperature variability and mortality: a multi-country study.
      • Zhao Q
      • Coelho MSZS
      • Li S
      • et al.
      Spatiotemporal and demographic variation in the association between temperature variability and hospitalizations in Brazil during 2000–2015: a nationwide time-series study.
      • Yang J
      • Zhou M
      • Li M
      • et al.
      Vulnerability to the impact of temperature variability on mortality in 31 major Chinese cities.
      • Ma C
      • Yang J
      • Nakayama SF
      • Honda Y
      The association between temperature variability and cause-specific mortality: evidence from 47 Japanese prefectures during 1972–2015.
      few have assessed the absolute mortality burden associated with temperature variability.
      We have seen an increased mortality burden attributable to hot temperatures,
      • Zhao Q
      • Guo Y
      • Ye T
      • et al.
      Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study.
      and it would be beneficial to explore how temperature variability-related mortality burden changes over time. Since the pre-industrial era, the global temperature has increased by more than 1°C.
      National Centers for Environmental information
      State of the climate. Global climate report—annual 2019.
      However, temperature variability was observed to vary in time and space without consistent temporal patterns.
      Australian GovernmentBureau of Meteorology
      Annual climate statement 2020.
      • Makowski K
      • Wild M
      • Ohmura A
      Diurnal temperature range over Europe between 1950 and 2005.
      • Sun X
      • Ren G
      • You Q
      • et al.
      Global diurnal temperature range (DTR) changes since 1901.
      The reasons for this temporal–spatial variation can be multifaceted. Dynamic temperature changes are highly correlated with long-wave radiation fluxes, which depend on both natural (eg, atmospheric circulation, cloud cover, and precipitation) and anthropogenic factors (eg, overexploitation and excessive grazing) that vary from region to region.
      • Guo F
      • Do V
      • Cooper R
      • et al.
      Trends of temperature variability: which variability and what health implications?.
      • Huntingford C
      • Jones PD
      • Livina VN
      • Lenton TM
      • Cox PM
      No increase in global temperature variability despite changing regional patterns.
      • Lee W
      • Kim Y
      • Sera F
      • et al.
      Projections of excess mortality related to diurnal temperature range under climate change scenarios: a multi-country modelling study.
      • Schär C
      • Vidale PL
      • Lüthi D
      • et al.
      The role of increasing temperature variability in European summer heatwaves.
      Therefore, it is necessary to understand the temporal trend in temperature variability-related mortality burden across the globe and to be able to make comparisons between regions during the same time window.
      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

      Data sources

      Daily death counts extracted from the MCC Collaborative Research Network database were used in this study. 750 cities across 43 countries or regions were included. International Classification of Diseases, 9th and 10th revision (ICD-9 and ICD-10) codes were used to identify causes of death. We extracted the data series on non-external causes of death (ICD-9: 0–799; ICD-10: A00–R99) or, if not available, all-cause mortality. Descriptive statistics by countries or regions are shown in the appendix (pp 3–4). Only 0·09% of all-cause death data were missing (appendix pp 5–6).
      Daily 1-h Tmax and Tmin data at 0·5° × 0·5° latitude–longitude resolution during 1979–2019 were collected from the Climate Prediction Centre Global Temperature data provided by the National Oceanic and Atmospheric Administration Physical Sciences Laboratory. The dataset originated from Tmin and Tmax data from 6000–7000 stations across the globe and is interpolated using the Shepard algorithm with orographic consideration to develop gridded data.
      • Yang Z
      • Zhang J
      Dataset of high temperature extremes over the major land areas of the Belt and Road for 1979–2018.
      The daily mean temperature was found by calculating the mean Tmin and Tmax. Temperature variability was calculated as the SD of the daily Tmin and Tmax for the current day (lag0) and lag days (eg, lag1, lag2, … lag7).
      • Guo Y
      • Gasparrini A
      • Armstrong BG
      • et al.
      Temperature variability and mortality: a multi-country study.
      For example, temperature variability over lag 0–3 days was calculated as follows:
      SD(Tmax-lag0,Tmax-lag0,Tmax-lag1,Tmax-lag1,Tmax-lag2,Tmax-lag2,Tmax-lag3,Tmax-lag3)


      Data on the global gross domestic product (GDP) and population in 0·5° grid between 1980 and 2020 by 10 years were obtained from the Global Carbon Project.
      • Murakami D
      • Yamagata Y
      Estimation of gridded population and GDP scenarios with spatially explicit statistical downscaling.
      GDP and population data were linearly interpolated over time to generate values for each year. GDP per capita was calculated by dividing the GDP by the population. All GDP per capita data were adjusted to 2010 $US.
      We obtained country-specific mortality rates for each year from the World Bank. For each year, the mean daily deaths for each grid cell were computed as the product of the grid-specific population and annual mortality rate of the country where the grid cell was located, divided by the number of days in a year. Mortality rates were assumed to be identical across all grid cells in the same country, which is widely used in the Global Burden of Diseases, Injuries, and Risk Factors Study.
      • Fitzmaurice C
      • Abate D
      • Abbasi N
      • et al.
      Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the Global Burden of Disease Study.
      GBD 2015 Chronic Respiratory Disease Collaborators
      Global, regional, and national deaths, prevalence, disability-adjusted life years, and years lived with disability for chronic obstructive pulmonary disease and asthma, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015.
      GBD 2015 Mortality and Causes of Death Collaborators
      Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980-2015: a systematic analysis for the Global Burden of Disease Study 2015.
      • Rudd KE
      • Johnson SC
      • Agesa KM
      • et al.
      Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the Global Burden of Disease Study.
      Ethics approval was not required for our analysis of aggregate anonymised data from the MCC Collaborative Research Network database.

      Statistical analysis

      Using a seasonal-trend decomposition procedure based on locally weighted smoothing (STL), we decomposed time-series data for temperature variability into seasonal, trend, and remainder components.
      • Cleveland RB
      • Cleveland WS
      • McRae JE
      • Terpenning I
      STL: a seasonal-trend decomposition.
      We applied the STL method to each grid cell to decompose the time-series data of temperature variability and extract the long-term trend. The global trend of temperature variability was then obtained by calculating the mean long-term trends across all grid cells.
      A three-stage approach established in previous research was applied to quantify the global temperature variability-related mortality burden at a spatial resolution of 0·5° × 0·5°.
      • Zhao Q
      • Guo Y
      • Ye T
      • et al.
      Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study.
      • Gasparrini A
      • Guo Y
      • Sera F
      • et al.
      Projections of temperature-related excess mortality under climate change scenarios.
      Briefly, we firstly obtained location-specific temperature variability–mortality associations based a on daily time series of 750 locations from the MCC Collaborative Research Network. Then, a multivariable meta-regression model was built 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.
      In the first stage, a generalised linear regression model with a quasi-Poisson family was applied in each location to obtain location-specific effect estimates for temperature variability–mortality association. The equation was as follows:
      • Ver Hoef JM
      • Boveng PL
      Quasi-Poisson vs. negative binomial regression: how should we model overdispersed count data?.
      Yit~Poisson(μ;θ)E(Yit)=exp(αi+βiTVit+cb(Tempit,lag=21)+ns(Timeit,df=7/year)+γiDOWit)VAR(Yit)=θμ


      where Yit denotes daily death count in location i on day t; αi represents the intercept in location i; βi and γI represent the coefficients in location i; and TVit stands for the linear function of temperature variability.
      • Guo Y
      • Gasparrini A
      • Armstrong BG
      • et al.
      Temperature variability and mortality: a multi-country study.
      cb(Tempit,lag = 21), built by distributed-lag nonlinear models (DLNMs), is a two-dimensional, parameterised cross-basis function of daily mean temperature. One dimension is for the space of temperature, featuring the non-linear effect of temperature, with a natural cubic spline function with three internal knots placed at the 25th, 50th, and 75th percentile of the location-specific temperature distribution. The other dimension is for the space of lag, featuring the delayed effect of temperature over 21 days of lag, with a natural cubic spline function with two internal knots placed at equally spaced values in the log scale, plus intercept; ns(Timeit,df = 7/year) is a natural cubic spline for time with seven degrees of freedom per year. The number of degrees of freedom determines the flexibility of the spline function. We chose seven degrees of freedom per year to provide adequate control for long-term trend and seasonality. DOWit stands for the day of the week coded as a categorical variable. VAR(Yit) and μ denote the variance and expectation of Yit and θ is an overdispersion parameter. The association between temperature variability and mortality was presented as the relative risk (RR) with 95% CI associated with each 1°C increase in temperature variability. The percentage change in mortality with an IQR increase in temperature variability was also computed.
      In the second stage, a multivariable meta-regression model was built to quantify the relationship between the location-specific effect estimates obtained from the first stage and a set of independent location-specific explanatory variables from each location. We identified five explanatory variables that were documented in previous studies to contribute to the heterogeneity of location-specific effect estimates, including continents, five climate groups of Köppen climate classification, GDP per capita, the yearly average of daily mean temperature, and the range of daily mean temperature.
      • Zhao Q
      • Guo Y
      • Ye T
      • et al.
      Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study.
      • Gasparrini A
      • Guo Y
      • Hashizume M
      • et al.
      Mortality risk attributable to high and low ambient temperature: a multicountry observational study.
      All five variables should have global data at the grid cell level. Mid-year GDP per capita (the middle year of the study period for each location) was calculated to reflect the mean GDP per capita for each location. The performance of multivariable meta-regression models was checked by the I2 statistic. The final model showed an I2 of 22·67% (appendix p 7). The coefficients of five explanatory variables were extracted from the constructed model and used in the third stage to estimate the temperature variability–mortality association at the grid cell level.
      In the third stage, the fitted meta-regression model obtained in the second stage with five grid-specific explanatory variables was used to estimate the temperature variability–mortality association between 2000 and 2019 at the grid cell level.
      We calculated the daily excess deaths associated with temperature variability in each grid cell using the following equation:
      RRit=exp(βper °C increase×TVit)EDit=(RRit-1)×Di


      where RRit is the RR of grid cell i on day; βper 1°C increase is the grid-specific association; TVit is the temperature variability of grid cell i on day t; EDit stands for the excess deaths in grid cell i on day t; and Di is the daily deaths in each grid cell.
      The total number of excess deaths was computed as the sum of daily excess deaths for each year and the entire study period at the global, regional, and national levels. The percentage excess in mortality was calculated by the ratio of excess deaths to total deaths. The average annual percentage excess over 20 years was further computed. Annual excess deaths per 100 000 residents (excess death rate) were also presented. For each region or continent, we calculated the percentage change per decade in both percentage excess in mortality and excess death rate, using a linear regression model considering a Gaussian distribution of percentage excess and excess death rate on the log scale. The 95% CI of percentage change per decade was obtained based on 1000 bootstrap replicates. To make our results easier to follow, we applied the length of exposure of 8 days (temperature variability 0–7) in the main analyses. Results for other lengths of exposure (from temperature variability 0–1 to temperature variability 0–6 and from temperature variability 0–8 to temperature variability 0–10) were shown in the sensitivity analyses.
      Several other sensitivity analyses were done to test the robustness of our results, as follows: extending the maximum lag periods of mean temperature from 21 days to 24 days and 28 days; using alternative degrees of freedom values for time trend (from seven degrees of freedom per year to six degrees of freedom and eight degrees of freedom per year) and lag days of mean temperature (from four degrees of freedom to five degrees of freedom and six degrees of freedom); and controlling the potential effect of relative humidity using a natural cubic spline with three degrees of freedom. A detailed description of sensitivity analyses is shown in the appendix (p 2). The significance of the difference in results between primary analyses and sensitivity analyses was tested using a fixed-effect meta-regression model. Additionally, as we used the counter-factual scenario of no variation in the main analyses, excess deaths represent those that would not have occurred if temperature variability never exceeded 0°C. Considering that temperature variability is less likely to be 0, we also calculated the excess deaths under the counterfactual scenario of the grid-specific minimum temperature variability, by excluding the excess deaths associated with temperature variability ranging from 0 to minimum value, to assess the mortality burden in the more stringent criteria.
      R (version 3.6.2) was used for all analyses. The R packages dlnm (version 2.4.2), mixmeta (version 1.0.8), and stR were used to perform DLNMs, meta-regression models, and seasonal-trend decomposition, respectively.

      Role of the funding source

      The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report.

      Results

      The mean annual temperature variability between 2000 and 2019 is shown in figure 1A. Globally, a large variation in temperature variability was observed. Several regions were identified to have higher temperature variability, such as North America, southern Africa, and northern Africa. The global mean temperature variability was 6·0°C (SD 1·3) in 2000 and 6·2°C (1·3) in 2019 (appendix p 8). After seasonal-trend decomposition, a rising long-term trend in temperature variability was found across the globe (figure 1B; appendix p 9). Among all regions, Australia and New Zealand had the largest increase in annual temperature variability (appendix p 9).
      Figure thumbnail gr1
      Figure 1Mean annual temperature variability at a spatial resolution of 0·5° × 0·5° (A) and the long-term trend in annual temperature variability after seasonal-trend decomposition (B) across the globe from 2000–19
      Triangles in A represent the 750 locations used in the first-stage analysis.
      In general, each IQR increase in temperature variability was associated with a mean 0·7% increase in mortality across all grid cells, with a median value of 0·6% (IQR 0·3–1·0; figure 2A). The country-specific temperature variability related-mortality risks are shown in the appendix (pp 10–14). The geographical variation was observed globally. South Asia had the highest mortality risk associated with temperature variability (figure 2A). Hotspot areas with the biggest contribution to excess deaths were recognised in most parts of south and east Asia (figure 2B). A higher percentage excess in mortality was observed in most of west Asia, the south of middle Asia, and the north of south Asia (figure 2C). The junction of western Africa and central Africa had the highest excess death rate (per 100 000 residents; figure 2D). The changing nature per decade of the percentage excess is shown in figure 2E. The percentage excess on the southeast coast of Australia increased markedly, along with separate areas in the north of western Asia. The excess death rates (per 100 000 residents) were shown to increase in east Asia, the south of North America, and the southeast coast of Australia (figure 2F).
      Figure thumbnail gr2
      Figure 2Percentage change in mortality associated with an IQR (for each grid cell) increase in temperature variability (A), mean annual excess deaths (B), percentage excess in mortality (C), excess deaths per 100 000 residents (D), change in percentage excess in mortality per decade (E), and change in excess deaths per 100 000 residents per decade (F) due to temperature variability in 2000–19 at a spatial resolution of 0·5° × 0·5°
      The scale in E represents change in percentage points, not percentage change.
      From 2000–19, globally, 1 753 392 (95% CI 1 159 901 to 2 357 718) excess deaths was associated with temperature variability per year (table), accounting for 3·4% (95% CI 2·2 to 4·6) of the total deaths and 26 (17 to 35) excess deaths per 100 000 residents (figure 3; appendix p 15). The three leading continents in terms of percentage excess in mortality were Asia (4·7%), Oceania (3·2%), and the Americas (2·7%; figure 3A; appendix p 15). Southern Asia had the highest excess death rate (39 per 100 000 residents, 95% CI 29 to 48) among all regions, whereas the lowest value was observed for other regions in Oceania (nine per 100 000 residents, –4 to 22; figure 3B; appendix p 15). In addition to the region, climate zones contributed to the variation in excess mortality (appendix p 16). Dry climates had the highest percentage excess in mortality (6·0, 95% CI 4·6 to 7·5).
      TablePercentage excess in mortality and excess deaths per 100 000 residents in 2000 and 2019 and percentage change per decade over 2000–19 by continent and region
      Mean annual excess deathsExcess in mortality (%)Excess deaths, per 100 000 residents
      20002019Change per decade (%)
      Percentage change per decade was estimated based on a linear regression model considering a Gaussian distribution of percentage excess and excess death rate on the log scale. The 95% CI of percentage change per decade was obtained based on 1000 bootstrap replicates.
      20002019Change per decade (%)
      Percentage change per decade was estimated based on a linear regression model considering a Gaussian distribution of percentage excess and excess death rate on the log scale. The 95% CI of percentage change per decade was obtained based on 1000 bootstrap replicates.
      Global1 753 392 (1 159 901 to 2 357 718)3·2 (2·1 to 4·3)3·5 (2·3 to 4·7)4·6 (3·7 to 5·3)27·7 (18·2 to 37·4)24·6 (16·4 to 32·9)−7·5 (−9·0 to −5·8)
      Americas160 207 (86 968 to 235 078)2·6 (1·4 to 3·9)2·7 (1·5 to 4·0)3·0 (1·6 to 4·6)18·2 (9·8 to 26·7)18·1 (9·9 to 26·4)1·0 (−1·5 to 4·5)
      Northern America86 097 (50 232 to 122 855)3·1 (1·8 to 4·4)3·1 (1·8 to 4·5)2·8 (1·3 to 4·4)25·9 (15·2 to 36·9)26·0 (15·1 to 37·1)1·4 (−1·3 to 4·0)
      Latin America and the Caribbean74 110 (36 736 to 112 222)2·3 (1·1 to 3·5)2·3 (1·2 to 3·5)3·3 (1·8 to 4·8)13·6 (6·6 to 20·6)13·4 (6·8 to 20·2)0·7 (−2·1 to 4·2)
      Europe127 890 (64 611 to 191 908)1·5 (0·8 to 2·2)1·7 (0·9 to 2·4)4·4 (2·2 to 5·6)17·7 (8·9 to 26·6)17·4 (9·1 to 25·8)−3·8 (−8·2 to −1·2)
      Northern Europe13 236 (8036 to 18 475)1·4 (0·8 to 1·9)1·6 (1·0 to 2·2)4·2 (0·6 to 6·5)14·2 (8·6 to 19·7)13·2 (8·1 to 18·3)−7·5 (−13·1 to −2·8)
      Eastern Europe55 396 (18 825 to 92 378)1·3 (0·4 to 2·1)1·4 (0·5 to 2·3)2·6 (0·7 to 4·0)18·9 (6·5 to 31·4)17·7 (6·4 to 29·2)−8·3 (−11·4 to −5·5)
      Western Europe29 868 (17 884 to 42 046)1·6 (0·9 to 2·2)1·8 (1·1 to 2·5)4·5 (−0·2 to 6·9)15·0 (9·0 to 21·2)17·3 (10·4 to 24·3)5·3 (−0·1 to 9·5)
      Southern Europe29 390 (19 867 to 39 009)2·1 (1·4 to 2·8)2·2 (1·5 to 2·9)2·9 (0·8 to 5·5)20·7 (13·9 to 27·5)19·7 (13·3 to 26·2)−2·0 (−6·7 to 3·1)
      Africa170 094 (50 522 to 291 300)1·7 (0·5 to 2·9)1·8 (0·5 to 3·0)3·3 (1·9 to 4·6)22·2 (6·3 to 38·4)12·9 (3·9 to 22·0)−31·3 (−32·7 to −29·6)
      Northern Africa42 326 (28 057 to 56 784)3·6 (2·4 to 4·8)3·6 (2·4 to 4·8)0·7 (−1·1 to 2·3)23·6 (15·6 to 31·7)19·0 (12·6 to 25·5)−12·2 (−13·4 to −10·7)
      Sub-Saharan Africa127 768 (22 465 to 234 516)1·4 (0·2 to 2·6)1·4 (0·2 to 2·7)0·2 (−1·1 to 1·6)21·9 (3·8 to 40·1)11·4 (1·8 to 21·2)−37·3 (−39·0 to −35·6)
      Asia1 288 284 (955 630 to 1 627 654)4·5 (3·4 to 5·7)4·7 (3·5 to 6·0)2·1 (1·0 to 3·0)33·1 (24·6 to 41·7)30·5 (22·5 to 38·6)−5·3 (−6·9 to −3·6)
      Southeast Asia113 360 (74 167 to 153 030)3·0 (1·9 to 4·0)3·3 (2·2 to 4·5)2·8 (0·6 to 4·8)20·7 (13·7 to 27·8)19·3 (12·6 to 26·1)−9·8 (−13·2 to −4·8)
      Western Asia63 169 (51 417 to 75 071)5·5 (4·5 to 6·5)5·7 (4·7 to 6·8)3·2 (1·7 to 5·4)32·2 (26·3 to 38·2)25·7 (20·9 to 30·6)−10·9 (−13·9 to −7·8)
      Central Asia19 959 (14 820 to 25 190)5·1 (3·8 to 6·4)5·1 (3·8 to 6·4)0·7 (−0·5 to 2·5)37·6 (28·0 to 47·3)28·2 (20·9 to 35·6)−14·9 (−17·1 to −13·4)
      Southern Asia646 213 (483 012 to 812 162)5·2 (3·9 to 6·6)5·2 (3·9 to 6·6)1·1 (−0·4 to 2·4)44·0 (32·9 to 55·4)34·8 (26·0 to 43·8)−12·5 (−15·2 to −10·2)
      Eastern Asia445 583 (332 214 to 562 201)4·1 (3·1 to 5·1)4·5 (3·3 to 5·7)3·1 (1·5 to 4·4)26·9 (20·3 to 33·6)30·7 (22·8 to 38·9)6·2 (4·7 to 7·2)
      Oceania6917 (2171 to 11 778)2·9 (0·7 to 5·2)3·6 (1·2 to 6·1)8·7 (6·1 to 11·6)21·1 (5·4 to 37·3)21·5 (7·3 to 36·2)−4·4 (−9·5 to 1·9)
      Australia and New Zealand6184 (2478 to 9980)3·4 (1·3 to 5·4)4·3 (1·7 to 6·9)7·3 (4·3 to 10·4)22·7 (9·1 to 36·6)26·3 (10·6 to 42·4)0·7 (−3·7 to 6·7)
      Other regions in Oceania733 (−308 to 1798)1·8 (−0·8 to 4·4)1·3 (−0·6 to 3·3)−5·8 (−18·4 to 5·6)15·9 (−6·9 to 39·6)6·9 (−2·9 to 17·0)−37·4 (−52·5 to −23·9)
      All regions in the table are defined according to UN Statistics Division (M49) regional groupings. Other regions in Oceania are defined as all areas outside Australia and New Zealand in Oceania. Corresponding scatter plots from 2000 and 2019 are shown in the appendix (pp 17–20).
      * Percentage change per decade was estimated based on a linear regression model considering a Gaussian distribution of percentage excess and excess death rate on the log scale. The 95% CI of percentage change per decade was obtained based on 1000 bootstrap replicates.
      Figure thumbnail gr3
      Figure 3Average annual percentage excess in mortality (A) and excess deaths per 100 000 residents (B) due to temperature variability in 2000–19 by continent and region
      The global percentage excess in mortality increased from 3·2% (95% CI 2·1 to 4·3) to 3·5% (2·3 to 4·7) between 2000 and 2019, representing an increased rate of 4·6% (3·7 to 5·3) per decade (table). Australia and New Zealand generated the largest increase in percentage excess, increasing from 3·4% (1·3 to 5·4) in 2000 to 4·3% (1·7 to 6·9) in 2019, representing an increased rate of 7·3% (4·3 to 10·4) per decade. The largest decline occurred in other regions in Oceania, with a decreased rate of 5·8% per decade (–18·4 to 5·6; table). Corresponding scatter plots from 2000 to 2019 are shown in the appendix (pp 17–20).
      Lists of the top 20 countries ranked by temperature variability-related mortality burden in both 2000 and 2019 included many of the same countries, but the order changed (appendix p 21; figure 4). Among the top ten countries in percentage excess in 2019, four were listed in the current world bank high-income economies, including Saudi Arabia (first), Kuwait (third), United Arab Emirates (fourth), and Qatar (eighth) (figure 4A). Compared with the percentage excess, excess death rates (per 100 000 residents) decreased markedly during the 20 years studied (figure 4B).
      Figure thumbnail gr4
      Figure 4Leading 20 countries for percentage excess in mortality (A) and excess deaths per 100 000 residents (B) in 2000 and 2019
      In the sensitivity analyses, the mortality burden associated with temperature variability decreased with shorter exposure to temperature variability (appendix p 22). After changing the model parameters, our results changed slightly (appendix pp 23–24). When the counterfactual scenario of grid-specific minimum temperature variability was applied, the percentage excess was 2·2% (95% CI 1·5–2·9), nearly two-thirds of that under the counterfactual scenario of zero temperature variability (appendix p 25).

      Discussion

      To our knowledge, this is the largest and first study to use global gridded observation data at a spatial resolution of 0·5° × 0·5° to systematically estimate the global burden of mortality associated with temperature variability and explore its temporal trend over 20 years. From 2000 to 2019, the daily mean value of temperature variability generally increased. A considerable number of deaths were associated with temperature variability per year, causing a substantial mortality burden worldwide. An increasing trend of the percentage excess in mortality was observed during the 20 years studied.
      Consistent with previous studies,
      • Yang J
      • Zhou M
      • Li M
      • et al.
      Vulnerability to the impact of temperature variability on mortality in 31 major Chinese cities.
      • Ma C
      • Yang J
      • Nakayama SF
      • Honda Y
      The association between temperature variability and cause-specific mortality: evidence from 47 Japanese prefectures during 1972–2015.
      • Hu K
      • Guo Y
      • Yang X
      • et al.
      Temperature variability and mortality in rural and urban areas in Zhejiang province, China: an application of a spatiotemporal index.
      • Cheng J
      • Xu Z
      • Bambrick H
      • Su H
      • Tong S
      • Hu W
      Impacts of heat, cold, and temperature variability on mortality in Australia, 2000-2009.
      • Tian Y
      • Liu H
      • Si Y
      • et al.
      Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: a national time-series study.
      we observed an increased mortality risk associated with temperature variability, accounting for a substantial mortality burden. The percentage change in mortality associated with an IQR increase in temperature variability ranged from 0 to 2% for most grid cells, which is similar to previous results based on 12 countries.
      • Guo Y
      • Gasparrini A
      • Armstrong BG
      • et al.
      Temperature variability and mortality: a multi-country study.
      The physiological mechanisms underlying this association might relate to thermal adjustment to temperatures through physiological and behavioural responses that are impeded by unstable weather over a short period of time.
      • Cheng J
      • Xu Z
      • Bambrick H
      • Su H
      • Tong S
      • Hu W
      The mortality burden of hourly temperature variability in five capital cities, Australia: time-series and meta-regression analysis.
      • Garrett AT
      • Goosens NG
      • Rehrer NJ
      • Patterson MJ
      • Cotter JD
      Induction and decay of short-term heat acclimation.
      During these processes, multiple organs can be involved (eg, respiratory, circulatory, and immune systems) by affecting heart rate, blood viscosity, fibrinogen, platelet count, arterial blood pressure, and oxygen uptake.
      • Lee W
      • Kim Y
      • Honda Y
      • Kim H
      Association between diurnal temperature range and mortality modified by temperature in Japan, 1972–2015: investigation of spatial and temporal patterns for 12 cause-specific deaths.
      • Martinez-Nicolas A
      • Meyer M
      • Hunkler S
      • et al.
      Daytime variation in ambient temperature affects skin temperatures and blood pressure: ambulatory winter/summer comparison in healthy young women.
      • Ballester F
      • Corella D
      • Pérez-Hoyos S
      • Sáez M
      • Hervás A
      Mortality as a function of temperature. A study in Valencia, Spain, 1991–1993.
      Although the biological mechanisms have not been fully elucidated, they imply a difficult process of thermal adjustment to temperature variability.
      To protect human health against temperature variability, proactive countermeasures such as warning systems, community-level responses, and instructions for self-protection are necessary. Many policies have been developed to cope with the threat of climate-related extreme events, for example, warning systems for heatwaves and air pollution.
      • Brennan M
      • O'Shea PM
      • Mulkerrin EC
      Preventative strategies and interventions to improve outcomes during heatwaves.
      • Lowe D
      • Ebi KL
      • Forsberg B
      Heatwave early warning systems and adaptation advice to reduce human health consequences of heatwaves.
      However, policies and strategies rarely exist to effectively cope with the adverse health impacts of temperature variability. Previous investigations reported an estimate of 7·6% for attributable mortality caused by ambient air pollution.
      • Cohen AJ
      • Brauer M
      • Burnett R
      • et al.
      Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015.
      Investigators who separated the hot and cold impact from non-optimum temperatures suggested an excess death ratios of 8·5% for cold-related temperatures and 0·9% for hot-related temperatures.
      • Zhao Q
      • Guo Y
      • Ye T
      • et al.
      Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study.
      As suggested by our findings, temperature variability has similar impacts to air pollution and non-optimum temperatures on global mortality. More attention should be paid to the health impacts of temperature variability. One solution is to develop early warning systems of temperature variability, as personal protection behaviours are highly correlated with risk perception. For example, outdoor workers can prepare enough clothes in advance with a timely warning to protect them from sudden temperature changes. Development of guidance on self-protection (eg, stay indoors, take clothes, and take care of vulnerable populations, such as children and older people) with community social programmes will be of great benefit to help people understand what they need to do when temperature fluctuates over a short time period. In the long run, measures to reduce the impact of climate change (eg, clean energy and greenhouse gas emission reduction) should be promoted to fundamentally mitigate global warming, as well as the increasing trend of temperature variability, although these measures might take time to implement and have an impact.
      • Wang F
      • Harindintwali JD
      • Yuan Z
      • et al.
      Technologies and perspectives for achieving carbon neutrality.
      Regions with a higher percentage excess in mortality due to temperature variability (eg, Asia, Australia and New Zealand, and northern Africa) are of great importance to contribute to coordinated actions for health. Some countries, especially developing countries, will suffer disproportionately more from the adverse effects of global climate change, which could be a potential driver for international inequality.
      • King AD
      • Harrington LJ
      The inequality of climate change from 1·5 to 2 C of global warming.
      In this study, we observed a small but significant increasing trend in both temperature variability and temperature-related percentage excess in mortality globally. The reduction in the excess death rate from 2000 to 2019 might be largely due to a decrease in mortality rate, whereas almost all regions showed an upward trend in percentage excess in mortality associated with temperature variability, indicating a persistent impact in the past two decades. Although few studies focused on the temporal trends of temperature variability-related mortality burden, investigations can be made through the relevant assessment of similar indicators. For example, a multicountry study of 20 countries or regions projected that a 1·4–10·3% increase in excess deaths attributable to the DTR will happen by the end of this century; the study inferred a more pronounced mortality burden in the future due to more unstable weather than in the past, although there might be an adaptation to climate change benefiting from socioeconomic development and investment in public health.
      • Lee W
      • Kim Y
      • Sera F
      • et al.
      Projections of excess mortality related to diurnal temperature range under climate change scenarios: a multi-country modelling study.
      • Todd N
      • Valleron A-J
      Space–time covariation of mortality with temperature: a systematic study of deaths in France, 1968–2009.
      • Åström DO
      • Tornevi A
      • Ebi KL
      • Rocklöv J
      • Forsberg B
      Evolution of minimum mortality temperature in Stockholm, Sweden, 1901–2009.
      More targeted policies should be implemented to avoid the negative health impacts of temperature variability, especially for regions with a higher increasing rate of temperature variability-related mortality burden (eg, Australia and New Zealand, northern Europe, Latin America and the Caribbean, and western Asia).
      This study has several strengths. First, to our knowledge, this is the first and largest study to systematically explore the mortality burden associated with temperature variability on a global scale. Compared with previous studies that were restricted to single or several countries,
      • Ma C
      • Yang J
      • Nakayama SF
      • Honda Y
      The association between temperature variability and cause-specific mortality: evidence from 47 Japanese prefectures during 1972–2015.
      • Tian Y
      • Liu H
      • Si Y
      • et al.
      Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: a national time-series study.
      • Zhang Y
      • Xiang Q
      • Yu C
      • et al.
      Mortality risk and burden associated with temperature variability in China, United Kingdom and United States: comparative analysis of daily and hourly exposure metrics.
      this study offers a finer spatial view of the mortality burden associated with temperature variability, which can provide new clues on geographical variations and allow within-country comparisons. Second, this study benefits from global gridded population and climate data. To minimise potential exposure misclassification from aggregating individual exposure to location or country level (aggregation bias), we used exposure data in a 0·5° grid to produce better countrywide and global estimates. Finally, we considered spatiotemporal trends over a 20-year period of fast climate change. The findings of this study provide a better understanding of how temperature variability has affected human health amid inevitable warming trends and gradual acclimatisation to climate change.
      This study also has some limitations. We used country-specific mortality rates rather than grid-specific mortality rates because of insufficient data. The assumption of an identical mortality rate across grid cells within the same country is widely used in the Global Burden of Diseases, Injuries, and Risk Factors Study,
      • Fitzmaurice C
      • Abate D
      • Abbasi N
      • et al.
      Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the Global Burden of Disease Study.
      GBD 2015 Chronic Respiratory Disease Collaborators
      Global, regional, and national deaths, prevalence, disability-adjusted life years, and years lived with disability for chronic obstructive pulmonary disease and asthma, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015.
      GBD 2015 Mortality and Causes of Death Collaborators
      Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980-2015: a systematic analysis for the Global Burden of Disease Study 2015.
      • Rudd KE
      • Johnson SC
      • Agesa KM
      • et al.
      Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the Global Burden of Disease Study.
      and should not have a substantial effect on our estimations at country, region, and global levels. However, this fact limits our ability to identify variation in temperature variation-related mortality burden within countries and these limitations should be supplemented in future studies by collection of mortality data at a finer level. Grid cell-specific data (eg, temperature, population, and GDP) applied in this study also introduced uncertainties, as interpolation and downscaling produced prediction errors.
      Although previous studies suggested variation in susceptibility to temperature variability across age, sex, and causes of mortality,
      • Tian Y
      • Liu H
      • Si Y
      • et al.
      Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: a national time-series study.
      • Wu Y
      • Xu R
      • Wen B
      • et al.
      Temperature variability and asthma hospitalisation in Brazil, 2000–2015: a nationwide case-crossover study.
      our study did not characterise these differences because of a lack of age-specific, sex-specific, and cause-specific mortality data at both the grid cell and country level. Future research could complement the evidence provided in this study if relevant data are available. Several predictors that could explain the heterogeneity in the grid-specific temperature variability related-mortality associations were used to estimate these associations. However, we must acknowledge that there could be unexplained heterogeneity contributed by both the paucity of grid-specific data and unknown factors. Further studies are warranted to provide more precise estimates of this association. Owing to the nature of time-series designs, we could not investigate the causal relationship between temperature variability and mortality. Grid cell-specific estimations were based on the temperature variability related-mortality association, but not causation. Finally, MCC data include only 43 countries or regions and have limited information on countries located in the Sahara desert. This paucity of data might affect the accuracy of effect estimates. Although we used Köppen climate classification as one of the predictors in the model, future studies are needed to further explore the association between temperature variability and mortality in the desert area.
      In conclusion, this study highlights the substantial mortality burden associated with temperature variability. This burden had a complex pattern of variation globally and a slightly increasing temporal trend over the past two decades. Considering climate change, our findings could assist in raising public awareness and improving the understanding of the impacts of temperature variability on health.
      Contributors
      YG, AG, MH, and BAr set up the collaborative network. YG, SL, YW, and QZ designed the study. YG, SL, QZ, and AG developed the statistical methods. YW, SL, and BW took the lead in drafting the manuscript and interpreting the results. All authors provided the data and contributed to the interpretation of the results and the submitted version of the manuscript. YG, SL, YW, and QZ accessed and verified the data. All authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.

      Data sharing

      Data were collected within the Multi-Country Multi-City (MCC) Collaborative Research Network under a data sharing agreement and cannot be made publicly available. Researchers can refer to MCC participants, who are listed as coauthors of this Article, for information on accessing the data for each country.

      Declaration of interests

      We declare no competing interests.
      Acknowledgments
      This study was supported by the Australian Research Council (DP210102076) and the Australian National Health and Medical Research Council (APP2000581). YW was supported by the China Scholarship Council (number 202006010044). SL was supported by an Emerging Leader Fellowship of the Australian National Health and Medical Research Council (number APP2009866). QZ was supported by the Program of Qilu Young Scholars of Shandong University, Jinan, China. BW was supported by the China Scholarship Council (number 202006010043). JK and AU were supported by the Czech Science Foundation (project number 20–28560S). NS was supported by the National Institute of Environmental Health Sciences-funded HERCULES Center (P30ES019776). S-CP and YLG were supported by the Ministry of Science and Technology (Taiwan; MOST 109–2621-M-002–021). YH was supported by the Environment Research and Technology Development Fund (JPMEERF15S11412) of the Environmental Restoration and Conservation Agency. MdSZSC and PHNS were supported by the São Paulo Research Foundation (FAPESP). ST was supported by the Science and Technology Commission of Shanghai Municipality (grant number 18411951600). HO and EI were supported by the Estonian Ministry of Education and Research (IUT34–17). JM was supported by a fellowship of Fundação para a Ciência e a Tecnlogia (SFRH/BPD/115112/2016). AG and FS 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). AS, SR, and FdD were supported by the EU's Horizon 2020 project, Exhaustion (grant ID 820655). VH was supported by the Spanish Ministry of Economy, Industry and Competitiveness (grant ID PCIN-2017–046). AT was supported by MCIN/AEI/10.13039/501100011033 (grant CEX2018-000794-S). YG was supported by the Career Development Fellowship (number APP1163693) and Leader Fellowship (number APP2008813) of the Australian National Health and Medical Research Council. Statistics South Africa kindly provided the mortality data, but had no other role in the study. This Article is published in memory of Simona Fratianni, who helped to contribute the data for Romania.
      Editorial note: the Lancet Group takes a neutral position with respect to territorial claims in published maps and institutional affiliations.

      Supplementary Material

      References

      1. 1.
        • Denman KL
        • Brasseur G
        • Chidthaisong A
        • et al.
        Couplings between changes in the climate system and biogeochemistry.
        in: Soloman SD Qin M Massing Z Climate change 2007: the physical science basis. Cambridge University Press, Cambridge2007: 499-587
      2. 2.
        • Meehl GA
        • Tebaldi C
        More intense, more frequent, and longer lasting heat waves in the 21st century.
        Science. 2004; 305: 994-997
      3. 3.
        • McBean G
        Climate change and extreme weather: a basis for action.
        Nat Hazards. 2004; 31: 177-190
      4. 4.
        • Stocker TF
        • Qin D
        • Plattner G-K
        • et al.
        Climate change 2013: the physical science basis. Working group I contribution to the fifth assessment report of the intergovernmental panel on climate change.
      5. 5.
        • Murray CJ
        • Aravkin AY
        • Zheng P
        • et al.
        Global burden of 87 risk factors in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019.
        Lancet. 2020; 396: 1223-1249
      6. 6.
        • Watts N
        • Amann M
        • Arnell N
        • et al.
        The 2020 report of The Lancet Countdown on health and climate change: responding to converging crises.
        Lancet. 2021; 397: 129-170
      7. 7.
        • Shi L
        • Kloog I
        • Zanobetti A
        • Liu P
        • Schwartz JD
        Impacts of temperature and its variability on mortality in New England.
        Nat Clim Chang. 2015; 5: 988-991
      8. 8.
        • Cheng J
        • Xu Z
        • Bambrick H
        • Su H
        • Tong S
        • Hu W
        The mortality burden of hourly temperature variability in five capital cities, Australia: time-series and meta-regression analysis.
        Environ Int. 2017; 109: 10-19
      9. 9.
        • Lee W
        • Kim Y
        • Sera F
        • et al.
        Projections of excess mortality related to diurnal temperature range under climate change scenarios: a multi-country modelling study.
        Lancet Planet Health. 2020; 4: e512-e521
      10. 10.
        • Zhan Z
        • Zhao Y
        • Pang S
        • Zhong X
        • Wu C
        • Ding Z
        Temperature change between neighboring days and mortality in United States: a nationwide study.
        Sci Total Environ. 2017; 584-585: 1152-1161
      11. 11.
        • Guo Y
        • Gasparrini A
        • Armstrong BG
        • et al.
        Temperature variability and mortality: a multi-country study.
        Environ Health Perspect. 2016; 124: 1554-1559
      12. 12.
        • Zhao Q
        • Coelho MSZS
        • Li S
        • et al.
        Spatiotemporal and demographic variation in the association between temperature variability and hospitalizations in Brazil during 2000–2015: a nationwide time-series study.
        Environ Int. 2018; 120: 345-353
      13. 13.
        • Yang J
        • Zhou M
        • Li M
        • et al.
        Vulnerability to the impact of temperature variability on mortality in 31 major Chinese cities.
        Environ Pollut. 2018; 239: 631-637
      14. 14.
        • Ma C
        • Yang J
        • Nakayama SF
        • Honda Y
        The association between temperature variability and cause-specific mortality: evidence from 47 Japanese prefectures during 1972–2015.
        Environ Int. 2019; 127: 125-133
      15. 15.
        • Tian Y
        • Liu H
        • Si Y
        • et al.
        Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: a national time-series study.
        PLoS Med. 2019; 16e1002738
      16. 16.
        • Zhao Q
        • Guo Y
        • Ye T
        • et al.
        Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study.
        Lancet Planet Health. 2021; 5: e415-e425
      17. 17.
        • National Centers for Environmental information
        State of the climate. Global climate report—annual 2019.
        https://www.ncdc.noaa.gov/sotc/global/201913
        Date: 2020
        Date accessed: September 27, 2021
      18. 18.
        • Australian Government
        • Bureau of Meteorology
        Annual climate statement 2020.
        http://www.bom.gov.au/climate/current/annual/aus/#tabs=About
        Date: 2021
        Date accessed: September 27, 2021
      19. 19.
        • Makowski K
        • Wild M
        • Ohmura A
        Diurnal temperature range over Europe between 1950 and 2005.
        Atmos Chem Phys. 2008; 8: 6483-6498
      20. 20.
        • Sun X
        • Ren G
        • You Q
        • et al.
        Global diurnal temperature range (DTR) changes since 1901.
        Clim Dyn. 2019; 52: 3343-3356
      21. 21.
        • Guo F
        • Do V
        • Cooper R
        • et al.
        Trends of temperature variability: which variability and what health implications?.
        Sci Total Environ. 2021; 768144487
      22. 22.
        • Huntingford C
        • Jones PD
        • Livina VN
        • Lenton TM
        • Cox PM
        No increase in global temperature variability despite changing regional patterns.
        Nature. 2013; 500: 327-330
      23. 23.
        • Lee W
        • Kim Y
        • Sera F
        • et al.
        Projections of excess mortality related to diurnal temperature range under climate change scenarios: a multi-country modelling study.
        Lancet Planet Health. 2020; 4: e512-e521
      24. 24.
        • Schär C
        • Vidale PL
        • Lüthi D
        • et al.
        The role of increasing temperature variability in European summer heatwaves.
        Nature. 2004; 427: 332-336
      25. 25.
        • Yang Z
        • Zhang J
        Dataset of high temperature extremes over the major land areas of the Belt and Road for 1979–2018.
        Big Earth Data. 2020; 4: 128-141
      26. 26.
        • Murakami D
        • Yamagata Y
        Estimation of gridded population and GDP scenarios with spatially explicit statistical downscaling.
        Sustainability (Basel). 2019; 112106
      27. 27.
        • Fitzmaurice C
        • Abate D
        • Abbasi N
        • et al.
        Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the Global Burden of Disease Study.
        JAMA Oncol. 2019; 5: 1749-1768
      28. 28.
        • GBD 2015 Chronic Respiratory Disease Collaborators
        Global, regional, and national deaths, prevalence, disability-adjusted life years, and years lived with disability for chronic obstructive pulmonary disease and asthma, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015.
        Lancet Respir Med. 2017; 5: 691-706
      29. 29.
        • GBD 2015 Mortality and Causes of Death Collaborators
        Global, regional, and national life expectancy, all-cause mortality, and cause-specific mortality for 249 causes of death, 1980-2015: a systematic analysis for the Global Burden of Disease Study 2015.
        Lancet. 2016; 388: 1459-1544
      30. 30.
        • Rudd KE
        • Johnson SC
        • Agesa KM
        • et al.
        Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the Global Burden of Disease Study.
        Lancet. 2020; 395: 200-211
      31. 31.
        • Cleveland RB
        • Cleveland WS
        • McRae JE
        • Terpenning I
        STL: a seasonal-trend decomposition.
        J Off Stat. 1990; 6: 3-73
      32. 32.
        • Gasparrini A
        • Guo Y
        • Sera F
        • et al.
        Projections of temperature-related excess mortality under climate change scenarios.
        Lancet Planet Health. 2017; 1: e360-e367
      33. 33.
        • Ver Hoef JM
        • Boveng PL
        Quasi-Poisson vs. negative binomial regression: how should we model overdispersed count data?.
        Ecology. 2007; 88: 2766-2772
      34. 34.
        • Gasparrini A
        • Guo Y
        • Hashizume M
        • et al.
        Mortality risk attributable to high and low ambient temperature: a multicountry observational study.
        Lancet. 2015; 386: 369-375
      35. 35.
        • Hu K
        • Guo Y
        • Yang X
        • et al.
        Temperature variability and mortality in rural and urban areas in Zhejiang province, China: an application of a spatiotemporal index.
        Sci Total Environ. 2019; 647: 1044-1051
      36. 36.
        • Cheng J
        • Xu Z
        • Bambrick H
        • Su H
        • Tong S
        • Hu W
        Impacts of heat, cold, and temperature variability on mortality in Australia, 2000-2009.
        Sci Total Environ. 2019; 651: 2558-2565
      37. 37.
        • Tian Y
        • Liu H
        • Si Y
        • et al.
        Association between temperature variability and daily hospital admissions for cause-specific cardiovascular disease in urban China: a national time-series study.
        PLoS Med. 2019; 16e1002738
      38. 38.
        • Garrett AT
        • Goosens NG
        • Rehrer NJ
        • Patterson MJ
        • Cotter JD
        Induction and decay of short-term heat acclimation.
        Eur J Appl Physiol. 2009; 107: 659-670
      39. 39.
        • Lee W
        • Kim Y
        • Honda Y
        • Kim H
        Association between diurnal temperature range and mortality modified by temperature in Japan, 1972–2015: investigation of spatial and temporal patterns for 12 cause-specific deaths.
        Environ Int. 2018; 119: 379-387
      40. 40.
        • Martinez-Nicolas A
        • Meyer M
        • Hunkler S
        • et al.
        Daytime variation in ambient temperature affects skin temperatures and blood pressure: ambulatory winter/summer comparison in healthy young women.
        Physiol Behav. 2015; 149: 203-211
      41. 41.
        • Ballester F
        • Corella D
        • Pérez-Hoyos S
        • Sáez M
        • Hervás A
        Mortality as a function of temperature. A study in Valencia, Spain, 1991–1993.
        Int J Epidemiol. 1997; 26: 551-561
      42. 42.
        • Brennan M
        • O'Shea PM
        • Mulkerrin EC
        Preventative strategies and interventions to improve outcomes during heatwaves.
        Age Ageing. 2020; 49: 729-732
      43. 43.
        • Lowe D
        • Ebi KL
        • Forsberg B
        Heatwave early warning systems and adaptation advice to reduce human health consequences of heatwaves.
        Int J Environ Res Public Health. 2011; 8: 4623-4648
      44. 44.
        • Cohen AJ
        • Brauer M
        • Burnett R
        • et al.
        Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015.
        Lancet. 2017; 389: 1907-1918
      45. 45.
        • Wang F
        • Harindintwali JD
        • Yuan Z
        • et al.
        Technologies and perspectives for achieving carbon neutrality.
        Innovation (N Y). 2021; 2100180
      46. 46.
        • King AD
        • Harrington LJ
        The inequality of climate change from 1·5 to 2 C of global warming.
        Geophys Res Lett. 2018; 45: 5030-5033
      47. 47.
        • Todd N
        • Valleron A-J
        Space–time covariation of mortality with temperature: a systematic study of deaths in France, 1968–2009.
        Environ Health Perspect. 2015; 123: 659-664
      48. 48.
        • Åström DO
        • Tornevi A
        • Ebi KL
        • Rocklöv J
        • Forsberg B
        Evolution of minimum mortality temperature in Stockholm, Sweden, 1901–2009.
        Environ Health Perspect. 2016; 124: 740-744
      49. 49.
        • Zhang Y
        • Xiang Q
        • Yu C
        • et al.
        Mortality risk and burden associated with temperature variability in China, United Kingdom and United States: comparative analysis of daily and hourly exposure metrics.
        Environ Res. 2019; 179108771
      50. 50.
        • Wu Y
        • Xu R
        • Wen B
        • et al.
        Temperature variability and asthma hospitalisation in Brazil, 2000–2015: a nationwide case-crossover study.
        Thorax. 2021; 76: 962-969