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.

Ambient fine particulate matter and daily mortality: a comparative analysis of observed and estimated exposure in 347 cities

BACKGROUND

Model-estimated air pollution exposure products have been widely used in epidemiological studies to assess the health risks of particulate matter with diameters of ≤2.5 µm (PM2.5). However, few studies have assessed the disparities in health effects between model-estimated and station-observed PM2.5 exposures.

METHODS

We collected daily all-cause, respiratory and cardiovascular mortality data in 347 cities across 15 countries and regions worldwide based on the Multi-City Multi-Country collaborative research network. The station-observed PM2.5 data were obtained from official monitoring stations. The model-estimated global PM2.5 product was developed using a machine-learning approach. The associations between daily exposure to PM2.5 and mortality were evaluated using a two-stage analytical approach.

RESULTS

We included 15.8 million all-cause, 1.5 million respiratory and 4.5 million cardiovascular deaths from 2000 to 2018. Short-term exposure to PM2.5 was associated with a relative risk increase (RRI) of mortality from both station-observed and model-estimated exposures. Every 10-μg/m3 increase in the 2-day moving average PM2.5 was associated with overall RRIs of 0.67% (95% CI: 0.49 to 0.85), 0.68% (95% CI: -0.03 to 1.39) and 0.45% (95% CI: 0.08 to 0.82) for all-cause, respiratory, and cardiovascular mortality based on station-observed PM2.5 and RRIs of 0.87% (95% CI: 0.68 to 1.06), 0.81% (95% CI: 0.08 to 1.55) and 0.71% (95% CI: 0.32 to 1.09) based on model-estimated exposure, respectively.

CONCLUSIONS

Mortality risks associated with daily PM2.5 exposure were consistent for both station-observed and model-estimated exposures, suggesting the reliability and potential applicability of the global PM2.5 product in epidemiological studies.

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.

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.

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.

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.

Joint effect of heat and air pollution on mortality in 620 cities of 36 countries

BACKGROUND

The epidemiological evidence on the interaction between heat and ambient air pollution on mortality is still inconsistent.

OBJECTIVES

To investigate the interaction between heat and ambient air pollution on daily mortality in a large dataset of 620 cities from 36 countries.

METHODS

We used daily data on all-cause mortality, air temperature, particulate matter ≤ 10 μm (PM10), PM ≤ 2.5 μm (PM2.5), nitrogen dioxide (NO2), and ozone (O3) from 620 cities in 36 countries in the period 1995-2020. We restricted the analysis to the six consecutive warmest months in each city. City-specific data were analysed with over-dispersed Poisson regression models, followed by a multilevel random-effects meta-analysis. The joint association between air temperature and air pollutants was modelled with product terms between non-linear functions for air temperature and linear functions for air pollutants.

RESULTS

We analyzed 22,630,598 deaths. An increase in mean temperature from the 75th to the 99th percentile of city-specific distributions was associated with an average 8.9 % (95 % confidence interval: 7.1 %, 10.7 %) mortality increment, ranging between 5.3 % (3.8 %, 6.9 %) and 12.8 % (8.7 %, 17.0 %), when daily PM10 was equal to 10 or 90 μg/m3, respectively. Corresponding estimates when daily O3 concentrations were 40 or 160 μg/m3 were 2.9 % (1.1 %, 4.7 %) and 12.5 % (6.9 %, 18.5 %), respectively. Similarly, a 10 μg/m3 increment in PM10 was associated with a 0.54 % (0.10 %, 0.98 %) and 1.21 % (0.69 %, 1.72 %) increase in mortality when daily air temperature was set to the 1st and 99th city-specific percentiles, respectively. Corresponding mortality estimate for O3 across these temperature percentiles were 0.00 % (-0.44 %, 0.44 %) and 0.53 % (0.38 %, 0.68 %). Similar effect modification results, although slightly weaker, were found for PM2.5 and NO2.

CONCLUSIONS

Suggestive evidence of effect modification between air temperature and air pollutants on mortality during the warm period was found in a global dataset of 620 cities.

Mortality risks associated with floods in 761 communities worldwide: time series study

OBJECTIVE

To evaluate lag-response associations and effect modifications of exposure to floods with risks of all cause, cardiovascular, and respiratory mortality on a global scale.

DESIGN

Time series study.

SETTING

761 communities in 35 countries or territories with at least one flood event during the study period.

PARTICIPANTS

Multi-Country Multi-City Collaborative Research Network database, Australian Cause of Death Unit Record File, New Zealand Integrated Data Infrastructure, and the International Network for the Demographic Evaluation of Populations and their Health Network database.

MAIN OUTCOME MEASURES

The main outcome was daily counts of deaths. An estimation for the lag-response association between flood and daily mortality risk was modelled, and the relative risks over the lag period were cumulated to calculate overall effects. Attributable fractions of mortality due to floods were further calculated. A quasi-Poisson model with a distributed lag non-linear function was used to examine how daily death risk was associated with flooded days in each community, and then the community specific associations were pooled using random effects multivariate meta-analyses. Flooded days were defined as days from the start date to the end date of flood events.

RESULTS

A total of 47.6 million all cause deaths, 11.1 million cardiovascular deaths, and 4.9 million respiratory deaths were analysed. Over the 761 communities, mortality risks increased and persisted for up to 60 days (50 days for cardiovascular mortality) after a flooded day. The cumulative relative risks for all cause, cardiovascular, and respiratory mortality were 1.021 (95% confidence interval 1.006 to 1.036), 1.026 (1.005 to 1.047), and 1.049 (1.008 to 1.092), respectively. The associations varied across countries or territories and regions. The flood-mortality associations appeared to be modified by climate type and were stronger in low income countries and in populations with a low human development index or high proportion of older people. In communities impacted by flood, up to 0.10% of all cause deaths, 0.18% of cardiovascular deaths, and 0.41% of respiratory deaths were attributed to floods.

CONCLUSIONS

This study found that the risks of all cause, cardiovascular, and respiratory mortality increased for up to 60 days after exposure to flood and the associations could vary by local climate type, socioeconomic status, and older age.

Rapid increase in the risk of heat-related mortality

Heat-related mortality has been identified as one of the key climate extremes posing a risk to human health. Current research focuses largely on how heat mortality increases with mean global temperature rise, but it is unclear how much climate change will increase the frequency and severity of extreme summer seasons with high impact on human health. In this probabilistic analysis, we combined empirical heat-mortality relationships for 748 locations from 47 countries with climate model large ensemble data to identify probable past and future highly impactful summer seasons. Across most locations, heat mortality counts of a 1-in-100 year season in the climate of 2000 would be expected once every ten to twenty years in the climate of 2020. These return periods are projected to further shorten under warming levels of 1.5 °C and 2 °C, where heat-mortality extremes of the past climate will eventually become commonplace if no adaptation occurs. Our findings highlight the urgent need for strong mitigation and adaptation to reduce impacts on human lives.

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.

Heat-related cardiorespiratory mortality: Effect modification by air pollution across 482 cities from 24 countries

BACKGROUND

Evidence on the potential interactive effects of heat and ambient air pollution on cause-specific mortality is inconclusive and limited to selected locations.

OBJECTIVES

We investigated the effects of heat on cardiovascular and respiratory mortality and its modification by air pollution during summer months (six consecutive hottest months) in 482 locations across 24 countries.

METHODS

Location-specific daily death counts and exposure data (e.g., particulate matter with diameters ≤ 2.5 µm [PM2.5]) were obtained from 2000 to 2018. We used location-specific confounder-adjusted Quasi-Poisson regression with a tensor product between air temperature and the air pollutant. We extracted heat effects at low, medium, and high levels of pollutants, defined as the 5th, 50th, and 95th percentile of the location-specific pollutant concentrations. Country-specific and overall estimates were derived using a random-effects multilevel meta-analytical model.

RESULTS

Heat was associated with increased cardiorespiratory mortality. Moreover, the heat effects were modified by elevated levels of all air pollutants in most locations, with stronger effects for respiratory than cardiovascular mortality. For example, the percent increase in respiratory mortality per increase in the 2-day average summer temperature from the 75th to the 99th percentile was 7.7% (95% Confidence Interval [CI] 7.6-7.7), 11.3% (95%CI 11.2-11.3), and 14.3% (95% CI 14.1-14.5) at low, medium, and high levels of PM2.5, respectively. Similarly, cardiovascular mortality increased by 1.6 (95%CI 1.5-1.6), 5.1 (95%CI 5.1-5.2), and 8.7 (95%CI 8.7-8.8) at low, medium, and high levels of O3, respectively.

DISCUSSION

We observed considerable modification of the heat effects on cardiovascular and respiratory mortality by elevated levels of air pollutants. Therefore, mitigation measures following the new WHO Air Quality Guidelines are crucial to enhance better health and promote sustainable development.

Effect modification of greenness on the association between heat and mortality: A multi-city multi-country study

BACKGROUND

Identifying how greenspace impacts the temperature-mortality relationship in urban environments is crucial, especially given climate change and rapid urbanization. However, the effect modification of greenspace on heat-related mortality has been typically focused on a localized area or single country. This study examined the heat-mortality relationship among different greenspace levels in a global setting.

METHODS

We collected daily ambient temperature and mortality data for 452 locations in 24 countries and used Enhanced Vegetation Index (EVI) as the greenspace measurement. We used distributed lag non-linear model to estimate the heat-mortality relationship in each city and the estimates were pooled adjusting for city-specific average temperature, city-specific temperature range, city-specific population density, and gross domestic product (GDP). The effect modification of greenspace was evaluated by comparing the heat-related mortality risk for different greenspace groups (low, medium, and high), which were divided into terciles among 452 locations.

FINDINGS

Cities with high greenspace value had the lowest heat-mortality relative risk of 1·19 (95% CI: 1·13, 1·25), while the heat-related relative risk was 1·46 (95% CI: 1·31, 1·62) for cities with low greenspace when comparing the 99^th temperature and the minimum mortality temperature. A 20% increase of greenspace is associated with a 9·02% (95% CI: 8·88, 9·16) decrease in the heat-related attributable fraction, and if this association is causal (which is not within the scope of this study to assess), such a reduction could save approximately 933 excess deaths per year in 24 countries.

INTERPRETATION

Our findings can inform communities on the potential health benefits of greenspaces in the urban environment and mitigation measures regarding the impacts of climate change.

FUNDING

This publication was developed under Assistance Agreement No. RD83587101 awarded by the U.S. Environmental Protection Agency to Yale University. It has not been formally reviewed by EPA. The views expressed in this document are solely those of the authors and do not necessarily reflect those of the Agency. EPA does not endorse any products or commercial services mentioned in this publication. Research reported in this publication was also supported by the National Institute on Minority Health and Health Disparities of the National Institutes of Health under Award Number R01MD012769. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Also, this work has been supported by the National Research Foundation of Korea (2021R1A6A3A03038675), Medical Research Council-UK (MR/V034162/1 and MR/R013349/1), Natural Environment Research Council UK (Grant ID: NE/R009384/1), Academy of Finland (Grant ID: 310372), European Union's Horizon 2020 Project Exhaustion (Grant ID: 820655 and 874990), Czech Science Foundation (22-24920S), Emory University's NIEHS-funded HERCULES Center (Grant ID: P30ES019776), and Grant CEX2018-000794-S funded by MCIN/AEI/ 10.13039/501100011033 The funders had no role in the design, data collection, analysis, interpretation of results, manuscript writing, or decision to publication.

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

BACKGROUND

Increased mortality risk is associated with short-term temperature variability. However, to our knowledge, there has been no comprehensive assessment of the temperature variability-related mortality burden worldwide. In this study, using data from the MCC Collaborative Research Network, we first explored the association between temperature variability and mortality across 43 countries or regions. Then, to provide a more comprehensive picture of the global burden of mortality associated with temperature variability, global gridded temperature data with a resolution of 0·5° × 0·5° were used to assess the temperature variability-related mortality burden at the global, regional, and national levels. Furthermore, temporal trends in temperature variability-related mortality burden were also explored from 2000-19.

METHODS

In this modelling study, we applied a three-stage meta-analytical approach to assess the global temperature variability-related mortality burden at a spatial resolution of 0·5° × 0·5° from 2000-19. Temperature variability was calculated as the SD of the average of the same and previous days' minimum and maximum temperatures. We first obtained location-specific temperature variability related-mortality associations based on a daily time series of 750 locations from the Multi-country Multi-city Collaborative Research Network. We subsequently constructed a multivariable meta-regression model with five predictors to estimate grid-specific temperature variability related-mortality associations across the globe. Finally, percentage excess in mortality and excess mortality rate were calculated to quantify the temperature variability-related mortality burden and to further explore its temporal trend over two decades.

FINDINGS

An increasing trend in temperature variability was identified at the global level from 2000 to 2019. Globally, 1 753 392 deaths (95% CI 1 159 901-2 357 718) were associated with temperature variability per year, accounting for 3·4% (2·2-4·6) of all deaths. Most of Asia, Australia, and New Zealand were observed to have a higher percentage excess in mortality than the global mean. Globally, the percentage excess in mortality increased by about 4·6% (3·7-5·3) per decade. The largest increase occurred in Australia and New Zealand (7·3%, 95% CI 4·3-10·4), followed by Europe (4·4%, 2·2-5·6) and Africa (3·3, 1·9-4·6).

INTERPRETATION

Globally, a substantial mortality burden was associated with temperature variability, showing geographical heterogeneity and a slightly increasing temporal trend. Our findings could assist in raising public awareness and improving the understanding of the health impacts of temperature variability.

FUNDING

Australian Research Council, Australian National Health & Medical Research Council.

Comparison of weather station and climate reanalysis data for modelling temperature-related mortality

Epidemiological analyses of health risks associated with non-optimal temperature are traditionally based on ground observations from weather stations that offer limited spatial and temporal coverage. Climate reanalysis represents an alternative option that provide complete spatio-temporal exposure coverage, and yet are to be systematically explored for their suitability in assessing temperature-related health risks at a global scale. Here we provide the first comprehensive analysis over multiple regions to assess the suitability of the most recent generation of reanalysis datasets for health impact assessments and evaluate their comparative performance against traditional station-based data. Our findings show that reanalysis temperature from the last ERA5 products generally compare well to station observations, with similar non-optimal temperature-related risk estimates. However, the analysis offers some indication of lower performance in tropical regions, with a likely underestimation of heat-related excess mortality. Reanalysis data represent a valid alternative source of exposure variables in epidemiological analyses of temperature-related risk.

Fluctuating temperature modifies heat-mortality association in the globe

Studies have investigated the effects of heat and temperature variability (TV) on mortality. However, few assessed whether TV modifies the heat-mortality association. Data on daily temperature and mortality in the warm season were collected from 717 locations across 36 countries. TV was calculated as the standard deviation of the average of the same and previous days’ minimum and maximum temperatures. We used location-specific quasi-Poisson regression models with an interaction term between the cross-basis term for mean temperature and quartiles of TV to obtain heat-mortality associations under each quartile of TV, and then pooled estimates at the country, regional, and global levels. Results show the increased risk in heat-related mortality with increments in TV, accounting for 0.70% (95% confidence interval [CI]: −0.33 to 1.69), 1.34% (95% CI: −0.14 to 2.73), 1.99% (95% CI: 0.29–3.57), and 2.73% (95% CI: 0.76–4.50) of total deaths for Q1–Q4 (first quartile–fourth quartile) of TV. The modification effects of TV varied geographically. Central Europe had the highest attributable fractions (AFs), corresponding to 7.68% (95% CI: 5.25–9.89) of total deaths for Q4 of TV, while the lowest AFs were observed in North America, with the values for Q4 of 1.74% (95% CI: −0.09 to 3.39). TV had a significant modification effect on the heat-mortality association, causing a higher heat-related mortality burden with increments of TV. Implementing targeted strategies against heat exposure and fluctuant temperatures simultaneously would benefit public health.

•

Increased temperature variability (TV) poses a greater mortality risk due to heat

•

TV has a more profound modification effect on extreme heat-mortality association

•

Strategies against heat and TV simultaneously would benefit public health

Mortality risk attributable to wildfire-related PM2·5 pollution: a global time series study in 749 locations

BACKGROUND

Many regions of the world are now facing more frequent and unprecedentedly large wildfires. However, the association between wildfire-related PM_2·5 and mortality has not been well characterised. We aimed to comprehensively assess the association between short-term exposure to wildfire-related PM_2·5 and mortality across various regions of the world.

METHODS

For this time series study, data on daily counts of deaths for all causes, cardiovascular causes, and respiratory causes were collected from 749 cities in 43 countries and regions during 2000-16. Daily concentrations of wildfire-related PM_2·5 were estimated using the three-dimensional chemical transport model GEOS-Chem at a 0·25° × 0·25° resolution. The association between wildfire-related PM_2·5 exposure and mortality was examined using a quasi-Poisson time series model in each city considering both the current-day and lag effects, and the effect estimates were then pooled using a random-effects meta-analysis. Based on these pooled effect estimates, the population attributable fraction and relative risk (RR) of annual mortality due to acute wildfire-related PM_2·5 exposure was calculated.

FINDINGS

65·6 million all-cause deaths, 15·1 million cardiovascular deaths, and 6·8 million respiratory deaths were included in our analyses. The pooled RRs of mortality associated with each 10 μg/m³ increase in the 3-day moving average (lag 0-2 days) of wildfire-related PM_2·5 exposure were 1·019 (95% CI 1·016-1·022) for all-cause mortality, 1·017 (1·012-1·021) for cardiovascular mortality, and 1·019 (1·013-1·025) for respiratory mortality. Overall, 0·62% (95% CI 0·48-0·75) of all-cause deaths, 0·55% (0·43-0·67) of cardiovascular deaths, and 0·64% (0·50-0·78) of respiratory deaths were annually attributable to the acute impacts of wildfire-related PM_2·5 exposure during the study period.

INTERPRETATION

Short-term exposure to wildfire-related PM_2·5 was associated with increased risk of mortality. Urgent action is needed to reduce health risks from the increasing wildfires.

FUNDING

Australian Research Council, Australian National Health & Medical Research Council.

Mortality risk attributable to wildfire-related PM2·5 pollution: a global time series study in 749 locations

BACKGROUND

Many regions of the world are now facing more frequent and unprecedentedly large wildfires. However, the association between wildfire-related PM2·5 and mortality has not been well characterised. We aimed to comprehensively assess the association between short-term exposure to wildfire-related PM2·5 and mortality across various regions of the world.

METHODS

For this time series study, data on daily counts of deaths for all causes, cardiovascular causes, and respiratory causes were collected from 749 cities in 43 countries and regions during 2000-16. Daily concentrations of wildfire-related PM2·5 were estimated using the three-dimensional chemical transport model GEOS-Chem at a 0·25° × 0·25° resolution. The association between wildfire-related PM2·5 exposure and mortality was examined using a quasi-Poisson time series model in each city considering both the current-day and lag effects, and the effect estimates were then pooled using a random-effects meta-analysis. Based on these pooled effect estimates, the population attributable fraction and relative risk (RR) of annual mortality due to acute wildfire-related PM2·5 exposure was calculated.

FINDINGS

65·6 million all-cause deaths, 15·1 million cardiovascular deaths, and 6·8 million respiratory deaths were included in our analyses. The pooled RRs of mortality associated with each 10 μg/m3 increase in the 3-day moving average (lag 0-2 days) of wildfire-related PM2·5 exposure were 1·019 (95% CI 1·016-1·022) for all-cause mortality, 1·017 (1·012-1·021) for cardiovascular mortality, and 1·019 (1·013-1·025) for respiratory mortality. Overall, 0·62% (95% CI 0·48-0·75) of all-cause deaths, 0·55% (0·43-0·67) of cardiovascular deaths, and 0·64% (0·50-0·78) of respiratory deaths were annually attributable to the acute impacts of wildfire-related PM2·5 exposure during the study period.

INTERPRETATION

Short-term exposure to wildfire-related PM2·5 was associated with increased risk of mortality. Urgent action is needed to reduce health risks from the increasing wildfires.

FUNDING

Australian Research Council, Australian National Health & Medical Research Council.

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

BACKGROUND

Exposure to cold or hot temperatures is associated with premature deaths. We aimed to evaluate the global, regional, and national mortality burden associated with non-optimal ambient temperatures.

METHODS

In this modelling study, we collected time-series data on mortality and ambient temperatures from 750 locations in 43 countries and five meta-predictors at a grid size of 0·5° × 0·5° across the globe. A three-stage analysis strategy was used. First, the temperature-mortality association was fitted for each location by use of a time-series regression. Second, a multivariate meta-regression model was built between location-specific estimates and meta-predictors. Finally, the grid-specific temperature-mortality association between 2000 and 2019 was predicted by use of the fitted meta-regression and the grid-specific meta-predictors. Excess deaths due to non-optimal temperatures, the ratio between annual excess deaths and all deaths of a year (the excess death ratio), and the death rate per 100 000 residents were then calculated for each grid across the world. Grids were divided according to regional groupings of the UN Statistics Division.

FINDINGS

Globally, 5 083 173 deaths (95% empirical CI [eCI] 4 087 967-5 965 520) were associated with non-optimal temperatures per year, accounting for 9·43% (95% eCI 7·58-11·07) of all deaths (8·52% [6·19-10·47] were cold-related and 0·91% [0·56-1·36] were heat-related). There were 74 temperature-related excess deaths per 100 000 residents (95% eCI 60-87). The mortality burden varied geographically. Of all excess deaths, 2 617 322 (51·49%) occurred in Asia. Eastern Europe had the highest heat-related excess death rate and Sub-Saharan Africa had the highest cold-related excess death rate. From 2000-03 to 2016-19, the global cold-related excess death ratio changed by -0·51 percentage points (95% eCI -0·61 to -0·42) and the global heat-related excess death ratio increased by 0·21 percentage points (0·13-0·31), leading to a net reduction in the overall ratio. The largest decline in overall excess death ratio occurred in South-eastern Asia, whereas excess death ratio fluctuated in Southern Asia and Europe.

INTERPRETATION

Non-optimal temperatures are associated with a substantial mortality burden, which varies spatiotemporally. Our findings will benefit international, national, and local communities in developing preparedness and prevention strategies to reduce weather-related impacts immediately and under climate change scenarios.

FUNDING

Australian Research Council and the Australian National Health and Medical Research Council.

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 burden of heat-related mortality attributable to recent human-induced climate change

Climate change affects human health; however, there have been no large-scale, systematic efforts to quantify the heat-related human health impacts that have already occurred due to climate change. Here, we use empirical data from 732 locations in 43 countries to estimate the mortality burdens associated with the additional heat exposure that has resulted from recent human-induced warming, during the period 1991-2018. Across all study countries, we find that 37.0% (range 20.5-76.3%) of warm-season heat-related deaths can be attributed to anthropogenic climate change and that increased mortality is evident on every continent. Burdens varied geographically but were of the order of dozens to hundreds of deaths per year in many locations. Our findings support the urgent need for more ambitious mitigation and adaptation strategies to minimize the public health impacts of climate change.

Short term associations of ambient nitrogen dioxide with daily total, cardiovascular, and respiratory mortality: multilocation analysis in 398 cities

OBJECTIVE

To evaluate the short term associations between nitrogen dioxide (NO₂) and total, cardiovascular, and respiratory mortality across multiple countries/regions worldwide, using a uniform analytical protocol.

DESIGN

Two stage, time series approach, with overdispersed generalised linear models and multilevel meta-analysis.

SETTING

398 cities in 22 low to high income countries/regions.

MAIN OUTCOME MEASURES

Daily deaths from total (62.8 million), cardiovascular (19.7 million), and respiratory (5.5 million) causes between 1973 and 2018.

RESULTS

On average, a 10 μg/m³ increase in NO₂ concentration on lag 1 day (previous day) was associated with 0.46% (95% confidence interval 0.36% to 0.57%), 0.37% (0.22% to 0.51%), and 0.47% (0.21% to 0.72%) increases in total, cardiovascular, and respiratory mortality, respectively. These associations remained robust after adjusting for co-pollutants (particulate matter with aerodynamic diameter ≤10 μm or ≤2.5 μm (PM₁₀ and PM_2.5, respectively), ozone, sulfur dioxide, and carbon monoxide). The pooled concentration-response curves for all three causes were almost linear without discernible thresholds. The proportion of deaths attributable to NO₂ concentration above the counterfactual zero level was 1.23% (95% confidence interval 0.96% to 1.51%) across the 398 cities.

CONCLUSIONS

This multilocation study provides key evidence on the independent and linear associations between short term exposure to NO₂ and increased risk of total, cardiovascular, and respiratory mortality, suggesting that health benefits would be achieved by tightening the guidelines and regulatory limits of NO₂.

Influence of air temperature and implemented veterinary measures on the incidence of human salmonellosis in the Czech Republic during 1998-2017

Background

The aim of our study was to analyse the influence of air temperature and implemented veterinary measures on salmonellosis incidence in the Czech Republic (CZ).

Methods

We conducted a descriptive analysis of salmonellosis as reported to the Czech national surveillance system during 1998–2017 and evaluated the influence of applied veterinary measures (started in January 2008) on salmonellosis incidence by comparing two 9-year periods (1998–2006, 2009–2017). Using a generalized additive model, we analysed association between monthly mean air temperature and log-transformed salmonellosis incidence over the entire twenty-year period.

Results

A total of 410,533 salmonellosis cases were reported during the study period in the CZ. Annual mean incidences of salmonellosis were 313.0/100,000 inhabitants before and 99.0/100,000 inhabitants after implementation of the veterinary measures. The time course of incidence was non-linear, with a sharp decline during 2006–2010. Significant association was found between disease incidence and air temperature. On average, the data indicated that within a common temperature range every 1 °C rise in air temperature contributed to a significant 6.2% increase in salmonellosis cases.

Conclusions

Significant non-linear effects of annual trend, within-year seasonality, and air temperature on the incidence of salmonellosis during 1998–2017 were found. Our study also demonstrates significant direct effect of preventive veterinary measures taken in poultry in reducing incidence of human salmonellosis in the CZ. The annual mean number of salmonellosis cases in the period after introducing the veterinary measures was only 32.5% of what it had been in the previous period.

Temporal changes in years of life lost associated with heat waves in the Czech Republic

Seniors constitute the population group generally most at risk of mortality due to heat stress. As life expectancy increases and health conditions of elderly people improve over time, vulnerability of the population to heat changes as well. We employed the years-of-life-lost (YLL) approach, considering life expectancy at the time of each death, to investigate how population ageing affects temporal changes in heat-related mortality in the Czech Republic. Using an updated gridded meteorological database, we identified heat waves during 1994-2017, and analysed temporal changes in their impacts on YLL and mortality. The mean impact of a heat-wave day on relative excess mortality and YLL had declined by approximately 2-3% per decade. That decline abated in the current decade, however, and the decreasing trend in mean excess mortality as well as YLL vanished when the short-term mortality displacement effect was considered. Moreover, the cumulative number of excess deaths and YLL during heat waves rose due to increasing frequency and intensity of heat waves during the examined period. The results show that in studies of temporal changes it is important to differentiate between mean effects of heat waves on mortality and the overall death burden associated with heat waves. Analysis of the average ratio of excess YLL/death per heat-wave day indicated that the major heat-vulnerable population group shifted towards older age (70+ years among males and 75+ years among females). Our findings highlight the importance of focusing heat-protection measures especially upon the elderly population, which is most heat-vulnerable and whose numbers are rising.

Short term association between ozone and mortality: global two stage time series study in 406 locations in 20 countries

OBJECTIVE

To assess short term mortality risks and excess mortality associated with exposure to ozone in several cities worldwide.

DESIGN

Two stage time series analysis.

SETTING

406 cities in 20 countries, with overlapping periods between 1985 and 2015, collected from the database of Multi-City Multi-Country Collaborative Research Network.

POPULATION

Deaths for all causes or for external causes only registered in each city within the study period. MAIN OUTCOME MEASURES: Daily total mortality (all or non-external causes only).

RESULTS

A total of 45 165 171 deaths were analysed in the 406 cities. On average, a 10 µg/m³ increase in ozone during the current and previous day was associated with an overall relative risk of mortality of 1.0018 (95% confidence interval 1.0012 to 1.0024). Some heterogeneity was found across countries, with estimates ranging from greater than 1.0020 in the United Kingdom, South Africa, Estonia, and Canada to less than 1.0008 in Mexico and Spain. Short term excess mortality in association with exposure to ozone higher than maximum background levels (70 µg/m³) was 0.26% (95% confidence interval 0.24% to 0.28%), corresponding to 8203 annual excess deaths (95% confidence interval 3525 to 12 840) across the 406 cities studied. The excess remained at 0.20% (0.18% to 0.22%) when restricting to days above the WHO guideline (100 µg/m³), corresponding to 6262 annual excess deaths (1413 to 11 065). Above more lenient thresholds for air quality standards in Europe, America, and China, excess mortality was 0.14%, 0.09%, and 0.05%, respectively.

CONCLUSIONS

Results suggest that ozone related mortality could be potentially reduced under stricter air quality standards. These findings have relevance for the implementation of efficient clean air interventions and mitigation strategies designed within national and international climate policies.

Predicted temperature-increase-induced global health burden and its regional variability

An increase in the global health burden of temperature was projected for 459 locations in 28 countries worldwide under four representative concentration pathway scenarios until 2099. We determined that the amount of temperature increase for each 100 ppm increase in global CO₂ concentrations is nearly constant, regardless of climate scenarios. The overall average temperature increase during 2010-2099 is largest in Canada (1.16 °C/100 ppm) and Finland (1.14 °C/100 ppm), while it is smallest in Ireland (0.62 °C/100 ppm) and Argentina (0.63 °C/100 ppm). In addition, for each 1 °C temperature increase, the amount of excess mortality is increased largely in tropical countries such as Vietnam (10.34%p/°C) and the Philippines (8.18%p/°C), while it is decreased in Ireland (-0.92%p/°C) and Australia (-0.32%p/°C). To understand the regional variability in temperature increase and mortality, we performed a regression-based modeling. We observed that the projected temperature increase is highly correlated with daily temperature range at the location and vulnerability to temperature increase is affected by health expenditure, and proportions of obese and elderly population.

The Role of Humidity in Associations of High Temperature with Mortality: A Multicountry, Multicity Study

BACKGROUND

There is strong experimental evidence that physiologic stress from high temperatures is greater if humidity is higher. However, heat indices developed to allow for this have not consistently predicted mortality better than dry-bulb temperature.

OBJECTIVES

We aimed to clarify the potential contribution of humidity an addition to temperature in predicting daily mortality in summer by using a large multicountry dataset.

METHODS

In 445 cities in 24 countries, we fit a time-series regression model for summer mortality with a distributed lag nonlinear model (DLNM) for temperature (up to lag 3) and supplemented this with a range of terms for relative humidity (RH) and its interaction with temperature. City-specific associations were summarized using meta-analytic techniques.

RESULTS

Adding a linear term for RH to the temperature term improved fit slightly, with an increase of 23% in RH (the 99th percentile anomaly) associated with a 1.1% [95% confidence interval (CI): 0.8, 1.3] decrease in mortality. Allowing curvature in the RH term or adding terms for interaction of RH with temperature did not improve the model fit. The humidity-related decreased risk was made up of a positive coefficient at lag 0 outweighed by negative coefficients at lags of 1-3 d. Key results were broadly robust to small model changes and replacing RH with absolute measures of humidity. Replacing temperature with apparent temperature, a metric combining humidity and temperature, reduced goodness of fit slightly.

DISCUSSION

The absence of a positive association of humidity with mortality in summer in this large multinational study is counter to expectations from physiologic studies, though consistent with previous epidemiologic studies finding little evidence for improved prediction by heat indices. The result that there was a small negative average association of humidity with mortality should be interpreted cautiously; the lag structure has unclear interpretation and suggests the need for future work to clarify. https://doi.org/10.1289/EHP5430.

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 (PM₁₀) and fine PM with an aerodynamic diameter of 2.5 μm or less (PM_2.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 PM₁₀ 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 PM_2.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 PM₁₀ and PM_2.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.).

The predictability of heat-related mortality in Prague, Czech Republic, during summer 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.

Temperature-related mortality impacts under and beyond Paris Agreement climate change scenarios

The Paris Agreement binds all nations to undertake ambitious efforts to combat climate change, with the commitment to Bhold warming well below 2 °C in global mean temperature (GMT), relative to pre-industrial levels, and to pursue efforts to limit warming to 1.5 °C". The 1.5 °C limit constitutes an ambitious goal for which greater evidence on its benefits for health would help guide policy and potentially increase the motivation for action. Here we contribute to this gap with an assessment on the potential health benefits, in terms of reductions in temperature-related mortality, derived from the compliance to the agreed temperature targets, compared to more extreme warming scenarios. We performed a multi-region analysis in 451 locations in 23 countries with different climate zones, and evaluated changes in heat and cold-related mortality under scenarios consistent with the Paris Agreement targets (1.5 and 2 °C) and more extreme GMT increases (3 and 4 °C), and under the assumption of no changes in demographic distribution and vulnerability. Our results suggest that limiting warming below 2 °C could prevent large increases in temperature-related mortality in most regions worldwide. The comparison between 1.5 and 2 °C is more complex and characterized by higher uncertainty, with geographical differences that indicate potential benefits limited to areas located in warmer climates, where direct climate change impacts will be more discernible.

Persistent multi-scale fluctuations shift European hydroclimate to its millennial boundaries

In recent years, there has been growing concern about the effect of global warming on water resources, especially at regional and continental scales. The last IPCC report on extremes states that there is medium confidence about an increase on European drought frequency during twentieth century. Here we use the Old World Drought Atlas palaeoclimatic reconstruction to show that when Europe’s hydroclimate is examined under a millennial, multi-scale perspective, a significant decrease in dryness can be observed since 1920 over most of central and northern Europe. On the contrary, in the south, drying conditions have prevailed, creating an intense north-to-south dipole. In both cases, hydroclimatic conditions have shifted to, and in some regions exceeded, their millennial boundaries, remaining at these extreme levels for the longest period of the 1000-year-long record.

In recent years, there has been an ongoing discussion about the hydroclimatic changes over Europe. Here, the authors show that since the beginning of the 20th century, hydroclimatic conditions have shifted to their millennial boundaries, remaining at these extreme levels for a period of unprecedented duration.

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.

Projections of temperature-related excess mortality under climate change scenarios

BACKGROUND

Climate change can directly affect human health by varying exposure to non-optimal outdoor temperature. However, evidence on this direct impact at a global scale is limited, mainly due to issues in modelling and projecting complex and highly heterogeneous epidemiological relationships across different populations and climates.

METHODS

We collected observed daily time series of mean temperature and mortality counts for all causes or non-external causes only, in periods ranging from Jan 1, 1984, to Dec 31, 2015, from various locations across the globe through the Multi-Country Multi-City Collaborative Research Network. We estimated temperature-mortality relationships through a two-stage time series design. We generated current and future daily mean temperature series under four scenarios of climate change, determined by varying trajectories of greenhouse gas emissions, using five general circulation models. We projected excess mortality for cold and heat and their net change in 1990-2099 under each scenario of climate change, assuming no adaptation or population changes.

FINDINGS

Our dataset comprised 451 locations in 23 countries across nine regions of the world, including 85 879 895 deaths. Results indicate, on average, a net increase in temperature-related excess mortality under high-emission scenarios, although with important geographical differences. In temperate areas such as northern Europe, east Asia, and Australia, the less intense warming and large decrease in cold-related excess would induce a null or marginally negative net effect, with the net change in 2090-99 compared with 2010-19 ranging from -1·2% (empirical 95% CI -3·6 to 1·4) in Australia to -0·1% (-2·1 to 1·6) in east Asia under the highest emission scenario, although the decreasing trends would reverse during the course of the century. Conversely, warmer regions, such as the central and southern parts of America or Europe, and especially southeast Asia, would experience a sharp surge in heat-related impacts and extremely large net increases, with the net change at the end of the century ranging from 3·0% (-3·0 to 9·3) in Central America to 12·7% (-4·7 to 28·1) in southeast Asia under the highest emission scenario. Most of the health effects directly due to temperature increase could be avoided under scenarios involving mitigation strategies to limit emissions and further warming of the planet.

INTERPRETATION

This study shows the negative health impacts of climate change that, under high-emission scenarios, would disproportionately affect warmer and poorer regions of the world. Comparison with lower emission scenarios emphasises the importance of mitigation policies for limiting global warming and reducing the associated health risks.

FUNDING

UK Medical Research Council.

Use of disdrometer data to evaluate the relationship of rainfall kinetic energy and intensity (KE-I)

Determination of rainfall kinetic energy (KE) is required to calculate erosivity, the ability of rainfall to detach soil particles and initiate erosion. Disdrometers can measure rainfall KE by measuring raindrop size and velocity. In the absence of such devices, KE is usually estimated with empirical equations that derive KE from measured rainfall intensity (I). We evaluated the performance of 14 different KE-I equations to estimate the 1min KE and event total KE, and compared these results with 821 observed rainfall events recorded by an optical disdrometer in the inner Ebro Basin, NE Spain. We also evaluated two sources of bias when using such relationships: bias from use of theoretical raindrop terminal velocities instead of measured values; and bias from time aggregation (recording rainfall intensity every 5, 10, 15, 30, and 60min). Empirical relationships performed well when complete events were considered (R(2)>0.90), but performed poorly for within-event variation (1min resolution). Also, several of the KE-I equations had large systematic biases. When raindrop size is known, estimation of terminal velocities by empirical laws led to overestimates of raindrop velocity and KE. Time aggregation led to large under-estimates of KE, although linear scaling successfully corrected for this bias.

Spatial Patterns of Heat-Related Cardiovascular Mortality in the Czech Republic

The study examines spatial patterns of effects of high temperature extremes on cardiovascular mortality in the Czech Republic at a district level during 1994–2009. Daily baseline mortality for each district was determined using a single location-stratified generalized additive model. Mean relative deviations of mortality from the baseline were calculated on days exceeding the 90th percentile of mean daily temperature in summer, and they were correlated with selected demographic, socioeconomic, and physical-environmental variables for the districts. Groups of districts with similar characteristics were identified according to socioeconomic status and urbanization level in order to provide a more general picture than possible on the district level. We evaluated lagged patterns of excess mortality after hot spell occurrences in: (i) urban areas vs. predominantly rural areas; and (ii) regions with different overall socioeconomic level. Our findings suggest that climatic conditions, altitude, and urbanization generally affect the spatial distribution of districts with the highest excess cardiovascular mortality, while socioeconomic status did not show a significant effect in the analysis across the Czech Republic as a whole. Only within deprived populations, socioeconomic status played a relevant role as well. After taking into account lagged effects of temperature on excess mortality, we found that the effect of hot spells was significant in highly urbanized regions, while most excess deaths in rural districts may be attributed to harvesting effects.

Contrasting patterns of hot spell effects on morbidity and mortality for cardiovascular diseases in the Czech Republic, 1994-2009

The study examines effects of hot spells on cardiovascular disease (CVD) morbidity and mortality in the population of the Czech Republic, with emphasis on differences between ischaemic heart disease (IHD) and cerebrovascular disease (CD) and between morbidity and mortality. Daily data on CVD morbidity (hospital admissions) and mortality over 1994-2009 were obtained from national hospitalization and mortality registers and standardized to account for long-term changes as well as seasonal and weekly cycles. Hot spells were defined as periods of at least two consecutive days with average daily air temperature anomalies above the 95% quantile during June to August. Relative deviations of mortality and morbidity from the baseline were evaluated. Hot spells were associated with excess mortality for all examined cardiovascular causes (CVD, IHD and CD). The increases were more pronounced for CD than IHD mortality in most population groups, mainly in males. In the younger population (0-64 years), however, significant excess mortality was observed for IHD while there was no excess mortality for CD. A short-term displacement effect was found to be much larger for mortality due to CD than IHD. Excess CVD mortality was not accompanied by increases in hospital admissions and below-expected-levels of morbidity prevailed during hot spells, particularly for IHD in the elderly. This suggests that out-of-hospital deaths represent a major part of excess CVD mortality during heat and that for in-hospital excess deaths CVD is a masked comorbid condition rather than the primary diagnosis responsible for hospitalization.

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.

Effects of sudden air pressure changes on hospital admissions for cardiovascular diseases in Prague, 1994-2009

Sudden weather changes have long been thought to be associated with negative impacts on human health, but relatively few studies have attempted to quantify these relationships. We use large 6-h changes in atmospheric pressure as a proxy for sudden weather changes and evaluate their association with hospital admissions for cardiovascular diseases (CVD). Winter and summer seasons and positive and negative pressure changes are analysed separately, using data for the city of Prague (population 1.2 million) over a 16-year period (1994-2009). We found that sudden pressure drops in winter are associated with significant rise in hospital admissions. Increased CVD morbidity was observed neither for pressure drops in summer nor pressure increases in any season. Analysis of synoptic weather maps shows that large pressure drops in winter are associated with strong zonal flow and rapidly moving low-pressure systems with centres over northern Europe and atmospheric fronts affecting western and central Europe. Analysis of links between passages of strong atmospheric fronts and hospital admissions, however, shows that the links disappear if weather changes are characterised by frontal passages. Sudden pressure drops in winter are associated also with significant excess CVD mortality. As climate models project strengthening of zonal circulation in winter and increased frequency of windstorms, the negative effects of such weather phenomena and their possible changes in a warmer climate of the twenty-first century need to be better understood, particularly as their importance in inducing excess morbidity and mortality in winter may increase compared to cold spells.

Impacts of hot and cold spells differ for acute and chronic ischaemic heart diseases

Background

Many studies have reported associations between temperature extremes and cardiovascular mortality but little has been understood about differences in the effects on acute and chronic diseases. The present study examines hot and cold spell effects on ischaemic heart disease (IHD) mortality in the Czech Republic during 1994–2009, with emphasis upon differences in the effects on acute myocardial infarction (AMI) and chronic IHD.

Methods

We use analogous definitions for hot and cold spells based on quantiles of daily average temperature anomalies, thus allowing for comparison of results for summer hot spells and winter cold spells. Daily mortality data were standardised to account for the long-term trend and the seasonal and weekly cycles. Periods when the data were affected by epidemics of influenza and other acute respiratory infections were removed from the analysis.

Results

Both hot and cold spells were associated with excess IHD mortality. For hot spells, chronic IHD was responsible for most IHD excess deaths in both male and female populations, and the impacts were much more pronounced in the 65+ years age group. The excess mortality from AMI was much lower compared to chronic IHD mortality during hot spells. For cold spells, by contrast, the relative excess IHD mortality was most pronounced in the younger age group (0–64 years), and we found different pattern for chronic IHD and AMI, with larger effects on AMI.

Conclusions

The findings show that while excess deaths due to IHD during hot spells are mainly of persons with chronic diseases whose health had already been compromised, cardiovascular changes induced by cold stress may result in deaths from acute coronary events rather than chronic IHD, and this effect is important also in the younger population. This suggests that the most vulnerable population groups as well as the most affected cardiovascular diseases differ between hot and cold spells, which needs to be taken into account when designing and implementing preventive actions.

Comparison of UTCI with other thermal indices in the assessment of heat and cold effects on cardiovascular mortality in the Czech Republic

We compare the recently developed Universal Thermal Climate Index (UTCI) with other thermal indices in analysing heat- and cold-related effects on cardiovascular (CVD) mortality in two different (urban and rural) regions in the Czech Republic during the 16-year period from 1994–2009. Excess mortality is represented by the number of deaths above expected daily values, the latter being adjusted for long-term changes, annual and weekly cycles, and epidemics of influenza/acute respiratory infections. Air temperature, UTCI, Apparent Temperature (AT) and Physiologically Equivalent Temperature (PET) are applied to identify days with heat and cold stress. We found similar heat effects on CVD mortality for air temperature and the examined thermal indices. Responses of CVD mortality to cold effects as characterised by different indices were much more varied. Particularly important is the finding that air temperature provides a weak cold effect in comparison with the thermal indices in both regions, so its application—still widespread in epidemiological studies—may underestimate the magnitude of cold-related mortality. These findings are important when possible climate change effects on heat- and cold-related mortality are estimated. AT and PET appear to be more universal predictors of heat- and cold- related mortality than UTCI when both urban and rural environments are of concern. UTCI tends to select windy rather than freezing days in winter, though these show little effect on mortality in the urban population. By contrast, significant cold-related mortality in the rural region if UTCI is used shows potential for UTCI to become a useful tool in cold exposure assessments.

Relationships between sudden weather changes in summer and mortality in the Czech Republic, 1986-2005

The study examines the relationship between sudden changes in weather conditions in summer, represented by (1) sudden air temperature changes, (2) sudden atmospheric pressure changes, and (3) passages of strong atmospheric fronts; and variations in daily mortality in the population of the Czech Republic. The events are selected from data covering 1986-2005 and compared with the database of daily excess all-cause mortality for the whole population and persons aged 70 years and above. Relative deviations of mortality, i.e., ratios of the excess mortality to the expected number of deaths, were averaged over the selected events for days D-2 (2 days before a change) up to D+7 (7 days after), and their statistical significance was tested by means of the Monte Carlo method. We find that the periods around weather changes are associated with pronounced patterns in mortality: a significant increase in mortality is found after large temperature increases and on days of large pressure drops; a decrease in mortality (partly due to a harvesting effect) occurs after large temperature drops, pressure increases, and passages of strong cold fronts. The relationship to variations in excess mortality is better expressed for sudden air temperature/pressure changes than for passages of atmospheric fronts. The mortality effects are usually more pronounced in the age group 70 years and above. The impacts associated with large negative changes of pressure are statistically independent of the effects of temperature; the corresponding dummy variable is found to be a significant predictor in the ARIMA model for relative deviations of mortality. This suggests that sudden weather changes should be tested also in time series models for predicting excess mortality as they may enhance their performance.

[Effects of sudden air temperature and pressure changes on mortality in the Czech Republic]

We have developed an algorithm for identifying sudden changes in air pressure and temperature over the Czech Republic. Such events were retrieved from the data covering in 1986-2005 and were matched with the daily numbers of all-cause deaths and deaths due to cardiovascular diseases from the national database, separately for the whole population and that aged 70 years and over. Excess daily mortality was determined by calculating deviations of the observed number of deaths from the expected number of deaths for each day in the respective groups. The relative deviation of the mortality the mean was calculated as the ratio of the excess mortality to the expected number of deaths. We used 3-hour air pressure data from 10 meteorological stations and hourly air temperature data from 9 stations representative of the Czech Republic. Pressure changes were evaluated on time scales of 3, 6 and 12 hours, separately for summer and winter time. Temperature changes were evaluated on a 24-hour time scale, separately for summer and winter season. Events characterized by pressure or temperature changes above the critical threshold and recorded within 24 hours at more than 50% of meteorological stations were retrieved. The critical thresholds were defined separately for each station using quantiles of distributions of air pressure and temperature changes. Relative mortality deviations for days D-2 (2 days before the change) to D+7 (7 days after the change) were averaged over the retrieved events. Statistical significance of the mean relative deviation was tested using the Monte Carlo method. Increased mortality followed large temperature increases and large pressure drops both in summer and winter months. Decreased mortality was observed after large pressure increases and large temperature drops in summer. Mortality variations are usually more pronounced in the population aged 70 years and over, and cardiovascular diseases account for most deaths after sudden temperature changes.

Classifications of atmospheric circulation patterns: recent advances and applications

We review recent advances in classifications of circulation patterns as a specific research area within synoptic climatology. The review starts with a general description of goals of classification and the historical development in the field. We put circulation classifications into a broader context within climatology and systematize the varied methodologies and approaches. We characterize three basic groups of classifications: subjective (also called manual), mixed (hybrid), and objective (computer-assisted, automated). The roles of cluster analysis and principal component analysis in the classification process are clarified. Several recent methodological developments in circulation classifications are identified and briefly described: the introduction of nonlinear methods, objectivization of subjective catalogs, efforts to optimize classifications, the need for intercomparisons of classifications, and the progress toward an optimum, if possible unified, classification method. Among the recent tendencies in the applications of circulation classifications, we mention a more extensive use in climate studies, both of past, present, and future climates, innovative applications in the ensemble forecasting, increasing variety of synoptic-climatological investigations, and steps above from the troposphere. After introducing the international activity within the field of circulation classifications, the COST733 Action, we briefly describe outputs of the inventory of classifications in Europe, which was carried out within the Action. Approaches to the evaluation of classifications and their mutual comparisons are also reviewed. A considerable part of the review is devoted to three examples of applications of circulation classifications: in historical climatology, in analyses of recent climate variations, and in analyses of outputs from global climate models.

Decreased impacts of the 2003 heat waves on mortality in the Czech Republic: an improved response?

The paper examines impacts on mortality of heat waves in 2003, the hottest summer on record in the Czech Republic, and compares them with previous similar events. While most summer heat waves over the period since 1986 were associated with significantly elevated mortality, this was not the case for three out of the four heat waves in 2003. The relatively weak mortality response was particularly noteworthy for the most severe heat wave which occurred in the first 10 days of August 2003 and resulted in enormous excess mortality in some western European countries. A mortality displacement effect and short-term adaptation to heat contributed to the reduced mortality impacts of the heat waves that followed after previous relatively warm periods. However, the decreased mortality response of the 2003 heat waves compared to previous heat waves in the 1990s is also likely to have arisen from positive health-care and other socio-economic changes in the post-communist central European region over the past decade, as well as a better public awareness of heat-related risks due to enhanced media coverage and regular biometeorological forecast and warnings.

Mortality and displaced mortality during heat waves in the Czech Republic

The aims of this study were to assess impacts of hot summer periods on mortality in the Czech Republic and to quantify the size of the short-term displacement effect which resulted in lower than expected mortality after heat waves. The analysis covered the period 1982-2000 when several extraordinarily hot summers occurred in central Europe. Daily total all-cause mortality and mortality due to cardiovascular diseases (CVD) in the entire population of the Czech Republic (approximately 10 million inhabitants) were examined. The daily death counts were standardized to account for the long-term decline in mortality and the seasonal and weekly cycles. Heat-related mortality is better expressed if 1-day lag after temperature is considered compared to the unlagged relationship. With the 1-day lag, both excess total mortality and excess CVD mortality were positive during all 17 heat waves, and in 14 (12) heat waves the increase in total (CVD) mortality was statistically significant ( P=0.05). The mean relative rise in total mortality during heat waves was 13%. The response was greater in females than males and similar regardless of whether total or CVD mortality was used. The largest relative increases, exceeding 20% in both total and CVD mortality, were associated with heat waves which occurred in early summer (the first half of July 1984 and June 1994). The mortality displacement effect played an important role since mortality tended to be lower than expected after hot periods. The mean net mortality change due to heat waves was estimated to be about a 1% increase in the number of deaths. The large relative increases during some heat waves were particularly noteworthy since the study (in contrast to most analyses of the heat stress/mortality relationship) was not restricted to an urban area and/or an elderly population.

[High summer temperatures and mortality in the Czech Republic 1982-2000]

The heat-stress-related mortality, which is among main impacts of periods of high summer temperature on society, was reported in many European countries, but analyses focusing on central European population have been rare. Results of the analysis for the period of 1982-2000 in the Czech Republic indicate that heat stress leads to a considerably increased all-causes mortality and mortality due to cardiovascular diseases. Periods with the highest deviations of the daily number of deaths from a baseline (in all-year data) are influenza epidemics and heat waves; the distribution of days with the highest excess mortality in a year is clearly bimodal, showing a main peak in winter and a secondary one in summer. Summer days with a considerably increased mortality are almost entirely days with a positive temperature deviation from the seasonal course. Deviations of mortality from the baseline exceed 100 deaths daily (more than 30% relative increase) in heat wave peaks, and the excess total mortality during the severe 1994 heat waves was +456 deaths (+10.3%) for June 17 to 30, and 8 deaths (+12.3%) for July 24 to August 8. The relative increase in mortality due to cardiovascular diseases was even more pronounced. The mortality displacement effect played an important role, since it was estimated to account for as much as 52% of the total number of victims for the June 1994 heat wave and 48% for the July-August heat wave. People who would die soon without oppressive weather conditions make about half of the total number of deaths, which is a larger value compared to what other studies reported. The increased mortality is observed at maximum (average, minimum) daily temperatures higher than 25 degrees C (18 degrees C, 14 degrees C) and their anomalies from mean seasonal courses larger than 3 degrees C. The same values hold for both the total and cardiovascular mortality. The mortality response at high temperatures is more pronounced in females than males. Correlations between mortality and temperature variables (including heat index) are positive and statistically significant (P = 0.01), stronger for deviations of meteorological variables from seasonal courses than for raw values, in females than males, and for heat index than for any temperature variable and summer simmer index. The unlagged correlations are stronger than correlations with lags 1-3 days; positive values of correlation coefficients hold for lags 0 to 3 days only while at lags of 4 to 25 days, the link is negative (mostly statistically significant) which demonstrates the mortality displacement effect and its time extent.