Global, regional, and national burden of heatwave-related mortality from 1990 to 2019: A three-stage modelling study.

BACKGROUND: The regional disparity of heatwave-related mortality over a long period has not been sufficiently assessed across the globe, impeding the localisation of adaptation planning and risk management towards climate change. We quantified the global mortality burden associated with heatwaves at a spatial resolution of 0.5°×0.5° and the temporal change from 1990 to 2019. METHODS AND FINDINGS: We collected data on daily deaths and temperature from 750 locations of 43 countries or regions, and 5 meta-predictors in 0.5°×0.5° resolution across the world. Heatwaves were defined as location-specific daily mean temperature ≥95th percentiles of year-round temperature range with duration ≥2 days. We first estimated the location-specific heatwave-mortality association. Secondly, a multivariate meta-regression was fitted between location-specific associations and 5 meta-predictors, which was in the third stage used with grid cell-specific meta-predictors to predict grid cell-specific association. Heatwave-related excess deaths were calculated for each grid and aggregated. During 1990 to 2019, 0.94% (95% CI: 0.68-1.19) of deaths [i.e., 153,078 cases (95% eCI: 109,950-194,227)] per warm season were estimated to be from heatwaves, accounting for 236 (95% eCI: 170-300) deaths per 10 million residents. The ratio between heatwave-related excess deaths and all premature deaths per warm season remained relatively unchanged over the 30 years, while the number of heatwave-related excess deaths per 10 million residents per warm season declined by 7.2% per decade in comparison to the 30-year average. Locations with the highest heatwave-related death ratio and rate were in Southern and Eastern Europe or areas had polar and alpine climates, and/or their residents had high incomes. The temporal change of heatwave-related mortality burden showed geographic disparities, such that locations with tropical climate or low incomes were observed with the greatest decline. The main limitation of this study was the lack of data from certain regions, e.g., Arabian Peninsula and South Asia. CONCLUSIONS: Heatwaves were associated with substantial mortality burden that varied spatiotemporally over the globe in the past 30 years. The findings indicate the potential benefit of governmental actions to enhance health sector adaptation and resilience, accounting for inequalities across communities.

Gender inequalities in heat-related mortality in the Czech Republic.

It is acknowledged that climate change exacerbates social inequalities, and women have been reported as more vulnerable to heat than men in many studies in Europe, including the Czech Republic. This study aimed at investigating the associations between daily temperature and mortality in the Czech Republic in the light of a sex and gender perspective, taking into account other factors such as age and marital status. Daily mean temperature and individual mortality data recorded during the five warmest months of the year (from May to September) over the period 1995-2019 were used to fit a quasi-Poisson regression model, which included a distributed lag non-linear model (DLNM) to account for the delayed and non-linear effects of temperature on mortality. The heat-related mortality risks obtained in each population group were expressed in terms of risk at the 99th percentile of summer temperature relative to the minimum mortality temperature. Women were found generally more at risk to die because of heat than men, and the difference was larger among people over 85 years old. Risks among married people were lower than risks among single, divorced, and widowed people, while risks in divorced women were significantly higher than in divorced men. This is a novel finding which highlights the potential role of gender inequalities in heat-related mortality. Our study underlines the relevance of including a sex and gender dimension in the analysis of the impacts of heat on the population and advocates the development of gender-based adaptation policies to extreme heat.

Ambient Particulate Air Pollution and Daily Mortality in 652 Cities.

BACKGROUND: The systematic evaluation of the results of time-series studies of air pollution is challenged by differences in model specification and publication bias. METHODS: We evaluated the associations of inhalable particulate matter (PM) with an aerodynamic diameter of 10 µm or less (PM10) and fine PM with an aerodynamic diameter of 2.5 µm or less (PM2.5) with daily all-cause, cardiovascular, and respiratory mortality across multiple countries or regions. Daily data on mortality and air pollution were collected from 652 cities in 24 countries or regions. We used overdispersed generalized additive models with random-effects meta-analysis to investigate the associations. Two-pollutant models were fitted to test the robustness of the associations. Concentration-response curves from each city were pooled to allow global estimates to be derived. RESULTS: On average, an increase of 10 µg per cubic meter in the 2-day moving average of PM10 concentration, which represents the average over the current and previous day, was associated with increases of 0.44% (95% confidence interval [CI], 0.39 to 0.50) in daily all-cause mortality, 0.36% (95% CI, 0.30 to 0.43) in daily cardiovascular mortality, and 0.47% (95% CI, 0.35 to 0.58) in daily respiratory mortality. The corresponding increases in daily mortality for the same change in PM2.5 concentration were 0.68% (95% CI, 0.59 to 0.77), 0.55% (95% CI, 0.45 to 0.66), and 0.74% (95% CI, 0.53 to 0.95). These associations remained significant after adjustment for gaseous pollutants. Associations were stronger in locations with lower annual mean PM concentrations and higher annual mean temperatures. The pooled concentration-response curves showed a consistent increase in daily mortality with increasing PM concentration, with steeper slopes at lower PM concentrations. CONCLUSIONS: Our data show independent associations between short-term exposure to PM10 and PM2.5 and daily all-cause, cardiovascular, and respiratory mortality in more than 600 cities across the globe. These data reinforce the evidence of a link between mortality and PM concentration established in regional and local studies. (Funded by the National Natural Science Foundation of China and others.).

Evaluation of the ERA5 reanalysis-based Universal Thermal Climate Index on mortality data in Europe.

Air temperature has been the most commonly used exposure metric in assessing relationships between thermal stress and mortality. Lack of the high-quality meteorological station data necessary to adequately characterize the thermal environment has been one of the main limitations for the use of more complex thermal indices. Global climate reanalyses may provide an ideal platform to overcome this limitation and define complex heat and cold stress conditions anywhere in the world. In this study, we explored the potential of the Universal Thermal Climate Index (UTCI) based on ERA5 - the latest global climate reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF) - as a health-related tool. Employing a novel ERA5-based thermal comfort dataset ERA5-HEAT, we investigated the relationships between the UTCI and daily mortality data in 21 cities across 9 European countries. We used distributed lag nonlinear models to assess exposure-response relationships between mortality and thermal conditions in individual cities. We then employed meta-regression models to pool the results for each city into four groups according to climate zone. To evaluate the performance of ERA5-based UTCI, we compared its effects on mortality with those for the station-based UTCI data. In order to assess the additional effect of the UTCI, the performance of ERA5-and station-based air temperature (T) was evaluated. Whilst generally similar heat- and cold-effects were observed for the ERA5-and station-based data in most locations, the important role of wind in the UTCI appeared in the results. The largest difference between any two datasets was found in the Southern European group of cities, where the relative risk of mortality at the 1st percentile of daily mean temperature distribution (1.29 and 1.30 according to the ERA5 vs station data, respectively) considerably exceeded the one for the daily mean UTCI (1.19 vs 1.22). These differences were mainly due to the effect of wind in the cold tail of the UTCI distribution. The comparison of exposure-response relationships between ERA5-and station-based data shows that ERA5-based UTCI may be a useful tool for definition of life-threatening thermal conditions in locations where high-quality station data are not available.

The predictability of heat-related mortality in Prague, Czech Republic, during 2015-a comparison of selected thermal indices.

We compared selected thermal indices in their ability to predict heat-related mortality in Prague, Czech Republic, during the extraordinary summer 2015. Relatively, novel thermal indices-Universal Thermal Climate Index and Excess Heat Factor (EHF)-were compared with more traditional ones (apparent temperature, simplified wet-bulb globe temperature (WBGT), and physiologically equivalent temperature). The relationships between thermal indices and all-cause relative mortality deviations from the baseline (excess mortality) were estimated by generalized additive models for the extended summer season (May-September) during 1994-2014. The resulting models were applied to predict excess mortality in 2015 based on observed meteorology, and the mortality estimates by different indices were compared. Although all predictors showed a clear association between thermal conditions and excess mortality, we found important variability in their performance. The EHF formula performed best in estimating the intensity of heat waves and magnitude of heat-impacts on excess mortality on the most extreme days. Afternoon WBGT, on the other hand, was most precise in the selection of heat-alert days during the extended summer season, mainly due to a relatively small number of "false alerts" compared to other predictors. Since the main purpose of heat warning systems is identification of days with an increased risk of heat-related death rather than prediction of exact magnitude of the excess mortality, WBGT seemed to be a slightly favorable predictor for such a system.

Application of spatial synoptic classification in evaluating links between heat stress and cardiovascular mortality and morbidity in Prague, Czech Republic.

Spatial synoptic classification (SSC) is here first employed in assessing heat-related mortality and morbidity in Central Europe. It is applied for examining links between weather patterns and cardiovascular (CVD) mortality and morbidity in an extended summer season (16 May-15 September) during 1994-2009. As in previous studies, two SSC air masses (AMs)-dry tropical (DT) and moist tropical (MT)-are associated with significant excess CVD mortality in Prague, while effects on CVD hospital admissions are small and insignificant. Excess mortality for ischaemic heart diseases is more strongly associated with DT, while MT has adverse effect especially on cerebrovascular mortality. Links between the oppressive AMs and excess mortality relate also to conditions on previous days, as DT and MT occur in typical sequences. The highest CVD mortality deviations are found 1 day after a hot spell's onset, when temperature as well as frequency of the oppressive AMs are highest. Following this peak is typically DT- to MT-like weather transition, characterized by decrease in temperature and increase in humidity. The transition between upward (DT) and downward (MT) phases is associated with the largest excess CVD mortality, and the change contributes to the increased and more lagged effects on cerebrovascular mortality. The study highlights the importance of critically evaluating SSC's applicability and benefits within warning systems relative to other synoptic and epidemiological approaches. Only a subset of days with the oppressive AMs is associated with excess mortality, and regression models accounting for possible meteorological and other factors explain little of the mortality variance.

Biometeorology for cities.

Improvements in global sustainability, health, and equity will largely be determined by the extent to which cities are able to become more efficient, hospitable, and productive places. The development and evolution of urban areas has a significant impact on local and regional weather and climate, which subsequently affect people and other organisms that live in and near cities. Biometeorologists, researchers who study the impact of weather and climate on living creatures, are well positioned to help evaluate and anticipate the consequences of urbanization on the biosphere. Motivated by the 60th anniversary of the International Society of Biometeorology, we reviewed articles published in the Society's International Journal of Biometeorology over the period 1974-2017 to understand if and how biometeorologists have directed attention to urban areas. We found that interest in urban areas has rapidly accelerated; urban-oriented articles accounted for more than 20% of all articles published in the journal in the most recent decade. Urban-focused articles in the journal span five themes: measuring urban climate, theoretical foundations and models, human thermal comfort, human morbidity and mortality, and ecosystem impacts. Within these themes, articles published in the journal represent a sizeable share of the total academic literature. More explicit attention from urban biometeorologists publishing in the journal to low- and middle-income countries, indoor environments, animals, and the impacts of climate change on human health would help ensure that the distinctive perspectives of biometeorology reach the places, people, and processes that are the foci of global sustainability, health, and equity goals.

Heat- and cold-stress effects on cardiovascular mortality and morbidity among urban and rural populations in the Czech Republic.

Several studies have examined the relationship of high and low air temperatures to cardiovascular mortality in the Czech Republic. Much less is understood about heat-/cold-related cardiovascular morbidity and possible regional differences. This paper compares the effects of warm and cold days on excess mortality and morbidity for cardiovascular diseases (CVDs) in the city of Prague and a rural region of southern Bohemia during 1994-2009. Population size and age structure are similar in the two regions. The results are evaluated for selected population groups (men and women). Excess mortality (number of deaths) and morbidity (number of hospital admissions) were determined as differences between observed and expected daily values, the latter being adjusted for long-term changes, annual and weekly cycles, and epidemics of influenza/acute respiratory infections. Generally higher relative excess CVD mortality on warm days than on cold days was identified in both regions. In contrast to mortality, weak excess CVD morbidity was observed for both warm and cold days. Different responses of individual CVDs to heat versus cold stress may be caused by the different nature of each CVD and different physiological processes induced by heat or cold stress. The slight differences between Prague and southern Bohemia in response to heat versus cold stress suggest the possible influence of environmental and socioeconomic factors such as the effects of urban heat island and exposure to air pollution, lifestyle differences, and divergence in population structure, which may result in differing vulnerability of urban versus rural population to temperature extremes.

Global, regional, and national burden of mortality associated with cold spells during 2000-19: a three-stage modelling study.

BACKGROUND: Exposure to cold spells is associated with mortality. However, little is known about the global mortality burden of cold spells. METHODS: A three-stage meta-analytical method was used to estimate the global mortality burden associated with cold spells by means of a time series dataset of 1960 locations across 59 countries (or regions). First, we fitted the location-specific, cold spell-related mortality associations using a quasi-Poisson regression with a distributed lag non-linear model with a lag period of up to 21 days. Second, we built a multivariate meta-regression model between location-specific associations and seven predictors. Finally, we predicted the global grid-specific cold spell-related mortality associations during 2000-19 using the fitted meta-regression model and the yearly grid-specific meta-predictors. We calculated the annual excess deaths, excess death ratio (excess deaths per 1000 deaths), and excess death rate (excess deaths per 100 000 population) due to cold spells for each grid across the world. FINDINGS: Globally, 205 932 (95% empirical CI [eCI] 162 692-250 337) excess deaths, representing 3·81 (95% eCI 2·93-4·71) excess deaths per 1000 deaths (excess death ratio), and 3·03 (2·33-3·75) excess deaths per 100 000 population (excess death rate) were associated with cold spells per year between 2000 and 2019. The annual average global excess death ratio in 2016-19 increased by 0·12 percentage points and the excess death rate in 2016-19 increased by 0·18 percentage points, compared with those in 2000-03. The mortality burden varied geographically. The excess death ratio and rate were highest in Europe, whereas these indicators were lowest in Africa. Temperate climates had higher excess death ratio and rate associated with cold spells than other climate zones. INTERPRETATION: Cold spells are associated with substantial mortality burden around the world with geographically varying patterns. Although the number of cold spells has on average been decreasing since year 2000, the public health threat of cold spells remains substantial. The findings indicate an urgency of taking local and regional measures to protect the public from the mortality burdens of cold spells. FUNDING: Australian Research Council, Australian National Health and Medical Research Council, EU's Horizon 2020 Project Exhaustion.

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.

All-cause, cardiovascular, and respiratory mortality and wildfire-related ozone: a multicountry two-stage time series analysis.

BACKGROUND: Wildfire activity is an important source of tropospheric ozone (O3) pollution. However, no study to date has systematically examined the associations of wildfire-related O3 exposure with mortality globally. METHODS: We did a multicountry two-stage time series analysis. From the Multi-City Multi-Country (MCC) Collaborative Research Network, data on daily all-cause, cardiovascular, and respiratory deaths were obtained from 749 locations in 43 countries or areas, representing overlapping periods from Jan 1, 2000, to Dec 31, 2016. We estimated the daily concentration of wildfire-related O3 in study locations using a chemical transport model, and then calibrated and downscaled O3 estimates to a resolution of 0·25°â€ˆ× 0·25° (approximately 28 km2 at the equator). Using a random-effects meta-analysis, we examined the associations of short-term wildfire-related O3 exposure (lag period of 0-2 days) with daily mortality, first at the location level and then pooled at the country, regional, and global levels. Annual excess mortality fraction in each location attributable to wildfire-related O3 was calculated with pooled effect estimates and used to obtain excess mortality fractions at country, regional, and global levels. FINDINGS: Between 2000 and 2016, the highest maximum daily wildfire-related O3 concentrations (≥30 µg/m3) were observed in locations in South America, central America, and southeastern Asia, and the country of South Africa. Across all locations, an increase of 1 µg/m3 in the mean daily concentration of wildfire-related O3 during lag 0-2 days was associated with increases of 0·55% (95% CI 0·29 to 0·80) in daily all-cause mortality, 0·44% (-0·10 to 0·99) in daily cardiovascular mortality, and 0·82% (0·18 to 1·47) in daily respiratory mortality. The associations of daily mortality rates with wildfire-related O3 exposure showed substantial geographical heterogeneity at the country and regional levels. Across all locations, estimated annual excess mortality fractions of 0·58% (95% CI 0·31 to 0·85; 31 606 deaths [95% CI 17 038 to 46 027]) for all-cause mortality, 0·41% (-0·10 to 0·91; 5249 [-1244 to 11 620]) for cardiovascular mortality, and 0·86% (0·18 to 1·51; 4657 [999 to 8206]) for respiratory mortality were attributable to short-term exposure to wildfire-related O3. INTERPRETATION: In this study, we observed an increase in all-cause and respiratory mortality associated with short-term wildfire-related O3 exposure. Effective risk and smoke management strategies should be implemented to protect the public from the impacts of wildfires. FUNDING: Australian Research Council and the Australian National Health and Medical Research Council.

Regional variation in the role of humidity on city-level heat-related mortality.

The rising humid heat is regarded as a severe threat to human survivability, but the proper integration of humid heat into heat-health alerts is still being explored. Using state-of-the-art epidemiological and climatological datasets, we examined the association between multiple heat stress indicators (HSIs) and daily human mortality in 739 cities worldwide. Notable differences were observed in the long-term trends and timing of heat events detected by HSIs. Air temperature (Tair) predicts heat-related mortality well in cities with a robust negative Tair-relative humidity correlation (CT-RH). However, in cities with near-zero or weak positive CT-RH, HSIs considering humidity provide enhanced predictive power compared to Tair. Furthermore, the magnitude and timing of heat-related mortality measured by HSIs could differ largely from those associated with Tair in many cities. Our findings provide important insights into specific regions where humans are vulnerable to humid heat and can facilitate the further enhancement of heat-health alert systems.

Temporal variations in the short-term effects of ambient air pollution on cardiovascular and respiratory mortality: a pooled analysis of 380 urban areas over a 22-year period.

BACKGROUND: Ambient air pollution, including particulate matter (such as PM10 and PM2·5) and nitrogen dioxide (NO2), has been linked to increases in mortality. Whether populations' vulnerability to these pollutants has changed over time is unclear, and studies on this topic do not include multicountry analysis. We evaluated whether changes in exposure to air pollutants were associated with changes in mortality effect estimates over time. METHODS: We extracted cause-specific mortality and air pollution data collected between 1995 and 2016 from the Multi-Country Multi-City (MCC) Collaborative Research Network database. We applied a two-stage approach to analyse the short-term effects of NO2, PM10, and PM2·5 on cause-specific mortality using city-specific time series regression analyses and multilevel random-effects meta-analysis. We assessed changes over time using a longitudinal meta-regression with time as a linear fixed term and explored potential sources of heterogeneity and two-pollutant models. FINDINGS: Over 21·6 million cardiovascular and 7·7 million respiratory deaths in 380 cities across 24 countries over the study period were included in the analysis. All three air pollutants showed decreasing concentrations over time. The pooled results suggested no significant temporal change in the effect estimates per unit exposure of PM10, PM2·5, or NO2 and mortality. However, the risk of cardiovascular mortality increased from 0·37% (95% CI -0·05 to 0·80) in 1998 to 0·85% (0·55 to 1·16) in 2012 with a 10 µg/m3 increase in PM2·5. Two-pollutant models generally showed similar results to single-pollutant models for PM fractions and indicated temporal differences for NO2. INTERPRETATION: Although air pollution levels decreased during the study period, the effect sizes per unit increase in air pollution concentration have not changed. This observation might be due to the composition, toxicity, and sources of air pollution, as well as other factors, such as socioeconomic determinants or changes in population distribution and susceptibility. FUNDING: None.

Seasonality of mortality under climate change: a multicountry projection study.

BACKGROUND: Climate change can directly impact temperature-related excess deaths and might subsequently change the seasonal variation in mortality. In this study, we aimed to provide a systematic and comprehensive assessment of potential future changes in the seasonal variation, or seasonality, of mortality across different climate zones. METHODS: In this modelling study, we collected daily time series of mean temperature and mortality (all causes or non-external causes only) via the Multi-Country Multi-City Collaborative (MCC) Research Network. These data were collected during overlapping periods, spanning from Jan 1, 1969 to Dec 31, 2020. We projected daily mortality from Jan 1, 2000 to Dec 31, 2099, under four climate change scenarios corresponding to increasing emissions (Shared Socioeconomic Pathways [SSP] scenarios SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5). We compared the seasonality in projected mortality between decades by its shape, timings (the day-of-year) of minimum (trough) and maximum (peak) mortality, and sizes (peak-to-trough ratio and attributable fraction). Attributable fraction was used to measure the burden of seasonality of mortality. The results were summarised by climate zones. FINDINGS: The MCC dataset included 126 809 537 deaths from 707 locations within 43 countries or areas. After excluding the only two polar locations (both high-altitude locations in Peru) from climatic zone assessments, we analysed 126 766 164 deaths in 705 locations aggregated in four climate zones (tropical, arid, temperate, and continental). From the 2000s to the 2090s, our projections showed an increase in mortality during the warm seasons and a decrease in mortality during the cold seasons, albeit with mortality remaining high during the cold seasons, under all four SSP scenarios in the arid, temperate, and continental zones. The magnitude of this changing pattern was more pronounced under the high-emission scenarios (SSP3-7.0 and SSP5-8.5), substantially altering the shape of seasonality of mortality and, under the highest emission scenario (SSP5-8.5), shifting the mortality peak from cold seasons to warm seasons in arid, temperate, and continental zones, and increasing the size of seasonality in all zones except the arid zone by the end of the century. In the 2090s compared with the 2000s, the change in peak-to-trough ratio (relative scale) ranged from 0·96 to 1·11, and the change in attributable fraction ranged from 0·002% to 0·06% under the SSP5-8.5 (highest emission) scenario. INTERPRETATION: A warming climate can substantially change the seasonality of mortality in the future. Our projections suggest that health-care systems should consider preparing for a potentially increased demand during warm seasons and sustained high demand during cold seasons, particularly in regions characterised by arid, temperate, and continental climates. FUNDING: The Environment Research and Technology Development Fund of the Environmental Restoration and Conservation Agency, provided by the Ministry of the Environment of Japan.

Association of holidays and the day of the week with suicide risk: multicounty, two stage, time series study.

OBJECTIVES: To assess the short term temporal variations in suicide risk related to the day of the week and national holidays in multiple countries. DESIGN: Multicountry, two stage, time series design. SETTING: Data from 740 locations in 26 countries and territories, with overlapping periods between 1971 and 2019, collected from the Multi-city Multi-country Collaborative Research Network database. PARTICIPANTS: All suicides were registered in these locations during the study period (overall 1 701 286 cases). MAIN OUTCOME MEASURES: Daily suicide mortality. RESULTS: Mondays had peak suicide risk during weekdays (Monday-Friday) across all countries, with relative risks (reference: Wednesday) ranging from 1.02 (95% confidence interval (CI) 0.95 to 1.10) in Costa Rica to 1.17 (1.09 to 1.25) in Chile. Suicide risks were lowest on Saturdays or Sundays in many countries in North America, Asia, and Europe. However, the risk increased during weekends in South and Central American countries, Finland, and South Africa. Additionally, evidence suggested strong increases in suicide risk on New Year's day in most countries with relative risks ranging from 0.93 (95% CI 0.75 to 1.14) in Japan to 1.93 (1.31 to 2.85) in Chile, whereas the evidence on Christmas day was weak. Suicide risk was associated with a weak decrease on other national holidays, except for Central and South American countries, where the risk generally increased one or two days after these holidays. CONCLUSIONS: Suicide risk was highest on Mondays and increased on New Year's day in most countries. However, the risk of suicide on weekends and Christmas varied by country and territory. The results of this study can help to better understand the short term variations in suicide risks and define suicide prevention action plans and awareness campaigns.

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

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

Ozone-related acute excess mortality projected to increase in the absence of climate and air quality controls consistent with the Paris Agreement.

Short-term exposure to ground-level ozone in cities is associated with increased mortality and is expected to worsen with climate and emission changes. However, no study has yet comprehensively assessed future ozone-related acute mortality across diverse geographic areas, various climate scenarios, and using CMIP6 multi-model ensembles, limiting our knowledge on future changes in global ozone-related acute mortality and our ability to design targeted health policies. Here, we combine CMIP6 simulations and epidemiological data from 406 cities in 20 countries or regions. We find that ozone-related deaths in 406 cities will increase by 45 to 6,200 deaths/year between 2010 and 2014 and between 2050 and 2054, with attributable fractions increasing in all climate scenarios (from 0.17% to 0.22% total deaths), except the single scenario consistent with the Paris Climate Agreement (declines from 0.17% to 0.15% total deaths). These findings stress the need for more stringent air quality regulations, as current standards in many countries are inadequate.

Comparison for the effects of different components of temperature variability on mortality: A multi-country time-series study.

BACKGROUND: Temperature variability (TV) is associated with increased mortality risk. However, it is still unknown whether intra-day or inter-day TV has different effects. OBJECTIVES: We aimed to assess the association of intra-day TV and inter-day TV with all-cause, cardiovascular, and respiratory mortality. METHODS: We collected data on total, cardiovascular, and respiratory mortality and meteorology from 758 locations in 47 countries or regions from 1972 to 2020. We defined inter-day TV as the standard deviation (SD) of daily mean temperatures across the lag interval, and intra-day TV as the average SD of minimum and maximum temperatures on each day. In the first stage, inter-day and intra-day TVs were modelled simultaneously in the quasi-Poisson time-series model for each location. In the second stage, a multi-level analysis was used to pool the location-specific estimates. RESULTS: Overall, the mortality risk due to each interquartile range [IQR] increase was higher for intra-day TV than for inter-day TV. The risk increased by 0.59% (95% confidence interval [CI]: 0.53, 0.65) for all-cause mortality, 0.64% (95% CI: 0.56, 0.73) for cardiovascular mortality, and 0.65% (95% CI: 0.49, 0.80) for respiratory mortality per IQR increase in intra-day TV0-7 (0.9 °C). An IQR increase in inter-day TV0-7 (1.6 °C) was associated with 0.22% (95% CI: 0.18, 0.26) increase in all-cause mortality, 0.44% (95% CI: 0.37, 0.50) increase in cardiovascular mortality, and 0.31% (95% CI: 0.21, 0.41) increase in respiratory mortality. The proportion of all-cause deaths attributable to intra-day TV0-7 and inter-day TV0-7 was 1.45% and 0.35%, respectively. The mortality risks varied by lag interval, climate area, season, and climate type. CONCLUSIONS: Our results indicated that intra-day TV may explain the main part of the mortality risk related to TV and suggested that comprehensive evaluations should be proposed in more countries to help protect human health.

Optimal heat stress metric for modelling heat-related mortality varies from country to country.

Combined heat and humidity is frequently described as the main driver of human heat-related mortality, more so than dry-bulb temperature alone. While based on physiological thinking, this assumption has not been robustly supported by epidemiological evidence. By performing the first systematic comparison of eight heat stress metrics (i.e., temperature combined with humidity and other climate variables) with warm-season mortality, in 604 locations over 39 countries, we find that the optimal metric for modelling mortality varies from country to country. Temperature metrics with no or little humidity modification associates best with mortality in ~40% of the studied countries. Apparent temperature (combined temperature, humidity and wind speed) dominates in another 40% of countries. There is no obvious climate grouping in these results. We recommend, where possible, that researchers use the optimal metric for each country. However, dry-bulb temperature performs similarly to humidity-based heat stress metrics in estimating heat-related mortality in present-day climate.

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

BACKGROUND: Heat and cold are established environmental risk factors for human health. However, mapping the related health burden is a difficult task due to the complexity of the associations and the differences in vulnerability and demographic distributions. In this study, we did a comprehensive mortality impact assessment due to heat and cold in European urban areas, considering geographical differences and age-specific risks. METHODS: We included urban areas across Europe between Jan 1, 2000, and Dec 12, 2019, using the Urban Audit dataset of Eurostat and adults aged 20 years and older living in these areas. Data were extracted from Eurostat, the Multi-country Multi-city Collaborative Research Network, Moderate Resolution Imaging Spectroradiometer, and Copernicus. We applied a three-stage method to estimate risks of temperature continuously across the age and space dimensions, identifying patterns of vulnerability on the basis of city-specific characteristics and demographic structures. These risks were used to derive minimum mortality temperatures and related percentiles and raw and standardised excess mortality rates for heat and cold aggregated at various geographical levels. FINDINGS: Across the 854 urban areas in Europe, we estimated an annual excess of 203 620 (empirical 95% CI 180 882-224 613) deaths attributed to cold and 20 173 (17 261-22 934) attributed to heat. These corresponded to age-standardised rates of 129 (empirical 95% CI 114-142) and 13 (11-14) deaths per 100 000 person-years. Results differed across Europe and age groups, with the highest effects in eastern European cities for both cold and heat. INTERPRETATION: Maps of mortality risks and excess deaths indicate geographical differences, such as a north-south gradient and increased vulnerability in eastern Europe, as well as local variations due to urban characteristics. The modelling framework and results are crucial for the design of national and local health and climate policies and for projecting the effects of cold and heat under future climatic and socioeconomic scenarios. FUNDING: Medical Research Council of UK, the Natural Environment Research Council UK, the EU's Horizon 2020, and the EU's Joint Research Center.