The associations between air pollution, especially particulate matter (PM), and adverse human health effects have been well documented (
1-10). PM is associated with decreased respiratory function, aggravation of existing respiratory and cardiovascular conditions, altered defense mechanisms, and even premature death. The most susceptible populations include those with preexisting respiratory or cardiovascular conditions, asthmatics, children, and the elderly (
8,11).
To date, few epidemiology studies have used PM measures other than size-segregated mass as the exposure metric. Schwartz et al. (12) looked at episodes of high coarse particle concentration in Spokane, Washington, and found that windblown dust episodes were not associated with increased mortality. In the Harvard Six Cities Study, Schwartz et al. (3) found a significant association between nonaccidental mortality and particulate matter
2.5 µm in aerodynamic diameter (PM2.5) and sulfur. They did not find a significant association with particulate matter
10 µm in aerodynamic diameter (PM10) or the coarse fraction of PM [PMCF (PM10 minus PM2.5)]. In contrast, Ostro et al. (13) found that PM10 dominated by coarse particles was associated with an increase in mortality in the Coachella Valley in California. The differences in the results from these two studies may be due to the particulate composition as well as the difference in the amount of PMCF. In the eastern United States, PM2.5 is dominated by sulfates (34%), whereas in the western and central United States it is dominated by organic carbon (OC) from motor vehicles and vegetative burning (39%) (14). The average PM2.5/PM10 ratio for the Six Cities Study (3) was 0.6 (based on the 50th percentiles) as compared to a ratio of 0.3 for Phoenix, Arizona (15).
The goal of the present study was to evaluate the associations between daily air pollution and total nonaccidental and cardiovascular mortality in Phoenix. Phoenix is an arid southwestern city with a population of approximately 1 million residents (16). It is an interesting location because of its large proportion of elderly people (9.7% of the population is > 65 years of age). The elderly are more susceptible to air pollution than the general public (2). The primary sources of PM in Phoenix are motor vehicles, paved road dust, and vegetative burning (15).
This study focused on the effects of air pollution on cardiovascular mortality for several reasons. First, the association between air pollutants and cardiovascular mortality has been consistent in previous studies (1,2,9,17). Second, a study in Baltimore, Maryland, found that heart rate variability was associated with PM2.5 in elderly subjects with cardiovascular conditions (18). Finally, in this study cardiovascular mortality had the largest sample size, accounting for 45% of the total nonaccidental deaths in the study region (based on zip codes). This may be reflective of the increased size of Phoenix's elderly population, which is more prone to cardiovascular disease.
A unique aspect of this study is that our pollution data include daily information not only on traditional gaseous pollutants, but also on PM in various size fractions and the chemical composition of PM2.5. From 1995 to 1997, the U.S. Environmental Protection Agency (EPA) National Exposure Research Laboratory (NERL) operated a comprehensive monitoring platform in Phoenix. They collected daily PM2.5 samples and subsequently analyzed them for various chemical components of PM. This provided an opportunity to examine more specific metrics for PM than simply mass, as well as an opportunity to identify selected chemical components of PM that are associated with mortality.
In addition to PM, this study also evaluated the association between total nonaccidental and cardiovascular mortality and other measured air pollutants: carbon monoxide, nitrogen dioxide, sulfur dioxide, and ozone. These EPA criteria pollutants are also associated with mortality (7,17,19,20).
Study area and data. Mortality data for all of Maricopa County from 1995 to 1997 were obtained from the Arizona Center for Health Statistics in Phoenix. Death certificate data included residence zip code and the primary cause of death as identified by the International Classification of Diseases, Ninth Revision (ICD-9, World Health Organization, Geneva). Only the deaths of residents in the zip codes located near the air pollution platform were included in this study. This zip-code region was recommended by the Arizona Department of Environmental Quality (Phoenix, AZ). We evaluated total nonaccidental mortality (ICD-9 codes < 800) and cardiovascular mortality (ICD-9 codes 390-448.9) in this study. Summary statistics for the mortality outcomes are presented in Table 1.
We obtained PM2.5, PM10, PMCF, and PM2.5 chemical composition data from the EPA NERL platform in central Phoenix. Chemical composition was only available for PM2.5. The monitoring platform is approx-imately 10 km west-northwest of downtown Phoenix at a state and local air monitoring station. Standard meteorologic parameters such as wind speed and direction, temperature, and relative humidity were continuously measured. The average temperature in Phoenix from 1995 to 1997 was 23.7 ± 8.1°C. The average relative humidity was 32 ± 15%.
NERL investigators made hourly PM2.5 and PM10 measurements each day using two colocated tapered element oscillation microbalance (TEOM) monitors (Rupprecht & Patasnick Co., Albany, NY). The TEOM-PM10 was fitted with an EPA-approved federal reference method PM10 impactor inlet (model 246b; Andersen Instruments, Smyrna, GA). The TEOM-PM2.5 was fitted with a PM2.5 cyclone inlet (University Research Glassware, Chapel Hill, NC). The PM2.5 cyclone on the TEOM was replaced with a well-impactor ninety-six (WINS) inlet on 20 December 1996. The WINS inlet has a sharper cut point as compared to the cyclone. We averaged the hourly concentrations to create a 24-hr average (0700-0700 hr), and we calculated the concentration of coarse fraction (TEOM PMCF) as TEOM PM10 minus TEOM PM2.5.
NERL investigators collected the daily gravimetric integrated 24-hr (starting at 0700) fine particle filter samples using a dual fine particle sequential sampler (DFPSS; University Research Glassware). The DFPSS was fitted with a cyclone that was identical to the cyclone on the TEOM-PM2.5. The DFPSS collected daily samples on both Teflon and quartz filters. The Teflon filter was used for mass and elemental analysis, whereas the quartz filter was used for carbon analysis. In addition to the DFPSS, NERL investigators operated a dichotomous sampler (Andersen Instruments, Inc.) every third day beginning 17 June 1996. Both the PM2.5 and PMCF samples were collected on Teflon filters. The investigators measured elemental concentrations at the EPA (Research Triangle Park, NC) with energy dispersive X-ray fluorescence. OC and elemental carbon (EC) were measured by Sunset Laboratory (Forest Grove, OR) using thermal optical transmittance (21).
PM and gaseous pollutant concentrations (range and mean ± SD) from 1995 to 1997 are presented in Table 2. We obtained gaseous criteria pollutant data for CO, NO2, O3, and SO2 from the EPA Aerometric Information Retrieval System (AIRS) database (22) for residential sites in the Phoenix region. We averaged CO values over four monitoring sites and we averaged NO2 over two sites. Only one residential monitoring site was available for SO2. We averaged the hourly averages for CO, NO2, and SO2 over 24 hr from 0700 to 0700. We used the maximum hourly O3 (O3 max) concentration in the same 24-hr period in the analysis.
The PM2.5 constituents that we evaluated for effects on mortality were sulfur, zinc, lead, soil-corrected potassium (KS) (23), OC, EC, total carbon (TC), and reconstructed soil. Soil was reconstructed by summing the oxides of Al, Si, Ca, Fe, and Ti using the formula recommended by Malm et al. (24). We also considered PM2.5 that was corrected for soil content (nonsoil PM2.5 = PM2.5 - reconstructed soil). Table 3 presents the percent of the total mass of PM2.5 accounted for by each component. The elements aluminum, silicon, calcium, titanium, and iron were not evaluated separately in the mortality analysis because they are the major elemental components of soil.
Statistical analysis. In our zip-code regions, we analyzed a total of 9,276 nonaccidental deaths from 1995 to 1997. Poisson regression was used to evaluate the association between the air pollutant exposure variables and the mortality outcomes (2,5).
We used Poisson regression because mortality data are discrete counts and death is a rare event. Poisson regression assumes the variance is equal to the mean. When the variance exceeds the mean, the variance is overdispersed. We adjusted standard errors for overdispersion; however, the amount of overdispersion was small. The overdispersion parameter was 1.05 and 1.00 for nonaccidental and cardiovascular mortality, respectively. We calculated all relative risks (RRs) for an interquartile increase (25th-75th percentile) in pollutant concentration.
The effect of air pollution on mortality is small and can be influenced by confounders. Therefore, base models for total mortality and cardiovascular mortality were constructed by adjusting for day of the week with indicator variables, and time trends, temperature, and relative humidity with smoothing functions (25). We determined degrees of freedom (df) for the function used to smooth time trend by minimizing autocorrelation as well as the Akaike information criterion (AIC) (26). We chose the df and lag for the smoothing functions for temperature and relative humidity to minimize the AIC. The base model for total mortality used indicator variables for day of the week, 10 df for time trends, 2 df for temperature with 2 days lag, and 2 df for relative humidity with 0 days lag. The base model for cardiovascular mortality used indicator variables for the day of the week, 10 df for time trends, 2 df for temperature with 1 day lag, and 2 df for relative humidity with 0 days lag.
We included continuous daily data from 1995 to 1997 (1,097 days) in the study. Each day was coded and included in the model to adjust for time trends. Little autocorrelation was observed after adjusting for day of week, time trends, temperature, and relative humidity. The autocorrelation for days 1-25 for both total and cardiovascular mortality were within the 95% confidence interval for an independent series.
We evaluated air pollution exposure variables by adding them individually as linear terms to the base model. The air pollution exposure metrics that were evaluated in this analysis included CO, NO2, O3, SO2, TEOM PM10, TEOM PM2.5, TEOM PMCF, PM2.5 (DFPSS), S, Zn, Pb, soil, KS, nonsoil PM, OC, EC, and TC. Lag days ranging from 0 to 4 were investigated. We evaluated the assumption of a linear relationship using a smooth function. This assumption was met if a straight line could be placed within the 95% confidence intervals (CIs). A p-value < 0.05 associated with the pollution exposure variable was considered significant. We conducted Poisson regression analyses using S-PLUS 4 (Mathsoft, Inc., Seattle, WA)
Factor analysis. We conducted a factor analysis on the daily concentrations of the chemical components of PM2.5 from samples collected by the DFPSS (Al, Si, S Ca, Fe, Zn, Mn, Pb, Br, KS, OC, and EC). The analysis also included the daily averages of the gaseous species emitted by combustion sources (CO, NO2, and SO2). Factor analysis is a technique used to explain the correlations between variables in terms of underlying factors that are not directly measurable. Each factor is a linear combination of the original variables and all such factors are orthogonal to each other. The factors were extracted using principal component analysis with a varimax rotation. We conducted factor analysis using SAS (SAS Institute Inc, Cary, NC). We used the resultant factor scores as surrogate exposure variables in predicting mortality outcomes with the Poisson regression model. Each factor was evaluated in a single source model. However, because the factor scores formed a set of orthogonal variables, we performed a separate regression analysis with all of the scores included in one multifactor model.
We also conducted a factor analysis on the daily concentrations of the chemical components of PMCF from samples collected by the dicot (Al, Si, Cl, S, K Ca, Mn, Fe, Zn, Br, Pb, Sr, Cu, and Rb). We did not use the scores from this analysis in the time-series analysis because the sampling period started in June 1996 and samples were only collected every third day.
Table 4 shows the correlation coefficients between PM, gaseous pollutants, temperature, and relative humidity for Phoenix in 1995-1997. PM2.5 (obtained from the DFPSS) was highly correlated with CO (r = 0.85) and NO2 (r = 0.79), but less so with SO2 (r =0.43). PM2.5 from the DFPSS was highly correlated with that measured with the TEOM (r = 0.93). Table 5 shows the correlation coefficients between selected chemical composition components of PM2.5 and the other air pollutants. TEOM PM10 was correlated with fine soil (r = 0.72), OC (r = 0.58), EC (r = 0.58), and TC (r = 0.59). TEOM PM2.5 was highly correlated with OC (r = 0.89), EC (r = 0.84), TC (r = 0.90), and to a lesser extent with Zn (r = 0.61), Pb (r = 0.67), and KS (r = 0.59). The high correlation coefficients between carbon and PM2.5 indicate that the majority of the variation in PM2.5 is due to combustion products. PMCF was correlated with soil (r = 0.66).
OC and EC concentrations follow a seasonal pattern--they are high in the colder months and low in the warmer months. This pattern is due to increased combustion emissions from space heating and the decreased mixing height during the winter months. Particulate sulfur concentrations peak in the warmer months. Soil concentration also follows a seasonal trend, with higher concentrations in the spring and fall. Measured soil concentrations decreased after 20 December 1996 because of the use of the WINS inlet.
Summaries of the RR between the exposure variables and both total and cardiovascular mortality are presented in Tables 6 and 7, respectively. Because of space limitations, we only present statistically significant (p < 0.05) and marginally significant (p < 0.10) results in the tables, although models were run using all of the pollutants listed in Table 5. Tables of all of the nonsignificant results are available from the authors by request. We evaluated the associations between total and cardiovascular mortality and the gaseous pollutants, PM mass metrics, and PM composition metrics using single-pollutant models.
IQR values were corrected on 22 March 2000.
We found significant associations between both mortality outcomes and selected gaseous air pollutants. CO and NO2 were positively associated with total mortality at 0- and 1-day lags. There was evidence of a weak association with SO2 at 0 days lag (p < 0.10). We found several strong associations with cardiovascular mortality. Cardiovascular mortality was positively associated with CO (0-4 days lag). This was the most consistent association because the association was significant for all 5 lag days. Statistically significant associations (p < 0.05) were also evident with NO2 on lag days 1 and 4, although the association was weaker on lag days 2 and 3. Cardiovascular mortality was also associated with SO2 (lag days 2, 3, and 4).
We also found significant associations between the mortality outcomes and particulate mass. The associations between PM10 and total mortality, and between PMCF and total mortality, were marginal (p < 0.10). Total mortality was not significantly associated with PM2.5; however, the RR was 1.02 (CI, 1.00-1.05). All PM mass metrics were associated with an excess risk of cardiovascular death. The strongest associations were with PM2.5 (TEOM), followed by PM10 and PMCF. PM2.5 adjusted for soil content (nonsoil PM2.5) was also related with cardiovascular mortality with 1 day lag (p < 0.10). Table 7 lists all of the statistically significant associations with cardiovascular mortality. Cardiovascular mortality showed a more consistent association with particulate mass concentrations than total mortality. We further investigated the associations between the mortality outcomes and PM by evaluating the association between the mortality outcomes and the PM composition. The PM2.5 composition data analysis revealed that EC and TC were significantly associated with cardiovascular mortality (1 day lag). Weaker associations were also evident with OC at 1 and 3 days lag and TC at 3 days lag. KS had a significant positive association with cardiovascular mortality (3 day lag).
We also found that soil, S, and Pb were negatively associated with total mortality. That is, these exposure variables were associated with a decrease in excess deaths.
We further evaluated the associations between the mortality outcomes and sources of both particulate and gas-phase pollutants using the scores from a factor analysis in place of the individual pollutant concentrations. The results from the analysis with five factors are presented in Table 8. Factor 1 probably represents the influence of motor vehicle exhaust and resuspended road dust with high loadings (loading > 0.5) on Mn, Fe, Zn, Pb, OC, EC, CO, and NO2. Factor 2 represents soil with high loadings on Al, Si, and Fe. Factor 3 represents vegetative burning with high loadings on OC and KS. Factor 4 represents a local source of SO2 with a high loading on SO2. Factor 5 represents predominately regional sulfate with a high loading on S. The RRs associated with an interquartile range increase in each factor are presented in Table 9. Total mortality had both a positive and a negative association with the factor representing regional sulfate, positive on lag day 0 (same day) and negative on lag day 3. The factor representing SO2 had a negative association with total mortality. We also found a significant negative association for fine soil on lag days 1 and 2, and a nearly significant negative association on lag days 3 and 4. Cardiovascular disease was significantly associated with the factors representing motor vehicles (lag day 1) and vegetative burning (lag day 3). Regression analysis with all of the factors included in a multisource model produced similar results.
Table 10 presents the results from the factor analysis on the daily concentrations of the chemical components of PMCF from samples collected by the dichotomous sampler. Factor 1 represents soil with high loadings on Al, Si, K, Ca, Mn, Fe, Sr, and Rb. Factor 2 represents a source of coarse fraction metals with high loadings on Zn, Pb, and Cu. Factor 3 represents a marine influence with a high loading on Cl. These three factors explain 91.8% of the variance in the PMCF data.
Sensitivity analysis. As a sensitivity analysis, we analyzed temperature as a cofactor rather than a confounder. That is, we evaluated the effects of temperature on mortality as an independent variable rather than adjusting for it in the model as a confounding variable. We evaluated the significance of temperature after adjusting for day of the week, time trends, and relative humidity. For total and cardiovascular mortality, we found that temperature was not associated with excess deaths. Temperature was not correlated with either PM10 (r = -0.08) or PM2.5 (r = -0.25). A second analysis examined the effect of extreme temperatures. If the average daily temperature was greater than or equal to the 95th percentile (35.4°C), we assigned a 1 to the predictor variable; otherwise we assigned a 0. We did not find an association between extreme temperature and total mortality. However, with cardiovascular mortality, extreme temperature was associated with excess deaths at 0 and 2 days lag (p < 0.1). To further assess the importance of the high temperature days to our analysis, we evaluated the association between PM2.5 and cardiovascular mortality after excluding the days when the temperature was above the 95th percentile. The effect of eliminating the high temperature days was negligible. The RR for cardiovascular mortality associated with PM2.5 (1 day lag) including all days was the same as that excluding the hottest days (RR = 1.06; CI, 1.02-1.10).
We also conducted a sensitivity analysis with relative humidity as a cofactor, with the model controlling for time trends and temperature. As a cofactor, relative humidity was not associated with either total mortality or cardiovascular mortality. To further assess the effects of extreme relative humidity, we eliminated the driest days (relative humidity < 25th percentile) from the data. We then found that the coarse fraction was no longer associated with total mortality. The association between cardiovascular mortality and coarse fraction was statistically significant (p < 0.05) on the concurrent day, but nonsignificant with 1-day lag.
We also used dew point rather than relative humidity in the base model. Controlling for dew point rather than relative humidity did not alter our results. We obtained similar regression coefficients.
To assess the effect of replacing the PM2.5 cyclone on the TEOM with the WINS, we evaluated the association between soil and total mortality from 1 January 1995 to 31 December 1996 and from 1 January 1997 to 31 December 1997. The latter period represented the WINS inlet measurements. The association between soil and mortality was not significant for the cyclone measurements alone. Analysis with only the WINS data revealed that the association between soil and mortality was positive and significant at 0 days lag, but not significant for any of the other days.
We estimated soil-related potassium using a correction ratio = K/Si (23). We then reevaluated the RR for cardiovascular mortality associated with KS using KS calculated from three slightly different values of K/Si. This correction ratio is dependent on where the soil was obtained: PM2.5 paved road dust (K/Si = 1.85/13.69), an agricultural field (K/Si = 1.98/14.35), or Phoenix desert soil (K/Si = 1.89/14.00) (27). We found similar RRs for cardiovascular mortality associated with KS when we used any of these three approaches. In contrast, total potassium was not associated with either total or cardiovascular mortality.
To our knowledge this is the first time-series analysis that has looked at the association between PM chemical composition and mortality and the association between the underlying factors influencing that composition and mortality. Ozkaynak and Thurston (28) combined source apportionment and epidemiologic methods to assess the effects of air pollution on mortality. However, their study was a cross-sectional analysis rather than a time-series analysis. The present study found significant associations between air pollutants and total nonaccidental and cardiovascular mortality. The association between PM10 and cardiovascular mortality is consistent with previous studies. Zmirou et al. (17) reported an RR for cardiovascular mortality from a 50-µg/m3 increase in PM10 (RR = 1.04) in a study of air pollution in 10 large Western European cities. Pope et al. (5) found an association between respiratory disease death and cardiovascular deaths with PM10 in Utah. Schwartz (2) also found that on high-pollution days (increased total suspended particulates) there was an increased risk of death from cardiovascular disease (RR = 1.09) in Philadelphia, Pennsylvania, and Birmingham, Alabama (1). Furthermore, Anderson et al. (29) found that black smoke was associated with a 0.58% increase in cardiovascular deaths in London.
The association between PM2.5 and cardiovascular mortality is similar to that of Schwartz et al. (3), who found that a 10-µg/m3 increase in PM2.5 was associated with a 1.5% increase in total mortality and 2.1% increase in mortality from ischemic heart disease in a study of six eastern U.S. cities. In contrast to Schwartz et al. (3), the present study also found a significant association between PMCF and total and cardiovascular mortality. Although Schwartz et al. (3) did not find a significant association between coarse fraction and mortality when the results from all six cities were combined, there was an association in Steubenville, Ohio, alone. Such observed differences may have been due to differences in regional coarse fraction composition. In Spokane, Schwartz et al. (12) also found no association between coarse particle concentration and total mortality. However, that study only looked at high episodes of coarse particle concentrations resulting from dust storms. Our findings are in agreement with Ostro et al. (13), who found a significant association between daily PM10 dominated by coarse particles and mortality.
We investigated the possibility that PMCF was a surrogate for dryness by eliminating the days with humidity less than the 25th percentile. Although the association with total mortality was no longer significant, we found a significant association with cardiovascular mortality.
The reason for the negative association between soil and total mortality is unclear. One possible explanation for this observation is related to the fact that the PM2.5 cyclone on the DFPSS was replaced with a WINS on 20 December 1996. The sharper cut point reduced the amount of soil intrusion into the PM2.5 sample, which could produce soil data that are essentially different between the 1995-1996 period and 1997. To assess this effect, we eliminated all 1997 soil data and reevaluated the RR for total mortality. After removing the WINS data, the association between reconstructed soil and total mortality was not significant. These observations are similar to that of Ozkaynak and Thurston (28), who in a study of the association between U.S. mortality rates and particle pollution levels in 1980 found that soil was the least significant predictor of mortality. We also evaluated the association between soil and mortality with only soil data obtained with the WINS. The association was positive and significant (p < 0.05) on the concurrent day, but not significant on any other lag days. However, this observation may be due to the low number of days used to evaluate the association between WINS PM2.5 soil and total mortality (n = 377).
With respect to the elemental components of PM2.5, we found that EC was significantly associated with cardiovascular mortality. EC is found in combustion-derived particles, most notably diesel exhaust (21). We found that KS [potassium from vegetative burning (23)] was also associated with cardiovascular mortality.
We found several associations that are potentially spurious. The associations with these variables were found with only total mortality and not with cardiovascular mortality. Lead was negatively associated with total mortality at lag day 3, although this may be reflective of the moderate correlation between Pb and soil in Phoenix (r = 0.49). Pb may have accumulated in the soil or in road dust from the past use of leaded gasoline. S was also negatively associated with total mortality on lag day 3. At present, the reason for the negative association with S is unclear. However, S accounts for a relatively small percentage of the mass of PM2.5 (15%). The significant negative associations between total mortality and Pb and S were not consistent with the lack of association between these exposure variables and cardiovascular mortality.
For the gaseous species, we found that total nonaccidental mortality and cardiovascular mortality were strongly associated with CO and NO2. These observations are similar to those of Burnett et al. (30), who found associations between CO and NO2 and
total nonaccidental mortality in Toronto, Canada. Burnett et al. (30) also found that cardiac mortality was associated with CO. CO exacerbates cardiac conditions (10). CO concentrations are also associated with hospital admissions for cardiovascular disease (31). In Phoenix the primary sources of CO and NO2 are motor vehicles.
The association between SO2 and cardiovascular mortality was similar to that of Zmirou et al. (17), who also found that an increase in SO2 was associated with an increase in cardiovascular deaths (RR = 1.04). In addition, Zmirou et al. (17) found weak but significant association between 1-hr maximum O3 concentrations and cardiovascular mortality (RR = 1.02). Hoek et al. (32) also found an association between total mortality and O3 in the Netherlands. We found no significant associations with O3.
The present study demonstrated the use of factor analysis in an epidemiologic study. Using factor analysis, we were able to identify those underlying factors of measured air pollution composition variability that were associated with excess mortality. Poisson regression with factor scores as exposure variables revealed that combustion-related pollutants associated with motor vehicles and vegetative burning as well as fine particulate SO4 concentrations were significantly associated with cardiovascular mortality. The soil factor, however, was associated with fewer than expected total deaths. These results are consistent with our time-series results for individual pollutants, specifically CO, NO2, KS, EC, OC, and reconstructed soil. It is interesting to note that the factor representing S was significantly associated with cardiovascular mortality, whereas S alone in an individual pollutant model was not associated with cardiovascular mortality. This may be reflective of the contribution of Pb and Br to the S factor.
A unique aspect of this study was the use of the chemical composition data of PM2.5 Using such data, we found positive associations between cardiovascular mortality and KS, OC, and EC as well as the more traditionally measured pollutants CO, NO2, SO2, PM10, PM2.5, and PMCF. Significant associations were also found with factors associated with incomplete combustion products and particulate S compounds. A limitation of this study is that the factor analysis results are only in terms of the variance explained by each factor, rather than in terms of the quantitative contribution from a specific source category. Although methods are available to include quantitative source apportionments in a time-series framework (33), such an analysis is beyond the scope of this initial investigation.