- Open Open Peer Review
Risk estimates of mortality attributed to low concentrations of ambient fine particulate matter in the Canadian community health survey cohort
Environmental Healthvolume 15, Article number: 18 (2016)
Understanding the shape of the relationship between long-term exposure to ambient fine particulate matter (PM2.5) concentrations and health risks is critical for health impact and risk assessment. Studies evaluating the health risks of exposure to low concentrations of PM2.5 are limited. Further, many existing studies lack individual-level information on potentially important behavioural confounding factors.
A prospective cohort study was conducted among a subset of participants in a cohort that linked respondents of the Canadian Community Health Survey to mortality (n = 299,500) with satellite-derived ambient PM2.5 estimates. Participants enrolled between 2000 and 2008 were followed to date of death or December 31, 2011. Cox proportional hazards models were used to estimate hazard ratios (HRs) for mortality attributed to PM2.5 exposure, adjusted for individual-level and contextual covariates, including smoking behaviour and body mass index (BMI).
Approximately 26,300 non-accidental deaths, of which 32.5 % were due to circulatory disease and 9.1 % were due to respiratory disease, occurred during the follow-up period. Ambient PM2.5 exposures were relatively low (mean = 6.3 μg/m3), yet each 10 μg/m3 increase in exposure was associated with increased risks of non-accidental (HR = 1.26; 95 % CI: 1.19-1.34), circulatory disease (HR = 1.19; 95 % CI: 1.07–1.31), and respiratory disease mortality (HR = 1.52; 95 % CI: 1.26–1.84) in fully adjusted models. Higher hazard ratios were observed for respiratory mortality among respondents who never smoked (HR = 1.97; 95 % CI: 1.24–3.13 vs. HR = 1.45; 95 % CI: 1.17–1.79 for ever smokers), and among obese (BMI ≥ 30) respondents (HR = 1.76; 95 % CI: 1.15-2.69 vs. HR = 1.41; 95 % CI: 1.04–1.91 for normal weight respondents), though differences between groups were not statistically significant. A threshold analysis for non-accidental mortality estimated a threshold concentration of 0 μg/m3 (+95 % CI = 4.5 μg/m3).
Increased risks of non-accidental, circulatory, and respiratory mortality were observed even at very low concentrations of ambient PM2.5. HRs were generally greater than most literature values, and adjusting for behavioural covariates served to reduce HR estimates slightly.
Ambient fine particulate air pollution (PM2.5) is known to contribute to cardiovascular and respiratory morbidity, and is recognized as an important contributor to global disease burden . Recent estimates from Global Burden of Disease suggest that ambient air pollution was responsible for nearly 2.9 million deaths per year in 2013 . While ambient PM2.5 concentrations in Canada are generally below national and international guidelines, analyses from the 1991 Canadian Census Health and Environment Cohort (CanCHEC) suggest that long-term exposure to PM2.5 in Canada (mean = 8.9 μg/m3) may contribute to non-accidental and cardiovascular mortality . However, that study did not include individual-level information on potentially important confounding factors such as smoking and obesity and applied an indirect approach to adjust for these and other factors [3, 4]. Analysis of the United States Agricultural Health Study (AHS) cohort also supports an association between cardiovascular mortality and long-term exposure to low concentrations of ambient PM2.5 (mean = 9.2 μg/m3) . Moreover, a recent meta-analysis of studies conducted in North America and internationally supports an association between long-term exposure to PM2.5 and mortality, with the strongest association observed for cardiovascular mortality .
The WHO PM2.5 guideline of 10 μg/m3 was based on the lower end of the exposure distribution in previous studies , though there are few studies that have evaluated concentration-response associations at very low exposures. The Global Burden of Disease 2010 study  developed a mortality risk model for PM2.5 over the global range of concentrations. This model incorporated a counterfactual uncertainty distribution, below which no excess risk was assumed, and was specified by a uniform distribution between 5.8 μg/m3 and 8.8 μg/m3. This uncertainty distribution was selected based on the lack of empirical evidence of any statistical association between ambient PM2.5 and mortality below their counterfactual distribution. These concentrations represent the 48.9th and 79.9th percentiles of the exposure distribution in this study, respectively. Therefore, it is of interest to examine the shape of the concentration-mortality association at these very low concentrations, as well as the statistical strength of evidence for such an association.
In this study, we examine the relationship between long-term exposure to ambient PM2.5 and non-accidental, respiratory, and cardiovascular mortality in the Canadian Community Health Survey (CCHS) cohort. Participants in this cross-sectional survey were enrolled across Canada between 2000 and 2008 and provided detailed individual-level information on potentially important confounding factors (e.g. smoking, obesity) that were not available for the previous analysis of PM2.5 and mortality in the CanCHEC study . As such, the primary aim of this study was to examine the relationship between very low concentrations of PM2.5 (mean = 6.3 μg/m3) and different causes of mortality in Canada and the impacts of adjusting for potential confounding factors. Finally, an improved, finer-scale, satellite-derived exposure model for PM2.5 (i.e., a 1 km2 grid) was used to reduce exposure misclassification.
The CCHS is a national, cross-sectional survey providing information about the health, behaviours, and health care use of the non-institutional Canadian population aged 12 or older. The survey excludes full-time members of the Canadian Armed Forces and residents of Indian reserves and certain remote areas. Exclusions represent less than 3 % of the target population of Canada . The annual component of the CCHS was conducted every two years from 2000/01 to 2007, after which the survey was conducted on an annual basis. The CCHS response rates are as follows: 84.7 % in Cycle 1.1 (2000/01), 80.7 % in Cycle 2.1 (2003), 78.9 % in Cycle 3.1 (2005), 77.6 % in 2007, and 75.0 % in 2008 . CCHS respondents were eligible for the CCHS-mortality cohort if they gave permission to share and link their information with other administrative datasets; 86.0 % of CCHS respondents agreed to the linkage.
The Canadian Mortality Database (CMDB) is a national database that contains all deaths registered in Canada since 1950. Deaths that occurred between January 1, 2000 and December 31, 2011 were eligible for linkage. The CMDB includes data on underlying cause of death and date of death.
The Historical Tax Summary File (HTSF) is a database of annual tax returns that represent all individuals who received a tax declaration in a given year. Tax years between 1996 and 2011 were eligible for linkage. The HTSF includes postal codes, names, and dates of death (if applicable).
The creation of the CCHS-Mortality Cohort was conducted in two steps. First, using a probabilistic linkage methodology based on the Fellegi-Sunter theory of record linkage , eligible CCHS respondents were linked to the HTSF (using date of birth, sex, name, and postal code), in order to capture these variables and date of death, as reported on tax files between 1996 and 2011. Approximately 85 % of eligible CCHS respondents were linked to the HTSF. Alternative postal codes and names were captured through this initial linkage and were used in the subsequent linkage to the CMDB, to improve linkage results. Second, all eligible CCHS respondents (regardless of whether they were linked to the HTSF) were also linked to the CMDB (which included underlying cause of death), using standard probabilistic linkage techniques (as described above) and followed for mortality from cohort entry (i.e., date of CCHS interview) to December 31, 2011.
A total of 457,300 eligible CCHS-mortality respondents were included, with 117,800 respondents in Cycle 1, 112,900 respondents in Cycle 2, 113,900 respondents in Cycle 3, and 112,700 respondents in 2007/08. CCHS respondents who were first linked to the HTSF had a greater probability to be linked to the CMDB since additional data in the HTSF (e.g., alternate postal codes, name, and date of death), were used in the probabilistic linkage. In order to reduce the probability of false-negative links, we excluded those CCHS respondents who were not linked to the HTSF (n = 69,300 respondents excluded) (Additional file 1).
Since the purpose of this analysis was to evaluate long-term effects of air pollution exposure, the study population was restricted to adults aged 25 to 90 years of age at enrollment (n = 72,000 respondents excluded). Adults older than 90 years of age were excluded from this study to ensure a sufficient sample size within all age strata. Similar to the CanCHEC study , immigrants living in Canada for less than 20 years (i.e., those who had arrived in Canada less than 20 years before the start date), were excluded from this study (n = 13,200 additional respondents excluded) for the following reasons. Immigrants are known to have better health and live longer than the Canadian-born population . Immigrants also more frequently live in areas of greater ambient air pollution (unpublished data), and their exposure to air pollution prior to living in Canada is largely unknown. Cause-specific mortality analyses among recent immigrants were also not meaningful due to small sample sizes in the CCHS cohort (i.e., < 250 deaths). Therefore, the use of a larger cohort would be necessary to examine the health effects of air pollution on recent immigrant populations. Finally, we excluded an additional 3,400 respondents who were not linked to air pollution estimates since they live beyond the boundaries of the air pollution models (i.e., in the northern Territories) (Additional file 1). The final analytical sample was 299,500 respondents (note slight inconsistencies due to rounding). All research using human data was carried out at Statistics Canada in accordance with the Statistics Act to meet standards of privacy and confidentiality associated with the internal use of survey data. The record linkage project was approved by the Executive Management Board at Statistics Canada (ref. num. 003–2015).
The place of residence of respondents at the date of entry into the Cohort was mapped in Geographic Information Systems (ArcGIS v.10; ESRI 2010) through the use of Statistics Canada’s Postal Code Conversion File plus (PCCF+) V.6B, which assigns geographic coordinates to postal codes based on a population-weighted random allocation algorithm . Respondent locations were then spatially linked to estimates from a surface layer of PM2.5 concentration derived by relating total column aerosol optical depth retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument to near-surface PM2.5 using the GEOS-Chem chemical transport model. Geographically weighted regression, which includes ground monitoring data and land use information, was subsequently applied to these estimates to produce average PM2.5 concentrations at a 0.01° × 0.01° (approximately 1 km2) resolution from 2004 to 2012 . These models included coverage for nearly all of mainland North America. These estimates were extended to 1998 to 2003 using the inter-annual variation of Boys et al. (2014) , who inferred global PM2.5 trends at 0.1° × 0.1° resolution using satellites from 1998 to 2012. Average PM2.5 levels were strongly correlated with ground-level observations in North America (R2 = 0.82, slope = 0.97; n = 1440) . Outliers that included PM2.5 values >20 ug/m3 were excluded from analysis (<1 % of respondents were excluded in this manner in any year). These outliers were likely due to inaccurate estimates of aerosol optical depth from satellite retrievals. For each year in the cohort, respondents were assigned a PM2.5 value corresponding to the mean of the three previous years to the follow-up year; therefore, exposure always preceded response. For example, for the follow-up year 2001, we assigned the mean PM2.5 estimates from 1998 to 2000.
Covariates and statistical methods
Standard Cox proportional hazards models  were used for survival analysis of non-accidental and cause-specific mortality within the cohort, from the date of interview for the CCHS to either the date of death recorded in the CMDB or the final date of the linkage project (i.e., 31 December, 2011). All models were stratified by sex and age (5-year intervals). Socioeconomic covariates included: immigrant status, visible minority status, Aboriginal status, and marital status, educational attainment, income adequacy quintile, and employment status (Table 1). Visible minority status was defined as in the Employment Equity Act, as “persons, other than Aboriginal peoples, who are non-Caucasian in race or non-white in colour” . Income adequacy quintiles were calculated based on the ratio of household income to the low-income cut-off for their household and community size. Low-income cut-offs represent families that spend more than 20 % of their income on food, shelter and clothing, and are adjusted for size of family and area of residence .
Neighbourhood socioeconomic status, including both social and material deprivation, contributes to increased risk of mortality in Canadian cities, although the presence of immigrants can reduce mortality risk . Ecological (contextual) covariates were derived from the long-form Canadian Census at the Census Division (CD) and Dissemination Area (DA) geographic scale, from the 2001 Census for respondents interviewed between 2000 and 2003, and the 2006 Census for respondents interviewed during or after 2004. Census Divisions are a subdivision of the provinces and territories that usually represent communities, regional districts, or several neighbouring municipalities, and range in size from several thousand to a few million persons . Dissemination Areas are the smallest geographical unit used by the Census and are delineated based on population counts based on the previous census, to target a population of 400–700 persons . There were 288 CDs and 54,623 DAs in Canada as of 2006 . These contextual covariates were then linked to individual respondents through a common geographic identifier (i.e., a numeric code identifying the DA or CD). For each CD and DA, the proportion of recent immigrants (<5 years residency in Canada), educational attainment (the proportion of persons aged 15 years or older who had not graduated from high school) and low income (the proportion of persons below the low-income cut-off) were derived for both Census years . The proportion of recent immigrants in a region may provide a health benefit in the form of social inclusion if the resident is a member of a unified community, though it also may represent social deprivation, since recent immigrants also include persons of very low SES upon arrival in Canada (e.g., refugees or temporary workers). The other two ecological covariates (educational attainment and low income) provide a more direct estimate of neighbourhood socioeconomic status. Although broader geographic scales such as Census Tracts (CTs) are more often used to derive neighbourhood contextual variables , CTs were not available for rural areas. Neighbourhood covariates were therefore calculated by taking the difference between CD and DA estimates. It was expected a priori that the ecological covariates would attenuate risk estimates as in previous work on CanCHEC .
In addition to the socioeconomic and ecological covariates, this study included four health status/behavioural covariates. Body Mass Index (BMI) was derived from the self-reported height and weight of respondents, and adjusted using correction factors that were developed for the CCHS to account for self-reporting bias in BMI data . The International Standard Classification was used to categorize Body Mass Index , with obesity subdivided into two categories (i.e., BMI 30 – 34.9 and BMI ≥ 35) to further differentiate health risks among obese persons within the study. Smoking behaviour was categorized as never, former, or current smokers. Detailed data on smoking behaviour (e.g., number of cigarettes smoked per day) were available only for daily smokers (ca. 21.3 % of respondents) and were therefore not included. Fruit and vegetable daily consumption and alcohol consumption were also included, as in previous studies  (Table 1).
Survival models were examined in a sequential manner by adding all of the socioeconomic covariates in a single model, then adding in the ecological covariates to the socioeconomic models, and finally by adding the behavioural covariates to create fully adjusted models for non-accidental mortality (ICD-10 codes A-R) and mortality attributed to circulatory disease (ICD-10: I00–I99, with and without diabetes, E10–E14), including the subgroups of ischemic heart disease (ICD-10: I20–I25), and cerebrovascular disease (ICD-10: I60–I69). We also considered models for mortality due to respiratory disease (ICD-10: J00–J99), also including chronic obstructive pulmonary disease (ICD-10: J19–J46), and lung cancer (ICD-10: C33–C34). We also examined a model of socioeconomic and behavioural covariates, excluding ecological covariates. We added groups of variables in this manner to specifically examine the influence of including the behavioural variables to a model which included both socioeconomic and ecological variables, as were available in previous cohort studies in Canada . Effect modification by sex, smoking behaviour (ever smoked vs. never smoked), BMI (obese: BMI ≥ 30 and obese II: BMI ≥ 35 vs. normal weight: BMI = 18.5–25), fruit and vegetable consumption (<5 servings vs. ≥ 5 servings), alcohol consumption (regular drinker vs. occasional/never/former drinker), and age (<75 years vs. ≥ 75 years) were also evaluated in separate Cox proportional hazards models, and Cochran’s Q-statistic heterogeneity tests were used to evaluate significant differences in HRs among groups . These covariates were chosen for effect modification analysis due to known physiological differences between these groups of respondents, and interest in previous studies .
To examine the shape of the relationship between non-accidental mortality hazard ratio (HR) and air pollution exposure, we fitted spline-based HR curves using the smoothing method in the R package “smoothHR” on the fully adjusted model . The package uses a combination of AIC and BIC to determine the optimal degrees of freedom to use in the model . We also estimated the PM2.5 threshold concentration (T) by fitting Cox proportional hazards models to a series of newly defined PM2.5 based variables of the form: PM2.5 (T) = PM2.5 – T; if PMi > T and 0 otherwise, for T = 1 to 10. Our estimate of T is the concentration corresponding to the largest (−2) log-likelihood value (−2LL) obtained from the Cox model. Ninety-five percent confidence intervals on T were based on changes in -2LL of 3.84 units.
All descriptive statistics reported from the survey were rounded to the nearest hundred for institutional confidentiality reasons.
A total of 299,500 respondents were included in the study after excluding respondents who were not linked to a tax file, respondents who were not within the 25 to 90 year age range and were not recent immigrants (i.e., < 20 years in Canada), and respondents who were not linked to air pollution estimates. Respondents were followed for mortality for up to 12 years after cohort entry (mean follow-up period (± SD) was 7.6 ± 2.7 years). The mean exposure (± SD) of respondents to PM2.5 estimated from the 3-year moving average was 6.3 ± 2.5 μg/m3. The PM2.5 person-year exposure percentiles within the final study cohort were: minimum: 1.0 μg/m3, 5th: 3.0 μg/m3, 25th: 4.2 μg/m3, median: 5.9 μg/m3, 75th: 8.3 μg/m3, 95th: 11.3 μg/m3, and maximum: 13.0 μg/m3. In large cities (metropolitan pop. > 1 million), PM2.5 estimates were generally greater than in surrounding areas, and there were areas of the downtown core exceeding 8 ug/m3 in all of these cities (Fig. 1). Mean PM2.5 exposure increased incrementally by decreasing income quintile and was highest for respondents in the poorest income quintile (Table 1). PM2.5 exposure was also greatest for the most highly educated respondents (Table 1). Obese respondents were exposed to less air pollution than those of increasingly lower weight classes, with the greatest exposure among respondents classified as underweight (Table 1). Hazard ratios for non-accidental mortality were calculated for all variables and ecological covariates (Table 1). Among ecological covariates for DAs and CDs, the proportion of recent immigrants, high school graduates and low income families were positively correlated with average PM2.5 air pollution exposure (Table 2). The proportion of recent immigrants was protective for mortality at the broader landscape level (i.e., the CT), though increased HRs at the neighbourhood scale (i.e., DAs) (Table 1). Associations between all combinations of the covariates are provided in Additional file 2.
Separate Cox proportional hazards models were run for all covariates in the fully adjusted models. Immigrant status, greater educational attainment, higher income, being overweight or obese (type I), and increasing consumption of fruits and vegetables were all associated with a lower risk of non-accidental mortality (Table 1). Aboriginal status, being unmarried, being underweight or obese (type II), not employed, smoking, and not regularly drinking alcohol were associated with a greater risk of non-accidental mortality (Table 1).
Covariates were added in a stepwise manner to a Cox proportional hazards model for non-accidental mortality to assess their contribution to the model (Table 3). In general, the addition of socioeconomic covariates improved the model fit, and resulted in a significantly increased HR from the unadjusted model (Table 3; Cochran’s Q = 4.29; p = 0.04). The additional of behavioural covariates to the socioeconomic model reduced HRs somewhat, though not significantly (Table 3; Cochran’s Q = 0.23, p = 0.63). The addition of ecological covariates to the socioeconomic model, particularly the percentage of recent immigrants and high school graduates, also improved model fit and significantly increased HRs (Table 3; Cochran’s Q = 27.30, p < 0.01). The addition of behavioural covariates to create a fully adjusted model also improved model fit, though the HRs declined non-significantly from the second adjusted model (Table 3; Cochran’s Q = 2.41, p = 0.12).
Table 4 presents the HRs and 95 % CI for Cox proportional hazard models for non-accidental mortality and mortality due to circulatory or respiratory causes. In the fully adjusted model, HR estimates for non-accidental mortality were 1.26 (95 % C.I.: 1.19–1.34) per 10 μg/m3 increase in ambient PM2.5. The strongest association was observed for respiratory disease mortality, with an HR of 1.52 (95 % C.I.: 1.26–1.84) per 10 μg/m3 increase in ambient PM2.5. In fully adjusted models, HRs were significantly greater than one for all causes of death except cerebrovascular disease and lung cancer, though the HRs were significant in the models that did not include behavioural covariates (Table 4). For all causes of death, HRs were greater in the fully adjusted model than in the unadjusted model, though were reduced after adding behavioural covariates (Table 4).
The results of effect modification by sex, age, BMI (i.e., obese vs. normal weight), fruit and vegetable consumption (i.e., < 5 or ≥ 5 daily servings), smoking (i.e., ever smoked vs. never smoked) and alcohol consumption are presented in Table 5. In a fully adjusted model, the HR for non-accidental mortality among men was 1.34 (95 % C.I.: 1.24–1.46) per 10 μg/m3 increase in ambient PM2.5 and was significantly greater than that of women (Cochran’s Q; Table 5). The HRs for circulatory and respiratory disease mortality among men were also greater than among women, though the differences in HRs were not statistically significant (Cochran’s Q; Table 5). None of the other comparisons among groups were statistically significant (Table 5). However, the HR of respiratory disease mortality was particularly high among never smokers (HR: 1.97; 95 % CI: 1.23–3.13 per 10 μg/m3 increase in PM2.5) and among obese respondents (HR = 1.76, 95 % CI: 1.15–2.69 per 10 μg/m3 increase in PM2.5) (Table 5).
We fitted a nonparametric smoothing (spline) to examine the shape of the association between exposure and non-accidental mortality within the fully adjusted model. The relationship between the logarithm of the hazard function and PM2.5 is presented in Fig. 2 in addition to its 95 % confidence intervals. We specified a reference concentration of 1 μg/m3 which forces the predicted log-hazard function to equal 0 at the reference level. The smoothed curve generally increased with increasing concentration, however the confidence intervals are relatively wide making it difficult to speculate on a specific shape of the concentration-mortality association based on this graphical evidence. Our estimate of the threshold concentration was 0 μg/m3 with an upper 95 % CI value of 4.5 μg/m3.
Within our cohort, exposure to PM2.5 assigned at baseline was associated with an increased risk of non-accidental mortality and mortality due to circulatory and respiratory disease. Risks for all causes of death examined were greatest after adjusting for socioeconomic and ecological covariates, though were reduced after adjusting for smoking, alcohol consumption, BMI, and fruit/vegetable consumption. The largest hazard ratios per 10 μg/m3 increase in PM2.5 were observed for respiratory mortality compared to the other cause-specific estimates. Elevated risk was observed for respiratory mortality associated with air pollution among obese respondents and never-smokers, though the differences between these and reference groups were not statistically significant. We also examined the shape of the exposure-response curve, and although the lowest measured concentration of PM2.5 was 1 μg/m3, we found no lower threshold for response. Although this finding is potentially informative for burden assessment, it is worth noting that we did not distinguish between anthropogenic and natural sources of PM2.5 in this study.
This study adds to previous work in Canada, which has a generally lower mean PM2.5 exposure than other countries, by providing direct adjustments for behavioural covariates (i.e., smoking and obesity) that are known contributors to mortality. This study used similar methodology to a previous study in Canada, the Canadian Census Health and Environment Cohort (CanCHEC) , but was unable to directly evaluate the role of behavioural covariates. In general, our HR estimates for non-accidental mortality (HR = 1.26; 95 % CI: 1.19–1.34) were greater than those in CanCHEC (HR = 1.15; 95 % CI: 1.13–1.16; Cochran’s Q = 9.3, p < 0.01), though our estimates for circulatory death were similar (CCHS HR = 1.19; 95 % CI: 1.07–1.31; CanCHEC HR = 1.16; 95 % C.I.: 1.13–1.18; Cochran’s Q = 0.1, p = 0.8) (all units per 10 μg/m3 increase in PM2.5) .
The fact that we found stronger associations between mortality and PM2.5 here than were observed in the CanCHEC study  might be due to improvements in estimates of PM2.5. A new PM2.5 model developed at a much finer scale (1 km2 grid rather than 10 km2 grid) allowed respondents to be assigned more accurate, finer-scale estimates exposure to of PM2.5. This improved exposure model may have a particularly strong effect on respondents who live in mid-sized cities (e.g., Calgary, Edmonton) that would otherwise have been assigned a lower, regional (i.e., rural) average (Fig. 1). However, this improvement is expected to be limited somewhat by the limitations of location error in geocoding residences based on postal code, as well as respondent mobility throughout the study area, resulting in differences in personal exposure. Another strength of this study was that it assigned exposures to respondents in the three years preceding death, thereby ensuring that exposure always preceded health effects rather than being assigned concurrently. This method also takes long-term variation of exposure into account.
In our study, HR estimates increased after the addition of ecological covariates, which differs from the earlier CanCHEC study, in which the addition of ecological covariates served to decrease the HR estimates . As described earlier, the ecological covariates used here were derived for smaller areas than in the CanCHEC study due to the absence of Census Tracts in rural areas. The methodological differences in deriving ecological covariates, particularly at a finer scale (i.e., assigning DA-level covariates rather than CT-level covariates), may also be one of the primary reasons why differences in HR estimates were observed between this study and CanCHEC, since fine scale covariates would be more spatially variable and covariates would more accurately reflect local socioeconomic conditions. Indeed, when the ecological covariates were removed from the Cox models of non-accidental mortality, the otherwise fully adjusted model provided an HR = 1.085 (Table 3), which is more consistent with the fully adjusted models in CanCHEC . Ecological covariates included in this study were all positively correlated with PM2.5 (Table 2). Given the much greater PM2.5 exposure in urban environments, this association for recent immigrants and persons of high educational attainment is possibly due to a higher population of both in cities. The correlation with PM2.5 was weaker for the proportion of low income families, which was consistent with the similar proportions of regional-adjusted low-income families in rural and urban environments .
Our HRs for non-accidental mortality were greater than those reported for all-cause mortality in other, international studies that had considered the same behavioural covariates, though were generally similar when ecological covariates were excluded from our estimates [6, 23–25]. For example, the American Cancer Society study, which included 1.2 million adults in the United States, estimated an HR for all-cause mortality of 1.06 per 10 μg/m3 increase in PM2.5 (95 % C.I.: 1.02–1.11) after controlling for behavioural covariates, though that study did not include ecological covariates . Similarly, a global, pooled meta-analysis estimated an all-cause mortality HR of 1.06 (95 % C.I.: 1.04–1.08) per 10 μg/m3 increase in PM2.5 . However, our results with ecological covariates were not significantly different from those of a large meta-analysis of European studies, where the pooled HR estimate for natural-cause mortality adjusted for socioeconomic and behavioural covariates (though not large-scale socioeconomic covariates) was 1.09 per 5 μg/m3 increase in PM2.5 (95 % C.I.: 1.03–1.14) . This estimate was not significantly lower than in our study (Cochran’s Q = 0.8, p = 0.4), where HR = 1.12 (95 % C.I.: 1.09–1.16) when scaled to a 5 μg/m3 increase in PM2.5.
Hazard ratio estimates for mortality due to circulatory disease (i.e., HR = 1.19) were generally consistent with those reported in the international literature, including the Harvard Six Cities study extended follow-up, which reported a HR of 1.28 per 10 μg/m3 increase in PM2.5 [6, 27], and a study in the U.K., which reported an HR of 1.05 per 1.9 μg/m3 increase in PM2.5 after adjustment for sex, age, BMI, and smoking (our study: HR = 1.03, 95 % CI: 1.01–1.05 when scaled to a 1.9 μg/m3 increase) . However, our estimate was much greater than that reported from a study in Rome (HR = 1.06, 95 % CI: 1.04–1.08), which adjusted for some individual and area-based socioeconomic covariates , and the Dutch Environmental Longitudinal Study (DUELS), which reported an HR of 1.09 (95 % CI: 1.06–1.12) per 10 μg/m3 increase in PM2.5 .
Our hazard ratio estimates for respiratory disease (HR = 1.52) were generally greater than those in the literature, though literature estimates for HRs vary among studies. For example, one study in Rome that used area-based socioeconomic covariates identified a non-significant HR of 1.03 for respiratory disease . On the other hand, the California Teachers Study identified an HR for respiratory mortality of 1.21 , and the Dutch cohort (DUELS) estimated an HR of 1.18 , which were similar to our HR estimate for respiratory mortality prior to adjustment for ecological covariates (HR = 1.21). Another study in the UK reported an HR of 1.17 (95 % CI: 1.12–1.22) per 1.9 μg/m3 increase in PM2.5 . Our HR estimate after adjustment for ecological covariates was lower than this study (HR = 1.08, 95 % CI: 1.04–1.12) when scaled to a 1.9 μg/m3 increase in PM2.5.
Our study also evaluated the role of effect modification by sex, age and behavioural covariates, and found a significantly greater HR estimate for non-accidental mortality among men than women. In a pooled European analysis of multiple cohorts, HRs were elevated among men but not women . Our results are overall similar, although our generally greater HR estimates for non-accidental mortality might explain why HR was significant for both men and women. Men also had a greater HR than did women for circulatory disease mortality (though the differences were not significant), similar to the AHS cohort . This finding was inconsistent with the results of a small (n = 3,239) cohort of white, non-smoking adults, where the relative risk of coronary heart disease mortality was elevated among women but not among men in a fully adjusted model . Observed differences might be, at least in part, explained by relatively small cohort sizes.
Our HR estimates for non-accidental and circulatory mortality among obese and normal weight groups were not significantly different. Effect modification of cardiovascular mortality by obesity had previously been evaluated elsewhere in two all-female cohorts. One study identified a significantly greater HR with increasing BMI, with an HR for obese women of 1.35 (95 % C.I.: 1.12–1.64 per 10 μg/m3 increase in PM2.5) . The other study did not test differences statistically among groups but did report an HR of 1.99 (95 % C.I.: 1.23–3.22 per 10 μg/m3 increase in PM2.5) for obese women . The ACS also reported a greater HR among obese men . In our study, obese respondents also had a high risk of respiratory mortality (HR = 1.76; 95 % C.I.: 1.15–2.69), though possibly due to a small number of deaths (n = 500), the HR estimate was not significantly different from the normal weight population.
In our study, persons who had never smoked had a qualitatively greater risk of non-accidental and circulatory mortality from fine particulate exposure than those who had smoked, though the difference between groups was non-significant (Table 5). This finding was consistent with the literature, where a marginally greater risk of cardiovascular mortality was observed among never smokers than among current or former smokers [6, 23, 33, 34]. In a Dutch cohort, respiratory mortality was qualitatively greater among current smokers than never smokers , a finding that was not consistent with our study.
There were several limitations with our study that may contribute to uncertainty in our estimates. The cohort was chosen because of the inclusion of various behavioural covariates, but it is generally much smaller than that of CanCHEC, which used the Census of population (i.e., 20 % of the population of Canada) . Mean estimates of PM2.5 in Canada are generally lower than in other study countries , and the effect size is relatively small, requiring a large sample size to have adequate power for HR estimation. As a result, in our study the 95 % CIs were very wide in comparison to other studies , and we were also unable to adequately assess the shape of the concentration-response curves for other causes of death. It is also worth mentioning that our study relied on self-reported estimates for BMI and smoking. Although we were able to mathematically adjust BMI for self-reporting error based on measured BMI from another survey, it is possible that estimates of smoking may underrepresent actual smoking rates. Additionally, the follow-up period in our study was relatively short, particularly for respondents who entered the cohort in the final survey year (i.e., 2008, with a maximum of 4 years of follow-up). However, respondents entering the cohort in the first year of survey and who had remained in the cohort for the entire period were followed for a maximum of 12 years, which is comparable to the mean follow-up period (i.e., 12.6 years) in a review of other cohorts examining the same relationship . The limitation of having a short follow-up period was mitigated somewhat by considering exposures that preceded the event.
In creating the cohort, 69,300 CCHS respondents were excluded since they were not linked to the HTSF (tax) file and were therefore not candidates for the probabilistic linkage. The excluded population were those who did not file a tax return, and the characteristics of this population differed somewhat from the cohort. In general, the excluded population was younger and had a lower educational attainment than the final cohort. Therefore, the cohort might be slightly biased towards higher educational attainment and those active in the labour market, though these same characteristics were used for adjustment in survival models.
Estimates of PM2.5 exposure were assigned at baseline at the person’s place of residence. Accuracy in geocoding residences was limited by the program PCCF+, which assigns residences to postal code representative points. The size of postal codes is relatively small (i.e. typically a few city blocks) in urban centres; therefore the PCCF+ program is highly accurate within these areas. However, estimates of PM2.5 exposure in rural areas are less likely to have been assigned accurately since postal code areas can be quite large. We performed a sensitivity analysis that considered only cohort members that lived within urban areas (i.e., Census Metropolitan Areas), and despite exposures being much greater in urban areas, results were not significantly different than those reported above (HR = 1.19, 95 % CI: 1.11–1.27, Cochran’s Q = 1.71, p = 0.19). Given the short follow-up period, we also did not assess mobility in this study, making the assumption that respondents did not move. By not assigning air pollution exposures based on changes to residential history, it is expected that there would be some degree of exposure misclassification associated with this limitation. A previous study using CanCHEC considered the assignment of exposures at baseline vs. considering mobility during the follow-up period on mortality risk attributed to PM2.5. In general, there was very little difference in HR estimates (i.e., HR = 1.03, 95 % CI: 1.02–1.03 from baseline exposure, vs. HR = 1.04, 95 % CI: 1.03–1.04 for exposure considering mobility) . Although about 41 % of Canadians moved within the five-year period of 2001 to 2006 , the majority of moves were within cities or regions of similar PM2.5 exposures (not published). To assess this limitation, we ran a sensitivity analysis where we included only persons who had at least 3 years of residence in the same postal code. HRs for non-accidental mortality were similar to those for the entire cohort (HR = 1.28, 95 % CI: 1.19–1.37).
Finally, the cohort was developed based on a probabilistic linkage methodology to assign deaths to CCHS members. We attempted to reduce the potential for linkage error by limiting our cohort to persons linked to a tax file, since mortality rates among cohort members not linked to a tax file were substantially lower due to fewer elements of respondent data that could be used for linkage.
In general, this study documented an association between non-accidental, circulatory, and respiratory mortality and fine particular matter in a cohort adjusted for socioeconomic, ecological, and behavioural covariates and exposed to a relatively low exposure distribution (mean = 6.3 μg/m3). Although our CI were wide in the concentration-response curve, an increased risk of mortality was observed even at very low concentrations of PM2.5 (Fig. 2), at values lower than the WHO guideline of 10 μg/m3 . Further studies on a larger cohort are needed to evaluate the shape of the concentration-response curve at these lower concentrations of PM2.5. We also updated the results of previous Canadian studies by using an improved, finer-scale exposure model to assign PM2.5 estimates to cohort members, which may have, in part, caused observed increases in HR estimates relative to CanCHEC . Finally, this study indicates that the addition of fine-scale behavioural covariates serves to reduce the HR estimates compared to the otherwise fully adjusted survival models.
United States Agricultural Health Study
Body Mass Index
Canadian Census Health and Environment Cohort
Canadian Community Health Survey
Canadian Mortality Database
Historical Tax Summary File
(−2) Log-likelihood ratio
moderate resolution imaging spectroradiometer
Postal Code Conversion File Plus
- PM2.5 :
Fine particulate matter
World Health Organization
Lin SS, Vos T, Flaxman AD, Danaei G, Shibuya K, Adair-Rohani H, et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: A systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2224–60.
GBD 2013 Risk Factors Collaborators. Global, regional and national comparative risk assessment of 79 behavioural, environmental/occupational and metabolic risks or clusters of risks in 188 countries 1990–2013: a systematic analysis for the GBD 2013. Lancet. 2015; in press.
Crouse DL, Peters PA, van Donkelaar A, Goldberg MS, Villeneuve PJ, Brion O, et al. Risk of nonaccidental and cardiovascular mortality in relation to long-term exposure to low concentrations of fine particulate matter: A Canadian national-level cohort study. Environ Health Persp. 2012;120:708–14.
Shin HH, Cakmak S, Brian O, Villeneuve P, Turner MC, Goldberg MS, et al. Indirect adjustment for multiple missing variables applicable to environmental epidemiology. Environ Res. 2014;134:482–7.
Weichenthal S, Villeneuve PJ, Burnett RT, van Donkelaar A, Martin RV, Jones RR, et al. Long-term exposure to fine particulate matter: association with nonaccidental and cardiovascular mortality in the Agricultural Health Study Cohort. Environ Health Persp. 2014;112:609–15.
Hoek G, Krishnan RM, Beelen R, Peters A, Ostro B, Brunekreef B, et al. Long-term air pollution exposure and cardio-respiratory mortality: a review. Environ Health. 2013;12:43.
Canada S. Canadian Community Health Survey (CCHS) annual component: User guide, 2007–2008 microdata files. Ottawa, ON, Canada: Statistics Canada; 2009.
Fellegi IP, Sunter AB. A theory for record linkage. J Am Stat Assoc. 1969;64:1183–210.
Wilkins R, Tjepkema M, Mustard C, Choinière R. The Canadian census mortality follow-up study, 1991 through 2001. Health Rep. 2008;19:25–43.
Wilkins R, Peters PA. PCCF + Version 5K* User’s Guide. Automated geographic coding based on the Statistics Canada Postal Code Conversion Files including Postal Codes through May 2011. Statistics Canada: Ottawa, ON, Canada, 2012. Catalogue no. 82F0086-XDB.
van Donkelaar A, Martin RV, Spurr RJD, Burnett RT. High-resolution satellite-derived PM2.5 from optimal estimation and geographically weighted regression over North America. Environ Sci Technol. 2015; in press (doi:10.1021/acs.est.5b02076).
Boys BL, Martin RV, van Donkelaar A, MacDonell RJ, Hsu NC, Cooper MJ, et al. Fifteen-year global time series of satellite-derived fine particulate matter. Environ Sci Technol. 2014;48:11109–18.
Cox DR. Regression models and life tables. J Royal Stat Soc B. 1972;20:187–220.
Statistics Canada. Census Dictionary 2006. Statistics Canada: Ottawa, ON, Canada, 2010. Catalogue no. 92-566-X.
Ross NA, Oliver LN, Villeneuve PJ. The contribution of neighbourhood material and social deprivation to survival: a 22-year follow-up of more than 500,000 Canadians. Int J Environ Res Public Health. 2013;10:1378–91.
Crouse DL, Peters PA, Villeneuve PJ, Proux M-O, Shin HH, Goldberg MS, et al. Within- and between-city contrasts in nitrogen dioxide and mortality in 10 Canadian cities; a subset of the Canadian Census Health and Environment Cohort (CanCHEC). J Expo Sci Env Epi. 2014. doi:10.1038/jes.2014.89.
Connor Gorber S, Shields M, Tremblay MS, McDowell I. The feasibility of establishing correction factors to adjust self-reported estimates of obesity in the Canadian Community Health Survey. Health Rep. 2008;19:71–82.
World Health Organization. BMI Classification: The international classification of adult underweight, overweight, and obesity according to BMI. http://apps.who.int/bmi/index.jsp?introPage=intro_3.html (2006). Accessed 28 Jul 2015.
Chen H, Burnett RT, Kwong JC, Villeneuve PJ, Goldberg MS, Brook RD, et al. Spatial associations between ambient fine particulate matter and incident hypertension. Circulation. 2013;129:562–9.
Conover W. Practical Nonparametric Statistics. 3rd ed. New York: Wiley; 1999.
Meira-Machado L, Cadarso-Suárez C, Gude F, Araújo A. smoothHR: An R package for pointwise nonparametric estimation of hazard ratio curves of continuous predictors. Comput Math Methods Med. 2013; doi:10.1155/2013/745742.
Fortin, M. A comparison of rural and urban workers living in low-income. Rural and Small Town Canada Analysis Bulletin (Statistics Canada). 2008;7:1–18. Cat. no. 21-006-XIE.
Pope CA, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, et al. Lung cancer, cardiopulmonary mortality and long-term exposure to fine particulate air pollution. JAMA. 2002;287:1132–41.
Cesaroni G, Badaloni C, Gariazzo C, Staffogia M, Sozzi R, Davoli M, et al. Long-term exposure to urban air pollution and mortality in a cohort of more than a million adults in Rome. Environ Health Perspect. 2013;121:324–31.
Fischer PH, Marra M, Ameling CB, Hoek G, Beelen R, de Hoogh K, et al. Air pollution and mortality in seven million adults: The Dutch Environmental Longitudinal Study (DUELS). Environ Health Perspect. 2015; doi: 10.1289/ehp.1408254.
Beelen R, Raaschou-Nielsen O, Stafoggia M, Jovanovic Andersen Z, Weinmayr G, Hoffmann B, et al. Effects of long-term exposure to air pollution on natural-cause mortality: an analysis of 22 European cohorts within the multicentre ESCAPE project. Lancet. 2014;383:785–95.
Laden F, Schwartz J, Speizer FE, Dockery DW. Reduction in fine particulate air pollution and mortality: extended follow-up of the Harvard Six Cities Study. Am J Resp Crit Care. 2006;173:667–72.
Carey IM, Atkinson RW, Kent AJ, van Staa T, Cook DG, Anderson HR. Mortality associations with long-term exposure to outdoor air pollution in a national English cohort. Am J Respir Crit Care Med. 2013;187:1226–33.
Ostro B, Lipsett M, Reynolds P, Goldberg D, Hertz A, Garcia C, et al. Long term exposure to constituents of fine particulate air pollution and mortality: results from the California teachers study. Environ Health Persp. 2010;118:363–9.
Chen LH, Knutsen SF, Shavlik D, Beeson WL, Petersen F, Ghamsary M, et al. The association between fatal coronary heart disease and ambient particulate air pollution: are females at greater risk? Environ Health Persp. 2005;113:1723–9.
Miller KA, Siscovick DS, Sheppard L, Shepherd K, Sullivan JH, Anderson GL, et al. Long-term exposure to air pollution and incidence of cardiovascular events in women. New Engl J Med. 2007;356:447–58.
Puett RC, Schwartz J, Hart JE, Yanosky JD, Speizer FE, Suh H, et al. Chronic particulate exposure, mortality, and coronary heart disease in the nurses’ health study. Am J Epidemiol. 2008;168:1161–8.
Krewski D, Jerrett M, Burnett RT, Ma R, Hughes E, Shi Y, et al. 2009. Extended follow-up and spatial analysis of the American Cancer Society study linking particulate air pollution and mortality. Res Rep Health Eff Inst. 2009;140:5–136.
Beelen R, Hoek G, van Den Brandt PA, Goldbohm RA, Fisher P, Shouten LJ, et al. Long-term effects of traffic-related air pollution on mortality in a Dutch cohort (NLCS-AIR study). Environ Health Persp. 2008;116:196–202.
Crouse DL, Peters PA, Hystad P, Brook JR, van Donkelaar A, Martin RV, et al. Ambient PM2.5, O3, and NO2 exposures and associations with mortality over 16 years of follow-up in the Canadian Census Health and Environment Cohort (CanCHEC). Environ Health Persp. 2015;123:1180–6.
Statistics Canada. 2006 Census data products: 2006 census trends. http://www12.statcan.ca/census-recensement/2006/dp-pd/92-596/P2-2.cfm?Lang=eng&T=CSD&LINE_ID=701&TOPIC_ID=700 (2010). Accessed 12 Dec 2015.
The authors wish to acknowledge the contributions of Masoud Nasari (Health Canada) for his assistance with R software and Saeeda Khan (Statistics Canada) for her assistance with the PCCF+ program.
The authors do not declare any competing interests.
LP designed the study, linked the cohort to the air pollution models, conducted the statistical and GIS analyses, and drafted the manuscript. MT participated in the study design, the cohort linkage, and provided feedback on the cohort. DLC participated in the study design, and provided assistance with statistical techniques and covariate preparation. SW conceived of the study and participated in the study design and in drafting the manuscript. AVD and RVM developed the air pollution models and provided feedback on their use in an epidemiologic context. MB and HC provided early feedback to improve the analyses and the presentation of results. RTB coordinated the study and participated in the study design and the statistical analysis. All authors actively edited the manuscript and approved the final manuscript.