Low concentrations of fine particle air pollution and mortality in the Canadian Community Health Survey cohort
Environmental Health volume 18, Article number: 84 (2019)
Approximately 2.9 million deaths are attributed to ambient fine particle air pollution around the world each year (PM2.5). In general, cohort studies of mortality and outdoor PM2.5 concentrations have limited information on individuals exposed to low levels of PM2.5 as well as covariates such as smoking behaviours, alcohol consumption, and diet which may confound relationships with mortality. This study provides an updated and extended analysis of the Canadian Community Health Survey-Mortality cohort: a population-based cohort with detailed PM2.5 exposure data and information on a number of important individual-level behavioural risk factors. We also used this rich dataset to provide insight into the shape of the concentration-response curve for mortality at low levels of PM2.5.
Respondents to the Canadian Community Health Survey from 2000 to 2012 were linked by postal code history from 1981 to 2016 to high resolution PM2.5 exposure estimates, and mortality incidence to 2016. Cox proportional hazard models were used to estimate the relationship between non-accidental mortality and ambient PM2.5 concentrations (measured as a three-year average with a one-year lag) adjusted for socio-economic, behavioural, and time-varying contextual covariates.
In total, 50,700 deaths from non-accidental causes occurred in the cohort over the follow-up period. Annual average ambient PM2.5 concentrations were low (i.e. 5.9 μg/m3, s.d. 2.0) and each 10 μg/m3 increase in exposure was associated with an increase in non-accidental mortality (HR = 1.11; 95% CI 1.04–1.18). Adjustment for behavioural covariates did not materially change this relationship. We estimated a supra-linear concentration-response curve extending to concentrations below 2 μg/m3 using a shape constrained health impact function. Mortality risks associated with exposure to PM2.5 were increased for males, those under age 65, and non-immigrants. Hazard ratios for PM2.5 and mortality were attenuated when gaseous pollutants were included in models.
Outdoor PM2.5 concentrations were associated with non-accidental mortality and adjusting for individual-level behavioural covariates did not materially change this relationship. The concentration-response curve was supra-linear with increased mortality risks extending to low outdoor PM2.5 concentrations.
Exposure to ambient fine particle air pollution (PM2.5) is responsible for an estimated 2.9 million deaths annually and 83 million disability-adjusted life years lost , with several large epidemiological cohort studies linking long-term exposure to PM2.5 to all-cause and cause-specific mortality [2,3,4]. Even in settings with relatively low concentrations of air pollution, such as Canada, the relationships persist [5, 6]. Despite these findings, there remain two key areas of potential bias and uncertainty that past work has been unable to address simultaneously. The first is the inability to directly adjust for individual-level behavioural risk factors associated with chronic disease mortality, such as smoking, diet, and exercise, or health measures such as body mass index; various indirect methods for adjustment have been applied elsewhere [7, 8]. The second regards the shape of the concentration-response curve for PM2.5 and mortality. This issue has become increasingly pertinent as clean air regulations have succeeded in reducing PM2.5 concentrations across North America and elsewhere, and thus a better understanding of the shape of the PM2.5-mortality associations at low concentrations are required for cost-benefit assessments of future reduction efforts.
The purpose of this study was to provide an updated and extended analysis of the Canadian Community Health Survey-Mortality cohort  including : additional years of follow-up to 2016 ; improvements in the resolution of PM2.5 exposure (approximately 1km2 grid) ; annual residential history from 1981 to 2016 for all cohort members from a linkage to postal code records ; time-varying contextual covariates ; inclusion of immigrants to Canada, and  an improved linkage between survey respondents and death records. We examine the shape of the concentration-response curves using a Shape Constrained Health Impact Function (SCHIF)  and perform sensitivity analyses.
The Canadian Community Health Survey (CCHS) is a national cross-sectional survey of the Canadian population that collects information related to health status, health care utilization, and health determinants. From 2000 to 2007 the survey was administered every 2 years to approximately 130,000 respondents; from 2007 onwards, data has been collected on an ongoing basis from 65,000 respondents per year and released annually with response rates declining over time (Fig. 1) [11,12,13,14,15,16]. The CCHS data are sampled from approximately 98% of the Canadian population aged 12 and older living in private dwellings within the 115 Heath Regions covering all provinces and territories. Individuals living on Indian Reserves and on Crown Land, institutional residents, full-time members of the Canadian Forces, and residents of certain remote regions are excluded.
Consent to record linkage and data sharing was obtained at time of survey (Fig. 1) and only those CCHS respondents who agreed were linked to death records and residential history through Statistics Canada’s Social Data Linkage Environment (SDLE) . The linkage was approved by Statistics Canada  and is governed by the Directive on Microdata Linkage . The linkage occurred within the Derived Record Depository, a highly secure linkage environment comprised of a national dynamic relational database of basic personal identifiers. Survey and administrative data are linked to the Derived Record Depository using a SAS-based generalized record linkage software (G-link) that supports deterministic and probabilistic linkage based on the mathematical theory of record linkage developed at Statistics Canada . Mortality linkage to the Derived Record Depository between 2000 and 2016 was 99.8% . A list of linked unique individuals is created through linkages that are deterministic (matching records based on unique identifiers) and probabilistic (matching records based on non-unique identifiers such as names, sex, date of birth, and postal code and estimating the likelihood that records are referring to the same entity). For the CCHS cycles considered, there was a linkage rate to the Derived Record Depository of 95.2% and a false error rate for the CCHS to SDLE linkage of 0.4% .
There were 666,000 CCHS respondents who agreed to data linkage (Fig. 1), reduced to 540,900 after excluding subjects with death dates prior to survey response dates (i.e. either the death record or linkage must be incorrect) or who were below the age of 25 or above 89 at time of survey as they are less likely to reside at the same postal code as their income tax mailing address . The CCHS to SDLE linkage rates across key indicators were consistently high, ranging from 94.4% for the 20–29 age group to 96.2% for the 80–89 age group, 95.5 and 95.3% for males and females respectively, and by province/territory from 91.8% for the Yukon to 96.7% for Newfoundland and Labrador .
After linkage, we stacked the CCHS cycles into one data file. We standardized the variable categorizations when discrepancies between cycles existed. The covariates (listed, along with categorizations in Table 1) included socio-economic, behavioural, and contextual measures. More information about the definitions and classifications of these variables can be found elsewhere . Provincially standardized deciles were calculated according to the distribution of residents in each province and the ratio of their total household income to the low-income cut-off for their corresponding household and community size. As this measure excluded subjects living in territories, we took the mean income within each decile and used these as cut-offs to categorize subjects living in territories by income into deciles. Once all subjects were placed in deciles, we merged groups to create quintiles.
Postal code history was complete from 1981 to 2016 for 35.0% of respondents and 12.6% of respondents had no postal code history. There were gaps in postal code histories for 52.4% of respondents, which is to be expected, as taxes may not have been filed for various reasons (e.g immigration, death, or age). We imputed complete or partial postal codes only when bookended by postal codes with sufficient similarity before and after the gap . For example, if a postal code in 2008 was K1A 0T6 and then K1A 0K9 in 2012, a partial postal code of K1A 0** would be imputed for the four missing years from 2009 to 2011. We did not impute postal codes if a gap existed at the beginning or end of the follow-up period or after a person’s death; full or partial postal codes (two to five digits) were imputed for 1.5% of person-years.
We organised the cohort into a person-year file with each year of exposure (1981–2016) per person representing a row of data. Subsequently we excluded specific person-years  once they turned age 90 during follow-up , if the person had immigrated to Canada less than 10 years prior to survey interview , if there was no postal code , if the postal code could not be linked to air pollution or contextual covariates , if the PM2.5 three-year moving average with a one-year lag was calculated by fewer than 2 years of exposure data, or  if the person-year was before survey interview date or after a person’s death (Fig. 1). We excluded recent immigrants to Canada (10 years or less) since they have spent the majority of their lives outside of Canada with unknown exposure, and this time exceeds the number of years in Canada where exposure can be estimated.
Exposure file and analytical file
The task of linking contextual covariates and air pollution values to the cohort required the creation of a master list of postal codes with their respective points of latitude and longitude and census geography. We produced this list from Statistics Canada’s June 2017 Postal Code Conversion File and the two previous versions (August 2015 and May 2011) to ensure coverage of retired postal codes [25,26,27]. Since census geography does not align with postal code locations, a single postal code can have multiple points of latitude and longitude. Each can represent the centroid of a blockface (i.e. a street block), dissemination block (i.e. an area bounded on all sides by roads), or dissemination area (i.e. adjacent dissemination blocks that collectively contain 400 to 700 persons) within a postal code.
We developed and used annual exposure estimates of PM2.5 from 1998 to 2012 by relating satellite retrievals of aerosol optical depth (AOD) to near-surface PM2.5 concentrations using the geophysically-based relationship simulated by a chemical transport model . Ground monitoring data from the National Air Pollution Surveillance (NAPS) network were then incorporated, along with other North American-based measurements, to constrain these PM2.5 estimates with geographically weighted regression. The resulting ambient PM2.5 surface provided estimates for North America at about a 1km2 resolution . Spatial variation from this surface was used with simulated PM2.5 and consistently constrained with local ground-based monitors to extend our PM2.5 coverage to 2015 .
The ambient warm season daily-maximum eight-hour average O3 surfaces were developed by Environment and Climate Change Canada for 2002–2015 using chemical transport modelling informed by surface observations as hourly estimates from 2002 to 2015 [30,31,32]. Estimates of NO2 were created using a national land use regression model (LUR) informed by on satellite-derived NO2 (10 km resolution), distances to highways and major roads, and roadway kernel density gradients . Ozone and NO2 values were back-casted to obtain exposures for 1981–2015 using ground monitoring data from the Canadian National Air Pollution Surveillance program. Annual adjustment factors were calculated at a census division level from the ratio of observed concentration to the values in the surface for the reference year (see Pinault et al. for more detail ).
We linked postal codes to PM2.5 in ArcGIS Desktop 10.5.1 using the points of latitude and longitude from the master postal code list and the air pollutant surfaces. In cases where there were multiple points of latitude and longitude for a single urban postal code, we used equal weighting of the multiple air pollutant values to provide a singular value. In rural communities, we took the population-weighted average of the values associated with duplicate postal codes. We used population-weighing to average multiple values to create inputs for partial postal codes (2 to 5 digits).
Contextual covariates were available at various census geographies and we merged these to individual person-years via postal codes (as described below). We created historic measures when possible to reflect neighbourhood-level changes over the time.
Regions of Canada that share air quality characteristics and movement patterns have been defined by the Air Quality Management System (AQMS) as six distinct airsheds . By subdividing the country into large geographic areas, adjustment for the broad spatial variation in mortality rates can be performed [9, 34]. We assigned airshed to the cohort by postal code. We used a population weighted mode in cases where there were multiple points of latitude and longitude for a single postal code.
We developed a historic community size variable to account for different sizes of metropolitan regions and changes in population over time, classifying Census Metropolitan Areas (CMAs: major urban core, 100,000+ residents) and Census Agglomerations (CAs: smaller urban cores, 10,000+ residents) by population counts . Since CMA/CAs cover large areas that can include farmland near the urban-rural fringe and residential enclaves of commuters to the city, we created a measure that accounts for differences in urban form within these CMA/CAs. We used population density measures (1991–2016) and frequencies for different modes of transportation at the neighbourhood level (1996–2016) to categorize census tracts as active urban core, transit-reliant suburb, car-reliant suburb, exurban, and non-CMA/CA . Both CMA/CA size and urban form were attached to the postal code list via census geography before merging with the cohort. In cases where there were multiple points of latitude and longitude representing a postal code, we used a population-weighted mode to assign categories.
The Canadian Marginalization Index (Can-MARG) is a measure of community-level marginalization comprised of four factors: material deprivation (e.g. proportion of people living in dwellings in need of repair), residential instability (e.g. proportion of people who live in a dwelling that they do not own) dependency (e.g. proportion of seniors and youth compared to those who are not), and ethnic concentration (e.g. proportion of recent immigrants and self-reported visible minorities) . We used historic census tract-level Can-MARG values in CMA/CAs, and a population-weighted aggregation of the dissemination area-level Can-MARG values at the census subdivision level in rural areas outside of CMA/CAs that are not covered by census tracts. We assigned Can-MARG values to points of latitude and longitude before quintiles were assigned.
We calculated for each individual and year of follow-up a three-year moving average for PM2.5 with a one-year lag, (e.g. the exposure in 2002 is the average of exposures in 1999, 2000, and 2001).
We ran standard Cox proportional hazard models to assess the relationship between PM2.5 exposure and non-accidental death (ICD-10 codes A to R) from survey interview year to the end of follow-up period or year of death. We started model building with a baseline hazard function for PM2.5 stratified by five-year age groups, sex, and survey cycle to ensure that respondents within these strata would be broadly comparable. We calculated new hazard ratios for models that included each socio-economic and behavioural covariate individually. We included covariates in the partially-adjusted model if the log difference between the new hazard ratio and the baseline was more than 10%. Subsequently, we added contextual covariates individually to the partially-adjusted model and included them in the final model using the same criteria (comparing to the partially-adjusted model that included socio-economic and behavioural covariates). All covariates considered for inclusion in the final model and the associated hazard ratios are found in Table 1.
We examined the shape of the association between PM2.5 and mortality with a SCHIF . This method is based on a construction of several transformations of concentration and fitting the transformed variable in a Cox model, estimating the log-hazard ratio for a unit change in the transformed variable and its standard error. An ensemble of all models examined was then constructed using a weighted average of the predicted log-hazard ratio and any concentration, with weights defined by the AIC of each model. The transformations are variations on a sigmoidal function which yields supra-linear, near linear, and sub-linear shapes.
We examined effect modification by select socio-economic and behavioural covariates, and by high- and low-exposure groups to the combined oxidant capacity of NO2 and O3 (henceforth: OX) which is calculated as the redox-weighted oxidant capacity  i.e. a weighted average using redox potentials as the weights (Oxwt = [(1.07 V × NO2) + (2.075 V × O3)]/3.145 V) (Table 4). We examined multiple pollutant models to investigate whether the inclusion of other common pollutants (NO2, O3, and OX) in the model may modify the PM2.5-mortality relationship [5, 39].
There were 4,452,700 person-years in the cohort after exclusion criteria were applied (Fig. 1) from 452,700 unique individuals. Entry into the cohort and length of the follow-up period varied by survey cycle, with the first cohort having up to 15 years of follow-up. For those who died, the average follow-up period was 5.1 years (s.d. 3.4); it was 6.5 years (s.d. 4.1) for those who survived the follow-up period. There were 50,700 non-accidental deaths. Of these, there were 7900 deaths from ischemic heart disease, 2800 from cerebrovascular disease, and 4300 from other cardiovascular diseases; 900 from pneumonia, 2800 from COPD, and 1100 from other respiratory diseases; 5500 from lung cancer, 1300 from colon cancer, 1300 from breast cancer, 1100 from pancreatic cancer, and 9900 from all other cancers. Further, there were 1700 deaths from diabetes, 3900 deaths from neuropsychiatric conditions, 2200 from digestive diseases, 1100 from genitourinary diseases and 3000 from all other non-accidental causes.
Exposure to PM2.5 was higher in women, more recent immigrants, and non-Indigenous people. Being single, university educated, and in the poorest income quintile were also associated with higher exposures (Table 1). We observed higher exposure to PM2.5 in people living in the largest CMAs and in the East Central airshed (which includes Toronto and Montreal). The distribution of exposure estimates for PM2.5, NO2, O3, and OX is found in Table 2.
The cohort was generally representative of the Canadian population, as seen through their mortality rates by subgroup (Table 1). Immigrants and non-Indigenous people had lower mortality rates compared to their counterparts. Being married, holding a university degree, and being employed were associated with a lower risk of mortality. As expected, there were clear trends in mortality risk with income, education, and immigrant status.
The unadjusted model had a hazard ratio of 0.96 (95% CI 0.92–1.00) which increased to 1.11 (95% CI 1.04–1.18) when adjusted by the socio-economic, behavioural, and contextual covariates that met the inclusion threshold (Table 3). All covariates except for body mass index (BMI), employment status, and urban form met the criteria and were included in the final model. When we added the behavioural covariates to a model that included only socio-economic covariates the hazard ratio increased from 1.05 (95% CI 1.00–1.09) to 1.09 (95% CI 1.05–1.15). Conversely, when we added the behavioural covariates to a model that included both the socio-economic and contextual covariates, they lowered the PM2.5 hazard ratio from 1.13 (95% CI 1.06–1.21) to 1.11 (95% CI 1.04–1.18).
The SCHIF characterisation of the PM2.5-mortality association (for all cohort members) displayed a supra-linear shape that rises in a steeper fashion compared to the standard log-linear model prediction for lower concentrations and changes in a more moderate manner for higher levels (Fig. 2). Note that the SCHIF displays wider uncertainty intervals compared to the log-linear model at low concentrations, in part due to the additional variation associated with model shape, a feature captured by the SCHIF but not the log-linear model. We observed a positive and statistically significant (p < 0.05) association between PM2.5 and non-accidental mortality for all concentrations examined as indicated by the SCHIF hazard ratio predictions.
We assessed effect modification within the PM2.5-mortality relationship by separating the cohort by age, sex, immigrant status (i.e. immigrants who had been in Canada for 10 or more years vs. non-immigrants), and educational attainment, and comparing resulting hazard ratios with Cochrane’s Q (Table 4). The hazard ratio was 4% higher for males (1.13 95% CI 1.03–1.23) than females (1.09 95% CI 0.99–1.19). When contrasted by age, the hazard ratio was 9% lower for those aged 75 years or more (1.04 95% CI 0.94–1.16) compared to those aged 65–74 (1.13 95% CI 1.01–1.27) and 10% lower compared to those aged 65 or less (1.14 95% CI 1.01–1.29). The hazard ratio for non-immigrants was higher than that of the final model (1.14 95% CI 1.07–1.23) and the immigrant group had a null hazard ratio (0.98 95% CI 0.83–1.16). The hazard ratio for those without a high school diploma (1.08 95% CI 098–1.19) was lower than those who graduated from high school (1.14 95% CI 1.04–1.24). The Cochrane’s Q p-values did not indicate that the above hazard ratios were significantly different between subgroups. We repeated the effect modification analyses for behavioural covariates. There was no significant difference between those who consumed fewer than five servings of fruits and vegetables per day compared to those who consumed five or more (1.10 95% CI 1.01–1.20 vs. 1.16 95% CI 1.04–1.30) although the hazard ratio was higher for those who consumed more fruits and vegetables. We found that hazard ratios were higher for regular drinkers (1.18 95% CI 1.09–1.28) and daily or occasional smokers (1.13 95% CI 0.99–1.27) compared to never or former drinkers (1.01 95% CI 0.90–1.12) or never or former smokers (1.11 95% CI 1.03–1.20), with a significant difference found between those who do and do not consume alcohol (p < 0.05). The HRs produced for each subgroup were pooled (Table 4), resulting in HRs that were similar to the full cohort final model. The high- and low- Ox groups had significantly different PM2.5-mortality hazard ratios. The inclusion of other pollutants (O3, NO2, and OX) attenuated the PM2.5 hazard ratios and produced confidence intervals that include a null value, with the greatest reduction seen in the model that included PM2.5, NO2, and O3 (1.00 95% CI 0.98–1.02, 1.03 95% CI 1.01–1.05, 1.05 95% CI 1.03–1.07 respectively) (Table 5).
Using a cohort comprised of several cycles of a health survey with up to a 15-year follow-up period and high resolution exposure estimates, we found that exposure to PM2.5 was associated with an 11% increase in non-accidental mortality per 10 μg/m3 after extensive adjustment for socio-economic, behavioural, and contextual covariates.
The hazard ratio for the full cohort was similar to that of the Nurse’s Health Study in the United States (1.13 95% CI 1.05–1.22) that adjusted for individual-level socio-economic and behavioural covariates  and a cohort from England (1.13 95% CI 1.00–1.25) that controlled for smoking, BMI, income, age, and sex . Burnett and colleagues  report hazard ratio estimates for a 10 μg/m3 change in long-term exposure to PM2.5 and non-accidental mortality in 41 cohorts conducted globally, 36 of which included adjustment for behavioural risk factors. The pooled hazard ratio among these 36 cohorts was 1.09 (95% CI 1.05–1.12), a value similar to that observed in our current study (1.11 95% CI 1.04–1.18). A version of the 2001 CanCHEC census-based cohort produced a hazard ratio that is similar to this work (1.09 95% CI 1.07–1.11) .
The impact of individual-level behavioural risk factors on the PM2.5-mortality association was assessed to address a common critique of many large administrative cohort studies examining the air pollution-mortality relationship. The inclusion of behavioural covariates to a model including socio-economic and ecological covariates lowered the PM2.5 hazard ratio 2% (from 1.13 to 1.11). This modest change in the hazard ratio can be interpreted to indicate that the behavioural covariates were being adequately controlled for by the socio-economic and ecological covariates in the established relationship between PM2.5 exposure and non-accidental mortality. This finding is similar to the previous CCHS cohort analysis and analysis of a Medicare-based cohort; both reported that adjustment for behavioural covariates had a minimal effect on hazard ratios [3, 7]. There is evidence (Tables 3 and 4) for a small increase in risk of PM2.5-related mortality in occasional or regular drinkers but this may be masked by null effects from the inclusion of other behavioural covariates (fruit and vegetable consumption, smoking behaviours) and this confounding is likely the result of the spatial distribution of drinking behaviours, with binge drinkers having the largest mortality risk but lower PM2.5 exposures. This study, through its inclusion of multiple covariates and an explicit a priori analysis approach for model building therefore provides the most extensive evidence to date that, in the Canadian context, missing data on behavioural risk factors for mortality have a minimal confounding bias on the PM2.5-mortality association.
The increase in the PM2.5 hazard ratio with the addition of the ecological covariates was largely driven by the addition of airsheds. Not only do these airsheds characterize broad air movement patterns, they also capture areas with similar composition of PM2.5 (e.g., proportion of PM2.5 composed of nitrate is highest in the Prairie airshed, whereas the Southern Atlantic airshed is composed of a notably higher proportion of black carbon) . They also delineate general socio-cultural groups with distinct mortality risk factors beyond those captured by the typical socioeconomic census variables included in our survival models. The three airsheds with the largest hazard ratios, along with high material deprivation, all have the lowest levels of air pollutants which would account for the negative confounding effect observed in Table 3. Further, the largest airshed (East Central) contains both Toronto and Montreal, the two largest CMAs in Canada and significant population hubs. High PM2.5 exposure and related mortality are largely driven by the population of Toronto (21% of the national population in 2006) where the mean PM2.5 exposure is 9.33 μg/m3 whereas the mean in the rest of the country is 7.68 μg/m3 . These results are consistent with a descriptive analysis of PM2.5 exposure in 2006 long-form census respondents . Although urban areas are the most common residence for both high income and highly educated Canadians, rural residences are more common among the high income earners than university graduates (i.e within the highest income quintile, 73.7% urban vs 26.3% urban fringe or rural; among those who are university educated, 82.6% urban vs. 17.3% urban fringe or rural). The greater tendency for high-income Canadians to live in rural areas is consistent with the findings in this paper. As a result, PM2.5 exposure by income categories is a slightly more linear pattern than education in both of these studies.
We estimated the shape of the concentration-response (CR) function for the PM2.5-mortality association. A slight supra-linear association (Fig. 2) was found, with a steep CR function at the low to median PM2.5 range which levelled off slightly after approximately 10 μg/m3. The SCHIF hazard ratio predictions indicated a positive and significant association between PM2.5 and non-accidental mortality for all concentrations, suggesting risks to concentrations below 2 μg/m3. Previous work using a CCHS-based cohort used a spline-based procedure and found that the shape of the relationship between non-accidental mortality and PM2.5 was supra-linear in shape with a threshold of 4.5 μg/m3, but was limited due to wide confidence intervals . A study in China using a SCHIF function found non-linear relationships for multiple causes of death . Such a relationship, when applied in a health impact framework, as in the Global Burden of Disease [45, 46] and in the recent Global Exposure Mortality Model  suggest benefits both from reducing PM2.5 concentrations areas with the highest concentrations and from continuing to reduce them in relatively cleaner areas, including Canada, where it is estimated that the entire population now lives in areas with ambient PM2.5 concentrations below the current WHO Guideline . Worldwide it is estimated that small absolute reductions under 3 μg/m3 could prevent hundreds of thousands of deaths in areas that comparatively have low levels of PM2.5 .
The risk of non-accidental mortality from exposure to PM2.5 was 4% higher in males over females (males 1.13, females 1.09), a pattern that has emerged in similar work. The hazard ratios from the current study are more aligned with the ESCAPE European pooled cohort (males 1.14 95% CI 1.04–1.24; females 0.99 95% CI 0.92–1.07)  albeit with a higher hazard ratio for women when compared to the previous version of CCHS-based cohort (males 1.34 95% CI 1.24–1.46; females 1.18 95% CI 1.09–1.28) . Hazard ratios were lowest for members of the cohort aged 75 and older (1.04) and were similar for those aged 65 and under (1.14) and 65 to 75 (1.13); this is similar to the European study which found that risk decreases with age (< 60 years 1.16 95% CI 1.00–1.34; 60–75 years 1.10 95% CI 1.00–1.20; ≥75 years 1.03 95% CI 0.95–1.11). When we divided the cohort into immigrants (in Canada for 10 years or more) and non-immigrants, the PM2.5-mortality association increased for non-immigrants and was null among the immigrant population. This result is consistent with prior Canadian census-based cohort studies  and is possibly the result of what is termed the “healthy immigrant effect” [49,50,51,52,53], likely intensified by the preferential settlement of immigrants into the largest cities which have higher PM2.5 exposure. The hazard ratio for high school graduates (1.14) was higher than for those without a diploma (1.08) which is to be expected given that the latter is more likely to live in rural areas , and have a mean PM2.5 exposure that is lower than other educational groups (Table 1) . We examined effect modification by behavioural covariates (i.e., fruit and vegetable consumption, smoking behaviour, and alcohol consumption) and found significant difference in the resulting hazard ratios only in the case of alcohol consumption. Effect modification analyses on the ESCAPE cohort also found no effect modification by fruit and vegetable consumption or smoking behaviour, but did not consider alcohol consumption .
The multiple pollutant models indicated that the relationship between non-accidental mortality and PM2.5 exposure are attenuated when we included other pollutants (NO2, O3, and OX) in the models. These findings indicate both that PM2.5 is associated with mortality and that the inclusion of gaseous co-pollutants, Ox in particular, may better characterize the biologically active aspects of PM2.5 and the overall air pollution mixture compared to the PM2.5 mass concentration . Weichenthal et al. looked at the effect modification of oxidant gases on PM2.5 more specifically and found that spatial variations in Ox concentrations may act as surrogates for the presence or absence of harmful air pollutant mixtures that enhance PM2.5 toxicity . We examined the PM2.5-mortality association in both low- and high- Ox person-years and found a 24% difference in risk. Our findings support these previous studies using different longitudinal Canadian cohorts and that knowledge of interactions between PM2.5 and oxidant gases leading to adverse health will improve risk management activities and public health.
We performed this analysis on an extended and updated version of a cohort described in a previous study by Pinault et al.  with improvements to the exposure assessment and linkage to death, postal code history, environmental exposures, and contextual covariates. While some of the results are comparable to the previous cohort (e.g. socio-economic + behavioural covariate models are within a 1% margin), there are differences in the covariates included in the final models and the resulting hazard ratios. This is not unexpected since the contextual covariates addressing area-level marginalization in the two studies were created differently (area-level proportions of specific variables vs. a principle component analysis which resulted in four factors), and measured at different geographical units (census divisions vs. census tracts and census subdivisions). Another difference is that the updated cohort and current work includes immigrants who have lived in Canada for ten or more years whereas the previous work only included those who had been in Canada for 20 or more years. This newly included group of semi-recent immigrants (10–19 years in Canada) have substantially lower hazard ratios of mortality compared to the non-immigrant population (Table 1). Their inclusion in the current study acts to reduce the overall PM2.5 hazard ratio (Table 4).
This large, national cohort is an extension and improvement to the previous CCHS-Mortality cohort, with an updated linkage and extended follow-up period for mortality and postal code history which now spans 36 years (1981 to 2016). More broadly the cohort has many strengths, including the fine resolution of the PM2.5 estimates (1km2), the ability to incorporate mobility across the follow-up years, an explicit a priori model building strategy, the inclusion of multiple time-varying contextual covariates to address spatial, neighbourhood- and city-level characteristics, and most uniquely the behavioural covariates such as smoking behaviours, alcohol consumption, diet, and exercise to control for health behaviours related to mortality that are not typically found on cohorts of this size.
This cohort and the analysis are limited by the data available. First, postal code history was derived from tax and administrative data. Historical postal codes reflect the mailing address as reported on a tax return and not necessarily a person’s residence; in 92.9% of cases the postal code reflects the person’s residence at time of survey . Similarly, outdoor ambient levels of PM2.5 at a person’s residence may not reflect their actual exposure. Sensitivity analysis performed with the 2001 CanCHEC found that finer scale resolution (1km2) estimates of PM2.5 resulted in lower AIC values and higher hazard ratios in the PM2.5-mortality model for non-accidental death compared to a 10km2 or 5km2 grid indicating that exposure estimates that are more specific to a person’s residence are appropriate . Gaps in postal code history are imputed under the assumption that the person did not leave the country or community during that time. In assigning contextual covariates by postal code, misclassification may occur from taking the mode or mean when estimating a single value to represent multiple points of latitude and longitude for a single postal code. Second, in contrast to the CanCHEC cohorts (Pappin AJ, Crouse DL, Christidis T, Pinault LL, Tjepkema M, Erickson A, Brauer M, Weichenthal S, van Donkelaar A, Martin RV, Brook J, Hystad P. Burnett RT. Associations between low levels of fine particulate matter andmortality within Canadian cohorts. Environ Health Persp., under review), this cohort does not completely represent the full Canadian population; the Canadian Community Health Survey is not a census of the population and survey weights were not used in this analysis. Further, in creating this cohort persons were removed if they did not consent to data linkage or if they could not be linked to the SDLE. The CCHS over-samples rural communities  which results in a disproportionate sample in areas with low levels of PM2.5 and higher rates of mortality. The sampling framework and un-weighted analysis likely caused the null unadjusted hazard ratio which became positive as covariates were added to the model to address confounding. These results are consistent with the Agricultural Health study which examined non-accidental death related to PM2.5 in rural communities in two American states (Iowa and North Carolina) and found a protective hazard ratio in minimally and fully adjusted models . Regardless, the protective unadjusted hazard ratio should not come as a surprise as contextual and socio-economic covariates are included in models because we know that they are related to both PM2.5 and mortality and can act as confounders (see Table 1 for the mortality Hazard Ratios by individual covariates). Given that these factors covary with both mortality and PM2.5 their inclusion in the models is crucial. We suggest that the unadjusted model is not reflective of the PM2.5-mortality relationship and that the direction or magnitude should not be over-interpreted. Third, although this cohort includes behavioural covariates these are self-reported and in some cases there are missing responses. To avoid introducing bias into the cohort, we used dummy variables to code missing information rather than excluding non-respondents outright. Finally, the cohort itself is limited by follow-up and some persons have as few as 4 years of follow-up (with a maximum follow-up of 15 years).
We provided an update to the Canadian Community Health Survey-Mortality cohort, with a new linkage of the survey respondents to death records, inclusion of additional survey cycles, an extension of the annual residential history and mortality follow-up period, a finer scale of air pollution exposure, time-varying contextual covariates, and the inclusion of immigrants who have lived in Canada for 10–20 years (rather than only those who have been in Canada for 20+ years). The risk of non-accidental mortality from ambient PM2.5 was found even at low levels although the hazard ratio was attenuated in models that included other pollutants (NO2, O3, and OX). The PM2.5-mortality association displayed a supra-linear concentration-response curve. The inclusion of behavioural covariates that could confound the PM2.5-mortality association (fruit and vegetable consumption, leisure exercise frequency, alcohol consumption, and smoking behaviours) did not appear to impact hazard ratios. Hazard ratios were higher for males, those aged 65 or less, and non-immigrants.
Availability of data and materials
The datasets generated and analysed in this study are not publicly available due to privacy and confidentiality standards stated in the Statistics Act.
(− 2) Log likelihood
Akaike information criterion
Aerosol optical depth
Air quality management system
Body mass index
Canadian Census Health and Environment Cohort
Canadian Marginalization Index
Canadian Community Health Survey
Census Metropolitan Area
Environmental Protection Agency
Health Effects Institute
International Classification of Diseases
National Air Pollution Surveillance
- OX :
- PM2.5 :
Fine particle matter
Parts per billion
Schwarz Bayesian Criterion
Shape constrained health impact function
Social data linkage environment
GBD 2017 Risk Factor Collaborators. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392:1923–94.
Beelen R, Raaschou-Nielsen O, Stafoggia M, Andersen ZJ, Weinmayr G, Hoffmann B, Wolf K, Samoli E, Fischer P, Nieuwenhuijsen M. 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.
Di Q, Wang Y, Zanobetti A, Wang Y, Koutrakis P, Choirat C, Dominici F, Schwartz JD. Air pollution and mortality in the Medicare population. N Engl J Med. 2017;376:2513–22.
Pope CA III, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, Thurston GD. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. JAMA. 2002;287:1132–41.
Crouse DL, Peters PA, Hystad P, Brook JR, van Donkelaar A, Martin RV, Villeneuve PJ, Jerrett M, Goldberg MS, Pope CA III, Brauer M, Brook RD, Robichaud A, Menard R, Burnett RT. Ambient PM2.5, O3, and N2 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.
Pinault LL, Weichenthal S, Crouse DL, Brauer M, Erickson A, van Donkelaar A, Martin RV, Hystad P, Chen H, Finès P. Associations between fine particulate matter and mortality in the 2001 Canadian census health and environment cohort. Environ Res. 2017;159:406–15.
Erickson A, Brauer M, Christidis T, Pinault L, Crouse D, Donkelaar A, Weichenthal S, Pappin A, Tjepkema M, Martin R, Brook J, Hystad P, Burnett R. Evaluation of a method to indirectly adjust for unmeasured covariates in the association between fine particulate matter and mortality. Environ Res. 2019;175:108–116.
Shin HH, Cakmak S, Brion O, Villeneuve P, Turner MC, Goldberg MS, Jerrett M, Chen H, Crouse D, Peters P, Pope CA III. Indirect adjustment for multiple missing variables applicable to environmental epidemiology. Environ Res. 2014;134:482–7.
Pinault L, Tjepkema M, Crouse DL, Weichenthal S, van Donkelaar A, Martin RV, Brauer M, Chen H, Burnett RT. Risk estimates of mortality attributed to low concentrations of ambient fine particulate matter in the Canadian Community Health Survey cohort. Environ Health. 2016;15:1–15.
Nasari MM, Szyszkowicz M, Chen H, Crouse D, Turner MC, Jerrett M, Pope CA III, Hubbell B, Fann N, Cohen A. A class of non-linear exposure-response models suitable for health impact assessment applicable to large cohort studies of ambient air pollution. Air Qual Atmos Hlth. 2016;9:961–72.
Statistics Canada. Canadian Community Health Survey: Detailed information for 2001 (Cycle 1.1). 2007. http://www23.statcan.gc.ca/imdb/p2SV.pl?Function=getSurvey&Id=3359. Accessed 28 Jan 2019.
Statistics Canada. Canadian Community Health Survey (CCHS): Detailed information for 2003 (Cycle 2.1). 2007. http://www23.statcan.gc.ca/imdb/p2SV.pl?Function=getSurvey&Id=4995. Accessed 28 Jan 2019.
Statistics Canada. Canadian Community Health Survey (CCHS): Detailed information for 2005 (Cycle 3.1). 2007. http://www23.statcan.gc.ca/imdb/p2SV.pl?Function=getSurvey&Id=22642. Accessed 28 Jan 2019.
Statistics Canada. Canadian Community Health Survey (CCHS): Detailed information for 2007 (Cycle 4.1). 2008. http://www23.statcan.gc.ca/imdb/p2SV.pl?Function=getSurvey&Id=29539. Accessed 28 Jan 2019.
Statistics Canada. Canadian Community Health Survey (CCHS): detailed information for 2010. 2014. http://www23.statcan.gc.ca/imdb/p2SV.pl?Function=getSurvey&Id=81424. Accessed 28 Jan 2019.
Statistics Canada. Canadian Community Health Survey (CCHS): detailed information for 2012. 2012. http://www23.statcan.gc.ca/imdb/p2SV.pl?Function=getSurvey&Id=135927. Accessed 28 Jan 2019.
Statistics Canada. Social data linkage environment (SDLE). 2017. https://www.statcan.gc.ca/eng/sdle/index. Accessed 28 Jan 2019. Accessed 28 Jan 2019.
Statistics Canada. Linkage of the Canadian Community Health Survey (CCHS) to mortality, Cancer, hospital administrative files, and tax data (007–2018). 2019. https://www.statcan.gc.ca/eng/record/2018. Accessed 8 Feb 2019.
Statistics Canada. Directive on Microdata Linkage. 2017. https://www.statcan.gc.ca/eng/record/policy4-1. Accessed 8 Feb 2019.
Fellegi IP, Sunter AB. A theory for record linkage. J Am Stat Assoc. 1969;64:1183–210.
St-Jean H. SDLE production section. Social data linkage environment (SDLE) methodology report: linkage between the Canadian mortality database (CMDB 2015–16) and the SDLE derived record depository (version 18). Statistics Canada: Ottawa; 2018.
Judd P. SDLE production section. Social data linkage environment (SDLE) methodology report: external linkage between the Canadian Community Health Survey (2001 to 2014) and the SDLE derived record depository (version 5). Statistics Canada: Ottawa; 2017.
Bérard-Chagnon J. Comparison of Place of Residence between the T1 Family File and the Census: Evaluation using record linkage. Ottawa: Statistics Canada; 2017. Catalogue no. 91F0015M — No.13
Finès P, Pinault L, Tjepkema M. Imputing postal codes to analyze ecological variables in longitudinal cohorts: exposure to particulate matter in the Canadian Census Health and Environment Cohort Database. Ottawa: Statistics Canada; 2017. Catalogue no. 11–633-X — No. 006
Statistics Canada. Postal CodeOM Conversion File (PCCF). Ottawa: Statistics Canada; 2017. Catalogue no. 92–154-X
Statistics Canada. Postal CodeOM Conversion File (PCCF) Reference Guide August 2015. Ottawa: Statistics Canada; 2016. Catalogue no. 92–154-G
Statistics Canada. Postal CodeOM Conversion File (PCCF) November 2014. Ottawa: Statistics Canada; 2015. Catalogue no. 92–154-G
van Donkelaar A, Martin RV, Spurr RJ, Burnett RT. High-resolution satellite-derived PM2. 5 from optimal estimation and geographically weighted regression over North America. Environ Sci Technol. 2015;49:10482–91.
Meng J, Li C, Martin RV, van Donkelaar A, Hystad P, Brauer M. Estimated long-term (1981−2016) concentrations of ambient fine particulate matter across North America from chemical transport modeling, satellite remote sensing and ground-based measurements. Environ Sci Technol. 2019;53:5071–9.
Robichaud A, Ménard R, Zaïtseva Y, Anselmo D. Multi-pollutant surface objective analyses and mapping of air quality health index over North America. Air Qual Atmos Hlth. 2016;9:743–59.
Robichaud A, Ménard R. Multi-year objective analyses of warm season ground-level ozone and PM2.5 over North America using real-time observations and Canadian operational air quality models. Atmos Chem Phys. 2014;14:1769–800.
Environment and Climate Change Canada. CHRONOS_Ground-Level_O3_NA_2002.nc to CHRONOS_Ground-Level_O3_NA_2009.nc inclusive. Toronto: Air Quality Research Division; 2017. [generated July 2017]
Hystad P, Setton E, Cervantes A, Poplawski K, Deschenes S, Brauer M, van Donkelaar A, Lamsal L, Martin R, Jerrett M, Demers P. Creating national air pollution models for population exposure assessment in Canada. Environ Health Persp. 2011;119:1123–9.
Crouse DL, Philip S, Van Donkelaar A, Martin RV, Jessiman B, Peters PA, Weichenthal S, Brook JR, Hubbell B, Burnett RT. A new method to jointly estimate the mortality risk of long-term exposure to fine particulate matter and its components. Sci Rep. 2016;6:18916.
Statistics Canada. Illustrated Glossary. Ottawa: Statistics Canada; 2017. Catalogue no. 92–195-X
Gordon DL, Janzen M. Suburban nation? Estimating the size of Canada's suburban population. J Archit Plan Res. 2013;30:197–220.
Matheson FI, Dunn JR, Smith KL, Moineddin R, Glazier RH. Development of the Canadian marginalization index: a new tool for the study of inequality. C J Public Health. 2012;103(Suppl 2):12–6.
Bratsch SG. Standard electrode potentials and temperature coefficients in water at 298.15 K. J Phys Chem Ref Data. 1989;18:1–21.
Weichenthal S, Pinault LL, Burnett RT. Impact of oxidant gases on the relationship between outdoor fine particulate air pollution and nonaccidental, cardiovascular, and respiratory mortality. Sci Rep. 2017;7:16401.
Hart JE, Liao X, Hong B, Puett RC, Yanosky JD, Suh H, Kioumourtzoglou MA, Spiegelman D, Laden F. The association of long-term exposure to PM 2.5 on all-cause mortality in the nurses’ health study and the impact of measurement-error correction. Environ Health. 2015;14:38.
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 Resp Crit Care. 2013;187:1226–33.
Burnett R, Chen H, Szyszkowicz M, Fann N, Hubbell B, Pope CA III, Apte JS, Brauer M, Cohen A, Weichenthal S, Coggins J. Global estimates of mortality associated with long-term exposure to outdoor fine particulate matter. Proc Natl Acad Sci U S A. 2018;115:9592–7.
Pinault L, van Donkelaar A, Martin RV. Exposure to fine particulate matter air pollution in Canada. Health Rep. 2017;28:9.
Yin P, Brauer M, Cohen A, Burnett RT, Liu J, Liu Y, Liang R, Wang W, Qi J, Wang L, Zhou M. Long-term fine particulate matter exposure and nonaccidental and cause-specific mortality in a large national cohort of Chinese men. Environ Health Persp. 2017;125:117002.
Burnett RT, Pope CA III, Ezzati M, Olives C, Lim SS, Mehta S, Shin HH, Singh G, Hubbell B, Brauer M, Anderson HR. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environ Health Persp. 2014;122:397–403.
Cohen AJ, Brauer M, Burnett R, Anderson HR, Frostad J, Estep K, Balakrishnan K, Brunekreef B, Dandona L, Dandona R, Feigin V. Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the global burden of diseases study 2015. Lancet. 2017;389:1907–18.
Brauer M. PM2.5 air pollution, population exposed to levels exceeding WHO guideline value (% of total) for the Global Burden of Disease Study 2017 https://data.worldbank.org/indicator/EN.ATM.PM25.MC.ZS?end=2017&locations=CA&start=1990&view=chart. Accessed 24 May 2019.
Apte JS, Marshall JD, Cohen AJ, Brauer M. Addressing global mortality from ambient PM2. 5. Environ Sci Technol. 2015;49:8057–66.
Beiser M. The health of immigrants and refugees in Canada. C J Public Health. 2005;96(Suppl 2):30–44.
Ng E. The healthy immigrant effect and mortality rates. Health Rep. 2011;22:1–5.
Vang ZM, Sigouin J, Flenon A, Gagnon A. Are immigrants healthier than native-born Canadians? A systematic review of the healthy immigrant effect in Canada. Ethn Health. 2017;22:209–17.
Trovato F, Newbold KB. Understanding the healthy immigrant effect: evidence from Canada. In: Trovato F, editor. Migration, health and survival: international perspectives. Cheltenham: Edward Elgar Publishing Ltd; 2017. p. 15–30.
Omariba DW, Ross NA, Sanmartin C, Tu JV. Neighbourhood immigrant concentration and hospitalization: a multilevel analysis of cardiovascular-related admissions in Ontario using linked data. C J Public Health. 2014;105:404–11.
Crouse DL, Christidis T, Erickson A, Pinault L, Martin RV, Tjepkema M, Hystad P, Pappin AJ, Burnett RT, Brauer M, Weichenthal S. Evaluating the Sensitivity of PM2.5-Mortality Associations to the Spatial and Temporal Scale of Exposure Assessment at Low Particle Mass Concentrations. Epidemiology. Accepted.
Statistics Canada. Canadian Community Health Survey (CCHS): Annual component User guide 2010 and 2009-2010 Microdata files. http://www23.statcan.gc.ca/imdb-bmdi/pub/document/3226_D7_T9_V8-eng.pdf. Accessed 23 May 2019.
Weichenthal S, Villeneuve PJ, Burnett RT, van Donkelaar A, Martin RV, Jones RR, DellaValle CT, Sandler DP, Ward MH, Hoppin JA. Long-term exposure to fine particulate matter: association with nonaccidental and cardiovascular mortality in the agricultural health study cohort. Environ Health Persp. 2014;122:609–15.
The authors would like to acknowledge the contribution of Hong Chen for his code to run the SCHIF and the Canadian Urban Environmental Health Research Consortium (CANUE) for supplying an ozone-postal code linked file.
Research described in this article was conducted under contract to the Health Effects Institute (HEI), an organization jointly funded by the United States Environmental Protection Agency (EPA; Assistance Award No. R-82811201) and certain motor vehicle and engine manufacturers. The contents of this article do not necessarily reflect the views of HEI, or its sponsors, nor do they necessarily reflect the views and policies of the EPA or motor vehicle and engine manufacturers.
Ethics approval and consent to participate
Canadian Community Health Survey respondents who agreed to linkage of their survey responses with administrative and tax data were eligible. The linkage (007–2018) was approved by Statistics Canada’s senior management (https://www.statcan.gc.ca/eng/record/2018) and is governed by the Directive on Microdata Linkage (https://www.statcan.gc.ca/eng/record/ policy4–1).
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Christidis, T., Erickson, A.C., Pappin, A.J. et al. Low concentrations of fine particle air pollution and mortality in the Canadian Community Health Survey cohort. Environ Health 18, 84 (2019). https://doi.org/10.1186/s12940-019-0518-y