Spatial and Temporal Variation in PM_2.5 Chemical Composition in the United States for Health Effects Studies

Abstract

Numerous studies have shown that airborne particulate matter (PM) is associated with adverse health effects, including increased risk of premature mortality, hospital admissions, and higher rates of adverse respiratory health indicators in children [Pope and Dockery 2006; U.S. Environmental Protection Agency (EPA) 2006]. Although the health effects of airborne particles have been investigated vigorously for decades, uncertainty persists concerning those characteristics of PM that determine toxicity. To date, studies on the health impacts of PM exposure have used a variety of metrics for PM, including total suspended particles (TSP), coefficient of haze (COH), black smoke, British smoke, KM (a measure of particulate optical reflectance), and PM₁₀ and PM_2.5 (PM with an aerodynamic diameter of < 10 μm and < 2.5 μm, respectively). These indicators reflect PM mass in particular size ranges but not composition specifically. For effective control of particle pollution, information is needed on which sources contribute to the PM characteristics associated with health risk

A growing number of studies have investigated the health effects of PM_2.5, a PM indicator incorporated in the 1997 National Ambient Air Quality Standards (NAAQS) (e.g., Dominici et al. 2006; Franklin et al. 2007; Laden et al. 2006; Schwartz et al. 2002). This indicator was selected in the 1997 NAAQS because of well-established knowledge of the dosimetry of particles in this size range within the respiratory tract and epidemiologic evidence indicating adverse effects of PM_2.5 specifically. However, lacking evidence on the characteristics of PM in this size range that may determine toxicity, a general mass-based standard was promulgated. Characteristics of PM_2.5 that may be relevant to toxicity include metals, organic compounds adsorbed onto particles or forming particles themselves, biologic components, sulfate (SO₄²⁻), nitrate (NO₃⁻), acidity, and surface-adsorbed reactive gases such as ozone (O₃) [Health Effects Institute 2002; National Research Council (NRC) 2004]. Studies have associated several chemical components of PM_2.5 with mortality including iron (Fe), nickel (Ni), zinc (Zn) (Burnett et al. 2000), ammonium nitrate (Fairley 1999), elemental carbon (EC), organic carbon (OC), nitrates (Ostro et al. 2007), and sulfates (Burnett et al. 2000; Ostro 1995).

Recognizing the need for further research on PM characteristics and health, the U.S. EPA has established a national monitoring network for PM_2.5 that provides data on the chemical composition of PM (U.S. EPA 2006). As the data accumulate, they will foster epidemiologic studies designed to assess health risks associated with spatial and temporal variation in PM characteristics. In this article we report analyses of a database constructed from the U.S. EPA monitoring results for 52 PM_2.5 components in 187 counties in the continental United States for the period 2000–2005. We describe the spatial and temporal patterns of variation of PM_2.5 chemical components and identify components that might be evaluated in studies of PM_2.5 and human health effects.

For an individual chemical constituent to be a mediator of the risk associated with PM_2.5 total mass, the concentration of that component must co-vary with the more general mass variable used in epidemiologic research (i.e., PM_2.5 total mass); however, we recognize that multiple components may contribute to the risk and that components may interact. Other components that may be harmful to human health may not be related to the observed relationships between PM_2.5 and health. We provide descriptive analyses intended to identify candidate PM_2.5 components that meet this criterion of being correlated with PM_2.5 total mass and to summarize the spatial and temporal variation of such components. These candidate components should be explored further to determine whether they mediate the effect of PM_2.5 total mass and to investigate the underlying biological mechanism.

Methods

Database development

We developed a database of concentrations for 52 PM_2.5 components and PM_2.5 total mass for 187 continental U.S. counties for the period February 2000 to December 2005, based on data obtained from the U.S. EPA’s Office of Air Quality Planning and Standards (U.S. EPA 2006). Counties and components are listed in the Tables S1 and S2 of the Supplemental Material (http://www.ehponline.org/docs/2007/9621/suppl.pdf). Not all counties had the full complement of data for the entire time period. Although most monitors provide data every 6 days, the average frequency of measurement by monitor ranged from 3.1 to 11.9 days. We generated countywide estimates for each PM_2.5 component based on an analysis of the monitor or monitors within each county.

In developing and analyzing the data set, we needed to address several key issues, described below, for which we developed a protocol to combine the data to generate a countywide average.

Suspect data

The U.S. EPA coded some observations as problematic or unusual (e.g., “lab issues”). These individual observations, which included many extreme values, were omitted.

Noncontinental counties

We omitted nonmainland counties—that is, those in Hawaii and Alaska.

Co-location of monitors

Fifteen of 259 U.S. EPA monitoring sites (5.8%) had multiple monitors for duplicate sampling on the same day. Data from multiple monitors in the same site were treated as repeat measurements at the same site and were averaged to generate an overall observation at that location. After values from co-located monitors were averaged, county-level exposures were estimated as the average across monitors within the county.

Counties with little data

We omitted counties with data collected on PM_2.5 total mass or for any of the individual PM_2.5 components for only a brief period (< 6 months or < 30 observations). A total of 11.4% of the individual counties (2.7% of observation days) were omitted for this reason.

Check of unusual values

Observation days were omitted if the highest PM_2.5 value recorded was over three times higher than the second highest value. This criterion excluded only two observation days.

OC measurements require adjustment to correct for field blanks and to account for elements such as oxygen and hydrogen that are associated with OC to estimate organic matter. Organic carbon matter (OCM) was calculated as

where OCM = organic carbon matter; k = adjustment factor to account for noncarbon organic matter (1.4); OC_m = measured organic carbon; OC_b = organic carbon for blank filters. Blank filter correction values were based on U.S. EPA data (U.S. EPA 2006). Recent analysis has shown that OC blank values from 2001 to 2005 have increased for some samplers and decreased for others (Frank NH, unpublished data). We performed a sensitivity analysis using a data set with OCM estimated by sampler- and year-specific blank correction values (Frank NH, unpublished data). This alternative value for OCM is specified as OCM2.

Analysis

First we determined which PM_2.5 components contributed a substantial fraction to total PM_2.5 mass, for either the yearly average or any seasonal average. Seasons were defined based on 3-month periods (e.g., summer was defined as June–August). Second, we identified the components that co-vary day to day with total PM_2.5 mass. We calculated yearly and seasonal correlations between PM_2.5 total mass and each component and with the corresponding seasonally detrended time series,

, defined as:

where X^c_k,t = the concentration of component k at time t for county c,

is the 91-day moving average of the concentration of component k for county c centered at time t, and n^c_k,t = the number of days with observations for component k for county c for a 91-day moving average centered at time t.

Analysis of detrended data included only counties with more than one full year of observations.

Results

The original data set included 62,690 observation days across all sites (i.e., monitor-days of data), which dropped to 48,591 observations for 187 counties with the exclusions described earlier. Many counties had data for only a portion of the study period. The average number of observations per county for PM_2.5 total mass was 260 days (range, 41–676). Figure S1 in the Supplemental Material (http://www.ehponline.org/docs/2007/9621/suppl.pdf) shows the number of PM_2.5 observations per county for the study period. For other components, the average number of observations per county ranged from 248 days [sodium ion (Na⁺)] to 260 days (OC). Every county had data available for every season, with 26.3% of the PM_2.5 total mass data from summer, 26.5% from autumn, 22.7% from winter, and 24.6% from spring. Figure 1 shows average PM_2.5 levels for the study period (2000–2005) for each county. Overall, PM_2.5 levels were higher in the eastern United States and California, and lowest in the central regions and Northwest. However, PM_2.5 concentrations had strong seasonal patterns that differed by region (Figure 2). On the west coast, levels peaked in winter and autumn, especially for northern California, whereas on the east coast higher levels were recorded for summer. Concentrations remained lower in the central United States throughout the year.

Open in a separate window

Figure 1

PM_2.5 average (μg/m³) for 187 U.S. counties, 2000–2005.

Open in a separate window

Figure 2

Seasonal PM_2.5 averages (μg/m³) for 187 U.S. counties, 2000–2005.

Table 1 provides summary statistics for each PM_2.5 component for the full year, and for summer and winter. Many components show strong seasonal patterns. For example, NO₃⁻, chlorine (Cl), Zn, Ni, and bromine (Br) are 3.6, 3.2, 1.5, 1.4, and 1.4 times higher in winter than in summer, respectively. Aluminum (Al), titanium (Ti), magnesium (Mg), silicon (Si), and SO₄²⁻ were 1.5 to > 2 times higher in summer than in winter. Other components did not show distinct seasonal patterns. The results in Table 1 and other summary measures that provide a national average conceal spatial heterogeneity on smaller spatial domains. Similarly, presentations of yearly values obscure seasonal differences.

If a PM_2.5 component contributes to the associations with risks for health outcomes observed in time-series studies based on PM_2.5 total mass, the component would be expected to exhibit strong day-to-day variation with PM_2.5 total mass. Many such components are likely to contribute substantially to PM_2.5 total mass. We first identified components that comprise the majority of overall PM_2.5 mass. Only 7 of the 52 components contributed ≥ 1% to the PM_2.5 total mass for the yearly average or any of the seasonal averages across all 187 counties. Those components (NH₄⁺, EC, OCM, NO₃⁻, Si, Na⁺, and SO₄²⁻) comprised 79–85% of the total PM_2.5 mass for the yearly or seasonal averages. Figure 3 shows the percentages of PM_2.5 contributed by these components for yearly, winter, and summer averages, for nationwide, eastern U.S., and western U.S. averages. SO₄²⁻ is a larger contributor in summer, whereas NO₃⁻ is a larger contributor in winter for both regions. Although Figure 3 presents results for the two U.S. regions, spatial heterogeneity can also exist within regions. Further, this analysis is limited to the components included in the database, and other components or chemical forms (e.g., ferric oxide) that were not measured could also have contributed ≥ 1% to total PM_2.5 mass.

Open in a separate window

Figure 3

Percent of PM_2.5 composition by component for yearly, winter, and summer averages, by region.

We also examined if any components contributed 1% or more to PM_2.5 within any individual county for either a yearly or seasonal average. Components meeting this criteria were Al, calcium (Ca), Cl, Fe, and potassium (K), which on average provide 0.18–0.62% of PM_2.5 total mass across the whole year, but in some cities contributed up to 5.4% for a given season. The contribution of these components to PM_2.5 total mass on average across all communities and the minimum and maximum values for any single community are provided in the Supplemental Material (Table S3; http://www.ehponline.org/docs/2007/9621/suppl.pdf) for yearly and seasonal averages.

Figures 4–7 map yearly and seasonal averages for the SO₄²⁻ and NO₃⁻ components. SO₄²⁻ PM_2.5 displays a strong east/west pattern (Figure 4). In the eastern United States, the SO₄²⁻ component of PM_2.5 typically peaks during summer (Figure 5). The NO₃⁻component of PM_2.5 shows a somewhat inverse pattern, with higher concentrations on the west coast, primarily in California (Figure 6). NO₃⁻ PM_2.5 also exhibits a north/south pattern, with higher levels in parts of the Northeast and decreasing levels towards the Southeast. This north/south gradient remains throughout all seasons (Figure 7), and highest concentrations in the eastern United States occur in winter. The western United States has the highest nitrate PM_2.5 concentrations during winter and autumn.

Open in a separate window

Figure 4

Sulfate PM_2.5 (μg/m³) averages for 187 U.S. counties, 2000–2005.

Open in a separate window

Figure 5

Seasonal sulfate PM_2.5 (μg/m³) averages for 187 U.S. counties, 2000–2005.

Open in a separate window

Figure 6

Nitrate PM_2.5 averages (μg/m³) for 187 U.S. counties, 2000–2005.

Open in a separate window

Figure 7

Seasonal nitrate PM_2.5 averages (μg/m³) for 187 U.S. counties, 2000–2005.

The Supplemental Material (Figures S2–S11; http://www.ehponline.org/docs/2007/9621/suppl.pdf) provides maps of yearly and seasonal averages for other key components: Na⁺, Si, EC, NH₄⁺, and OCM. Sodium ion PM_2.5 levels are higher in coastal regions, relating to sea salt (Figure S2), and do not exhibit a strong seasonal pattern outside of a moderate trend in the western United States (Figure S3). The highest overall Si levels were noted in Texas (Figure S4), with high concentrations in other areas by season (Figure S5). EC PM_2.5 showed spatial and temporal patterns similar to those of NO₃⁻, without the north/south gradient in the eastern United States (Figures S6 and S7). Both EC and NO₃⁻ PM_2.5 were higher in California, and peaked in winter and autumn. Because ammonium (NH₄⁺) is commonly observed in the forms of ammonium nitrate or ammonium sulfate, the ammonium component of PM_2.5 is correlated with SO₄²⁻ and NO₃⁻ components and consequently exhibits a mix of those components’ spatial and temporal patterns for yearly averages (Figure S8) and seasonal averages (Figure S9). OCM is higher on the west coast (Figures S10 and S11).

Table 2 provides correlations among yearly and seasonal averages for these seven components and all other components for which the correlation coefficient reaches ≥ 0.5. Additional correlations among PM_2.5 components are provided in the Supplemental Material (Table S4; http://www.ehponline.org/docs/2007/9621/suppl.pdf). These tables were created by first calculating the yearly and seasonal averages in each county for each component, and then calculating the correlation coefficients. NH₄⁺ was most strongly correlated with SO₄²⁻ and NO₃⁻, with a stronger relationship with SO₄²⁻ in summer (0.88) and NO₃⁻ in winter (0.86). EC concentrations covary with Fe concentrations in all seasons; with OCM in winter, summer, and autumn; and Ti in winter and autumn. The concentrations of OCM are correlated with the levels of K in winter and Ti in autumn. Si concentrations are associated with those of crustal elements including Ca. Na⁺ concentrations are most closely associated with levels of Cl, but less so in winter. The strongest correlations for SO₄²⁻ are with NH₄⁺.

Table 3 provides data on the correlations between day-to-day variations of the key component concentrations and of the total PM_2.5 mass for nationwide, eastern U.S., and western U.S. regions. Additional correlations are provided in the Supplemental Material in Table S5 (http://www.ehponline.org/docs/2007/9621/suppl.pdf). Components with the greatest contributions to total PM_2.5 mass also had the strongest temporal correlations with PM_2.5 total mass. The components typically co-varied with PM_2.5 total mass when they reached peak concentrations, such as summer for SO₄²⁻ and winter for NO₃⁻. Of the components not listed in Table 3, K was correlated with total PM_2.5 in winter (0.52) and Br in autumn (0.63), spring (0.55), and for yearly averages (0.56).

Using the seasonally detrended data for PM_2.5 mass and each component, on average across the 180 counties with ≥ 1 year of data, the following components were found to have strong day-to-day variation with PM_2.5 total mass: NH₄+ (average correlation 0.84); SO₄²⁻(0.78); OCM (0.68); NO₃⁻ (0.51); Br (0.51); and EC (0.51). The relationship between daily PM_2.5 and component concentrations varied by county. For these six components (NH₄⁺, SO₄²⁻, OCM, NO₃⁻, Br, and EC), the percentages of counties with correlation coefficients > 0.6 were 95, 90, 81, 34, 33, and 22%, respectively.

We applied an alternative method of adjustment to calculate OCM, discussed previously, using blank filter values specific to the year and type of sampler (Frank NH, unpublished data). The alternative measure, designated OCM2, provided comparable results to our original OCM measure. The correlation coefficient between OCM and OCM2 was 0.99 on average across all counties (range, 0.97–1.00). Both measures of organic carbon matter (OCM and OCM2) had similar values for yearly and seasonal concentrations (Supplemental Material, Table S6; http://www.ehponline.org/docs/2007/9621/suppl.pdf), the percentage of PM_2.5 total mass comprised of OCM by year or season (Table S7), and the correlation between OCM and PM_2.5 total mass, NH₄⁺, EC, NO₃⁻, Si, Na⁺, or SO₄²⁻, by year or season (Table S8).

Discussion

The PM_2.5 mixture varies strongly by region and by season, and the degree of spatial and temporal variability differs by component, which has implications for epidemiologic research on PM_2.5 characteristics. National studies have already demonstrated that the estimated short-term effects of PM₁₀ on mortality (Dominici et al. 2003; Peng et al. 2005) and of PM_2.5 on hospital admissions (Dominici et al. 2006) vary by season and by region, with the highest effect estimates for mortality and hospital admissions in the northeastern United States during summer. These regional and temporal differences may reflect variation in the PM_2.5 mixture and its sources.

These findings indicate the complexity of interpreting regional differences in the effect of PM_2.5 and of designing studies directed at characterizing effects of particular components. Because of variations in the PM_2.5 mixture, the risk associated with a particular component is assessed against a continually varying background of other pollutants. Our descriptive analyses of the new data on PM_2.5 components illustrate the challenge of testing hypotheses around specific components such as explaining observed seasonal and regional variation in the effect of PM_2.5 (Peng et al. 2005).

Many techniques are available to determine the sources of PM_2.5 components, including factor analysis, Gaussian plume modeling, and backward trajectory modeling; each has its own set of advantages and limitations (Hopke et al. 2006; Ito et al. 2006; Laden et al. 2000; Lapina and Paterson 2004; Mar et al. 2006; Paatero et al. 2003; Thurston et al. 2005). Methodologies that assign specific components or sets of components to sources face the challenge that any individual PM_2.5 component comes from a variety of sources. Table S9 in the Supplemental Material (http://www.ehponline.org/docs/2007/9621/suppl.pdf) lists some sources of the seven key components identified. Approaches to linking components to specific sources may be more suitable for localized studies in which dominant sources can be identified, for example in cases where detailed knowledge is available regarding sources for the region.

For studies based on national data or large regions, multiple sources of each component may complicate such efforts (Table S9). For example, Selenium (Se) and SO₄²⁻ had an overall correlation of 0.47 in this data set, and both can result from combustion of coal, oil, or biomass. However, Se emissions also come from smelters and coke production, and SO₄²⁻ emissions result from motor vehicles, incineration, electronics manufacturing, steel mills, and other sources. EC and Fe both come from traffic emissions, vegetative burning, oil combustion, and casting processes. Na⁺ and Cl co-vary because of a common origin in oceans; Na⁺, however, also comes from other sources, and Cl can result from combustion emissions from cooking, coal, automobiles, vegetation burning, and incinerators.

These findings suggest that the new data on PM_2.5 components may not lead to satisfactory, definitive source apportionment for national studies. In localized settings, source apportionment and related methodologies are more appropriate. Correlated concentrations and multiple sources complicate the identification of individual effects of various PM_2.5 components on a national scale. For example, SO₄²⁻ concentrations are associated with NH₄⁺ and Se concentrations. Therefore, a study identifying PM_2.5 sulfate as associated with adverse health impacts may be detecting effects of co-varying pollutants (e.g., Se, ammonium sulfate, or other components with similar sources for that region). Because every component in the data set has multiple and shared sources, no pairs of the components are perfectly correlated; the highest correlation of yearly averages (0.998) was for cerium and lanthanum. Therefore, methods using a single component or set of components as source surrogates [e.g., SO₄²⁻ and Se for coal combustion, vanadium (V) for oil combustion, EC for traf-fic] for national studies can be affected by some misclassification of contributing sources that might also vary by region and season. In local studies with less spatial heterogeneity of source profiles, methods such as source apportionment are more likely to be successful. Further, in studies of smaller regions, additional PM_2.5 chemical component data may be available, including concentrations of ammonium sulfate, rather than NH₄⁺ and SO₄²⁻ separately.

We found that of the 52 components considered, only seven contributed ≥ 1% to total PM_2.5 mass for the yearly average or any seasonal average. We also found that several of these seven components are also correlated with day-to-day variations in the PM_2.5 total mass. Results indicate that the strongest correlations with PM_2.5 total mass are NH₄+ (yearly, all seasons), OCM (especially winter), NO₃⁻(winter), and SO₄²⁻ (yearly, spring, autumn, and summer), with particularly strong correlations for NH₄+ or SO₄²⁻ in summer.

The observed health risks of PM_2.5 could be a function of the key components identified above; however, other explanations are also possible. These alternative scenarios include a component contributing < 1% to total PM_2.5 mass, but with concentrations below detection limits that co-vary with PM_2.5 total mass, or a component or set of components that co-vary with the identified key components. To gather evidence toward these alternative explanations, we evaluated which components co-varied with the seven key components as shown in Table 2, andwhich components co-varied with PM_2.5 detrended data, which identified six of the key components.

Limitations of these data and analyses include measurement error and detection limits, which may hinder identification of relationships among components or a component’s contribution to PM_2.5 total mass (Flanagan et al. 2006; Frank 2006; Schwab et al. 2006). These limitations may affect some chemical components more than others because of differing instrument abilities for detection and measurement. For example, the ratio of OCM to OC_m−OC_b [i.e., k (Equation 1)] can vary by site and season (Bae et al. 2006, Turpin and Lim 2001) although such specific adjustments to OC_m are currently not possible. In particular, k may be higher in rural settings than in urban settings. Even the levels of PM_2.5 total mass are subject to measurement error. Because not all possible PM_2.5 components were measured, the sum of measured PM_2.5 component concentrations was generally but not universally less than the total PM_2.5 mass. However, the sum of components can exceed PM_2.5 total mass because of negative artifacts such as loss of ammonium nitrate and other semivolatile organics (Frank 2006).

Because of these limitations, health risks could be associated with the true concentrations of a component or set of components that co-varies with PM_2.5 total mass, even if measured concentrations in this data set do not co-vary with PM_2.5 total mass because of measurement error. Further, we did not investigate the possibility that observed PM_2.5 health effects could result from a set of components with a collective concentration that co-varies with PM_2.5 total mass, although individual component concentrations do not.

Footnotes

References