Mortality effects of a copper smelter strike and reduced ambient sulfate particulate matter air pollution.

BACKGROUND: Numerous studies have reported associations between fine particulate and sulfur oxide air pollution and human mortality. Yet there continues to be concern that public policy efforts to improve air quality may not produce actual improvement in human health.

OBJECTIVES: This study retrospectively explored a natural experiment associated with a copper smelter strike from 15 July 1967 through the beginning of April 1968.

METHODS: In the 1960s, copper smelters accounted for approximately 90% of all sulfate emissions in the four Southwest states of New Mexico, Arizona, Utah, and Nevada. Over the 8.5-month strike period, a regional improvement in visibility accompanied an approximately 60% decrease in concentrations of suspended sulfate particles. We collected monthly mortality counts for 1960-1975 and analyzed them using Poisson regression models.

RESULTS: The strike-related estimated percent decrease in mortality was 2.5% (95% confidence interval, 1.1-4.0%), based on a Poisson regression model that controlled for time trends, mortality counts in bordering states, and nationwide mortality counts for influenza/pneumonia, cardiovascular, and other respiratory deaths.

CONCLUSIONS: These results contribute to the growing body of evidence that ambient sulfate particulate matter and related air pollutants are adversely associated with human health and that the reduction in this pollution can result in reduced mortality.

KEY WORDS: air pollution, copper smelter, mortality, natural experiment, particulate matter, sulfates. Environ Health Perspect 115:679-683 (2007). doi:10.1289/ehp.9762 available via http://dx.doi.org/ [Online 4 January 2007]

**********

Numerous epidemiologic studies have reported associations between ambient particulate matter (PM) air pollution and mortality (Pope and Dockery 2006). Most of these studies have evaluated associations between daily mortality counts and short-term (day-to-day) changes in pollution concentrations using either time-series (Aga et al. 2003; Analitis et al. 2006; Dominici et al. 2003; Samet et al. 2000; Schwartz 2000) or case-crossover (Schwartz 2004) study designs. Various population-based cross-sectional (Evans et al. 1984; Lave and Seskin 1970; Ozkaynak and Thurston 1987) and prospective cohort (Dockery et al. 1993; Jerrett et al. 2005; Laden et al. 2006; Pope et al. 2002; Pope et al. 1995) studies have also evaluated mortality associations with long-term (years to decades) differences in average concentrations across different metropolitan areas or communities. These studies found mortality risk to be consistently associated with concentrations of P[M.sub.2.5] (a common indicator of fine particulate matter including particles with an aerodynamic diameter [less than or equal to] 2.5 [micro]m cut point, which include most sulfate particles). Only a few studies of natural experiments have evaluated large-scale mortality effects associated with sharp, prolonged reductions in air pollution. Reductions in mortality accompanied substantial reductions in PM concentrations associated with the closure of a steel mill in Utah Valley, Utah (Pope et al. 1992), the banning of coal burning in Dublin, Ireland (Clancy et al. 2002), and a reduction in sulfur oxides through the imposition of fuel restrictions in Hong Kong (Hedley et al. 2002).

In the present study we explored an interesting natural experiment associated with a nationwide copper smelter strike from 15 July 1967 through early April 1968 (Cole 1969; Rowland and Greenspoon 1968). The strike had an especially large impact on four Southwest states with copper smelters: New Mexico, Arizona, Utah, and Nevada. In the 1960s, copper smelters accounted for approximately 90% of all sulfate emissions in these Southwest states. As reported by Trijonis and Yuan (1978) and Trijonis (1979), the 1967-1968 strike led to measurable regional reductions in suspended sulfate PM and increases in visibility throughout the Southwest states. Over the 8.5-month strike period, a regional improvement in visibility accompanied an approximately 60% (~ 2.5 [micro]g/[m.sup.3]) decrease in concentrations of suspended sulfate particles. The primary objective of this study is to retrospectively analyze the mortality data and evaluate the evidence of strike-related reductions in mortality.

Methods

Study area. The primary study area in this analysis included the four Southwest states most influenced by the copper smelter strike: New Mexico, Arizona, Utah, and Nevada (Figure 1). As reported and analyzed by Trijonis (1979), visibility data in these four states were collected from urban airports, suburban/nonurban airports, urban NASN (National Air Surveillance Network) sites, and nonurban NASN sites. Phoenix and Tucson, Arizona; Las Vegas, Nevada; Salt Lake City, Utah; and Albuquerque, New Mexico, were the population centers in the states most affected by the strike. Colorado was not included because of the lack of copper smelters, the distance of its population centers from the smelter locations, and the insignificant change in air quality in the Denver area during the strike period (Trijonis 1979). Despite El Paso's close proximity to smelters, Texas was not included because El Paso represented a very small portion of the state's population.

Mortality data. We collected mortality data from the National Center for Health Statistics' yearly mortality reports of the United States from 1960 to 1975 (National Center for Health Statistics 1964-1979). Monthly mortality counts for the four Southwest states, both individually and combined, were compiled. The United States was also divided into four additional regions (rest of United States, Eastern states, neighboring states, and bordering states) according to their proximity to the four Southwest states (Figure 2). Monthly mortality counts for each of these four regions were compiled. Nationwide monthly mortality counts for several specific causes of death were also compiled. Because International Classification of Diseases (ICD) codes changed in January 1968 from the 7th [World Health Organization (WHO) 1957] to the 8th Revision (WHO 1967), ICD-7 codes were used for 1960-1967 and ICD-8 codes were used for 1968-1975. Six nationwide cause-of-death mortality variables were created: a) influenza/pneumonia 1960-1967 (ICD-7 codes 480-483, 490-493); b) influenza/pneumonia 1968-1975 (ICD-8 codes 470-474, 480-486); c) cardiovascular disease 1960-1967 (ICD-7 codes 400-468, 331, 332); d) cardiovascular disease 1968-1975 (ICD-8 codes 390-458); e) respiratory system 1960-1967 (ICD-7 codes 470-527); and f) respiratory system 1968-1975 (ICD-8 codes 460-519).

Analysis. We plotted mortality counts for the four Southwest states over time and over mortality counts for Eastern states, neighboring states, and bordering states. We also analyzed variation in death counts using 12 different Poisson regression models for the four Southwest states as a region as well as for each Southwest state individually. All of the models controlled for time trends in mortality. Because of potential nonlinearities in the time trends, the models included spline smoothers of time with 1-3 degrees of freedom using generalized cross-validation to select the best fit. To further control for non-event-related regional variability in mortality, models 1-4 included total mortality counts from all other U.S. states, the Eastern U.S. states, neighboring states, and bordering states (Figure 2), respectively. Models 5-8 were the same as models 1-4 except that the two variables for nationwide monthly mortality counts for influenza/pneumonia (using ICD-7 and ICD-8 codes for relevant time period) were also added to the models to control for the effects of a nationwide epidemic or other regional epidemics. Models 9-12 most aggressively controlled for U.S. non-event-related mortality trends by including all six nationwide cause-of-death variables (influenza/pneumonia, cardiovascular disease, and other respiratory deaths, as discussed above).

Because the strike began in the middle of July and because there may be a lagged pollution effect on mortality, we created three strike period indicator variables to account for different possible lagged effects: a) July 1967-March 1968 (all strike months including the first month that was only partially effected by the strike), b) August 1967-March 1968 (only the full strike months), and c) August 1967-April 1968 (full strike months plus 1-month lag). Statistical analyses were conducted using R and SAS software estimating the cubic spline smooths using the "gam" function in the MGCV package (Wood 2003; Wood 2004) in R (R Development Core Team 2006) and the "spline" function using PROC GAM in SAS (release 9.1.3; SAS Institute, Inc., Cary, NC).

Results

Monthly mortality counts from the four Southwest states plotted over time are presented in Figure 3. Characteristic nonstationary Poissonian variability, time trend, and seasonal components of mortality counts are observable, but the effect of the strike period is small compared with other sources of variability. Figure 4 plots monthly mortality counts in the four Southwest states over monthly mortality counts for bordering states, which presumably have similar time trends, seasonality, and other regional similarities. Although the mortality counts in the four Southwest states were consistently below what would be expected based on the mortality trends of bordering states, the relative reductions in mortality were small.

Small but statistically significant strike-related decreases in mortality were observed based on the Poisson regression models. Figure 5 presents the estimates of the percent decrease in mortality [and 95% confidence intervals (CIs)] in the four Southwest states associated with the three different strike indicator periods for the 12 different Poisson regression models. Strike period mortality effect estimates indicated a decrease in mortality of approximately 1.5-4.0% depending on the model and strike period indicator used. For example the estimated percent decrease in mortality was 2.5% (95% CI, 1.1-4.0%), based on a model that used a strike period indicator for full strike months plus a 1-month lag (August 1967-April 1968) and that controlled for time trends, mortality counts in bordering states, and nationwide monthly mortality counts for influenza/pneumonia, cardiovascular, and other respiratory deaths. As can be seen in Figure 6, similar patterns in the effect estimates were observed for each of the four Southwest states individually.

Figure 7 presents plots of nationwide monthly mortality counts for influenza and pneumonia. The Hong Kong pandemic of 1968-1969 (Kilbourne 2006) along with other periods of elevated influenza/pneumonia, are observable. However, as demonstrated by the results presented in Figures 5 and 6, the strike period effect estimates were not highly sensitive to controlling for nation-wide monthly mortality counts for influenza/pneumonia, cardiovascular, and other respiratory deaths. The effect estimates were also not sensitive to whether the models controlled for mortality trends in the rest of the United States, eastern U.S. states, or neighboring states, but were consistently smaller when counts for the bordering states were included. Somewhat larger effect estimates were observed when the strike indicator did not include the first partial month of the strike and included a 1-month lag.

Discussion

The 8.5-month copper smelter strike in 1967-1968 provided a unique and interesting natural experiment to explore the mortality effects of air pollution. Relative to the natural experiment studies of Utah Valley (Pope et al. 1992) and Dublin (Clancy et al. 2002), the copper smelter strike resulted in smaller reductions in fine particulate pollution, but the reductions occurred over a much larger geographic region. This retrospective analysis of available mortality data suggests that the approximately 60% reduction in regional sulfate concentrations (~ 2.5 [micro]g/[m.sup.3]) associated with the copper smelter strike produced a small but measurable reduction in mortality (1.5-4.0%). The effect estimates were reasonably stable and robust to controlling for time trends, mortality trends in other areas of the United States, and nationwide monthly mortality counts for influenza/pneumonia, cardiovascular, and other respiratory deaths. The effect estimates were consistently reduced when mortality counts in bordering states were included in the models. Mortality counts from bordering states are more likely to have similar mortality trends and seasonality and other regional similarities, as well as common influences of regional communicable disease. Therefore, the use of mortality counts from bordering states as a control variable in the regression may be preferable to including mortality counts from more distant areas of the United States. However, because the air quality of the bordering states may also have been somewhat improved because of the smelter strike, using counts in the bordering states may result in overcontrolling and underestimating the mortality effect.

The estimated strike-related reduction in mortality of approximately 1.5-4.0% is consistent with what would be expected given the reduction in ambient concentrations of sulfate particles and estimated mortality effects from the relevant literature. Based on data presented by Trijonis and Yuan (1978), the average, seasonally adjusted, regional strike-related reduction in concentrations of sulfate particles was approximately 2.5 [micro]g/[m.sup.3]. Based on reported mortality reductions from the natural experiment studies of Utah Valley (Pope et al. 1992) and Dublin (Clancy et al. 2002), a 2.5 [micro]g/[m.sup.3] reduction in P[M.sub.2.5] concentrations would result in approximately a 0.8-1.1% decline in overall mortality, depending on assumptions about PM conversion factors and methods of extrapolation (Pope and Dockery 2006). The estimated percent change in mortality associated with a 2.5-[micro]g/[m.sup.3] change in P[M.sub.2.5] from long-term prospective cohort studies would equal approximately 1.5% based on nationwide American Cancer Society studies (Pope et al. 1995, 2002) and up to 4.0% based on the Harvard Six Cities (Dockery et al. 1993; Laden et al. 2006) and American Cancer Society, Los Angeles intrametropolitan area studies (Jerrett et al. 2005). In a recent U.S. Environmental Protection Agency (EPA) expert elicitation (U.S. EPA 2006) on the assessment of the concentration response relationship between P[M.sub.2.5] exposure and mortality, multiple experts suggest that a long-term 1-[micro]g/[m.sup.3] decrease in P[M.sub.2.5] would lead to a mean reduction in mortality of 0.4-2.0% with credible estimates ranging from 0 to 3.0%--comparable with effect estimates observed in this analysis.

There are several different reasons to expect different mortality effects from this natural experiment versus other relevant particulate air pollution studies. Given the ubiquitous regional exposure and the high propensity of sulfate particles to penetrate indoors (Sarnat et al. 2000, 2006; Suh et al. 1993), the smelter strike-related reduction in ambient air pollution may reflect relatively large reductions in personal exposures. Specific copper smelter-related sulfate particles and related co-pollutants may have different relative toxicities. Although detailed data on changes in levels of transition metals are not available, it is likely that there were important reductions in transition metals associated with the temporary closure of the copper smelters. There is evidence that fine particles with proportionally high levels of transition metals are relatively more toxic (Ghio 2004; Ghio and Cohen 2005). The strike may have also caused an exodus of smelter workers to neighboring regions, which may have artificially lowered the death counts in the four affected Southwest states during the strike period. Evidence that this is not the explanation for the strike-related reduction in mortality is provided by models that adjusted for mortality counts in bordering states. Adjustment for bordering state mortality counts reduced the strike effect estimates, whereas the opposite would be expected if the workers left the affected states to increase disproportionately the population of the bordering states.

A limitation of this study is that other events or factors that occurred approximately concurrent with the smelter strike may have caused the reduction in mortality in these four Southwest states and were not adequately controlled for by time trend variables, mortality counts in neighboring areas, and nationwide mortality counts for influenza/pneumonia, cardiovascular, and other respiratory deaths.

Another limitation of this analysis is the lack of spatial resolution with regard to the pollution and mortality data. For example, the results of this analysis are generally consistent with the interpretation that the copper smelter strike resulted in an 8.5-month reduction in fine sulfate PM, contributing to a concurrent reduction in mortality. However, one finding that is somewhat inconsistent with this interpretation is that similar strike-related mortality associations were observed in all four Southwest states, including New Mexico. However, based on analysis presented by Trijonis (1979), changes in visibility and sulfate concentrations in New Mexico were not as largely affected by the strike as the other three states. Another limitation of this analysis comes from the lack of available monthly data on cause of mortality for each state. These data would allow for an analysis to be conducted specifically on the effect of PM on cardiorespiratory related deaths.

Evidenced by the elevated effect estimate when the lagged month of April was included as part of the strike period indicator, this analysis suggests a lag effect of a few weeks to a month. There is no evidence that the lag effect lasted longer than 1 month.

In conclusion, the results of this study add to the growing body of evidence that ambient sulfate particulate matter and related air pollutants are associated with adverse human health effects and that the reduction in this pollution can reduce mortality. These results suggest that the approximately 60% reduction in regional sulfate concentrations (~ 2.5 [micro]g/[m.sup.3]) associated with the copper smelter strike produced improvements not just in visibility but also in small but measurable reductions in mortality.

CORRECTION

The concentration of sulfate particles noted in the manuscript originally published online as 1-1.5 [micro]g/[m.sup.3] has been corrected to ~ 2.5 [micro]g/[m.sup.3] in the introduction and the "Discussion." Also in the "Discussion," the decline in overall mortality originally was 0.2-0.7% and the estimated percent changes in mortality were 0.6-0.9% and 1.7-2.6%. They have been corrected here to 0.8-1.1%, 1.5%, and 4.0%, respectively.

REFERENCES

Aga E, Samoli E, Touloumi G, Anderson HR, Cadum E, Forsberg B, et al. 2003. Short-term effects of ambient particles on mortality in the elderly: results from 28 cities in the APHEA2 project. Eur Respir J 21 (40):28s-33s.

Analitis A, Katsouyanni K, Dimakopoulou K, Samoli E, Nikoloulopoulos AK, Petasakis Y, et al. 2006. Short-term effects of ambient particles on cardiovascular and respiratory mortality. Epidemiol 17:230-233.

Clancy L, Goodman P, Sinclair H, Dockery DW. 2002. Effect of air-pollution control on death rates in Dublin, Ireland: an intervention study. Lancet 360:1210-1214.

Cole JW. 1969. Bureau of Mines Minerals Yearbook: Metals, Minerals, and Fuels 1968, Copper. Washington, DC:U.S. Government Printing Office. Available: http://digital.library.wisc.edu/1711.dl/EcoNatRes.MinYB1968v1-2 [accessed 7 August 2006].

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:1753-1759.

Dominici F, Daniels M, McDermott A, Zeger SL, Samet J. 2003. Shape of the exposure-response relation and mortality displacement in the NMAPS database. In: Revised Analysis of Time-Series of Air Pollution and Health. Special Report. Boston:Health Effects Institute, 91-96.

Evans JS, Tosteson T, Kinney PL. 1984. Cross-sectional mortality studies and air pollution risk assessment. Environ Int 10:55-83.

Ghio AJ. 2004. Biological effects of Utah Valley ambient air particles in humans: a review. J Aeros Med 17:157-164.

Ghio AJ, Cohen MD. 2005. Disruption of iron homeostasis as a mechanism of biologic effect by ambient air pollution particles. Inhal Toxicol 17:709-716.

Hedley AJ, Wong CM, Thach TQ, Ma S, Lam TH, Anderson HR. 2002. Cardiorespiratory and all-cause mortality after restrictions on sulphur content of fuel in Hong Kong: an intervention study. Lancet 360:1646-1652.

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. Epidemiol 16:727-736.

Kilbourne E. 2006. Influenza pandemics of the 20th century. Emerg Infect Dis 12(1):9-14.

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):581-582.

Lave LB, Seskin EP. 1970. Air pollution and human health. Science 169:723-733.

National Center for Health Statistics. 1964-1979. Vital Statistics of the United States, 1960-1975. Vol II, Mortality, Part A. Washington, DC:National Center for Health Statistics.

Ozkaynak H, Thurston GD. 1987. Associations between 1980 U.S. mortality rates and alternative measures of airborne particle concentration. Risk Anal 7: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:1132-1141.

Pope CA III, Dockery DW. 2006. 2006 Critical review: health effects of fine particulate air pollution: lines that connect. Air Waste 56:709-742.

Pope CA III, Schwartz J, Ransom MR. 1992. Daily mortality and PM10 pollution in the Utah Valley. Arch Environ Health 47:211-217.

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:669-674.

R Development Core Team. 2006. R: A Language and Environment for Statistical Computing. Vienna, Austria:R Foundation for Statistical Computing. Available: http://www.R-project.org.

Rowland JA, Greenspoon GN. 1968. Bureau of Mines Minerals Yearbook: Metals, Minerals, and Fuels 1967, Copper. Washington, DC:U.S. Government Printing Office. Available: http://digital.library.wisc.edu/1711.dl/EcoNatRes.MinYB1967v1-2 [accessed 7 August 2006].

Samet JM, Dominici F, Curriero FC, Coursac I, Zeger SL. 2000. Fine particulate air pollution and mortality in 20 U.S. cities, 1987-1994. N Engl J Med 343:1742-1749.

Sarnat JA, Koutrakis P, Suh HH. 2000. Assessing the relationship between personal particulate and gaseous exposures of senior citizens living in Balimore, MD. Air Waste 50(7):1184-1198.

Sarnat SE, Coull BA, Schwartz J, Gold DR, Suh HH. 2006. Factors affecting the association between ambient concentrations and personal exposures to particles and gases. Environ Health Perspect 114:649-654.

Schwartz J. 2000. Assessing confounding, effect modification, and thresholds in the association between ambient particles and daily deaths. Environ Health Perspect 108:563-568.

Schwartz J. 2004. The effects of particulate air pollution on daily deaths: a multi-city case-crossover analysis. Occup Environ Med 61:956-961.

Suh HH, Koutrakis P, Spengler JD. 1993. Validation of personal exposure models for sulfate and aerosol strong acidity. Air Waste 43(6):845-850.

Trijonis J. 1979. Visibility in the southwest--an exploration of the historical data base. Atmos Environ 13:833-843.

Trijonis J, Yuan K. 1978. Visibility in the Southwest: An Exploration of the Historical Data Base. EPA 600/3-78-039,82-100. Washington, DC:U.S. Environmental Protection Agency.

U.S. EPA. 2006. Expanded Expert Judgment Assessment of the Concentration-Response Function Relationship between P[M.sub.2.5] Exposure and Mortality. Washington DC:U.S. Environmental Protection Agency. Available: http://www.epa.gov/ttn/ecas/regdata/Uncertainty/pm_ee_report.pdf [accessed 30 November 2006].

WHO. 1957. Manual of the International Statistical Classification of Diseases, Injuries, and Causes of Death. Based on the Recommendation of the 7th Revision Conference, 1955, and Adapted by the 9th World Health Assembly under the WHO Nomenclature Regulations. Geneva:World Health Organization.

WHO. 1967. International Classification of Diseases Manual of the International Statistical Classification of Diseases, Injuries, and Causes of Death. Based on the Recommendation of the 8th Revision Conference, 1965, and Adapted by the 19th World Health Assembly. Geneva:World Health Organization.

Wood SN. 2003. Thin-plate regression splines. J R Stat Soc [Ser B] 65(1):95-114.

Wood SN. 2004. Stable and efficient multiple smoothing parameter estimation for generalized additive models. J R Stat Soc 99:673-686.

C. Arden Pope III, Douglas L. Rodermund, and Matthew M. Gee

Department of Economics, Brigham Young University, Provo, Utah, USA

Address correspondence to C.A. Pope III, 142 FOB, Brigham Young University, Provo, UT 84602-2363 USA. Telephone: [phone]. Fax: [phone]. E-mail: [e-mail]

We thank W. Wilson and D. Dockery for directing us to the Trijonis visibility paper and for helping motivate this research project.

This research was supported in part by funds from the Mary Lou Fulton Professorship, Brigham Young University, Provo, Utah.

The authors declare they have no competing financial interests.

Received 26 September 2006; accepted 4 January 2007.