Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE): Phase 2.
Oct 12, 2022Research report (Health Effects Institute)
Links Between Air Pollution and Death Rates in Areas with Low Pollution: Phase 2
AI simplified
Abstract
In a cohort of 7.1 million Canadians, each 10-μg/m increase in outdoor fine particulate matter (PM) concentration is associated with a 8.4% increase in the risk of nonaccidental mortality.
- Positive associations between outdoor PM concentrations and nonaccidental mortality were consistently observed across all cohorts.
- Risk estimates indicated that exposure to PM concentrations as low as 2.5 μg/m may be linked to increased mortality.
- Specific causes of death, including ischemic heart disease, respiratory disease, cardiovascular disease, and diabetes, were associated with PM exposure.
- The concentration-response relationship for nonaccidental mortality flattened between PM concentrations of 5 to 9 μg/m before increasing again at higher levels.
- Sensitivity analyses suggested that the presence of other pollutants like ozone could attenuate the associations between PM and mortality.
AI simplified
INTRODUCTION: Mortality is associated with long-term exposure to fine particulate matter (particulate matter ≤2.5 μm in aerodynamic diameter; PM), although the magnitude and form of these associations remain poorly understood at lower concentrations. Knowledge gaps include the shape of concentration-response curves and the lowest levels of exposure at which increased risks are evident and the occurrence and extent of associations with specific causes of death. Here, we applied improved estimates of exposure to ambient PMto national population-based cohorts in Canada, including a stacked cohort of 7.1 million people who responded to census year 1991, 1996, or 2001. The characterization of the shape of the concentration-response relationship for nonaccidental mortality and several specific causes of death at low levels of exposure was the focus of the Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE) Phase 1 report. In the Phase 1 report we reported that associations between outdoor PMconcentrations and nonaccidental mortality were attenuated with the addition of ozone (O) or a measure of gaseous pollutant oxidant capacity (O), which was estimated from Oand nitrogen dioxide (NO) concentrations. This was motivated by our interests in understanding both the effects air pollutant mixtures may have on mortality and also the role of Oas a copollutant that shares common sources and precursor emissions with those of PM. In this Phase 2 report, we further explore the sensitivity of these associations with Oand O, evaluate sensitivity to other factors, such as regional variation, and present ambient PMconcentration-response relationships for specific causes of death. 2.5 2.5 2.5 3x3 2 3 2.5 3x2.5
METHODS: PMconcentrations were estimated at 1 kmspatial resolution across North America using remote sensing of aerosol optical depth (AOD) combined with chemical transport model (GEOS-Chem) simulations of the AOD:surface PMmass concentration relationship, land use information, and ground monitoring. These estimates were informed and further refined with collocated measurements of PMand AOD, including targeted measurements in areas of low PMconcentrations collected at five locations across Canada. Ground measurements of PMand total suspended particulate matter (TSP) mass concentrations from 1981 to 1999 were used to backcast remote-sensing-based estimates over that same time period, resulting in modeled annual surfaces from 1981 to 2016. 2.5 2.5 2.5 2.5 2.5 2
UNLABELLED: Annual exposures to PMwere then estimated for subjects in several national population-based Canadian cohorts using residential histories derived from annual postal code entries in income tax files. These cohorts included three census-based cohorts: the 1991 Canadian Census Health and Environment Cohort (CanCHEC; 2.5 million respondents), the 1996 CanCHEC (3 million respondents), the 2001 CanCHEC (3 million respondents), and a Stacked CanCHEC where duplicate records of respondents were excluded (Stacked CanCHEC; 7.1 million respondents). The Canadian Community Health Survey (CCHS) mortality cohort (mCCHS), derived from several pooled cycles of the CCHS (540,900 respondents), included additional individual information about health behaviors. Follow-up periods were completed to the end of 2016 for all cohorts. Cox proportional hazard ratios (HRs) were estimated for nonaccidental and other major causes of death using a 10-year moving average exposure and 1-year lag. All models were stratified by age, sex, immigrant status, and where appropriate, census year or survey cycle. Models were further adjusted for income adequacy quintile, visible minority status, Indigenous identity, educational attainment, labor-force status, marital status, occupation, and ecological covariates of community size, airshed, urban form, and four dimensions of the Canadian Marginalization Index (Can-Marg; instability, deprivation, dependency, and ethnic concentration). The mCCHS analyses were also adjusted for individual-level measures of smoking, alcohol consumption, fruit and vegetable consumption, body mass index (BMI), and exercise behavior. 2.5
UNLABELLED: In addition to linear models, the shape of the concentration-response function was investigated using restricted cubic splines (RCS). The number of knots were selected by minimizing the Bayesian Information Criterion (BIC). Two additional models were used to examine the association between nonaccidental mortality and PM. The first is the standard threshold model defined by a transformation of concentration equaling zero if the concentration was less than a specific threshold value and concentration minus the threshold value for concentrations above the threshold. The second additional model was an extension of the Shape Constrained Health Impact Function (SCHIF), the eSCHIF, which converts RCS predictions into functions potentially more suitable for use in health impact assessments. Given the RCS parameter estimates and their covariance matrix, 1,000 realizations of the RCS were simulated at concentrations from the minimum to the maximum concentration, by increments of 0.1 μg/m. An eSCHIF was then fit to each of these RCS realizations. Thus, 1,000 eSCHIF predictions and uncertainty intervals were determined at each concentration within the total range. 2.5 3
UNLABELLED: Sensitivity analyses were conducted to examine associations between PMand mortality when in the presence of, or stratified by tertile of, Oor O. Additionally, associations between PMand mortality were assessed for sensitivity to lower concentration thresholds, where person-years below a threshold value were assigned the mean exposure within that group. We also examined the sensitivity of the shape of the nonaccidental mortality-PMassociation to removal of person-years at or above 12 μg/m(the current U.S. National Ambient Air Quality Standard) and 10 μg/m(the current Canadian and former [2005] World Health Organization [WHO] guideline, and current WHO Interim Target-4). Finally, differences in the shapes of PM-mortality associations were assessed across broad geographic regions (airsheds) within Canada. 2.5 3x2.5 2.5 2.5 3 3
RESULTS: The refined PMexposure estimates demonstrated improved performance relative to estimates applied previously and in the MAPLE Phase 1 report, with slightly reduced errors, including at lower ranges of concentrations (e.g., for PM<10 μg/m). 2.5 2.5 3
UNLABELLED: Positive associations between outdoor PMconcentrations and nonaccidental mortality were consistently observed in all cohorts. In the Stacked CanCHEC analyses (1.3 million deaths), each 10-μg/mincrease in outdoor PMconcentration corresponded to an HR of 1.084 (95% confidence interval [CI]: 1.073 to 1.096) for nonaccidental mortality. For an interquartile range (IQR) increase in PMmass concentration of 4.16 μg/mand for a mean annual nonaccidental death rate of 92.8 per 10,000 persons (over the 1991-2016 period for cohort participants ages 25-90), this HR corresponds to an additional 31.62 deaths per 100,000 people, which is equivalent to an additional 7,848 deaths per year in Canada, based on the 2016 population. In RCS models, mean HR predictions increased from the minimum concentration of 2.5 μg/mto 4.5 μg/m, flattened from 4.5 μg/mto 8.0 μg/m, then increased for concentrations above 8.0 μg/m. The threshold model results reflected this pattern with -2 log-likelihood values being equal at 2.5 μg/mand 8.0 μg/m. However, mean threshold model predictions monotonically increased over the concentration range with the lower 95% CI equal to one from 2.5 μg/mto 8.0 μg/m. The RCS model was a superior predictor compared with any of the threshold models, including the linear model. 2.5 2.5 2.5 3 3 3 3 3 3 3 3 3 3 3
UNLABELLED: In the mCCHS cohort analyses inclusion of behavioral covariates did not substantially change the results for both linear and nonlinear models. We examined the sensitivity of the shape of the nonaccidental mortality-PMassociation to removal of person-years at or above the current U.S. and Canadian standards of 12 μg/mand 10 μg/m, respectively. In the full cohort and in both restricted cohorts, a steep increase was observed from the minimum concentration of 2.5 μg/mto 5 μg/m. For the full cohort and the <12 μg/mcohort the relationship flattened over the 5 to 9 μg/mrange and then increased above 9 μg/m. A similar increase was observed for the <10 μg/mcohort followed by a clear decline in the magnitude of predictions over the 5 to 9 μg/mrange and an increase above 9 μg/m. Together these results suggest that a positive association exists for concentrations >9 μg/mwith indications of adverse effects on mortality at concentrations as low as 2.5 μg/m. 2.5 3 3 3 3 3 3 3 3 3 3 3 3
UNLABELLED: Among the other causes of death examined, PMexposures were consistently associated with an increased hazard of mortality due to ischemic heart disease, respiratory disease, cardiovascular disease, and diabetes across all cohorts. Associations were observed in the Stacked CanCHEC but not in all other cohorts for cerebrovascular disease, pneumonia, and chronic obstructive pulmonary disease (COPD) mortality. No significant associations were observed between mortality and exposure to PMfor heart failure, lung cancer, and kidney failure. 2.5 2.5
UNLABELLED: In sensitivity analyses, the addition of Oand Oattenuated associations between PMand mortality. When analyses were stratified by tertiles of copollutants, associations between PMand mortality were only observed in the highest tertile of Oor O. Across broad regions of Canada, linear HR estimates and the shape of the eSCHIF varied substantially, possibly reflecting underlying differences in air pollutant mixtures not characterized by PMmass concentrations or the included gaseous pollutants. Sensitivity analyses to assess regional variation in population characteristics and access to healthcare indicated that the observed regional differences in concentration-mortality relationships, specifically the flattening of the concentration-mortality relationship over the 5 to 9 μg/mrange, was not likely related to variation in the makeup of the cohort or its access to healthcare, lending support to the potential role of spatially varying air pollutant mixtures not sufficiently characterized by PMmass concentrations. 3x2.5 2.5 3x2.5 2.5 3
CONCLUSIONS: In several large, national Canadian cohorts, including a cohort of 7.1 million unique census respondents, associations were observed between exposure to PMwith nonaccidental mortality and several specific causes of death. Associations with nonaccidental mortality were observed using the eSCHIF methodology at concentrations as low as 2.5 μg/m, and there was no clear evidence in the observed data of a lower threshold, below which PMwas not associated with nonaccidental mortality. 2.5 2.5 3
Related papers
Jan '20
Links Between Air Pollution and Death Rates in Areas with Low Pollution: Phase 1
cited by 46 papers
technical report
Sep '22
Long-Term Exposure to Low Levels of Air Pollutants and Its Links to Death and Illness in European Groups
cited by 90 papers
meta-analysis
Jan '13
Air pollution and death rates in multiple Latin American cities
cited by 67 papers
research support, non-u.s. gov't
Jun '09
Long-term traffic air pollution linked to lung and heart death rates in the Netherlands
cited by 125 papers
research support, u.s. gov't, non-p.h.s.
Oct '24
Long-term exposure to tiny outdoor particles and black carbon linked to mortality in Montreal and Toronto, Canada
cited by 4 papers
journal article
Aug '22
How Social Factors Affect Heart Disease Risk from Different Air Pollutants
cited by 14 papers
case reports