Long-Term Exposure to Constituents of Fine Particulate Air Pollution and Mortality: Results from the California Teachers Study
Bart Ostro1, Michael Lipsett2, Peggy Reynolds3, Debbie Goldberg3, Andrew Hertz3, Cynthia Garcia4, Katherine D. Henderson5, Leslie Bernstein5
1 California Environmental Protection Agency, Office of Environmental Health Hazard Assessment, Oakland, California, USA, 2 California Department of Public Health, Richmond, California, USA, 3 Northern California Cancer Center, Berkeley, California, USA, 4 California Air Resources Board, Sacramento, California, USA, 5 City of Hope, Duarte, California, USA
Environ Health Perspect 118:363-369 (2010). http://dx.doi.org/10.1289/ehp.0901181 [online 26 October 2009]
Background: Several studies have reported associations between long-term exposure to ambient fine particulate matter (PM) and cardiovascular mortality. However, the health impacts of long-term exposure to specific constituents of PM2.5 (PM with aerodynamic diameter ≤ 2.5 µm) have not been explored.
Methods: We used data from the California Teachers Study, a prospective cohort of active and former female public school professionals. We developed estimates of long-term exposures to PM2.5 and several of its constituents, including elemental carbon, organic carbon (OC), sulfates, nitrates, iron, potassium, silicon, and zinc. Monthly averages of exposure were created using pollution data from June 2002 through July 2007. We included participants whose residential addresses were within 8 and 30 km of a monitor collecting PM2.5 constituent data. Hazard ratios (HRs) were estimated for long-term exposure for mortality from all nontraumatic causes, cardiopulmonary disease, ischemic heart disease (IHD), and pulmonary disease.
Results: Approximately 45,000 women with 2,600 deaths lived within 30 km of a monitor. We observed associations of all-cause, cardiopulmonary, and IHD mortality with PM2.5 mass and each of its measured constituents, and between pulmonary mortality and several constituents. For example, for cardiopulmonary mortality, HRs for interquartile ranges of PM2.5, OC, and sulfates were 1.55 [95% confidence interval (CI), 1.431.69], 1.80 (95% CI, 1.681.93), and 1.79 (95% CI, 1.582.03), respectively. Subsequent analyses indicated that, of the constituents analyzed, OC and sulfates had the strongest associations with all four outcomes.
Conclusions: Long-term exposures to PM2.5 and several of its constituents were associated with increased risks of all-cause and cardiopulmonary mortality in this cohort. Constituents derived from combustion of fossil fuel (including diesel), as well as those of crustal origin, were associated with some of the greatest risks. These results provide additional evidence that reduction of ambient PM2.5 may provide significant public health benefits.
Key words: cardiopulmonary mortality, chronic exposure, cohort study, elemental carbon, fine particles, organic carbon, PM2.5, species, sulfates
Address correspondence to B. Ostro, California EPA, Office of Environmental Health Hazard Assessment, 1515 Clay St., 16th Floor, Oakland, CA 94612 USA. Telephone: (510) 622-3157. Fax: (510) 622-3210. E-mail: email@example.com
This work was supported in part by California ARB contract 03-313 and National Cancer Institute grant R01 CA077398.
The authors declare they have no competing financial interests.
R. Broadwin (Office of Environmental Health Hazard Assessment) and K. Turkiewicz and C. Taylor [California Air Resources Board (ARB)] provided assistance with PM2.5 species air quality data compilation.
Received 09 July 2009; accepted 26 October 2009; online 26 October 2009.
The opinions expressed in this article are solely those of the authors and do not represent the policy or position of the State of California, the California Environmental Protection Agency, or the California Department of Public Health.
Several cohort studies have provided evidence linking total and cardiovascular mortality to long-term exposure to fine particulate matter (PM) air pollution [PM ≤ 2.5 µm in diameter (PM2.5)] (Abbey et al. 1999; Dockery et al. 1993; Hoek et al. 2002; Laden et al. 2006; Miller et al. 2007; Pope et al. 1995, 2002, 2004). Ambient PM2.5 represents a heterogeneous mixture of solid and liquid particles generated by many sources, however, and there is little evidence for which constituents of PM2.5 are associated with the greatest risks. Among the available evidence from long-term exposure studies, Pope et al. (1995) and Dockery et al. (1993) both indicated that exposure to fine PM sulfate (SO4), likely generated from combusted fossil fuel, was associated with cardiopulmonary mortality.
The U.S. National Research Council has highlighted the importance of investigating characteristics and constituents of PM that contribute to their toxicity (National Research Council 2004). Data regarding differential PM2.5 constituent toxicity can have important implications for establishing ambient air quality standards, targeting control strategies, and enhancing the accuracy of health impact assessments. Routine collection of speciated PM2.5 data now provides opportunities to examine relationships between long-term exposure to specific PM constituents and morbidity and mortality.
We assessed whether PM2.5 constituents were associated with total and cardiopulmonary mortality among participants in the California Teachers Study (CTS), a prospective cohort study of female public school teachers and administrators, initiated in 1995. Within this cohort, the wealth of individual-level data collected allows for control of numerous potential confounders. Active smoking at cohort inception was very low (5%), and indoor occupational exposures among the cohort have generally been low, except to the extent that a small percentage of schools might be located near highways. Thus, residual confounding from smoking or occupational exposures is unlikely. Finally, the large sample size allowed us to restrict the study population to those within close proximity of sites that monitor PM2.5 species, thereby reducing potential impacts of pollution measurement error.
Materials and Methods
Study population. The CTS is a prospective cohort investigation of 133,479 current and former female public school professionals who completed baseline questionnaires mailed to enrollees in the California State Teachers Retirement System in 1995. Subsequent questionnaires were mailed to CTS participants in 1997 and 2000. Annual follow-up includes updating name and residential information of CTS members and outcome linkages. Participants’ ages varied from 22 to 104 years at enrollment, with a median of 54 years. The cohort is primarily white (86.7%) and born in the United States (93.6%). A full description of the CTS cohort is available elsewhere (Bernstein et al. 2002). All residential addresses from study enrollment forward were geocoded and linked with pollutant data to generate estimates of long-term exposure. Use of data involving human subjects in the CTS was approved by the Committee for the Protection of Human Subjects, California Health and Human Services Agency, and by the institutional review board at each participating institution.
Health data. CTS records are linked annually to government-maintained mortality and hospitalization files. Data on mortality were obtained from the California Department of Health Services (2008), the U.S. Social Security Administration death master file (U.S. Department of Commerce, 2008), and the National Death Index (Centers for Disease Control and Prevention 2008). Data on the underlying cause of death from August 2002 through July 2007 were used in our analyses of four mortality categories: all-cause mortality [excluding external causes; all International Classification of Diseases, 10th Revision (ICD-10) (World Health Organization 1993), codes except S through Z], cardiopulmonary mortality (ICD-10 codes I00I99 and J00J98), mortality due to ischemic heart disease (IHD; ICD-10 codes I20I25), and mortality due to pulmonary disease (ICD-10 codes J00J98). Following the practice of the National Center for Health Statistics, coding of cause of death for CTS participants is automated, adhering to World Health Organization rules by using the Automated Classification of Medical Entities system, which eliminates intercoder differences that arise with manual coding. We excluded women who lived outside of California at baseline or whose causes of death were unknown. We restricted our analysis to participants > 30 years of age.
We calculated person-months at risk based on the number of days between 1 August 2002 and a) a woman’s date of death, b) 31 July 2007, or c) the date of first non-California address for women who moved out of state before either a) or b) occurred.
Air pollution exposure estimates. PM2.5 and species were collected and analyzed by laboratories contracted by the U.S. Environmental Protection Agency (EPA) as part of the Speciation Trends Network (STN) (U.S. EPA 2008). The 24-hr averaged measurements were usually obtained on an every third- or sixth-day basis. Based on prior evidence of associations from time-series studies (Mar et al. 2000; Ostro et al. 2007) and from other epidemiologic or toxicologic studies, we examined PM2.5 mass and the following eight constituents: elemental carbon (EC), organic carbon (OC), SO4, nitrate (NO3), iron (Fe), potassium (K), silicon (Si), and zinc (Zn). Filters were analyzed by EPA staff for EC and OC using the total optical transmittance method; for SO4, NO3, and K using ion chromatography; and for trace elements using X-ray fluorescence. One monitor in each of the following eight counties collected data on PM2.5 and these constituents: Fresno, Kern, Los Angeles, Riverside, Santa Clara, San Diego, Sacramento, and Ventura. All monitors were operating as of 1 June 2002, and data collection for this study ended on 31 July 2007. After that date, changes in collection methods for EC and OC were initiated at several of the sites, rendering comparability with previous measures uncertain.
Each subject was assigned a monthly exposure value based on the monitor nearest her geocoded residential address. Months were included when at least 75% of the potential observations were available. For each individual and each pollutant, the values for all person-months of exposure were summed and then divided by the total months of exposure to create an average measure of overall long-term exposure. Two different exposure assessments were examined. First, we restricted the analyses to subjects whose residences were within 8 km of a monitor. Second, we restricted the analyses to subjects whose residences were within 30 km of a monitor. The smaller distance was selected to minimize potential exposure misclassification while providing a sample size at least as large as that of Harvard Six City (H6C) Study (Dockery et al. 1993), which was among the smallest of the major published air pollution cohort studies. The larger distance was used a) to examine the impact on the risk estimates of different distances from the monitor, and therefore different degrees of measurement error, and b) to increase the overall study power. Both buffer distances included some participants from counties adjacent to those in which monitors were located.
Covariates. For covariates, we selected most of the individual-level predictor variables for the regression analysis based on risk factors identified in previous observational studies of the effects of air pollution on cardiovascular and respiratory disease (Dockery et al. 1993; Jerrett et al. 2005; Pope et al. 2002). Sixteen individual-level covariates (a total of 47 terms) were included in the model: marital status; smoking status and pack-years of smoking; second-hand smoke exposure; body mass index (BMI); lifetime physical activity; alcohol consumption; average daily dietary intake of fat, fiber, and calories; menopausal status; hormone replacement therapy use; family history of myocardial infarction or stroke; and use of blood pressure medication or aspirin. Women with BMI < 16 or ≥ 54.9 (< 0.5% of the women) were excluded. In the Cox regression models, described below, the sample was stratified by age (divided into 2-year categories between ages 30 and 79, 3-year categories between ages 80 and 88, and one category for women ≥ 89 years of age) and by race/ethnicity (categorized into three groups: non-Hispanic white, all others, and unknown).
In addition to the individual-level covariates, six ecologic variables were added to the model. These variables, obtained at the census block group level from the 2000 Census (U.S. Census Bureau 2002), were included to consider contextual or neighborhood confounding. The variables, selected based on results from prior analysis (Jerrett et al. 2005), included income (median household income), income inequality (percent living in poverty), education (percent with bachelor’s degree and above), block group population, racial composition (percent black, white, and Hispanic), and unemployment (percent > 16 years of age who were unemployed).
Statistical methods. The statistical analyses were limited to participants who a) were living in California and were at least 30 years of age at baseline, b) lived at addresses that were successfully geocoded within either 8 or 30 km of one of the study monitors, c) had information available on all continuous variables used in the statistical models, and d) did not leave the state for more than four months during the study period.
We used Cox proportional hazards models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for each pollutantoutcome combination. Age at the start and end of follow-up (in days) was used to define time on study. We began the exposure period on 1 June 2002—2 months before the start of the cohort follow-up—to allow for analysis of potential impacts from prior exposures. In the analysis, each constituent was examined separately. HRs and the associated 95% CIs were scaled to the interquartile range (IQR), based on the average distributions for each pollutant for all women. Three different models were run, all adjusted for age and race. The first model included only a pollutant variable, each entered in a separate regression that was otherwise unadjusted. A second model included the full set of individual-level risk factors summarized above. The final model also included all of the contextual variables to test for the influence of local geographic impacts. Because there was little difference in the risk estimates for the unadjusted models versus those adjusted for individual covariates and those adjusted for both individual covariates and contextual risk factors, we present only the results for the fully adjusted models. However, the results of all three models using OC as the pollutant are presented in the Supplemental Material [see Table S-1 (doi:10.1289/ehp.0901181.S1 via http://dx.doi.org)]. As indicated above, we examined the risks for participants within either 8 or 30 km of the PM2.5 species monitors. Also, as a sensitivity analysis, we conducted a forward-selection regression procedure of the PM2.5 constituents to determine their relative importance. Because of the likelihood of multicollinearity, we constrained the selection to include only positive coefficients with p < 0.05. Several other sensitivity analyses were conducted using cardiopulmonary mortality as the end point. First, we examined several two-pollutant models focusing on those pollutants that appeared to be most strongly associated with mortality. Second, we examined the implications of using a full year of exposure data (rather than only 2 months) before the follow-up period. Thus, in this analysis, exposures began in June 2002 but cohort follow-up began in June 2003 (rather than August 2002 as in the base case). Both the exposure and cohort follow-up continued through July 2007. Finally, we examined the effects of using an exposure estimate consisting of only a single annual average for each pollutant, similar to some previous studies (Miller et al. 2007). In this case, the exposures consisted of a given pollutant’s annual average for the period June 2002 through May 2003 based on the monitor at the participants’ residence at the start of the follow-up period (June 2003). Proportional hazards analyses were conducted using SAS software (version 9.1; SAS Institute Inc., Cary, NC).
A total of 9,208 participants resided within 8 km of a monitor, and 52,226 within 30 km. About 5% of each group were < 30 years of age, 5% moved or died before August 2002, 3% were missing individual-level risk factors, and 1% were missing contextual risk factors. This left totals of 7,888 and 44,847 participants within 8 and 30 km of a monitor, respectively. Table 1 summarizes the pollutant variables used in the analyses and provides the mean, minimum, maximum, and IQR for each. Measures are based on individual-level estimates of exposure and do not represent the average readings for any specific monitor or set of monitors. The long-term mean PM2.5 concentrations for the participants within 8 km and 30 km were 17.0 and 17.5 µg/m3, respectively. For each buffer designation, OC constituted about one-third and NO3 about one-fourth of the total PM2.5 mass.
Table 2 summarizes some of the individual baseline characteristics of the CTS participants eligible for the analysis and living within the 8- and 30-km buffers, as well as their residential locations, by county. Individual characteristics (e.g., age, BMI, smoking status, marital status, menopausal status, hormone therapy use, and family history of heart disease) were similar for the two buffer configurations. Thirteen percent of all respondents within the 8-km buffers and 36% within the 30-km buffers of the species monitors resided in Los Angeles County.
Table 3 summarizes the correlations of the individual-level exposures to the constituents of PM2.5 based on the 8- and 30-km buffers around the species monitors. We observed relatively high correlations (r > 0.8) between PM2.5 and EC, NO3, and Zn. Among the constituents, we observed high correlations (r > 0.8) between EC and NO3, Fe, and Zn; between OC and K; and between Zn and NO3 and Fe.
The results of the regression analyses are presented in Tables 4 and 5 for the 8-km and 30-km buffers, respectively. For the 8-km buffer, 58% of all deaths were from cardiopulmonary disease versus 54% for the 30-km buffer. In addition, 6.8% of all women in the 8-km buffer and 5.8% of women in 30-km buffer died during the 5-year follow-up period. In the regression analysis, we entered each constituent separately into the model. Using a buffer of 8 km around each species monitor, we observed associations with all-cause, cardiopulmonary, and IHD mortality for PM2.5 mass and all of the constituents. For pulmonary mortality, we observed significant associations for OC, SO4, and Si but not PM2.5. We generally observed modestly higher HRs for cardiopulmonary versus all-cause mortality for PM2.5 mass and several constituents, whereas those for IHD were markedly greater. In addition, in micrograms per cubic meter, all of the constituents had greater HRs than did PM2.5 mass (data not shown).
Using the 30-km buffer around each PM2.5 species monitor (Table 5), we again observed associations between most of the constituents and all-cause, cardiopulmonary, and IHD mortality. EC was not associated with all-cause or cardiopulmonary mortality but was associated with IHD mortality. For pulmonary mortality, we observed associations for PM2.5 mass, OC, SO4, NO3, and Si. With few exceptions, the HRs were similar to or slightly less than those observed with the 8-km buffer.
Using the 8-km buffer, we conducted a forward-selection regression procedure to determine the relative importance of the constituents. For all four outcomes, only two constituents, OC followed by SO4, successfully entered the model. Forcing additional constituents into the model generated clear signs of multicollinearity, including widely varying (and high) HRs and negative β-coefficients. When PM2.5 was also included as a candidate variable for inclusion, OC and SO4 were still selected first in the model.
Table 6 summarizes the results of two-pollutant models for cardiopulmonary mortality using the 8-km buffer and including each of the other constituents in the model specification with either OC or SO4. In the models with OC plus another constituent, OC remained associated with cardiopulmonary mortality, with a fairly consistent HR. Likewise, in the models with SO4 plus another constituent, SO4 was consistently associated with cardiopulmonary mortality with an HR similar to models in which SO4 was the only pollutant.
We conducted an additional sensitivity analysis using exposures starting a full year before the beginning of the follow-up period. The exposures were calculated monthly from June 2002 through July 2007, whereas the cohort follow-up period was June 2003 through July 2007. The results were essentially unchanged from the basic model when the follow-up period began in August 2002, as summarized in Table 4 [see Supplemental Material, Table S-2 (doi:10.1289/ehp.0901181.S1)]. The two exceptions were reductions in HRs for OC (HR = 1.46; 95% CI, 1.241.71) and K (HR = 1.43; 95% CI, 1.131.82), with the latter becoming more consistent with the other estimates. In our final sensitivity analysis, we considered a single annual average (June 2002 through May 2003) as our exposure metric using the reported residence of each study participant for June 2003. We used a cohort follow-up period similar to the previous analysis (i.e., June 2003 through July 2007). As summarized in Figure 1, the HRs for cardiopulmonary mortality were reduced for each of the constituents (except for EC, for which the HR was similar), and we observed significant associations only for PM2.5, SO4, NO3,and Zn.
In this study we found strong and consistent associations between long-term exposure to PM2.5 mass, as well as several of its constituents, and all-cause, cardiopulmonary, IHD, and pulmonary mortality. Specifically, in single-pollutant models using the 8-km buffer, PM2.5 mass and all of the constituents were associated with all-cause, cardiopulmonary, and IHD mortality. For pulmonary mortality, we observed associations for OC, SO4, NO3, and Si but not for PM2.5 mass. When a buffer size of 30 km was used, we observed generally similar results except that EC was no longer associated with all-cause or cardiopulmonary mortality and PM2.5 was associated with pulmonary mortality. Forward-selection regression analysis indicated that long-term exposures to OC and SO4, in particular, had stronger associations with all four outcomes than did the other constituents or PM2.5 mass. Subsequent analysis of cardiopulmonary mortality using two-pollutant models that included either OC or SO4 plus one other constituent provided additional support for the importance of these two pollutants, because their HRs were robust to inclusion of other pollutants in the model.
We focused particular attention on attempting to improve the exposure assessment by using information on monthly residential history of the cohort and by limiting the sample to 8- or 30-km buffers around each monitor. Comparing results between the 8- and 30-km buffers, we observed lower HRs and a small reduction in the number of constituents associated with the mortality outcomes, which would be expected from increased exposure measurement error. We also conducted several sensitivity analyses of our exposure metric. The results differed little when we started the exposure period 2 months or 12 months before the initiation of cohort follow-up. However, we observed significant differences in the results when participants were all assigned a single annual average exposure based on their residential location in the first month of the follow-up period. It is not clear whether this result would hold for other data sets, but it merits further examination.
To our knowledge, only a few previous studies of long-term exposure have examined the effects of any constituents of PM2.5. In analyses of the cohorts from the ACS and the H6C Study, Pope et al. (1995, 2002) and Dockery et al. (1993) reported associations between SO4 and both all-cause and cardiopulmonary mortality. In addition, a cross-sectional study of mortality in the United States reported associations between metropolitan area-wide mortality rates and SO4 (Özkaynak and Thurston 1987). Using both the 8-km and 30-km buffer, we also found associations between SO4 and all four outcomes. In addition, in our forward stepwise selection regressions, SO4 was selected to be included in the model for all four outcomes. Thus, long-term exposure to SO4 is consistently associated with mortality even at the relatively low concentrations observed in this study (mean of 2.5 µg/m3 using the 30-km buffer) relative to those reported by Pope et al. (1995) and Dockery et al. (1993) (~ 11 and 8 µg/m3, respectively). Besides low-sulfur motor vehicle fuels and diesel-powered ships, there are relatively few sources of sulfur emissions in California. Therefore, SO4 makes up a much smaller proportion of PM2.5 in California than in the eastern or midwestern United States (Bell et al. 2007).
Regarding other constituents of PM2.5, a Dutch cohort study included an examination of relationships between mortality and long-term exposure to black smoke, a primary pollutant measured by reflectance on a Teflon filter, which has been shown to be highly correlated with EC (Beelen et al. 2008; Cyrys et al. 2003). In that study, associations were reported between black smoke and all-cause mortality and less so with cardiovascular mortality. In the CTS cohort using the 8-km buffer, we also observed associations for EC. However, we observed no association using the 30-km buffer, which may have been due in part to significant exposure misclassification. In general, researchers have observed exponential declines in EC concentrations with downwind distance from roadways. For instance, in Los Angeles EC drops to background levels within 100150 m downwind from busy roads (Zhu et al. 2002). Wind speed and direction can significantly affect EC concentrations; PM2.5 concentrations 150 m from the highway can be reduced by 50% depending on the wind direction (Hitchins et al. 2000).
We observed associations for OC with all four outcomes that we examined using both buffers. Moreover, using a forward-selection procedure, OC was the first pollutant selected for model inclusion for all-cause, cardiopulmonary, and IHD mortality (even when PM2.5 was included as a candidate pollutant). OC is both directly emitted and, after transformation through atmospheric chemical reactions, a secondary product of fuel combustion. Key sources in California include gasoline and diesel vehicles, residential wood combustion, agricultural and prescribed burning, and industrial combustion of fossil fuels. Diesel PM2.5 emissions consist of both OC and EC fractions, along with trace amounts of inorganic compounds (Abu-Allaban et al. 2004). The OC fraction of diesel exhaust contains heavy hydrocarbons, such as lubricating oils and non- and semivolatile polycyclic aromatic hydrocarbons. Although no previous epidemiologic studies of long-term exposure have included OC, there is evidence of serious health effects related to short-term exposures. For example, in studies of daily exposures in six California counties and in Phoenix, Arizona (Mar et al. 2000; Ostro et al. 2007), associations were reported between OC and cardiovascular mortality. In another study of short-term exposures, Metzger et al. (2004) reported an association between OC and emergency department visits for cardiovascular disease in Atlanta, Georgia. Possible pathophysiologic mechanisms linking OC with mortality, especially from IHD, may include oxidative stress and effects on blood pressure (Brook et al. 2004).
In single-pollutant models, we also observed associations of all-cause, cardiopulmonary, and IHD mortality with NO3, K, Fe, and Zn. K is often used as a marker of biomass combustion and residential wood burning, which tend to be a local and neighborhood problem in several parts of California (Lipsett et al. 1997). Because of high correlations with other pollutants in our study, however, it is difficult to determine whether long-term exposures to these constituents have independent effects on health. Previous toxicologic and epidemiologic studies, however, provide support for significant cardiovascular effects from exposure to all four of these constituents (Brook et al. 2004; Fairley 2003; Hoek 2003; Naeher et al. 2007; Ostro et al. 2007).
We observed strong associations between Si and several mortality categories, especially pulmonary mortality. Si is a crustal element that is a large component of soil and resuspended road dust. As such, it may be enriched by and serve as a surrogate for many toxic constituents found in road dust, including combustion-based material, brake dust, tire debris, and semivolatile compounds (Rogge et al. 1993). It may also serve as a general marker for proximity to traffic. Several studies of long-term exposure to traffic or traffic-based pollutants have reported associations with cardiopulmonary mortality (Beelen et al. 2008; Filleul et al. 2005; Jerrett et al. 2005) and with the development of cardiovascular disease (Hoffmann et al. 2006). Finally, a few studies have reported cardiovascular effects from direct exposure to Si. For example, Si has been found to be associated with a) ST-segment elevation, a marker for myocardial ischemia, in dogs; b) vasoconstriction of pulmonary arteries in rats; c) heart rate variability in humans; and d) cardiovascular mortality in time-series studies in Arizona and California (Batalha et al. 2002; Cavallari et al. 2008; Lipsett et al. 2006; Mar et al. 2000; Ostro et al. 2007; Wellenius et al. 2003). In contrast, in a study of six eastern and midwestern cities in the United States, the authors did not observe an association between mortality and daily exposure to Si (Laden et al. 2000). Thus, whereas it is likely that Si serves as a proxy either for the toxic constituents found in road dust or for exposures to traffic-related pollutants, the possibility of direct cardiotoxicity, although unlikely at the low concentrations observed in this investigation, cannot be ruled out.
Our findings of an association between long-term exposure to PM2.5 and all-cause and cardiovascular mortality in women are consistent with those reported in previous studies (Chen et al. 2005; Dockery et al. 1993; Eftim et al. 2008; Krewski et al. 2000; Miller et al. 2007; Pope et al. 1995, 2002, 2004). Several studies reported similar relative risks for females and males. However, our estimated HRs for the all-female teacher cohort were generally higher than results present for most other female cohorts, with the exception of the recent study of the Women’s Health Initiative (WHI) observational study (Miller et al. 2007). For cardiopulmonary mortality in females associated with a 10 µg/m3 change in PM2.5, HRs of 1.16 (95% CI, 1.081.27), 1.20 (95% CI, 1.091.32), and 1.76 (95% CI, 1.252.47) were reported for cohorts from the ACS, H6C Study, and WHI (for cardiovascular mortality), respectively (Dockery et al. 1993; Krewski et al. 2000; Miller et al. 2007; Pope et al. 1995). In contrast, extrapolating from the 30-km buffer results in the present study, a 10-µg/m3 change in PM2.5 is associated with an HR of 2.05 (95% CI, 1.802.36). There are several possible explanations for the higher estimate in this study. First, it may be attributable to efforts to improve the exposure assessment, because our study incorporated monthly residential history and small buffers around each pollution monitor. Second, the overall exposures may be greater because of the mild climate in California, where people are likely to have their windows open and spend more time outdoors. Finally, the greater relative risk estimates may be simply attributable to stochastic variability.
We observed particularly high relative risks from exposure to PM2.5 mass and several of its constituents for mortality from IHD. For example, using the 8-km buffer, the HRs associated with the IQRs for PM2.5, OC, and SO4 were 2.10, 2.02, and 1.82, respectively. Several other studies have reported associations for this outcome, as well (Jerrett et al. 2005; Miller et al. 2007; Pope et al. 2002).
This study is subject to several potential limitations. First, the study cohort may not be representative of the general population of California women or the adult female population of the United States. Likewise, the mix of pollutants in California is different from that observed in the U.S. Midwest and East Coast, because traffic is the dominant source of PM in California and there are few emissions from major industrial facilities such as coal-fired power plants. Third, we were limited to contemporaneous exposures for the PM2.5 species, which could bias the ultimate effect estimates if historical exposures are more important. Fourth, we based exposure assessment only linkage of nearest fixed-site monitors with residence; exposures derived from sources such as motor vehicle exhaust while commuting or indoor environments at homes and schools could not be accounted for in this data set. Fifth, exposure misclassification was likely lower within the smaller buffer radius, but this came at the cost of reduced statistical power. Few IHD and pulmonary events occurred among women within the 8-km buffer, which produced somewhat imprecise estimates. Finally, some caution is warranted in attributing a causal relationship to any single constituent because its effect estimate may be attributable to its intrinsic toxicity or its correlation with other, more toxic substances. In addition, the constituents may have differential measurement error due to differences in spatial variability and in indoor/outdoor penetration. Additional research using more years for this cohort, when available, or from other studies will be necessary to confirm the importance of the constituents associated with higher HRs in this investigation.
In summary, we analyzed a relatively homogeneous cohort of female teachers and administrators, coupled with spatial restrictions to reduce the effects of exposure misclassification. The results provide evidence for effects from long-term exposures to PM2.5 and several of its constituents on mortality. The constituents that appear to generate most of the risk are derived from combustion of fossil fuel (including diesel) and biomass, as well as from PM of crustal origin. Reduction of ambient PM2.5, particularly from fuel combustion, may provide significant public health benefits.
In an article published in Environmental Health Perspectives (Ostro et al. 2010), we analyzed the relationships of long-term exposure to fine particulate matter (≤ 2.5 µm in aerodynamic diameter; PM2.5) and its components with mortality in a cohort of > 100,000 active and retired female professionals participating in the California Teachers Study (CTS) cohort. We used a Cox proportional hazards model in which pollution exposure was measured as a continuous variable over the study period. Monthly average pollutant concentrations were obtained for each participant from measurements at the nearest PM2.5 monitor within either 8 or 30 km of her geocoded residential address. Each participant was assigned a single exposure value over the follow-up period, defined as the average pollutant concentration from the beginning of the observation period (1 June 2002) to the woman’s date of death, loss to follow-up, or study termination (31 July 2007). Thus, exposure assignment was dependent on the duration of follow-up for each participant.
In our article (Ostro et al. 2010), we reported associations of mortality from all causes, cardiopulmonary disease, and ischemic heart disease (IHD) with PM2.5 mass and several of its components. However, the estimated hazard ratios (HRs) were generally higher than those reported from previous cohort studies (Dockery et al. 1993; Eftim et al. 2008; Krewski et al. 2009; Laden et al. 2006; Pope et al. 1995). Part of this difference was likely due to the nature of the exposure assignment. Most previous cohort studies have assigned the same exposure period to all study subjects, regardless of when deaths occurred. Thus, estimated exposures for some study participants in several studies occurred after their deaths. In addition, exposures have usually been assigned to participants based on their residential address at enrollment only, without taking into account exposure changes that may have occurred throughout the study period or when participants relocated. Finally, many previous studies measured exposure for only a subset of the years during which the cohort was followed. In an effort to reduce these aspects of exposure misclassification, we estimated exposures beginning prior to the cohort follow-up period, continuing to the end of the study or until the participant died or relocated out of state, incorporating updated exposure assignments when the subjects moved.
Importantly, measured concentrations of several pollutants in California declined substantially from 2002 through 2007; annual average PM2.5, organic carbon (OC), and nitrates decreased by around 30% each. These marked decreases in ambient PM2.5 concentrations resulted in lower average exposure estimates for cohort members who survived to the end of our study. Thus, the exposure assigned to a participant who died at time t would tend to be greater for events occurring early in the observation period, compared with the long-term average exposures of the participants who comprised the remainder of the risk set (i.e., those who were still part of the cohort study at time t and who subsequently experienced lower ambient pollution levels).
We have reanalyzed the CTS data using time-dependent pollution metrics—in which the exposure estimates for everyone remaining alive in the risk set were recalculated at the time of each death—in order to compare their average exposures up to that time with that of the individual who had died. In this way, decedents and survivors comprising the risk set had similar periods of pollution exposure, without subsequent pollution trends influencing the surviving women’s exposure estimates.
As in our previous study (Ostro et al. 2010), we restricted the sample in this reanalysis to women living within 30 km of one of eight fixed-site monitors in the U.S. Environmental Protection Agency’s Speciation Trend Network (STN), resulting in a study population of almost 44,000 women. Residential addresses from study enrollment forward were geocoded and linked with monthly pollutant averages at the nearest STN monitor to generate estimates of long-term exposure. We also used the same set of individual and ecological covariates in a Cox proportional hazards model as was used in the original study. Pollutants entered separately into the model included PM2.5 mass, elemental carbon (EC), OC, sulfate, nitrate, iron, potassium, silicon, and zinc. We used data on primary cause of death from August 2002 through July 2007 to examine the relationships between pollutants and mortality from all causes and cardiopulmonary, pulmonary, and IHDs.
The results are summarized in Erratum Table 1, scaled to the interquartile range (IQR) for each pollutant. HRs were significantly attenuated from our previous results. No associations were observed between all-cause mortality and PM2.5 or its components. For cardiopulmonary mortality, we observed significant associations for PM2.5 mass, nitrate, sulfate, and silicon, with more modest associations for zinc. PM2.5 mass and all of its components were associated with mortality from IHD, whereas none of the pollutants was associated with pulmonary mortality. This Erratum Table 1 should replace Table 5 in our previous article (Ostro et al. 2010).
Compared with our previous results (Ostro et al. 2010), these updated PM2.5 HRs are more consistent with several other published estimates of mortality risks, which are scaled to an increment of 10 μg/m3 of long-term average PM2.5 and summarized in Erratum Table 2. For example, relative to our revised HR of 1.19 for cardiopulmonary disease, analogous HRs from previous studies include 1.09 (95% CI, 1.03–1.16) from the American Cancer Society–Cancer Prevention II (ACS) cohort (cardiopulmonary disease; Pope et al. 2004), 1.28 (95% CI, 1.13–1.44) from the Harvard Six Cities study (cardiovascular disease; Laden et al. 2006), and 1.10 (95% CI, 0.94–1.28) from the Los Angeles subcohort of the ACS study (cardiopulmonary disease; Jerrett et al. 2005). Much higher HRs were observed in the observational study of the Women’s Health Initiative cohort for cardiovascular and IHD mortality (Miller et al. 2007).
These revised results still support the existence of elevated risks of PM2.5-associated cardiopulmonary disease and IHD, and illustrate the importance of considering the impact of long-term pollution trends in modeling estimates of exposure.
The authors declare they have no actual or potential competing financial interests.
Bart Ostro California Office of Environmental Health Hazard Assessment, Oakland, California, E-mail: Bostro@Creal.cat
Peggy Reynolds, Debbie Goldberg, Andrew Hertz Cancer Prevention Institute of California, Berkeley, California
Richard T. Burnett, Hwashin Shin Health Canada, Ottawa, Ontario, Canada
Edward Hughes Edward Hughes Consulting, Ottawa, Ontario, Canada
Cynthia Garcia California Air Resources Board Sacramento, California
Katherine D. Henderson, Leslie Bernstein City of Hope, Duarte, California
Michael Lipsett California Department of Public Health, Richmond, California
Related EHP Articles
- Abbey DE, Nishino N, McDonnell WF, Burchette RJ, Knutsen SF, Lawrence Beeson W, et al. 1999. Long-term inhalable particles and other air pollutants related to mortality in nonsmokers. Am J Respir Crit Care Med 159(2):373–382.
- Abu-Allaban M, Rogers CF, Gertler AW. 2004. A quantitative description of vehicle exhaust particle size distributions in a highway tunnel. J Air Waste Manag Assoc 54(3):360–366.
- Batalha JR, Saldiva PH, Clarke RW, Coull BA, Stearns RC, Lawrence J, et al. 2002. Concentrated ambient air particles induce vasoconstriction of small pulmonary arteries in rats. Environ Health Perspect 110:1191–1197.
- Beelen R, Hoek G, van den Brandt PA, Goldbohm RA, Fischer P, Schouten LJ, et al. 2008. Long-term effects of traffic-related air pollution on mortality in a Dutch cohort (NLCS-AIR study). Environ Health Perspect 116:196–202.
- Bell ML, Dominici F, Ebisu K, Zeger SL, Samet JM. 2007. Spatial and temporal variation in PM2.5 chemical composition in the United States for health effects studies. Environ Health Perspect 115:989–995.
- Bernstein L, Allen M, Anton-Culver H, Deapen D, Horn-Ross PL, Peel D, et al. 2002. High breast cancer incidence rates among California teachers: results from the California Teachers Study (United States). Cancer Causes Control 13(7):625–635.
- Brook RD, Franklin B, Cascio W, Hong Y, Howard G, Lipsett M, et al. 2004. Air pollution and cardiovascular disease: a statement for healthcare professionals from the Expert Panel on Population and Prevention Science of the American Heart Association. Circulation 109(21):2655–2671.
- California Department of Public Health 2008. Birth, Death, Fetal Death, Still Birth & Marriage Certificates. Available: http://www.cdph.ca.gov/certlic/birthdeathmar/Pages/default.aspx [accessed 17 December 2008]
- Cavallari JM, Eisen EA, Fang SC, Schwartz J, Hauser R, Herrick RF, Christiani DC. 2008. PM2.5 metal exposures and nocturnal heart rate variability: a panel study of boilermaker construction workers. Environ Health 7:36.; doi: 10.1186/1476-069X-7-36
- Centers for Disease Control and Prevention 2008. National Death Index. Available: http://www.cdc.gov/nchs/ndi.htm [accessed 17 December 2008]
- Chen LH, Knutsen SF, Shavlik D, Beeson WL, Petersen F, Ghamsary M, et al. 2005. The association between fatal coronary heart disease and ambient particulate air pollution: are females at greater risk? Environ Health Perspect 113:1723–1729.
- Cyrys J, Heinrich J, Hoek G, Meliefste K, Lewne M, Gehring U, et al. 2003. Comparison between different traffic-related particle indicators: elemental carbon (EC), PM2.5 mass, and absorbance. J Expo Anal Environ Epidemiol 13(2):134–143.
- Dockery DW, Pope CA III, Xu X, Spengler JD, Ware JH, Fay ME, et al. 1993. An association between air pollution and mortality in six U.S. cities. N Engl J Med 329(24):1753–1759.
- Eftim SE, Samet JM, Janes H, McDermott A, Dominici F. 2008. Fine particulate matter and mortality: a comparison of the six cities and American Cancer Society cohorts with a Medicare cohort. Epidemiology 19(2):209–216.
- Fairley D. 2003. Mortality and air pollution for Santa Clara County, California, 19891996. In: Revised Analyses of Time-Series Studies of Air Pollution and Health. Boston:Health Effects Institute. pp. 97–106.
- Filleul L, Rondeau V, Vandentorren S, Le Moual N, Cantagrel A, Annesi-Maesano I, et al. 2005. Twenty five year mortality and air pollution: results from the French PAARC survey. Occup Environ Med 62(7):453–460.
- Hitchins J, Morawsaka L, Wolff R, Gilbert D.. 2000. Concentrations of submicrometre particles from vehicle emissions near a major road. Atmos Environ 34:51–59.
- Hoek G. 2003. Daily mortality and air pollution in the Netherlands. In: Revised Analyses of Time Series Studies of Air Pollution and Health. Boston:Health Effects Institute. pp. 133–142.
- Hoek G, Brunekreef B, Goldbohm S, Fischer P, van den Brandt PA. 2002. Association between mortality and indicators of traffic-related air pollution in the Netherlands: a cohort study. Lancet 360(9341):1203–1209.
- Hoffmann B, Moebus S, Stang A, Beck EM, Dragano N, Mohlenkamp S, et al. 2006. Residence close to high traffic and prevalence of coronary heart disease. Eur Heart J 27(22):2696–2702.
- Jerrett M, Burnett RT, Ma R, Pope CA III, Krewski D, Newbold KB, et al. 2005. Spatial analysis of air pollution and mortality in Los Angeles. Epidemiology 16(6):727–736.
- Krewski D, Burnett RT, Goldberg M, Hoover K, Siemiatycki J, Jerrett M, et al. 2000. Reanalysis of the Harvard Six Cities Study and the American Cancer Society Study of Particulate Air Pollution and Mortality: A Special Report of the Institute’s Particle Epidemiology Reanalysis Project. Boston, MA:Health Effects Institute.
- Laden F, Neas LM, Dockery DW, Schwartz J. 2000. Association of fine particulate matter from different sources with daily mortality in six U.S. cities. Environ Health Perspect 108:941–947.
- Laden F, Schwartz J, Speizer FE, Dockery DW. 2006. Reduction in fine particulate air pollution and mortality: extended follow-up of the Harvard Six Cities study. Am J Respir Crit Care Med 173(6):667–672.
- Lipsett M, Hurley S, Ostro B.. 1997. Air pollution and emergency room visits for asthma in Santa Clara County, California. Environ Health Perspect 105:216–222.
- Lipsett MJ, Tsai FC, Roger L, Woo M, Ostro BD. 2006. Coarse particles and heart rate variability among older adults with coronary artery disease in the Coachella Valley, California. Environ Health Perspect 114:1215–1220.
- Mar TF, Norris GA, Koenig JQ, Larson TV. 2000. Associations between air pollution and mortality in Phoenix, 19951997. Environ Health Perspect 108:347–353.
- Metzger KB, Tolbert PE, Klein M, Peel JL, Flanders WD, Todd K, et al. 2004. Ambient air pollution and cardiovascular emergency department visits. Epidemiology 15(1):46–56.
- Miller KA, Siscovick DS, Sheppard L, Shepherd K, Sullivan JH, Anderson GL, et al. 2007. Long-term exposure to air pollution and incidence of cardiovascular events in women. N Engl J Med 356(5):447–458.
- Naeher LP, Brauer M, Lipsett M, Zelikoff JT, Simpson CD, Koenig JQ, et al. 2007. Woodsmoke health effects: a review. Inhal Toxicol 19(1):67–106.
- National Research Council 2004. Research Priorities for Airborne Particulate Matter. IV. Continuing Research Progress. Washington, DC:National Academies Press.
- Ostro B, Feng WY, Broadwin R, Green S, Lipsett M. 2007. The effects of components of fine particulate air pollution on mortality in California: results from CALFINE. Environ Health Perspect 115:13–19.
- Ostro B, Lipsett M, Reynolds P, Goldberg D, Hertz A, Garcia C, et al. 2010. Long-term exposure to constituents of fine particulate air pollution and mortality: results from the California Teachers Study. Environ Health Perspect 118:363–369.
- Özkaynak H, Thurston GD. 1987. Associations between 1980 U.S. mortality rates and alternative measures of airborne particle concentration. Risk Anal 7(4):449–461.
- Pope CA III, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, et al. 2002. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. JAMA 287(9):1132–1141.
- Pope CA III, Burnett RT, Thurston GD, Thun MJ, Calle EE, Krewski D, et al. 2004. Cardiovascular mortality and long-term exposure to particulate air pollution: epidemiological evidence of general pathophysiological pathways of disease. Circulation 109(1):71–77.
- Pope CA III, Thun MJ, Namboodiri MM, Dockery DW, Evans JS, Speizer FE, et al. 1995. Particulate air pollution as a predictor of mortality in a prospective study of U.S. adults. Am J Respir Crit Care Med 151(3 pt 1):669–674.
- Puett RC, Hart JE, Yanosky JD, Paciorek C, Schwartz J, Suh H, et al. 2009. Chronic fine and coarse particulate exposure, mortality and coronary heart disease in the Nurses’ Health Study. Environ Health Perspect 117:1697–1701.
- Rogge WF, Hildemann LM, Mazurek MA, Cass GR. 1993. Sources of fine organic aerosol. 2. Noncatalyst and catalyst-equipped automobiles and heavy-duty diesel trucks. Environ Sci Technol 27(4):636–651.
- U.S. Census Bureau 2002. Census of Population and Housing, 2000: Summary File 3—California. Washington, DC:U.S. Census Bureau.
- U.S. Department of Commerce 2008. Social Security Administration’s Death Master File. Available: http://www.ntis.gov/products/ssa-dmf.aspx [accessed ]
- Wellenius GA, Coull BA, Godleski JJ, Koutrakis P, Okabe K, Savage ST, et al. 2003. Inhalation of concentrated ambient air particles exacerbates myocardial ischemia in conscious dogs. Environ Health Perspect 111:402–408.
- World Health Organization 1993. International Classification of Diseases, 10th Revision. Geneva:World Health Organization.
- Zhu Y, Hinds WC, Kim S, Shen S, Sioutas C. 2002. Study of ultrafine particles near a major highway with heavy-duty diesel traffic. Atmos Environ 36:4323–4335.
Come visit with EHP’s Science Editor, Jane Schroeder, at the 2016 ISEE Meeting in Rome, Italy, from 1–4 September 2016. This is a great opportunity to meet a member of our team, and to learn more about the journal. A number of our Associate Editors will be in attendance as well. Jane is also co-hosting a half-day writing and publishing workshop on 4 September, following the conference.
EHP is pleased to present the abstracts from the 28th annual meeting of the International Society for Environmental Epidemiology (ISEE), held in Rome, Italy, 1–4 September 2016, and hosted by the Department of Epidemiology Lazio Regional Health Service, ASL Roma 1, and the Italian Epidemiological Association. The focus of this year’s conference is current and future challenges in exposure assessment, study design, and data analyses.
Featured Children’s Health
Birgit Claus Henn, Adrienne S. Ettinger, Marianne R. Hopkins, Rebecca Jim, Chitra Amarasiriwardena, David C. Christiani, Brent A. Coull, David C. Bellinger, and Robert O. Wright
Diane Gilbert-Diamond, Jennifer A. Emond, Emily R. Baker, Susan A. Korrick, and Margaret R. Karagas
Benjamin B. Green, Margaret R. Karagas, Tracy Punshon, Brian P. Jackson, David J. Robbins, E. Andres Houseman, and Carmen J. Marsit