Rainfall events and daily mortality across 645 global locations: two stage time series analysis.

OBJECTIVE: To examine the associations between characteristics of daily rainfall (intensity, duration, and frequency) and all cause, cardiovascular, and respiratory mortality. DESIGN: Two stage time series analysis. SETTING: 645 locations across 34 countries or regions. POPULATION: Daily mortality data, comprising a total of 109 954 744 all cause, 31 164 161 cardiovascular, and 11 817 278 respiratory deaths from 1980 to 2020. MAIN OUTCOME MEASURE: Association between daily mortality and rainfall events with return periods (the expected average time between occurrences of an extreme event of a certain magnitude) of one year, two years, and five years, with a 14 day lag period. A continuous relative intensity index was used to generate intensity-response curves to estimate mortality risks at a global scale. RESULTS: During the study period, a total of 50 913 rainfall events with a one year return period, 8362 events with a two year return period, and 3301 events with a five year return period were identified. A day of extreme rainfall with a five year return period was significantly associated with increased daily all cause, cardiovascular, and respiratory mortality, with cumulative relative risks across 0-14 lag days of 1.08 (95% confidence interval 1.05 to 1.11), 1.05 (1.02 to 1.08), and 1.29 (1.19 to 1.39), respectively. Rainfall events with a two year return period were associated with respiratory mortality only, whereas no significant associations were found for events with a one year return period. Non-linear analysis revealed protective effects (relative risk <1) with moderate-heavy rainfall events, shifting to adverse effects (relative risk >1) with extreme intensities. Additionally, mortality risks from extreme rainfall events appeared to be modified by climate type, baseline variability in rainfall, and vegetation coverage, whereas the moderating effects of population density and income level were not significant. Locations with lower variability of baseline rainfall or scarce vegetation coverage showed higher risks. CONCLUSION: Daily rainfall intensity is associated with varying health effects, with extreme events linked to an increasing relative risk for all cause, cardiovascular, and respiratory mortality. The observed associations varied with local climate and urban infrastructure.

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.

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.

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.

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.

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

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

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.

Risk of ischemic and hemorrhagic stroke in relation to cold spells in four seasons.

BACKGROUND: Cold winter weather increases the risk of stroke, but the evidence is scarce on whether the risk increases during season-specific cold weather in the other seasons. The objective of our study was to test the hypothesis of an association between personal cold spells and different types of stroke in the season-specific context, and to formally assess effect modification by age and sex. METHODS: We conducted a case-crossover study of all 5396 confirmed 25-64 years old cases with stroke in the city of Kaunas, Lithuania, 2000-2015. We assigned to each case a one-week hazard period and 15 reference periods of the same calendar days of other study years. A personal cold day was defined for each case with a mean temperature below the fifth percentile of the frequency distribution of daily mean temperatures of the hazard and reference periods. Conditional logistic regression was applied to estimate odds ratios (OR) and 95% confidence intervals (95% CI) representing associations between time- and place-specific cold weather and stroke. RESULTS: There were positive associations between cold weather and stroke in Kaunas, with each additional cold day during the week before the stroke increases the risk by 3% (OR 1.03; 95% CI 1.00-1.07). The association was present for ischemic stroke (OR 1.05; 95% CI 1.01-1.09) but not hemorrhagic stroke (OR 0.98; 95% CI 0.91-1.06). In the summer, the risk of stroke increased by 8% (OR 1.08; 95% CI 1.00-1.16) per each additional cold day during the hazard period. Age and sex did not modify the effect. CONCLUSIONS: Our findings show that personal cold spells increase the risk of stroke, and this pertains to ischemic stroke specifically. Most importantly, cold weather in the summer season may be a previously unrecognized determinant of stroke.

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.

Associations Between Extreme Temperatures and Cardiovascular Cause-Specific Mortality: Results From 27 Countries.

BACKGROUND: Cardiovascular disease is the leading cause of death worldwide. Existing studies on the association between temperatures and cardiovascular deaths have been limited in geographic zones and have generally considered associations with total cardiovascular deaths rather than cause-specific cardiovascular deaths. METHODS: We used unified data collection protocols within the Multi-Country Multi-City Collaborative Network to assemble a database of daily counts of specific cardiovascular causes of death from 567 cities in 27 countries across 5 continents in overlapping periods ranging from 1979 to 2019. City-specific daily ambient temperatures were obtained from weather stations and climate reanalysis models. To investigate cardiovascular mortality associations with extreme hot and cold temperatures, we fit case-crossover models in each city and then used a mixed-effects meta-analytic framework to pool individual city estimates. Extreme temperature percentiles were compared with the minimum mortality temperature in each location. Excess deaths were calculated for a range of extreme temperature days. RESULTS: The analyses included deaths from any cardiovascular cause (32 154 935), ischemic heart disease (11 745 880), stroke (9 351 312), heart failure (3 673 723), and arrhythmia (670 859). At extreme temperature percentiles, heat (99th percentile) and cold (1st percentile) were associated with higher risk of dying from any cardiovascular cause, ischemic heart disease, stroke, and heart failure as compared to the minimum mortality temperature, which is the temperature associated with least mortality. Across a range of extreme temperatures, hot days (above 97.5th percentile) and cold days (below 2.5th percentile) accounted for 2.2 (95% empirical CI [eCI], 2.1-2.3) and 9.1 (95% eCI, 8.9-9.2) excess deaths for every 1000 cardiovascular deaths, respectively. Heart failure was associated with the highest excess deaths proportion from extreme hot and cold days with 2.6 (95% eCI, 2.4-2.8) and 12.8 (95% eCI, 12.2-13.1) for every 1000 heart failure deaths, respectively. CONCLUSIONS: Across a large, multinational sample, exposure to extreme hot and cold temperatures was associated with a greater risk of mortality from multiple common cardiovascular conditions. The intersections between extreme temperatures and cardiovascular health need to be thoroughly characterized in the present day-and especially under a changing climate.

Short-term exposure to ambient air pollution and individual emergency department visits for COVID-19: a case-crossover study in Canada.

BACKGROUND: Ambient air pollution is thought to contribute to increased risk of COVID-19, but the evidence is controversial. OBJECTIVE: To evaluate the associations between short-term variations in outdoor concentrations of ambient air pollution and COVID-19 emergency department (ED) visits. METHODS: We conducted a case-crossover study of 78 255 COVID-19 ED visits in Alberta and Ontario, Canada between 1 March 2020 and 31 March 2021. Daily air pollution data (ie, fine particulate matter with diameter less than 2.5 µm (PM2.5), nitrogen dioxide (NO2) and ozone were assigned to individual case of COVID-19 in 10 km × 10 km grid resolution. Conditional logistic regression was used to estimate associations between air pollution and ED visits for COVID-19. RESULTS: Cumulative ambient exposure over 0-3 days to PM2.5 (OR 1.010; 95% CI 1.004 to 1.015, per 6.2 µg/m3) and NO2 (OR 1.021; 95% CI 1.015 to 1.028, per 7.7 ppb) concentrations were associated with ED visits for COVID-19. We found that the association between PM2.5 and COVID-19 ED visits was stronger among those hospitalised following an ED visit, as a measure of disease severity, (OR 1.023; 95% CI 1.015 to 1.031) compared with those not hospitalised (OR 0.992; 95% CI 0.980 to 1.004) (p value for effect modification=0.04). CONCLUSIONS: We found associations between short-term exposure to ambient air pollutants and COVID-19 ED visits. Exposure to air pollution may also lead to more severe COVID-19 disease.

The association of COVID-19 incidence with temperature, humidity, and UV radiation - A global multi-city analysis.

BACKGROUND AND AIM: The associations between COVID-19 transmission and meteorological factors are scientifically debated. Several studies have been conducted worldwide, with inconsistent findings. However, often these studies had methodological issues, e.g., did not exclude important confounding factors, or had limited geographic or temporal resolution. Our aim was to quantify associations between temporal variations in COVID-19 incidence and meteorological variables globally. METHODS: We analysed data from 455 cities across 20 countries from 3 February to 31 October 2020. We used a time-series analysis that assumes a quasi-Poisson distribution of the cases and incorporates distributed lag non-linear modelling for the exposure associations at the city-level while considering effects of autocorrelation, long-term trends, and day of the week. The confounding by governmental measures was accounted for by incorporating the Oxford Governmental Stringency Index. The effects of daily mean air temperature, relative and absolute humidity, and UV radiation were estimated by applying a meta-regression of local estimates with multi-level random effects for location, country, and climatic zone. RESULTS: We found that air temperature and absolute humidity influenced the spread of COVID-19 over a lag period of 15 days. Pooling the estimates globally showed that overall low temperatures (7.5 °C compared to 17.0 °C) and low absolute humidity (6.0 g/m3 compared to 11.0 g/m3) were associated with higher COVID-19 incidence (RR temp =1.33 with 95%CI: 1.08; 1.64 and RR AH =1.33 with 95%CI: 1.12; 1.57). RH revealed no significant trend and for UV some evidence of a positive association was found. These results were robust to sensitivity analysis. However, the study results also emphasise the heterogeneity of these associations in different countries. CONCLUSION: Globally, our results suggest that comparatively low temperatures and low absolute humidity were associated with increased risks of COVID-19 incidence. However, this study underlines regional heterogeneity of weather-related effects on COVID-19 transmission.

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 99th 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.

Coarse Particulate Air Pollution and Daily Mortality: A Global Study in 205 Cities.

Rationale: The associations between ambient coarse particulate matter (PM2.5-10) and daily mortality are not fully understood on a global scale. Objectives: To evaluate the short-term associations between PM2.5-10 and total, cardiovascular, and respiratory mortality across multiple countries/regions worldwide. Methods: We collected daily mortality (total, cardiovascular, and respiratory) and air pollution data from 205 cities in 20 countries/regions. Concentrations of PM2.5-10 were computed as the difference between inhalable and fine PM. A two-stage time-series analytic approach was applied, with overdispersed generalized linear models and multilevel meta-analysis. We fitted two-pollutant models to test the independent effect of PM2.5-10 from copollutants (fine PM, nitrogen dioxide, sulfur dioxide, ozone, and carbon monoxide). Exposure-response relationship curves were pooled, and regional analyses were conducted. Measurements and Main Results: A 10 µg/m3 increase in PM2.5-10 concentration on lag 0-1 day was associated with increments of 0.51% (95% confidence interval [CI], 0.18%-0.84%), 0.43% (95% CI, 0.15%-0.71%), and 0.41% (95% CI, 0.06%-0.77%) in total, cardiovascular, and respiratory mortality, respectively. The associations varied by country and region. These associations were robust to adjustment by all copollutants in two-pollutant models, especially for PM2.5. The exposure-response curves for total, cardiovascular, and respiratory mortality were positive, with steeper slopes at lower exposure ranges and without discernible thresholds. Conclusions: This study provides novel global evidence on the robust and independent associations between short-term exposure to ambient PM2.5-10 and total, cardiovascular, and respiratory mortality, suggesting the need to establish a unique guideline or regulatory limit for daily concentrations of PM2.5-10.

Differential Mortality Risks Associated With PM2.5 Components: A Multi-Country, Multi-City Study.

BACKGROUND: The association between fine particulate matter (PM2.5) and mortality widely differs between as well as within countries. Differences in PM2.5 composition can play a role in modifying the effect estimates, but there is little evidence about which components have higher impacts on mortality. METHODS: We applied a 2-stage analysis on data collected from 210 locations in 16 countries. In the first stage, we estimated location-specific relative risks (RR) for mortality associated with daily total PM2.5 through time series regression analysis. We then pooled these estimates in a meta-regression model that included city-specific logratio-transformed proportions of seven PM2.5 components as well as meta-predictors derived from city-specific socio-economic and environmental indicators. RESULTS: We found associations between RR and several PM2.5 components. Increasing the ammonium (NH4+) proportion from 1% to 22%, while keeping a relative average proportion of other components, increased the RR from 1.0063 (95% confidence interval [95% CI] = 1.0030, 1.0097) to 1.0102 (95% CI = 1.0070, 1.0135). Conversely, an increase in nitrate (NO3-) from 1% to 71% resulted in a reduced RR, from 1.0100 (95% CI = 1.0067, 1.0133) to 1.0037 (95% CI = 0.9998, 1.0077). Differences in composition explained a substantial part of the heterogeneity in PM2.5 risk. CONCLUSIONS: These findings contribute to the identification of more hazardous emission sources. Further work is needed to understand the health impacts of PM2.5 components and sources given the overlapping sources and correlations among many components.