Title: Mortality Displacement in the Association of Ozone with Mortality: An Analysis of 48 US Cities Authors: Antonella Zanobetti and Joel Schwartz ONLINE DATA SUPPLEMENT
Additional Information on Materials and Methods and Results Study population We obtained individual mortality data for forty-eight cities in the United States from the National Center for Health Statistics (NCHS) for the years 1989 to 2000. The chosen cities were a convenience sample of cities for which we had previously extracted data, enhanced by additional cities chosen to ensure a reasonable spread of cities with respect to climatic patterns and geography. The mortality files provided information on the exact date of death, and the underlying cause of death. For this study we selected all-cause daily mortality excluding any deaths from accidental causes (ICD- code 9 th revision: 1-799, ICD-code 10 th revision: V01-Y98); we also examined any cardiovascular disease (ICD-9:390-429, ICD-10:I01-I60), any respiratory disease (ICD-9:460-519, ICD-10:J01-J99), and stroke (ICD-9: 436, ICD-10: I64). City characteristics such as population density, percentage of population 65 years of age and older in poverty status, and Percent of Households with Householder 65+ and income > $50000 were obtained from the 2000 United States census and we calculated the mean and variance of the daily summer (June through August) apparent temperature using data from the National Weather Service Surface Station (E1). We calculated the percentage of households with central air conditioning (AC) using data from the American Housing Survey of the US Census Bureau (E2). 1
Environmental Data We obtained ozone (8-hour) data from US Environmental Protection Agency s Air Quality System Technology Transfer Network (USEPA Technology Transfer Network, 2005). When multiple monitors were present in a city we estimated an average daily value for the city, and we accounted for the impact of occasional missing values using an algorithm previously described (E3, E4). We obtained local meteorological data such as mean temperature, and dew point temperature, from the United States Surface Airways and Airways Solar Radiation hourly data (1). We also computed apparent temperature (AT), defined as an individual s perceived air temperature given the humidity. AT was calculated with the following formula (E5, E6). AT = -2.653 + (0.994*Ta) + (0.0153*Td 2 ) where Ta is mean temperature and Td is dew point temperature. Statistical Methods In each city we computed daily counts of mortality; the mortality data was then merged to the air pollution and weather daily data, and we computed lags for ozone and temperature. Since the ozone effect is either restricted to the warm season, or at least is substantially different outside that period, we wanted to ensure that any changes in effect estimate reflect the effects of lags, and not of seasonal effect modification. Since we wanted to examine the effect up to the previous 20 days, we want to ensure that the ozone levels of the previous 20 days are part of the warmer season. Mortality declines from 2
January, but is still noticeably higher in April and early May than in August; and in some cities there are still cold temperatures in the month of April; these could cause seasonal confounding in the longer lags. Therefore we analyzed daily deaths that occurred in June- August, so that all the lags were still in the warm season. Temperature also may affect daily deaths with a lag, however a previous examination of the pattern in US cities showed mortality effects of warm temperatures disappeared within seven days (E7). Previous studies have generally not included this many lags of temperature in examining the effects of ozone, which might have allowed residual confounding. We investigated the association between daily concentrations of air pollutants and mortality using a two stage approach. First, we fit a time series analysis in each city separately and then we combined the city-specific results using a maximum likelihood meta-analytic approach. The time series of daily counts of mortality and ozone were investigated with a generalized linear model, with a quasi-poisson link function to account for overdispersion. In the model we controlled for season using natural splines with 2 degrees of freedom for each 3-month summer season, day of the week with indicator variables, and an unconstrained distributed lag for apparent temperature to take into account the effect of temperature today and the previous 7 days. Each lag of apparent temperature was treated as linear in these analyses, which were restricted to the warm season. We examined the dependence of daily deaths on 8-hr ozone concentrations on the day of death to use as a baseline analysis. We then refit models using ozone on the day of 3
death, and up to the previous 20 days using an unconstrained distributed lag model. If the pollution-related deaths are only being advanced by a few days to a few weeks we would see this harvesting effect expressed by a negative associations between air pollution and deaths several days to weeks subsequently. The effect of air pollution, net of any such short term rebound up to 20 days is the sum of the positive and negative effect estimates for all 21 days. Hence if this second estimate were smaller than the estimate using only lag 0 ozone, this would support the harvesting hypothesis. For Poisson regression the unconstrained distributed lag model can be written as: Log( E[ Y ]) = covariates + b Z + b Z +... + b Z (1) t 0 t 1 t- 1 q t- q Where Z t-p is the 8-hour ozone concentration p days prior to day t. The unconstrained distributed lag model is too noisy to provide any information about the shape of the effect versus lag, but it does give an unbiased estimate of the overall effect, computed as the sum of the ß j. Several methods have been used in order to study the shape of the distributed lag. One common approach is to constrain the β s to follow a flexible polynomial, such as cubic or quadratic. An extension of this idea, which aims for more flexibility than that afforded by a polynomial functions, is to subject the coefficients to non-parametric smoothing. As previously described (E8) we used a penalized quasi-likelihood to estimated the coefficient of the distributed lag. These models were fit in R (R Development Core Team (2007). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.r-project.org.). 4
In a second stage of the analysis, the city specific results were combined using the meta-regression technique of Berkey and coworkers (E9). To be conservative we report the results incorporating a random effect, whether or not there was a significant heterogeneity. The I 2 statistic is a generalization of the Χ 2 or Q test for heterogeneity, and express the proportion of variance explained. We used the I 2 statistic to assess the proportion of total variation in effect estimates that was due to between-cities heterogeneity (E10). We used the following formula: I 2 = [Q / (k - 1)] - 1 / [Q / (k - 1)], where Q is the Q-test for heterogeneity and k is the number of cities. We examined effect modification by city characteristics by entering them as predictor variables in the meta-regression. These included measures of socio-economic condition (Percent of Households with Householder 65+ and income > $50000, Percent population aged 65+ in poverty status), exposure related measures (Mean and Variance of summer apparent temperature), general social factors (population density), and percent population with central air conditioning (AC). 5
Results The cities analyzed in this study are listed in Table E1 together with the total population and the count of total and cause-specific summer time deaths per city. Figure E1 shows the city specific distribution of the daily all cause mortality during the June- August months; New York City, Los Angeles and Chicago are the biggest cities. Overall, we examined 1,614,124 deaths. Table E2 presents descriptive statistics of apparent temperature and ozone showing a wide range of ozone exposures and accompanying temperatures across the 48 cities. For example, the 75 th percentile of 8-hour ozone ranged from 19.8 ppb in Honolulu to 75.9 ppb in Los Angeles, while the 75 th percentile for apparent temperature ranged from 20.6 C in Spokane to 34.6 C in Houston. Figure 2 shows boxplots of the daily ozone distribution in each city, sorted by level of concentration of ozone. Figure E3 presents the city-specific distribution of apparent temperature during the summer months. 6
References E1. National Environmental Satellite, and a. I. S. Data. TD-3280 U.S. Surface Airways and Airways Solar Radiation Hourly. 2003. E2. US Census Bureau. American Housing Survey for the United States: 2001. U.S. Department of Commerce/U.S. Department of Housing and Urban Development, Washington, DC. 2002.http://www.census.gov/prod/2002pubs/h150-01.pdf. E3. Zanobetti A, Schwartz J, and Dockery DW. Airborne particles are a risk factor for hospital admissions for heart and lung disease. Environ Health Perspect 2000; 108(11):1071-7. E4. Schwartz J. The distributed lag between air pollution and daily deaths. Epidemiology 2000; 11(3):320-326. E5. Steadman RG The assessment of sultriness. Part II: Effects of wind, extra radiation and barometric pressure on apparent temperature. Journal of Applied Meteorology 1979; 18:874-885. E6. Kalkstein LS, Valimont KM. An evaluation of summer discomfort in the United States using a relative climatological index. Bulletin American Meteorological Society 1986; 67(7):842-848. E7. Braga AL, Zanobetti A, Schwartz J. The time course of weather-related deaths. Epidemiology 2001; 12(6):662-7. E8. Zanobetti A, Wand MP, Schwartz J, Ryan LM. Generalized additive distributed lag models: quantifying mortality displacement. Biostatistics 2000; 1(3):279-92. 7
E9. Berkey CS, Hoaglin DC, Antczak-Bouckoms A, Mosteller F, Colditz GA. Metaanalysis of multiple outcomes by regression with random effects. Stat Med 1998; 17(22):2537-50. E10. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med 2002; 21( ):1539-1558. 8
FIGURE LEGENDS Figure E1: City-specific distribution of the daily counts of all cause mortality for the June-August period for the years 1989-2000. Figure E2: City-specific distribution of the daily levels of concentrations of ozone during the June-August period for the years 1989-2000. Figure E3: City-specific distribution of the daily levels of apparent temperature during the June-August period for the years 1989-2000. 9
Table E1: city-specific descriptive: population, and counts of all causes and cause-specific deaths deaths during the June-August months Population Number of deaths by 1000 All cause CVD Respiratory Stroke City State by 1000 Albuquerque NM 557 11 5 2 1 Atlanta GA 2678 42 24 8 5 Austin TX 812 10 5 2 1 Baltimore MD 1405 44 24 11 6 Birmingham AL 992 23 12 5 3 Boston MA 2806 65 34 16 7 Boulder CO 291 3 1 1 0 Broward FL 1623 39 24 7 4 Canton OH 378 10 6 3 2 Charlotte NC 695 13 7 3 2 Chicago IL 5377 123 65 26 14 Cincinnati OH 845 25 14 6 3 Cleveland OH 1394 43 26 10 5 Colorado Springs CO 517 6 3 2 1 Columbus OH 1069 24 14 5 3 Dallas TX 2219 36 19 7 4 Denver CO 555 14 7 4 1 Detroit MI 2061 54 31 11 6 Greensborough NC 421 9 5 2 1 Honolulu HI 876 15 9 4 2 Houston TX 3401 55 28 11 7 Jersy city NJ 609 12 8 3 1 Kansas City KS 451 24 13 6 3 Los Angeles CA 9519 158 106 39 21 Miami FL 2253 52 34 12 5 Milwaukee WI 940 26 14 6 3 Nashville TN 570 18 9 4 2 New Haven CT 824 21 13 5 2 New Orleans LA 485 15 8 2 2 NYC NY 8008 177 150 28 12 Oklahoma City OK 660 19 11 4 2 Orlando FL 896 16 9 4 2 Philidelphia PA 1518 49 26 9 5 Phoenix AZ 3072 50 24 11 6 Pittsburg PA 1282 46 26 11 6 Provo/Orem UT 369 3 2 1 0 Sacramento CA 1223 22 13 6 3 St. Louis MO 1584 43 24 10 5 Salt Lake City UT 898 13 6 3 2 San Diego CA 2814 48 29 13 7 San Francisco CA 777 20 13 5 3 Seattle WA 1737 32 15 8 5 Spokane WA 418 10 5 3 1 Tampa FL 999 22 12 5 3 Terra Haute IN 106 4 2 1 0 Tulsa OK 563 15 8 3 2 Washington DC DC 762 21 12 4 2 Youngstown OH 370 10 6 3 1
Table E2: City-specific descriptive for apparent temperature and ozone during the June-August months Apparent temperature ( C) 8-hr Ozone (ppb) City State min 25th p mean 75th p max min 25th p mean 75th p max Albuquerque NM 13.0 21.7 23.1 24.8 29.9 23.6 45.8 51.3 57.2 74.9 Atlanta GA 15.0 27.5 29.0 31.1 35.8 7.3 43.3 58.8 73.2 142.0 Austin TX 22.5 31.5 32.5 34.0 41.6 8.5 29.8 42.3 53.3 102.5 Baltimore MD 11.5 23.3 26.3 29.7 39.7 1.9 43.7 56.2 68.2 115.2 Birmingham AL 16.6 28.0 29.5 31.8 37.4 13.0 37.2 51.1 63.5 104.6 Boston MA 8.9 19.3 22.6 26.0 37.0 6.1 32.1 44.9 55.8 105.8 Boulder CO 2.1 18.1 19.9 22.3 27.2 5.9 45.5 53.7 62.4 100.5 Broward FL 26.8 32.5 33.5 34.7 37.1 11.2 20.4 26.9 30.5 75.8 Canton OH 6.2 19.3 22.4 25.7 33.8 4.3 45.5 56.7 68.1 100.2 Charlotte NC 12.7 26.4 28.3 30.6 36.6 12.6 49.1 61.2 73.9 123.0 Chicago IL 7.8 20.1 23.5 27.1 39.4 1.9 31.4 41.8 50.7 96.8 Cincinnati OH 9.8 22.3 25.2 28.4 36.7 5.5 41.1 52.7 64.2 104.1 Cleveland OH 7.7 20.0 23.1 26.6 35.4 4.3 35.4 48.1 59.4 96.0 Colorado Springs CO 2.4 16.6 18.2 20.6 25.1 10.5 40.3 46.1 52.5 81.3 Columbus OH 9.2 21.6 24.6 27.8 36.8 6.3 43.9 54.8 66.6 108.2 Dallas TX 19.2 30.7 32.3 34.6 38.6 7.9 32.9 47.7 60.3 109.8 Denver CO 2.1 18.1 19.9 22.3 27.2 3.2 37.7 45.3 52.9 92.8 Detroit MI 7.3 19.8 23.1 26.4 36.5 4.4 33.0 45.0 56.2 102.4 Greensborough NC 10.5 24.9 26.8 29.4 35.5 9.0 47.0 58.2 68.9 102.9 Honolulu HI 26.2 28.8 29.8 30.7 33.9 2.0 10.3 15.1 19.8 34.3 Houston TX 22.3 31.7 32.9 34.6 37.4 7.1 25.6 42.7 55.9 118.6 Jersy city NJ 11.4 22.8 25.9 29.3 40.8 2.0 39.5 54.7 69.1 142.5 Kansas City KS 12.3 23.5 26.9 30.7 37.4 4.9 40.8 50.2 58.9 101.0 Los Angeles CA 12.8 19.4 21.0 22.3 29.3 7.3 47.5 62.8 75.4 146.2 Miami FL 26.8 32.5 33.5 34.7 37.1 11.2 20.4 26.9 30.5 75.8 Milwaukee WI 6.9 18.5 22.0 25.9 39.9 5.7 35.5 47.1 56.5 118.0 Nashville TN 15.0 26.6 28.5 31.2 37.2 6.0 36.8 48.5 60.3 102.1 New Haven CT 9.7 19.7 22.9 26.1 36.3 4.8 35.4 48.6 59.3 126.3 New Orleans LA 23.6 31.5 32.8 34.5 38.2 2.4 24.4 37.7 50.1 102.3 NYC NY 12.1 22.5 25.5 28.7 40.1 0.9 30.6 44.0 55.3 116.0 Oklahoma City OK 15.6 26.8 29.2 32.1 36.8 5.3 42.5 53.0 63.3 103.5 Orlando FL 24.7 30.8 31.9 33.1 36.3 8.1 29.2 41.6 53.1 100.4 Philidelphia PA 10.8 23.3 26.5 29.9 41.1 2.0 37.7 51.3 64.4 119.1 Phoenix AZ 17.6 31.4 33.5 36.2 40.5 20.5 48.3 56.2 64.4 90.3 Pittsburg PA 7.0 20.2 23.1 26.3 34.8 1.2 39.7 52.5 64.4 112.8 Provo/Orem UT 5.8 19.7 22.3 25.3 30.7 12.6 50.2 56.8 63.1 106.0 Sacramento CA 11.1 19.7 22.2 24.4 34.2 17.7 44.1 55.4 65.7 109.0 St. Louis MO 11.3 24.8 28.1 31.7 40.0 5.2 39.7 51.0 61.6 104.1 Salt Lake City UT 5.8 19.7 22.3 25.3 30.7 11.5 48.3 56.6 64.6 97.8 San Diego CA 13.3 19.6 21.4 23.1 31.1 17.7 37.9 45.2 50.8 90.8 San Francisco CA 9.7 14.4 15.8 17.0 28.0 3.1 15.0 19.9 23.8 65.6 Seattle WA 7.6 14.3 16.7 18.7 28.3 7.1 24.1 35.6 44.3 103.5 Spokane WA 4.2 14.0 17.2 20.6 29.7 11.8 36.8 44.3 51.4 77.4 Tampa FL 25.2 31.8 32.9 34.3 37.8 10.1 28.1 39.8 50.0 95.1 Terra Haute IN 9.4 21.8 25.0 28.3 37.0 6.4 43.1 54.5 65.0 107.1 Tulsa OK 15.4 27.4 30.1 33.4 37.2 11.4 44.4 54.9 64.5 111.0 Washington DC DC 11.9 24.8 27.7 30.9 39.1 2.5 41.3 53.7 65.8 121.9 Youngstown OH 5.9 18.5 21.6 25.1 35.3 2.0 37.9 51.3 65.1 106.1
FIG E1 New York City, NY Los Angeles, LA Chicago, IL Boston, MA Houston, TX Detroit, MI Miam i, FL Phoenix, AZ Philadelphia, PA San Diego, CA Pittsburgh, PA Baltimore, MD Saint Louis, MO Cleveland, OH Atlanta, GA Broward, FL Dallas, TX Seattle, WA Milwaukee, WI Cincinnati, OH Kansas City, KS Columbus, OH Birmingham, AL Tampa, FL Sacramento, CA Washington DC New Haven, CT San Francisco, CA Oklahoma City, OK Nashville, TN Orlando, FL New Orleans, LA Honolulu, HI Tulsa, OK Denver, CO Salt Lake City, UT Jersey City, NJ Charlotte, NC Albuquerque, NM Youngstown, OH Spokane, WA Canton, OH Aus tin, TX Greensborough, NC Colorado Springs, CO Terra Haute, IN Provo, UT Boulder, CO 0 20 40 60 80 100 120 140 160 180 200 All cause mortality
FIG E2 Charlotte, NC Los Angeles, LA Atlanta, GA Greensborough, NC Provo, UT Baltimore, MD Phoenix, AZ Canton, OH Salt Lake City, UT Washington DC Tulsa, OK Boulder, CO Columbus, OH Jersey City, NJ Terra Haute, IN Sacramento, CA Cincinnati, OH Oklahoma City, OK Pittsburgh, PA Albuquerque, NM Philadelphia, PA Saint Louis, MO Youngstown, OH Birmingham, AL Kansas City, KS Nashville, TN Cleveland, OH Colorado Springs, CO Denver, CO Dallas, TX Milwaukee, WI New Haven, CT San Diego, CA Spokane, WA Detroit, MI Boston, MA New York City, NY Chicago, IL Orlando, FL Aus tin, TX Tampa, FL New Orleans, LA Houston, TX Seattle, WA Miam i, FL Broward, FL San Francisco, CA Honolulu, HI 0 20 40 60 80 100 120 140 Ozone (ppb)
FIG E3 Phoenix, AZ Miam i, FL Broward, FL Houston, TX Dallas, TX New Orleans, LA Tampa, FL Aus tin, TX Orlando, FL Tulsa, OK Oklahoma City, OK Birmingham, AL Honolulu, HI Atlanta, GA Nashville, TN Charlotte, NC Saint Louis, MO Washington DC Greensborough, NC Kansas City, KS Philadelphia, PA Baltimore, MD Jersey City, NJ New York City, NY Cincinnati, OH Terra Haute, IN Columbus, OH Chicago, IL Salt Lake City, UT Provo, UT Albuquerque, NM Pittsburgh, PA Cleveland, OH Detroit, MI New Haven, CT Boston, MA Canton, OH Milwaukee, WI Sacramento, CA Youngstown, OH San Diego, CA Los Angeles, LA Denver, CO Boulder, CO Colorado Springs, CO Spokane, WA Seattle, WA San Francisco, CA 5 10 15 20 25 30 35 40 Apparent Temperature ( C)