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

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

Mental health is positively associated with biodiversity in Canadian cities.

Cities concentrate problems that affect human well-being and biodiversity. Exploring the link between mental health and biodiversity can inform more holistic public health and urban planning. Here we examined associations between bird and tree species diversity estimates from eBird community science datasets and national forest inventories with self-rated mental health metrics from the Canadian Community Health Survey. We linked data across 36 Canadian Metropolitan Areas from 2007-2022 at a postal code level. After controlling for covariates, we found that bird and tree species diversity were significantly positively related to good self-reported mental health. Living in a postal code with bird diversity one standard deviation higher than the mean increased reporting of good mental health by 6.64%. Postal codes with tree species richness one standard deviation more than the mean increased reporting of good mental health by 5.36%. Our results suggest that supporting healthy urban ecosystems may also benefit human well-being.

Estimating the heat-related mortality and morbidity burden in the province of Quebec, Canada.

BACKGROUND: As climate change increases the frequency and intensity of extreme heat events, there is an urgent need to quantify the heat-related health burden. However, most past studies have focussed on a single health outcome (mainly mortality) or on specific heatwaves, thus providing limited knowledge of the total pressure heat exerts on health services. OBJECTIVES: This study aims to quantify the heat-related mortality and morbidity burden for five different health outcomes including all-cause mortality, hospitalizations, emergency department (ED) visits, ambulance transports and calls to a health hotline, using the province of Quebec (Canada) as a case study. METHODS: A two-step statistical analysis was employed to estimate regional heat-health relationships using Distributed Lag Non-Linear Models (DLNM) and pooled estimates using a multivariate meta-regression. Heat burden was quantified by attributable fraction (AF) and attributable number (AN) for two temperature ranges: all heat (above the minimum mortality/morbidity temperature) and extreme heat (above the 95th percentile of temperature). RESULTS: Higher temperatures were associated with greater risk ratios for all health outcomes studied, but at different levels. Significant AF ranging from 2 to 3% for the all heat effect and 0.4-1.0% for extreme heat were found for all health outcomes, except for hospitalizations that had an AF of 0.1% for both heat exposures. The estimated burden of all heat (and extreme heat) every summer across the province was 470 (200) deaths, 225 (170) hospitalizations, 36 000 (6 200) ED visits, 7 200 (1 500) ambulance transports and 15 000 (3 300) calls to a health hotline, all figures significant. DISCUSSION: This new knowledge on the total heat load will help public health authorities to target appropriate actions to reduce its burden now and in the future. The proposed state-of-the-art framework can easily be applied to other regions also experiencing the adverse effects of extreme heat.

High-resolution spatial and spatiotemporal modelling of air pollution using fixed site and mobile monitoring in a Canadian city.

The development of high-resolution spatial and spatiotemporal models of air pollutants is essential for exposure science and epidemiological applications. While fixed-site sampling has conventionally provided input data for statistical predictive models, the evolving mobile monitoring method offers improved spatial resolution, ideal for measuring pollutants with high spatial variability such as ultrafine particles (UFP). The Quebec Air Pollution Exposure and Epidemiology (QAPEE) study measured and modelled the spatial and spatiotemporal distributions of understudied pollutants, such as UFPs, black carbon (BC), and brown carbon (BrC), along with fine particulate matter (PM2.5), nitrogen dioxide (NO2), and ozone (O3) in Quebec City, Canada. We conducted a combined fixed-site (NO2 and O3) and mobile monitoring (PM2.5, BC, BrC, and UFPs) campaign over 10-months. Mobile monitoring routes were monitored on a weekly basis between 8am-10am and designed using location/allocation modelling. Seasonal fixed-site sampling campaigns captured continuous 24-h measurements over two-week periods. Generalized Additive Models (GAMs), which combined data on pollution concentrations with spatial, temporal, and spatiotemporal predictor variables were used to model and predict concentration surfaces. Annual models for PM2.5, NO2, O3 as well as seven of the smallest size fractions in the UFP range, had high out of sample predictive accuracy (range r2: 0.54-0.86). Varying spatial patterns were observed across UFP size ranges measured as Particle Number Counts (PNC). The monthly spatiotemporal models for PM2.5 (r2 = 0.49), BC (r2 = 0.27), BrC (r2 = 0.29), and PNC (r2 = 0.49) had moderate or moderate-low out of sample predictive accuracy. We conducted a sensitivity analysis and found that the minimum number of 'n visits' (mobile monitoring sessions) required to model annually representative air pollution concentrations was between 24 and 32 visits dependent on the pollutant. This study provides a single source of exposure models for a comprehensive set of air pollutants in Quebec City, Canada. These exposure models will feed into epidemiological research on the health impacts of ambient UFPs and other pollutants.

Associations between paediatric obesity, chemical mixtures and environmental factors, in a national cross-sectional study of Canadian children.

BACKGROUND: Whilst single chemical exposures are suspected to be obesogenic, the combined role of chemical mixtures in paediatric obesity is not well understood. OBJECTIVES: We aimed to evaluate the potential associations between chemical mixtures and obesity in a population-based sample of Canadian children. METHODS: We ascertained biomonitoring and health data for children aged 3-11 from the cross-sectional Canadian Health Measures Survey from 2007 to 2019. Several chemicals of interest were measured in blood or urine and paediatric obesity was defined based on measured anthropometrics. Using quantile-based G computational analysis, we quantified the effects of three chemical mixtures selected a priori. Models were adjusted for sociodemographic and environmental factors identified through a directed acyclic graph. Results are presented through adjusted relative risks (RR) with 95% confidence intervals (95% CI). RESULTS: We included 9147 children. Of these, 24.1% were overweight or obese. Exposure to the mixture of bisphenol A, acrylamide, glycidamide, metals, parabens and arsenic increased the risk of childhood overweight or obesity by 45% (95% CI 1.09, 1.93), obesity by 109% (95% CI 1.27, 3.42) and central obesity by 82% (95% CI 1.30, 2.56). CONCLUSIONS: Our findings support the role of early childhood chemical exposures in paediatric obesity and the potential combined effects of chemicals.

Airborne Nanoparticle Concentrations Are Associated with Increased Mortality Risk in Canada's Two Largest Cities.

RATIONALE: Outdoor fine particulate air pollution (PM2.5) contributes to millions of deaths around the world each year, but much less is known about the long-term health impacts of other particulate air pollutants including ultrafine particles (a.k.a. nanoparticles) which are in the nanometer size range (<100 nm), widespread in urban environments, and not currently regulated. OBJECTIVES: Estimate the associations between long-term exposure to outdoor ultrafine particles and mortality. METHODS: Outdoor air pollution levels were linked to the residential addresses of a large, population-based cohort from 2001 - 2016. Associations between long-term exposure to outdoor ultrafine particles and nonaccidental and cause-specific mortality were estimated using Cox proportional hazards models. MEASUREMENTS: An increase in long-term exposure to outdoor ultrafine particles was associated with an increased risk of nonaccidental mortality (Hazard Ratio = 1. 073, 95% Confidence Interval = 1. 061, 1. 085) and cause-specific mortality, the strongest of which was respiratory mortality (Hazard Ratio = 1.174, 95% Confidence Interval = 1.130, 1.220). MAIN RESULTS: Long-term exposure to outdoor ultrafine particles was associated with increased risk of mortality. We estimated the mortality burden for outdoor ultrafine particles in Montreal and Toronto, Canada to be approximately 1100 additional nonaccidental deaths every year. Furthermore, we observed possible confounding by particle size which suggests that previous studies may have underestimated or missed important health risks associated with ultrafine particles. CONCLUSIONS: As outdoor ultrafine particles are not currently regulated, there is great potential for future regulatory interventions to improve population health by targeting these common outdoor air pollutants.

Extreme Temperatures and Stroke Mortality: Evidence From a Multi-Country Analysis.

BACKGROUND: Extreme temperatures contribute significantly to global mortality. While previous studies on temperature and stroke-specific outcomes presented conflicting results, these studies were predominantly limited to single-city or single-country analyses. Their findings are difficult to synthesize due to variations in methodologies and exposure definitions. METHODS: Within the Multi-Country Multi-City Network, we built a new mortality database for ischemic and hemorrhagic stroke. Applying a unified analysis protocol, we conducted a multinational case-crossover study on the relationship between extreme temperatures and stroke. In the first stage, we fitted a conditional quasi-Poisson regression for daily mortality counts with distributed lag nonlinear models for temperature exposure separately for each city. In the second stage, the cumulative risk from each city was pooled using mixed-effect meta-analyses, accounting for clustering of cities with similar features. We compared temperature-stroke associations across country-level gross domestic product per capita. We computed excess deaths in each city that are attributable to the 2.5% hottest and coldest of days based on each city's temperature distribution. RESULTS: We collected data for a total of 3 443 969 ischemic strokes and 2 454 267 hemorrhagic stroke deaths from 522 cities in 25 countries. For every 1000 ischemic stroke deaths, we found that extreme cold and hot days contributed 9.1 (95% empirical CI, 8.6-9.4) and 2.2 (95% empirical CI, 1.9-2.4) excess deaths, respectively. For every 1000 hemorrhagic stroke deaths, extreme cold and hot days contributed 11.2 (95% empirical CI, 10.9-11.4) and 0.7 (95% empirical CI, 0.5-0.8) excess deaths, respectively. We found that countries with low gross domestic product per capita were at higher risk of heat-related hemorrhagic stroke mortality than countries with high gross domestic product per capita (P=0.02). CONCLUSIONS: Both extreme cold and hot temperatures are associated with an increased risk of dying from ischemic and hemorrhagic strokes. As climate change continues to exacerbate these extreme temperatures, interventional strategies are needed to mitigate impacts on stroke mortality, particularly in low-income countries.

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.

Temperature frequency and mortality: Assessing adaptation to local temperature.

Assessing the association between temperature frequency and mortality can provide insights into human adaptation to local ambient temperatures. We collected daily time-series data on mortality and temperature from 757 locations in 47 countries/regions during 1979-2020. We used a two-stage time series design to assess the association between temperature frequency and all-cause mortality. The results were pooled at the national, regional, and global levels. We observed a consistent decrease in the risk of mortality as the normalized frequency of temperature increases across the globe. The average increase in mortality risk comparing the 10th to 100th percentile of normalized frequency was 13.03% (95% CI: 12.17-13.91), with substantial regional differences (from 4.56% in Australia and New Zealand to 33.06% in South Europe). The highest increase in mortality was observed for high-income countries (13.58%, 95% CI: 12.56-14.61), followed by lower-middle-income countries (12.34%, 95% CI: 9.27-15.51). This study observed a declining risk of mortality associated with higher temperature frequency. Our findings suggest that populations can adapt to their local climate with frequent exposure, with the adapting ability varying geographically due to differences in climatic and socioeconomic characteristics.

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.

Associations of neighborhood greenspace, and active living environments with autism spectrum disorders: A matched case-control study in Ontario, Canada.

BACKGROUND: Increasing evidence links early life residential exposure to natural urban environmental attributes and positive health outcomes in children. However, few studies have focused on their protective effects on the risk of autism spectrum disorder (ASD). The aim of this study was to investigate the associations of neighborhood greenspace, and active living environments during pregnancy with ASD in young children (≤6 years). METHODS: We conducted a population-based matched case-control study of singleton term births in Ontario, Canada for 2012-2016. The ASD and environmental data was generated using the Ontario Autism Spectrum Profile, the Better Outcomes Registry & Network Ontario, and Canadian Urban Environmental Health Research Consortium. We employed conditional logistic regressions to estimate the odds ratio (OR) between ASD and environmental factors characterizing selected greenspace metrics and neighborhoods conducive to active living (i.e., green view index (GVI), normalized difference vegetation index (NDVI), tree canopy, park proximity and active living environments index (ALE)). RESULTS: We linked 8643 mother-child pairs, including 1554 cases (18%). NDVI (OR 1.034, 0.944-1.024, per Inter Quartile Range [IQR] = 0.08), GVI (OR 1.025, 95% CI 0.953-1.087, per IQR = 9.45%), tree canopy (OR 0.992, 95% CI 0.903-1.089, per IQR = 6.24%) and the different categories of ALE were not associated with ASD in adjusted models for air pollution. In contrast, living closer to a park was protective (OR 0.888, 0.833-0.948, per 0.06 increase in park proximity index), when adjusted for air pollution. CONCLUSIONS: This study reported mixed findings showing both null and beneficial effects of green spaces and active living environments on ASD. Further investigations are warranted to elucidate the role of exposure to greenspaces and active living environments on the development of ASD.

Estimates of global mortality burden associated with short-term exposure to fine particulate matter (PM2·5).

BACKGROUND: The acute health effects of short-term (hours to days) exposure to fine particulate matter (PM2·5) have been well documented; however, the global mortality burden attributable to this exposure has not been estimated. We aimed to estimate the global, regional, and urban mortality burden associated with short-term exposure to PM2·5 and the spatiotemporal variations in this burden from 2000 to 2019. METHODS: We combined estimated global daily PM2·5 concentrations, annual population counts, country-level mortality rates, and epidemiologically derived exposure-response functions to estimate the mortality attributable to short-term PM2·5 exposure from 2000 to 2019, in the continental regions and in 13 189 urban centres worldwide at a spatial resolution of 0·1°â€ˆ× 0·1°. We tested the robustness of our mortality estimates with different theoretical minimum risk exposure levels, lag effects, and exposure-response functions. FINDINGS: Approximately 1 million (95% CI 690 000-1·3 million) premature deaths per year from 2000 to 2019 were attributable to short-term PM2·5 exposure, representing 2·08% (1·41-2·75) of total global deaths or 17 (11-22) premature deaths per 100 000 population. Annually, 0·23 million (0·15 million-0·30 million) deaths attributable to short-term PM2·5 exposure were in urban areas, constituting 22·74% of the total global deaths attributable to this cause and accounting for 2·30% (1·56-3·05) of total global deaths in urban areas. The sensitivity analyses showed that our worldwide estimates of mortality attributed to short-term PM2·5 exposure were robust. INTERPRETATION: Short-term exposure to PM2·5 contributes a substantial global mortality burden, particularly in Asia and Africa, as well as in global urban areas. Our results highlight the importance of mitigation strategies to reduce short-term exposure to air pollution and its adverse effects on human health. FUNDING: Australian Research Council and the Australian National Health and Medical Research Council.

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.

Residential surrounding greenness and the incidence of childhood asthma: Findings from a population-based cohort in Ontario, Canada.

Several epidemiological studies have investigated the possible role that living in areas with greater amounts of greenspace has on the incidence of childhood asthma. These findings have been inconsistent, and few studies explored the relevance of timing of exposure. We investigated the role of residential surrounding greenness on the risk of incident asthma using a population-based retrospective cohort study. We included 982,131 singleton births in Ontario, Canada between 2006 and 2013. Two measures of greenness, the Normalized Difference Vegetation Index (NDVI) and the Green View Index (GVI), were assigned to the residential histories of these infants from pregnancy through to 12 years of age. Longitudinally-based diagnoses of asthma were determined by using provincial administrative health data. The extended Cox hazards model was used to characterize associations between greenness measures and asthma (up to age 12 years) while adjusting for several risk factors. In a fully adjusted model, that included a term for traffic-related air pollution (NO2), we found no association between an interquartile range increase (0.08) of the NDVI during childhood and asthma incidence (HR = 0.99; 95 % CI = 0.99-1.01). In contrast, we found that an 0.08 increase in NDVI during childhood reduced the risk of asthma in children 7-12 years of age by 14 % (HR = 0.86, 95 % CI:0.79-0.95). Seasonal differences in the association between greenness and asthma were noted. Our findings suggest that residential proximity to greenness reduces the risk of asthma in children aged 7-12.

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.

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.

Tropical cyclone-specific mortality risks and the periods of concern: A multicountry time-series study.

BACKGROUND: More intense tropical cyclones (TCs) are expected in the future under a warming climate scenario, but little is known about their mortality effect pattern across countries and over decades. We aim to evaluate the TC-specific mortality risks, periods of concern (POC) and characterize the spatiotemporal pattern and exposure-response (ER) relationships on a multicountry scale. METHODS AND FINDINGS: Daily all-cause, cardiovascular, and respiratory mortality among the general population were collected from 494 locations in 18 countries or territories during 1980 to 2019. Daily TC exposures were defined when the maximum sustained windspeed associated with a TC was ≥34 knots using a parametric wind field model at a 0.5° × 0.5° resolution. We first estimated the TC-specific mortality risks and POC using an advanced flexible statistical framework of mixed Poisson model, accounting for the population changes, natural variation, seasonal and day of the week effects. Then, a mixed meta-regression model was used to pool the TC-specific mortality risks to estimate the overall and country-specific ER relationships of TC characteristics (windspeed, rainfall, and year) with mortality. Overall, 47.7 million all-cause, 15.5 million cardiovascular, and 4.9 million respiratory deaths and 382 TCs were included in our analyses. An overall average POC of around 20 days was observed for TC-related all-cause and cardiopulmonary mortality, with relatively longer POC for the United States of America, Brazil, and Taiwan (>30 days). The TC-specific relative risks (RR) varied substantially, ranging from 1.04 to 1.42, 1.07 to 1.77, and 1.12 to 1.92 among the top 100 TCs with highest RRs for all-cause, cardiovascular, and respiratory mortality, respectively. At country level, relatively higher TC-related mortality risks were observed in Guatemala, Brazil, and New Zealand for all-cause, cardiovascular, and respiratory mortality, respectively. We found an overall monotonically increasing and approximately linear ER curve of TC-related maximum sustained windspeed and cumulative rainfall with mortality, with heterogeneous patterns across countries and regions. The TC-related mortality risks were generally decreasing from 1980 to 2019, especially for the Philippines, Taiwan, and the USA, whereas potentially increasing trends in TC-related all-cause and cardiovascular mortality risks were observed for Japan. CONCLUSIONS: The TC mortality risks and POC varied greatly across TC events, locations, and countries. To minimize the TC-related health burdens, targeted strategies are particularly needed for different countries and regions, integrating epidemiological evidence on region-specific POC and ER curves that consider across-TC variability.

The relationship between hot temperatures and hospital admissions for psychosis in adults diagnosed with schizophrenia: A case-crossover study in Quebec, Canada.

INTRODUCTION: Some studies have found hot temperatures to be associated with exacerbations of schizophrenia, namely psychoses. As climate changes faster in Northern countries, our understanding of the association between temperature and hospital admissions (HA) for psychosis needs to be deepened. OBJECTIVES: 1) Among adults diagnosed with schizophrenia, measure the relationship between mean temperatures and HAs for psychosis during summer. 2) Determine the influence of individual and ecological characteristics on this relationship. METHODS: A cohort of adults diagnosed with schizophrenia (n = 30,649) was assembled using Quebec's Integrated Chronic Disease Surveillance System (QICDSS). The follow-up spanned summers from 2001 to 2019, using hospital data from the QICDSS and meteorological data from the National Aeronautics and Space Administration's (NASA) Daymet database. In four geographic regions of the province of Quebec, a conditional logistic regression was used for the case-crossover analysis of the relationship between mean temperatures (at lags up to 6 days) and HAs for psychosis using a distributed lag non-linear model (DLNM). The analyses were adjusted for relative humidity, stratified according to individual (age, sex, and comorbidities) and ecological (material and social deprivation index and exposure to green space) factors, and then pooled through a meta-regression. RESULTS: The statistical analyses revealed a statistically significant increase in HAs three days (lag 3) after elevated mean temperatures corresponding to the 90th percentile relative to a minimum morbidity temperature (MMT) (OR 1.040; 95% CI 1.008-1.074), while the cumulative effect over six days was not statistically significant (OR 1.052; 95% IC 0.993-1.114). Stratified analyses revealed non statistically significant gradients of increasing HAs relative to increasing material deprivation and decreasing green space levels. CONCLUSIONS: The statistical analyses conducted in this project showed the pattern of admissions for psychosis after hot days. This finding could be useful to better plan health services in a rapidly changing climate.