低暴露环境中死亡率与空气污染的关联(MAPLE):第2阶段。

M Brauer, J R Brook, T Christidis, Y Chu, D L Crouse, A Erickson, P Hystad, C Li, R V Martin, J Meng, A J Pappin, L L Pinault, M Tjepkema, A van Donkelaar, C Weagle, S Weichenthal, R T Burnett
{"title":"低暴露环境中死亡率与空气污染的关联(MAPLE):第2阶段。","authors":"M Brauer,&nbsp;J R Brook,&nbsp;T Christidis,&nbsp;Y Chu,&nbsp;D L Crouse,&nbsp;A Erickson,&nbsp;P Hystad,&nbsp;C Li,&nbsp;R V Martin,&nbsp;J Meng,&nbsp;A J Pappin,&nbsp;L L Pinault,&nbsp;M Tjepkema,&nbsp;A van Donkelaar,&nbsp;C Weagle,&nbsp;S Weichenthal,&nbsp;R T Burnett","doi":"","DOIUrl":null,"url":null,"abstract":"<p><strong>Introduction: </strong>Mortality is associated with long-term exposure to fine particulate matter (particulate matter ≤2.5 μm in aerodynamic diameter; PM<sub>2.5</sub>), 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 PM<sub>2.5</sub> to 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 PM<sub>2.5</sub> concentrations and nonaccidental mortality were attenuated with the addition of ozone (O<sub>3</sub>) or a measure of gaseous pollutant oxidant capacity (O<sub>x</sub>), which was estimated from O<sub>3</sub> and nitrogen dioxide (NO<sub>2</sub>) concentrations. This was motivated by our interests in understanding both the effects air pollutant mixtures may have on mortality and also the role of O<sub>3</sub> as a copollutant that shares common sources and precursor emissions with those of PM<sub>2.5</sub>. In this Phase 2 report, we further explore the sensitivity of these associations with O<sub>3</sub> and O<sub>x</sub>, evaluate sensitivity to other factors, such as regional variation, and present ambient PM<sub>2.5</sub> concentration-response relationships for specific causes of death.</p><p><strong>Methods: </strong>PM<sub>2.5</sub> concentrations were estimated at 1 km<sup>2</sup> spatial resolution across North America using remote sensing of aerosol optical depth (AOD) combined with chemical transport model (GEOS-Chem) simulations of the AOD:surface PM<sub>2.5</sub> mass concentration relationship, land use information, and ground monitoring. These estimates were informed and further refined with collocated measurements of PM<sub>2.5</sub> and AOD, including targeted measurements in areas of low PM<sub>2.5</sub> concentrations collected at five locations across Canada. Ground measurements of PM<sub>2.5</sub> and 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.</p><p><p>Annual exposures to PM<sub>2.5</sub> were 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.</p><p><p>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<sub>2.5</sub>. 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<sup>3</sup>. 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.</p><p><p>Sensitivity analyses were conducted to examine associations between PM<sub>2.5</sub> and mortality when in the presence of, or stratified by tertile of, O<sub>3</sub> or O<sub>x</sub>. Additionally, associations between PM<sub>2.5</sub> and 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-PM<sub>2.5</sub> association to removal of person-years at or above 12 μg/m<sup>3</sup> (the current U.S. National Ambient Air Quality Standard) and 10 μg/m<sup>3</sup> (the current Canadian and former [2005] World Health Organization [WHO] guideline, and current WHO Interim Target-4). Finally, differences in the shapes of PM<sub>2.5</sub>-mortality associations were assessed across broad geographic regions (airsheds) within Canada.</p><p><strong>Results: </strong>The refined PM<sub>2.5</sub> exposure 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<sub>2.5</sub> <10 μg/m<sup>3</sup>).</p><p><p>Positive associations between outdoor PM<sub>2.5</sub> concentrations and nonaccidental mortality were consistently observed in all cohorts. In the Stacked CanCHEC analyses (1.3 million deaths), each 10-μg/m<sup>3</sup> increase in outdoor PM<sub>2.5</sub> concentration 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 PM<sub>2.5</sub> mass concentration of 4.16 μg/m<sup>3</sup> and 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/m<sup>3</sup> to 4.5 μg/m<sup>3</sup>, flattened from 4.5 μg/m<sup>3</sup> to 8.0 μg/m<sup>3</sup>, then increased for concentrations above 8.0 μg/m<sup>3</sup>. The threshold model results reflected this pattern with -2 log-likelihood values being equal at 2.5 μg/m<sup>3</sup> and 8.0 μg/m<sup>3</sup>. However, mean threshold model predictions monotonically increased over the concentration range with the lower 95% CI equal to one from 2.5 μg/m<sup>3</sup> to 8.0 μg/m<sup>3</sup>. The RCS model was a superior predictor compared with any of the threshold models, including the linear model.</p><p><p>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-PM<sub>2.5</sub> association to removal of person-years at or above the current U.S. and Canadian standards of 12 μg/m<sup>3</sup> and 10 μg/m<sup>3</sup>, respectively. In the full cohort and in both restricted cohorts, a steep increase was observed from the minimum concentration of 2.5 μg/m<sup>3</sup> to 5 μg/m<sup>3</sup>. For the full cohort and the <12 μg/m<sup>3</sup> cohort the relationship flattened over the 5 to 9 μg/m<sup>3</sup> range and then increased above 9 μg/m<sup>3</sup>. A similar increase was observed for the <10 μg/m<sup>3</sup> cohort followed by a clear decline in the magnitude of predictions over the 5 to 9 μg/m<sup>3</sup> range and an increase above 9 μg/m<sup>3</sup>. Together these results suggest that a positive association exists for concentrations >9 μg/m<sup>3</sup> with indications of adverse effects on mortality at concentrations as low as 2.5 μg/m<sup>3</sup>.</p><p><p>Among the other causes of death examined, PM<sub>2.5</sub> exposures 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 PM<sub>2.5</sub> for heart failure, lung cancer, and kidney failure.</p><p><p>In sensitivity analyses, the addition of O<sub>3</sub> and O<sub>x</sub> attenuated associations between PM<sub>2.5</sub> and mortality. When analyses were stratified by tertiles of copollutants, associations between PM<sub>2.5</sub> and mortality were only observed in the highest tertile of O<sub>3</sub> or O<sub>x</sub>. 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 PM<sub>2.5</sub> mass 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/m<sup>3</sup> range, 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 PM<sub>2.5</sub> mass concentrations.</p><p><strong>Conclusions: </strong>In several large, national Canadian cohorts, including a cohort of 7.1 million unique census respondents, associations were observed between exposure to PM<sub>2.5</sub> with 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<sup>3</sup>, and there was no clear evidence in the observed data of a lower threshold, below which PM<sub>2.5</sub> was not associated with nonaccidental mortality.</p>","PeriodicalId":74687,"journal":{"name":"Research report (Health Effects Institute)","volume":" 212","pages":"1-91"},"PeriodicalIF":0.0000,"publicationDate":"2022-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9556709/pdf/hei-2022-212.pdf","citationCount":"0","resultStr":"{\"title\":\"Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE): Phase 2.\",\"authors\":\"M Brauer,&nbsp;J R Brook,&nbsp;T Christidis,&nbsp;Y Chu,&nbsp;D L Crouse,&nbsp;A Erickson,&nbsp;P Hystad,&nbsp;C Li,&nbsp;R V Martin,&nbsp;J Meng,&nbsp;A J Pappin,&nbsp;L L Pinault,&nbsp;M Tjepkema,&nbsp;A van Donkelaar,&nbsp;C Weagle,&nbsp;S Weichenthal,&nbsp;R T Burnett\",\"doi\":\"\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<p><strong>Introduction: </strong>Mortality is associated with long-term exposure to fine particulate matter (particulate matter ≤2.5 μm in aerodynamic diameter; PM<sub>2.5</sub>), 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 PM<sub>2.5</sub> to 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 PM<sub>2.5</sub> concentrations and nonaccidental mortality were attenuated with the addition of ozone (O<sub>3</sub>) or a measure of gaseous pollutant oxidant capacity (O<sub>x</sub>), which was estimated from O<sub>3</sub> and nitrogen dioxide (NO<sub>2</sub>) concentrations. This was motivated by our interests in understanding both the effects air pollutant mixtures may have on mortality and also the role of O<sub>3</sub> as a copollutant that shares common sources and precursor emissions with those of PM<sub>2.5</sub>. In this Phase 2 report, we further explore the sensitivity of these associations with O<sub>3</sub> and O<sub>x</sub>, evaluate sensitivity to other factors, such as regional variation, and present ambient PM<sub>2.5</sub> concentration-response relationships for specific causes of death.</p><p><strong>Methods: </strong>PM<sub>2.5</sub> concentrations were estimated at 1 km<sup>2</sup> spatial resolution across North America using remote sensing of aerosol optical depth (AOD) combined with chemical transport model (GEOS-Chem) simulations of the AOD:surface PM<sub>2.5</sub> mass concentration relationship, land use information, and ground monitoring. These estimates were informed and further refined with collocated measurements of PM<sub>2.5</sub> and AOD, including targeted measurements in areas of low PM<sub>2.5</sub> concentrations collected at five locations across Canada. Ground measurements of PM<sub>2.5</sub> and 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.</p><p><p>Annual exposures to PM<sub>2.5</sub> were 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.</p><p><p>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<sub>2.5</sub>. 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<sup>3</sup>. 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.</p><p><p>Sensitivity analyses were conducted to examine associations between PM<sub>2.5</sub> and mortality when in the presence of, or stratified by tertile of, O<sub>3</sub> or O<sub>x</sub>. Additionally, associations between PM<sub>2.5</sub> and 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-PM<sub>2.5</sub> association to removal of person-years at or above 12 μg/m<sup>3</sup> (the current U.S. National Ambient Air Quality Standard) and 10 μg/m<sup>3</sup> (the current Canadian and former [2005] World Health Organization [WHO] guideline, and current WHO Interim Target-4). Finally, differences in the shapes of PM<sub>2.5</sub>-mortality associations were assessed across broad geographic regions (airsheds) within Canada.</p><p><strong>Results: </strong>The refined PM<sub>2.5</sub> exposure 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<sub>2.5</sub> <10 μg/m<sup>3</sup>).</p><p><p>Positive associations between outdoor PM<sub>2.5</sub> concentrations and nonaccidental mortality were consistently observed in all cohorts. In the Stacked CanCHEC analyses (1.3 million deaths), each 10-μg/m<sup>3</sup> increase in outdoor PM<sub>2.5</sub> concentration 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 PM<sub>2.5</sub> mass concentration of 4.16 μg/m<sup>3</sup> and 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/m<sup>3</sup> to 4.5 μg/m<sup>3</sup>, flattened from 4.5 μg/m<sup>3</sup> to 8.0 μg/m<sup>3</sup>, then increased for concentrations above 8.0 μg/m<sup>3</sup>. The threshold model results reflected this pattern with -2 log-likelihood values being equal at 2.5 μg/m<sup>3</sup> and 8.0 μg/m<sup>3</sup>. However, mean threshold model predictions monotonically increased over the concentration range with the lower 95% CI equal to one from 2.5 μg/m<sup>3</sup> to 8.0 μg/m<sup>3</sup>. The RCS model was a superior predictor compared with any of the threshold models, including the linear model.</p><p><p>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-PM<sub>2.5</sub> association to removal of person-years at or above the current U.S. and Canadian standards of 12 μg/m<sup>3</sup> and 10 μg/m<sup>3</sup>, respectively. In the full cohort and in both restricted cohorts, a steep increase was observed from the minimum concentration of 2.5 μg/m<sup>3</sup> to 5 μg/m<sup>3</sup>. For the full cohort and the <12 μg/m<sup>3</sup> cohort the relationship flattened over the 5 to 9 μg/m<sup>3</sup> range and then increased above 9 μg/m<sup>3</sup>. A similar increase was observed for the <10 μg/m<sup>3</sup> cohort followed by a clear decline in the magnitude of predictions over the 5 to 9 μg/m<sup>3</sup> range and an increase above 9 μg/m<sup>3</sup>. Together these results suggest that a positive association exists for concentrations >9 μg/m<sup>3</sup> with indications of adverse effects on mortality at concentrations as low as 2.5 μg/m<sup>3</sup>.</p><p><p>Among the other causes of death examined, PM<sub>2.5</sub> exposures 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 PM<sub>2.5</sub> for heart failure, lung cancer, and kidney failure.</p><p><p>In sensitivity analyses, the addition of O<sub>3</sub> and O<sub>x</sub> attenuated associations between PM<sub>2.5</sub> and mortality. When analyses were stratified by tertiles of copollutants, associations between PM<sub>2.5</sub> and mortality were only observed in the highest tertile of O<sub>3</sub> or O<sub>x</sub>. 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 PM<sub>2.5</sub> mass 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/m<sup>3</sup> range, 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 PM<sub>2.5</sub> mass concentrations.</p><p><strong>Conclusions: </strong>In several large, national Canadian cohorts, including a cohort of 7.1 million unique census respondents, associations were observed between exposure to PM<sub>2.5</sub> with 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<sup>3</sup>, and there was no clear evidence in the observed data of a lower threshold, below which PM<sub>2.5</sub> was not associated with nonaccidental mortality.</p>\",\"PeriodicalId\":74687,\"journal\":{\"name\":\"Research report (Health Effects Institute)\",\"volume\":\" 212\",\"pages\":\"1-91\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2022-07-01\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9556709/pdf/hei-2022-212.pdf\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Research report (Health Effects Institute)\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Research report (Health Effects Institute)","FirstCategoryId":"1085","ListUrlMain":"","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

摘要

导言:死亡率与长期暴露于细颗粒物(空气动力学直径≤2.5 μm;PM2.5),尽管在较低浓度下,这些关联的程度和形式仍然知之甚少。知识差距包括浓度-反应曲线的形状和明显增加风险的最低接触水平,以及与特定死亡原因的关联的发生和程度。在这里,我们将环境PM2.5暴露的改进估计应用于加拿大以全国人口为基础的队列,包括1991年、1996年或2001年人口普查的710万人的队列。低暴露水平下非意外死亡率和几种特定死亡原因的浓度-反应关系的形状特征是低暴露环境中死亡率-空气污染关联(MAPLE)第1阶段报告的重点。在第一阶段报告中,我们报道了室外PM2.5浓度与非意外死亡率之间的关联随着臭氧(O3)的增加或气体污染物氧化能力(Ox)的测量而减弱,这是由O3和二氧化氮(NO2)浓度估计的。这是出于我们的兴趣,了解空气污染物混合物可能对死亡率的影响,以及O3作为一种与PM2.5具有共同来源和前体排放的共污染物的作用。在第二阶段报告中,我们进一步探讨了这些与O3和Ox相关的敏感性,评估了对其他因素(如区域差异)的敏感性,并提出了环境PM2.5浓度与特定死亡原因的响应关系。方法:利用气溶胶光学深度(AOD)遥感数据,结合化学输运模式(GEOS-Chem)模拟AOD、地表PM2.5质量浓度关系、土地利用信息和地面监测数据,估算北美地区1 km2空间分辨率下PM2.5浓度。这些估计是通过PM2.5和AOD的同时测量得到的,包括在加拿大五个地点收集的PM2.5浓度较低的地区的目标测量。研究人员利用1981年至1999年PM2.5和总悬浮颗粒物(TSP)质量浓度的地面测量数据,反演了同一时期基于遥感的估算数据,得出了1981年至2016年的年地表模型。然后,根据所得税档案中的年度邮政编码条目得出的居住历史,对加拿大几个以国家人口为基础的队列中受试者的PM2.5年暴露量进行了估计。这些队列包括三个基于人口普查的队列:1991年加拿大人口普查健康和环境队列(CanCHEC);250万受访者),1996年的CanCHEC(300万受访者),2001年的CanCHEC(300万受访者),以及排除了受访者重复记录的堆叠CanCHEC(堆叠CanCHEC;710万受访者)。加拿大社区健康调查(CCHS)死亡率队列(mCCHS)来自多个CCHS合并周期(540,900名受访者),包括有关健康行为的额外个人信息。所有队列的随访期均完成至2016年底。使用10年移动平均暴露和1年滞后估计非意外死亡和其他主要死亡原因的Cox比例风险比(HRs)。所有模型均按年龄、性别、移民身份以及适当的普查年份或调查周期进行分层。模型进一步调整了收入充足五分位数、可见少数民族地位、土著身份、受教育程度、劳动力状况、婚姻状况、职业,以及社区规模、空气密度、城市形态等生态协变量,以及加拿大边缘化指数(canmarg;不稳定、贫困、依赖和种族集中)。mCCHS分析还针对吸烟、饮酒、水果和蔬菜消费、身体质量指数(BMI)和运动行为的个人水平测量进行了调整。除了线性模型外,还使用限制三次样条(RCS)研究了浓度-响应函数的形状。通过最小化贝叶斯信息准则(BIC)来选择节点数量。另外两个模型被用来检验非意外死亡率和PM2.5之间的关系。第一种是标准阈值模型,如果浓度低于特定阈值,则浓度转换为零,如果浓度高于阈值,则浓度减去阈值。第二个附加模型是形状受限健康影响函数(SCHIF)的扩展,即eSCHIF,它将RCS预测转换为可能更适合用于健康影响评估的函数。 根据RCS参数估计值及其协方差矩阵,从最小浓度到最大浓度,以0.1 μg/m3的增量模拟了1000种RCS的实现。然后,一个eSCHIF适合于这些RCS实现中的每一个。因此,在总范围内的每个浓度下确定了1000个eSCHIF预测和不确定区间。在O3或Ox存在或按三分位数分层的情况下,进行敏感性分析以检查PM2.5与死亡率之间的关系。此外,评估PM2.5与死亡率之间的关系对较低浓度阈值的敏感性,其中低于阈值的人年被分配为该组内的平均暴露量。我们还研究了非意外死亡率- pm2.5相关性的形状对去除人年≥12 μg/m3(现行的美国国家环境空气质量标准)和10 μg/m3(现行的加拿大和前[2005]世界卫生组织指南,以及现行的世卫组织中期目标-4)的敏感性。最后,评估了加拿大境内不同地理区域(大气区)pm2.5死亡率关联形态的差异。结果:与先前和MAPLE第一阶段报告中应用的估计相比,改进的PM2.5暴露估计显示出更好的性能,误差略有减少,包括在较低的浓度范围内(例如PM2.5 3)。在所有队列中一致观察到室外PM2.5浓度与非意外死亡率之间的正相关。在堆叠CanCHEC分析(130万例死亡)中,室外PM2.5浓度每增加10 μg/m3,非意外死亡率的HR为1.084(95%可信区间[CI]: 1.073至1.096)。PM2.5质量浓度四分位数范围(IQR)增加4.16 μg/m3,年平均非意外死亡率为92.8 / 10,000人(1991-2016年期间25-90岁队列参与者),这一HR相当于每10万人中额外的31.62人死亡,相当于加拿大每年额外的7,848人死亡(基于2016年人口)。在RCS模型中,平均HR预测值从最低浓度2.5 μg/m3增加到4.5 μg/m3,从4.5 μg/m3降至8.0 μg/m3,然后在浓度高于8.0 μg/m3时增加。阈值模型结果反映了这一模式,在2.5 μg/m3和8.0 μg/m3时,-2对数似然值相等。然而,在浓度范围内,平均阈值模型预测单调增加,95% CI的下限为1,从2.5 μg/m3到8.0 μg/m3。与任何阈值模型(包括线性模型)相比,RCS模型是一个优越的预测器。在mCCHS队列分析中,纳入行为协变量并没有实质性地改变线性和非线性模型的结果。我们研究了非意外死亡率- pm2.5相关性的形状对去除人年的敏感性,这些人年分别达到或高于当前美国和加拿大的标准12 μg/m3和10 μg/m3。在全队列和两个受限队列中,观察到从最低浓度2.5 μg/m3急剧增加到5 μg/m3。对于全队列和3队列,在5 ~ 9 μg/m3范围内,关系趋于平缓,然后在9 μg/m3以上增加。在3队列中观察到类似的增加,随后在5至9 μg/m3范围内的预测幅度明显下降,在9 μg/m3以上的预测幅度增加。总之,这些结果表明,浓度>9 μg/m3与低至2.5 μg/m3的浓度对死亡率的不利影响呈正相关。在研究的其他死亡原因中,PM2.5暴露始终与缺血性心脏病、呼吸系统疾病、心血管疾病和糖尿病导致的死亡风险增加有关。在堆叠CanCHEC中观察到脑血管疾病、肺炎和慢性阻塞性肺疾病(COPD)死亡率的关联,但在所有其他队列中没有观察到。心力衰竭、肺癌和肾衰竭的死亡率与暴露于PM2.5之间没有显著关联。在敏感性分析中,O3和Ox的加入减弱了PM2.5与死亡率之间的关联。当按共污染物的分位数进行分层分析时,PM2.5与死亡率之间的关联仅在O3或Ox的最高分位数中观察到。在加拿大的广大地区,线性HR估计和eSCHIF的形状变化很大,可能反映了空气污染物混合物的潜在差异,而不是以PM2.5质量浓度或所包括的气态污染物为特征。 评估人口特征和获得医疗保健的区域差异的敏感性分析表明,观察到的浓度-死亡率关系的区域差异,特别是浓度-死亡率关系在5至9 μg/m3范围内趋于平缓,不太可能与队列组成或其获得医疗保健的机会的差异有关。这为PM2.5质量浓度未充分表征的空间变化的空气污染物混合物的潜在作用提供了支持。结论:在几个大型的加拿大国家队列中,包括710万人口普查受访者的队列,观察到PM2.5暴露与非意外死亡和几种特定死亡原因之间的关联。使用eSCHIF方法观察到PM2.5浓度低至2.5 μg/m3时与非意外死亡率的关联,在观察到的数据中没有明确的证据表明存在更低的阈值,低于该阈值PM2.5与非意外死亡率无关。
本文章由计算机程序翻译,如有差异,请以英文原文为准。

Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE): Phase 2.

Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE): Phase 2.

Mortality-Air Pollution Associations in Low Exposure Environments (MAPLE): Phase 2.

Introduction: Mortality is associated with long-term exposure to fine particulate matter (particulate matter ≤2.5 μm in aerodynamic diameter; PM2.5), 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 PM2.5 to 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 PM2.5 concentrations and nonaccidental mortality were attenuated with the addition of ozone (O3) or a measure of gaseous pollutant oxidant capacity (Ox), which was estimated from O3 and nitrogen dioxide (NO2) concentrations. This was motivated by our interests in understanding both the effects air pollutant mixtures may have on mortality and also the role of O3 as a copollutant that shares common sources and precursor emissions with those of PM2.5. In this Phase 2 report, we further explore the sensitivity of these associations with O3 and Ox, evaluate sensitivity to other factors, such as regional variation, and present ambient PM2.5 concentration-response relationships for specific causes of death.

Methods: PM2.5 concentrations were estimated at 1 km2 spatial resolution across North America using remote sensing of aerosol optical depth (AOD) combined with chemical transport model (GEOS-Chem) simulations of the AOD:surface PM2.5 mass concentration relationship, land use information, and ground monitoring. These estimates were informed and further refined with collocated measurements of PM2.5 and AOD, including targeted measurements in areas of low PM2.5 concentrations collected at five locations across Canada. Ground measurements of PM2.5 and 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.

Annual exposures to PM2.5 were 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.

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 PM2.5. 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/m3. 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.

Sensitivity analyses were conducted to examine associations between PM2.5 and mortality when in the presence of, or stratified by tertile of, O3 or Ox. Additionally, associations between PM2.5 and 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-PM2.5 association to removal of person-years at or above 12 μg/m3 (the current U.S. National Ambient Air Quality Standard) and 10 μg/m3 (the current Canadian and former [2005] World Health Organization [WHO] guideline, and current WHO Interim Target-4). Finally, differences in the shapes of PM2.5-mortality associations were assessed across broad geographic regions (airsheds) within Canada.

Results: The refined PM2.5 exposure 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 PM2.5 <10 μg/m3).

Positive associations between outdoor PM2.5 concentrations and nonaccidental mortality were consistently observed in all cohorts. In the Stacked CanCHEC analyses (1.3 million deaths), each 10-μg/m3 increase in outdoor PM2.5 concentration 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 PM2.5 mass concentration of 4.16 μg/m3 and 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/m3 to 4.5 μg/m3, flattened from 4.5 μg/m3 to 8.0 μg/m3, then increased for concentrations above 8.0 μg/m3. The threshold model results reflected this pattern with -2 log-likelihood values being equal at 2.5 μg/m3 and 8.0 μg/m3. However, mean threshold model predictions monotonically increased over the concentration range with the lower 95% CI equal to one from 2.5 μg/m3 to 8.0 μg/m3. The RCS model was a superior predictor compared with any of the threshold models, including the linear model.

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-PM2.5 association to removal of person-years at or above the current U.S. and Canadian standards of 12 μg/m3 and 10 μg/m3, respectively. In the full cohort and in both restricted cohorts, a steep increase was observed from the minimum concentration of 2.5 μg/m3 to 5 μg/m3. For the full cohort and the <12 μg/m3 cohort the relationship flattened over the 5 to 9 μg/m3 range and then increased above 9 μg/m3. A similar increase was observed for the <10 μg/m3 cohort followed by a clear decline in the magnitude of predictions over the 5 to 9 μg/m3 range and an increase above 9 μg/m3. Together these results suggest that a positive association exists for concentrations >9 μg/m3 with indications of adverse effects on mortality at concentrations as low as 2.5 μg/m3.

Among the other causes of death examined, PM2.5 exposures 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 PM2.5 for heart failure, lung cancer, and kidney failure.

In sensitivity analyses, the addition of O3 and Ox attenuated associations between PM2.5 and mortality. When analyses were stratified by tertiles of copollutants, associations between PM2.5 and mortality were only observed in the highest tertile of O3 or Ox. 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 PM2.5 mass 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/m3 range, 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 PM2.5 mass concentrations.

Conclusions: In several large, national Canadian cohorts, including a cohort of 7.1 million unique census respondents, associations were observed between exposure to PM2.5 with 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/m3, and there was no clear evidence in the observed data of a lower threshold, below which PM2.5 was not associated with nonaccidental mortality.

求助全文
通过发布文献求助,成功后即可免费获取论文全文。 去求助
来源期刊
自引率
0.00%
发文量
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信