Effect of Air Pollution on Daily Mortality in Hong Kong.
Time-series methods are widely used for assessment of short-term health effects of air pollution (1). Although limitations arise from ecologic fallacy (2) and the harvesting effect (3-5), time-series methods are more powerful and better able to characterize the population exposure effects than those based on geographic aggregations in cross-sectional studies (6). Also, methods to control for time-related confounding factors are well established (7). Daily time-series analysis is not applicable to the estimation of longer-term chronic exposure effects of air pollution (8), which are public health concerns. Daily time-series analysis may be better estimated from longitudinal studies, but it can be used to assess the potential health benefits of air quality intervention in terms of the number of hospital admissions and deaths avoidable if days with high concentrations (according to a chosen reference value) were eliminated, thus providing information to support the setting of air quality objectives (9,10).
To date, there is coherent evidence that air pollution has short-term effects on mortality (9,11-15), but the questions whether there are independent effects of a single pollutant to account for a health outcome under study and whether there are thresholds and linear or non-linear relationships are still not settled.
In the United States, particulates are regarded as the pollutants that account for most excess mortality due to air pollution (16), but in Europe several studies indicated a stronger association with sulfur dioxide (17). Some showed that it might be the sulfuric acid (18), acid aerosol (19), and mass concentration (20) associated with particulates that are responsible for the effects. Other studies showed that independent effects of individual pollutants cannot be identified in light of the complexity and variability of the air pollution mixtures to which people are exposed (21,22).
In this study we assessed the effects of air pollution on mortality outcomes and identified which pollutants would contribute most to the effect in Hong Kong, a subtropical city in the Asian Pacific rim. Patterns of exposure-response relationships for four criteria pollutants, nitrogen dioxide, sulfur dioxide, particulate matter [is less than] 10 [micro]m in aerodynamic diameter ([PM.sub.10]), and ozone were assessed during warm and cool seasons with a view to ascertaining their effects on the commonly used mortality outcomes.
Materials and Methods
Data. For the period 1995-1997, we obtained daily death counts for all nonaccidental causes [International Classification of Diseases, Revision 9 (ICD-9) [is less than] 800 (23)], respiratory disease (ICD-9 460-519), and cardiovascular disease (ICD-9 390-459) from the Census & Statistics Department (Hong Kong Special Administrative Region, People's Republic of China); meteorologic data (daily mean temperature and relative humidity) from the Hong Kong Observatory; and air pollutant concentrations (from two to seven monitoring stations) from the Environmental Protection Department. Daily means of 24hr concentrations of [NO.sub.2], [SO.sub.2], and [PM.sub.10] and 8-hr (900 hr-1700 hr) concentrations of [O.sub.3] were derived if they were non-missing. Daily concentrations were defined as non-missing if more than 17/24 hr concentrations and more than 5/8 hr concentrations were valid. According to the second phase guidelines of APHEA (Air Pollution and Health: a European Approach), non-missing daily means were first centered for each station i [i.e., individual daily concentrations ([X.sub.ij]) were subtracted by an annual station mean ([X.sub.i]) for each day j]. The centered data from all centers were then combined and added into the annual mean of all stations (X) to form [X'.sub.ij] = ([X.sub.ij] - [X.sub.i] + X). The daily (mean) concentrations of individual pollutants were computed for analysis by taking the mean of [X'.sub.ij] over all stations (24).
Statistical methods. We used Poisson regression with daily mortality counts as the dependent variable. To obtain a core model for each of the mortality outcomes for all ages, nonparametric smoothing (by means of the Loess function) terms for trend on days (1-1,096), seasonality, temperature, and humidity; and dummy variables for days of the week, holidays, and influenza epidemics [weeks with number of hospital admissions for influenza (ICD-9 487) in the upper quartile, which was on average over 8/week in 1995, 1996, and 1997, respectively] were fitted as the independent variables. In addition, we also considered the lag effects of temperature and humidity in building the core models. Residuals of each core model were examined to check whether there were discernible patterns and autocorrelation by means of residual plots and partial autocorrelation function plots, respectively (7). If necessary, both overdispersions and autocorrelations were further adjusted for the model using statistical procedures (7) implemented in S-Plus (MathSoft, Inc., Seattle, WA, USA). We paid special attention to ensure that there were no differences in the residuals between warm (April-September) and cool (other months) seasons.
We estimated concentrations of current day up to the previous 5 days for [O.sub.3] and up to 3 days for other pollutants and identified the best lagged day by a modified version of Akaike's Information Criterion (ACI) (25). The analysis was also performed using the Loess smoothing function to adjust for nonlinear effect of a copollutant. Differences in pollutant effects between seasons were assessed by a season-by-pollutant concentration interaction term in each model, and the effect estimates for cool and warm seasons were derived from the model with the interaction terms.
To perform the stratified analyses, we first obtained expected mortality counts ([Xi]) from the core model for all seasons. Poisson regression for the mortality outcomes (Y) was then fitted on pollutant concentrations (X) to obtain the log relative risk ([Beta]) estimate with offset on log([Xi]) (26) separately for warm and cool seasons. Offset is a computation procedure to treat log([Xi]) as a reference value and does not proceed to estimate a parameter for it in the Poisson regression log[E(Y)] = log([Xi]) + [Alpha] + [Beta]x (where [Alpha] is a parameter for the constant term). Exposure-response curves in warm and cool seasons were examined using generalized additive modeling (25).
Summary statistics. Summary statistics of mortality counts, air pollutant concentrations, and meteorologic measurements are presented in Table 1. There were more deaths, higher concentrations of pollutants (except for [SO.sub.2], which was about the same), and drier weather conditions in the cool season than in the warm season.
Table 1. Summary statistics of mortality outcome, air pollution levels, and meteorologic measures by season. No. (day) Mean SD Min Mortality counts Nonaccident (ICD-9: < 800) Warm 552 75.0 9.7 47 Cool 544 87.4 12.7 53 Cardiovascular (ICD-9: 390-456) Warm 552 19.5 4.7 8 Cool 544 26.2 6.4 12 Respiratory (ICD-9: 460-519) Warm 552 15.9 4.9 5 Cool 544 18.3 5.3 3 Air pollution concentrations ([micro]g/[m.sup.3]) [NO.sub.2] (24-hr) Warm 552 48.1 18.2 15.3 Cool 544 63.8 17.5 28.7 [S0.sub.2] (24-hr) Warm 550 18.3 13.0 1.9 Cool 544 17.2 11.6 1.1 [PM.sub.10] (24-hr) Warm 552 42.2 21.3 14.1 Cool 544 61.7 24.7 14.1 [0.sub.3] (8-hr) Warm 548 32.0 24.5 0 Cool 538 35.1 21.3 0 Meteorologic measurements: Temperature ([degrees] C) Warm 552 27.3 1.9 21.0 Cool 544 19.0 3.6 6.9 Humidity (%) Warm 552 80.7 7.4 46 Cool 544 74.7 12.4 31 [P.sub.10] Median [P.sub.90] Max Mortality counts Nonaccident (ICD-9: < 800) Warm 62 75 87 103 Cool 71 88 103 129 Cardiovascular (ICD-9: 390-456) Warm 14 19 26 35 Cool 18 26 35 53 Respiratory (ICD-9: 460-519) Warm 10 16 22 31 Cool 12 18 26 33 Air pollution concentrations ([micro]g/[m.sup.3]) [NO.sub.2] (24-hr) Warm 27.4 45.5 72.8 125.8 Cool 45.2 60.6 87.3 151.5 [S0.sub.2] (24-hr) Warm 5.9 15.0 35.3 83.6 Cool 6.4 14.4 30.8 90.1 [PM.sub.10] (24-hr) Warm 23.0 35.6 70.6 163.8 Cool 33.3 58.7 95.1 156.6 [0.sub.3] (8-hr) Warm 8.1 23.9 64.7 168.9 Cool 7.9 33.2 62.8 101.6 Meteorologic measurements: Temperature ([degrees] C) Warm 24.5 27.4 29.6 30.9 Cool 14.5 18.9 23.8 27.4 Humidity (%) Warm 73 80 91 97 Cool 58 76.5 89 95 Abbreviations: Max, maximum; Min, minimum; [P.sub.10], 10th percentile; [P.sub.90], 90th percentile.
Lag effects. In whole-year analysis, [NO.sub.2], [SO.sub.2], and [PM.sub.10] showed similar patterns for their effects on all the mortality outcomes in that the relative risks (RRs) increased from lag-day 0, were maximal at either lag-day 1 or lag-day 2, and declined to the lowest at lag-day 3. The RRs at the best lagged day (i.e., the day with minimum AIC) were significant in all three categories of deaths (p [is less than] 0.01) for [NO.sub.2] and [SO.sub.2], and was significant (p = 0.024) only for respiratory mortality for [PM.sub.10]. For [O.sub.3], the RRs were not significant (p [is greater than] 0.05) for any of the lagged days or for any mortality outcomes. With adjustment for autocorrelation, there was little change in the RRs and p-values. However, with adjustment for a copollutant, only the RR for [NO.sub.2] in cardiovascular mortality (p = 0.046) and for [SO.sub.2] in nonaccidental (p = 0.003) and cardiovascular (p = 0.023) mortality remained significant (Table 2).
Table 2. Relative risk (RR) and 95% confidence interval (Cl) of the best single lagged-day effects by linear extrapolation for a 10th-90th percentile change in pollutant concentration (1995-1997): whole year. Unadjusted Causes of mortality Lag day RR (95% Cl) p-Value [NO.sub.2] Nonaccident 1 1.04 (1.01-1.05) 0.001 Cardiovascular 2 1.06 (1.03-1.10) 0.001 Respiratory 0 1.08 (1.02-1.13) 0.003 [SO.sub.2] Nonaccident 1 1.03 (1.02-1.05) 0.000 Cardiovascular 1 1.05 (1.02-1.08) 0.001 Respiratory 0 1.04 (1.01-1.08) 0.010 [PM.sub.10] Nonaccident 1 1.02 (1.00-1.04) 0.102 Cardiovascular 2 1.03 (0.99-1.06) 0.165 Respiratory 1 1.06 (1.01-1.11) 0.024 [O.sub.3] Nonaccident 5 1.01 (0.99-1.03) 0.224 Cardiovascular 3 1.01 (0.98-1.05) 0.479 Respiratory 4 1.04 (1.00-1.08) 0.078 Autocorrelation adjusted Causes of mortality Lag day RR (95% CI) p-Value [NO.sub.2] Nonaccident 1 1.03 (1.01-1.05) 0.003 Cardiovascular 2 1.06 (1.02-1.10) 0.003 Respiratory 0 1.07 (1.02-1.12) 0.008 [SO.sub.2] Nonaccident 1 1.03 (1.01-1.05) 0.000 Cardiovascular 1 1.05 (1.02-1.08) 0.003 Respiratory 0 1.04 (1.01-1.07) 0.016 [PM.sub.10] Nonaccident 1 1.02 (1.00-1.04) 0.132 Cardiovascular 2 1.02 (0.99-1.06) 0.201 Respiratory 1 1.05 (1.01-1.10) 0.028 [O.sub.3] Nonaccident 5 1.01 (0.99-1.03) 0.226 Cardiovascular 3 1.01 (0.98-1.05) 0.426 Respiratory 4 1.03 (0.99-1.07) 0.145 Adjusted for copollutant Copol- Causes of mortality Lag day lutant(a) RR (95% Cl) p-Value [NO.sub.2] Nonaccident 1 [SO.sub.2] 1.00 (0.97-1.03) 0.896 Cardiovascular 2 [SO.sub.2] 1.04 (1.00-1.08) 0.046 Respiratory 0 [SO.sub.2] 1.05 (0.98-1.12) 0.168 [SO.sub.2] Nonaccident 1 [NO.sub.2] 1.03 (1.01-1.05) 0.003 Cardiovascular 1 [NO.sub.2] 1.04 (1.00-1.07) 0.023 Respiratory 0 [NO.sub.2] 1.02 (0.97-1.06) 0.450 [PM.sub.10] Nonaccident 1 [SO.sub.2] 0.99 (0.97-1.01) 0.397 Cardiovascular 2 [SO.sub.2] 0.98 (0.92-1.03) 0.363 Respiratory 1 [SO.sub.2] 1.04 (0.99-1.1O) 0.093 [O.sub.3] Nonaccident 5 [NO.sub.2] 1.01 (0.99-1.03) 0.288 Cardiovascular 3 [NO.sub.2] 1.00 (0.96-1.04) 0.997 Respiratory 4 [NO.sub.2] 1.03 (0.99-1.07) 0.163 (*) The copollutant that produced the least significant effect in the pollutant after adjustment.
Effects by seasons. During the warm season, there were no significant effects (p [is greater than] 0.1) for all pollutants for all mortality outcomes (Table 3). During the cool season, without adjustment for copollutants, a) all of the RRs at the best lagged days were significantly greater than unity in all the mortality outcomes (p [is less than] 0.015) for [NO.sub.2]; b) they were significant in nonaccidental and cardiovascular mortality (p [is less than] 0.002) for [SO.sub.2]; c) they were marginally significant in respiratory mortality (p = 0.054) for [PM.sub.10]; and d) they were all significant (p [is less than] 0.038) for [O.sub.3]. During the cool season with adjustment for copollutants, only the effects of [NO.sub.2] on cardiovascular mortality, [SO.sub.2] on nonaccident and cardiovascular mortality, and [O.sub.3] on nonaccident and respiratory mortality remained significant (p [is less than] 0.05). However, the between-season differences were statistically significant for [NO.sub.2] in cardiovascular mortality with and without adjustment for copollutants (p [is less than] 0.039) and for [O.sub.3] in all mortality outcomes (p [is less than] 0.044) without adjustment and in nonaccident and respiratory mortality (p [is less than] 0.032) with adjustment for copollutant.
Table 3. Relative risk (RR) and 95% confidence interval (Cl) of best single lagged day effects by linear extrapolation for a 10th-90th percentile change in pollutant concentration (1995-1997), without(a) and with(b) adjustment for a copollutant. Warm season Causes of mortality Copollutant RR (95% Cl) p-Value [NO.sub.2] Nonaccident -- 1.02 (0.99-1.05) 0.243 [SO.sub.2] 1.00 (0.97-1.04) 0.927 Cardiovascular -- 1.00 (0.94-1.06) 0.981 [SO.sub.2] 0.99 (0.94-1.05) 0.793 Respiratory -- 1.05 (0.99-1.13) 0.126 [SO.sub.2] 1.03 (0.95-1.12) 0.529 [SO.sub.2] Nonaccident -- 1.02 (0.99-1.04) 0.170 [NO.sub.2] 1.02 (0.99-1.04) 0.252 Cardiovascular -- 1.01 (0.97-1.05) 0.546 [NO.sub.2] 1.01 (0.97-1.06) 0.520 Respiratory -- 1.04 (0.99-1.09) 0.101 [NO.sub.2] 1.03 (0.97-1.08) 0.363 [PM.sub.10] Nonaccident -- 1.01 (0.98-1.04) 0.529 [SO.sub.2] 1.00 (0.96-1.03) 0.802 Cardiovascular -- 1.00 (0.94-1.06) 0.911 [NO.sub.2] 1.00 (0.91-1.10) 0.983 Respiratory -- 1.05 (0.98-1.12) 0.194 [NO.sub.2] 1.04 (0.96-1.12) 0.379 [O.sub.3] Nonaccident -- 0.99 (0.97-1.02) 0.609 [NO.sub.2] 0.99 (0.97-1.02) 0.537 Cardiovascular -- 0.98 (0.94-1.03) 0.485 [NO.sub.2] 0.98 (0.94-1.03) 0.513 Respiratory -- 0.99 (0.94-1.05) 0.750 [NO.sub.2] 0.99 (0.94-1.04) 0.710 Between Cool season season Causes of mortality Copollutant RR (95% Cl) p-Value p-value [NO.sub.2] Nonaccident -- 1.05 (1.02-1.08) 0.003 0.193 [SO.sub.2] 1.01 (0.97-1.05) 0.694 0.795 Cardiovascular -- 1.10 (1.05-1.16) 0.000 0.013 [SO.sub.2] 1.08 (1.02-1.14) 0.007 0.039 Respiratory -- 1.09 (1.02-1.16) 0.015 0.509 [SO.sub.2] 1.08 (0.98-1.19) 0.120 0.408 [SO.sub.2] Nonaccident -- 1.04 (1.02-1.07) 0.001 0.101 [NO.sub.2] 1.04 (1.00-1.07) 0.030 0.292 Cardiovascular -- 1.07 (1.02-1.11) 0.002 0.070 [NO.sub.2] 1.05 (1.00-1.09) 0.045 0.310 Respiratory -- 1.04 (1.00-1.09) 0.079 0.877 [NO.sub.2] 1.01 (0.94-1.08) 0.890 0.625 [PM.sub.10] Nonaccident -- 1.02 (0.99-1.05) 0.168 0.659 [SO.sub.2] 0.99 (0.96-1.02) 0.437 0.715 Cardiovascular -- 1.04 (0.99-1.09) 0.135 0.306 [NO.sub.2] 0.97 (0.91-1.03) 0.349 0.614 Respiratory -- 1.06 (1.00-1.13) 0.054 0.761 [NO.sub.2] 1.05 (0.99-1.12) 0.139 0.810 [O.sub.3] Nonaccident -- 1.04 (1.01-1.06) 0.012 0.026 [NO.sub.2] 1.03 (1.01-1.06) 0.021 0.032 Cardiovascular -- 1.05 (1.00-1.11) 0.038 0.044 [NO.sub.2] 1.04 (0.99-1.09) 0.150 0.132 Respiratory -- 1.08 (1.02-1.15) 0.011 0.027 [NO.sub.2] 1.08 (1.02-1.15) 0.013 0.030 Warm season, April-September; cool season, October-March. (a)<1> Estimated from core model + pollutant + season + pollutant x season. (b)Estimated from <1> + copollutant + copollutant x season.
Seasonal exposure-response relationships. Figures 1-4 show the exposure-response relationships for each pollutant for the three mortality outcomes at the best lagged day.
[NO.sub.2]. During the warm season, we observed no clear exposure-response relationships for the three outcomes for [NO.sub.2]. However, during the cool season there were observable linear exposure-response relationships throughout the concentration levels in nonaccidental mortality, but the curves showed positive and nonlinear relationships at concentrations higher than 80 [micro]g/[m.sup.3] in the other two outcomes.
[SO.sub.2]. During the warm season, no exposure-response relationships were observed when [SO.sub.2] was [is less than] 30 [micro]g/[m.sup.3], but there were some linear or nonlinear relationships above that concentration. During the cool season, we observed positive exposure-response relationships for concentrations of 0-40 [micro]g/[m.sup.3] [SO.sup.2].
[PM.sub10]. For [PM.sub.10], no clear exposure-response relationships were observed for the three outcomes in warm seasons, but in the cool season there was a positive exposure-response relationship for respiratory mortality for concentrations up to 80 [micro]g/[m.sup.3].
[O.sub.3]. There were no dear relationships for any of the three outcomes for [O.sub.3] during the warm season. However, during the cool season all of the mortality outcomes tended to increase with increasing concentrations.
All pollutant levels are high in Hong Kong. Although [SO.sub.2] has been reduced substantially due to government limits on the sulfur content of fuels in the early 1990s (2,7), the level of [SO.sub.2] in Hong Kong still ranks in the middle among more than 30 metropolitan cities in the world. The [SO.sub.2] level in Hong Kong is higher than those in Berlin, Germany; Boston, Massachusetts (USA); Brisbane, Australia; Kuala Lumpur, Malaysia; London, United Kingdom; and Paris, France (28).
The levels of [NO.sub.2] and [O.sub.3] have been increasing along with increasing vehicular traffic volume. Levels of [PM.sub.10], which is primarily related to the use of diesel engines, in Hong Kong are among the highest in the world: they are only lower than those in the most polluted cities such as Barcelona, Spain; Guangzhou, China; Manila, Republic of the Philippines; Mexico City, Mexico; Philadelphia, Pennsylvania (USA); Santiago, Chile; Shanghai, China; and Taipei, Taiwan.
In the present study the estimated effects of the pollutants on mortality reached a maximum at a lag of 1-2 days. These observations are consistent with those reported by Bremner et al. (29) in London: the effects increased from lag-day 0 to a maximum at lag-day 1 for [NO.sub.2], [SO.sub.2], and [PM.sub.10].
A major finding of this study is that [O.sub.3] had effects on all three mortality outcomes during the cool season, and the effects were greater than those in the warm season; this is unlike several other reports in which the effects were found in the warm season (9,12,29). This is consistent with our previous report on the effects of pollution on hospital admissions due to heart failure in subjects [is greater than or equal to] 65 years of age (30). The effects of the other oxidant pollutants ([NO.sub.2] and [SO.sub.2]) were also significant for all of the mortality outcomes in the cool season but not in the warm season. In Athens, Greece, effects of [SO.sub.2] on all causes of nonaccidental mortality were also observed in the cool season (31), but in London, the effects for [NO.sub.2] and [SO.sub.2] were observed in the warm season (12). When the data from five western European cities and four central European cities were combined, [SO.sub.2] also showed slightly stronger effects during the warm season than during the cool season (32).
In Hong Kong in the cool season, air pollutant levels were higher ([NO.sub.2], 64 vs. 48; [PM.sub.10], 62 vs. 42; [O.sub.3], 35 vs. 32 [micro]g/[m.sup.3]) than those in the warm season, except [SO.sub.2], which was slightly lower (17 vs. 18). Because pollutants were correlated (r = 0.54-0.72 between [NO.sub.2], [SO.sub.2], and [PM.sub.10] during the cool season), greater effects observed during cool weather may be due to other pollutants that were also at higher levels during the cool season. The cool season in Hong Kong is drier (humidity 75% vs. 81%), less cloudy (63% vs. 72%), and less variable, so people are more likely to go outdoors and open the windows, thus being exposed to higher levels of air pollution. In contrast, during the warm season (temperatures of 25 [degrees] C-30 [degrees] C and humidity of 73o/o-91% between 10th to 90th percentiles) people usually use air-conditioning, thus reducing the risks of outdoor ambient air pollution exposure.
Another major finding in this study is the positive exposure-response relationships for [NO.sub.2] and [SO.sub.2] and all the outcomes during the cool season. There were no thresholds, and the effects showed an inverted "J" shape at higher concentrations. At very high concentrations, the risks of mortality could be reduced possibly because vulnerable subjects may have died before the concentration had reached the maximum levels (4). During the warm season, we observed no consistent positive or negative relationships for all the pollutants. In Hang Kong, there are greater variations in weather conditions in the warm season, when heavy rain, rain storms, and typhoons are common. These factors, in addition to the frequent use of air-conditioning, would prevent the actual exposure-response relationships between air pollution and mortality from being readily observable.
In the absence of an observed linear exposure-response relationship, generalized additive modelling (GAM) could be used to examine whether there are any other forms of relationships. Instead of obtaining a single parameter for the effect, GAM is fitted to obtain a parameter at each point of the independent variable after applying some smoothing function to the data. The fitted values (presented as deviation from an overall mean), along with values of the independent variable, produce an exposure-response plot. It is useful, as demonstrated in this study, in the interpretation of results of daily time-series studies for health effects of air pollution.
Morris and Naumova (33) reported synergistic effects of carbon monoxide and lower temperatures on hospital admissions due to congestive heart failure in Chicago, Illinois (USA). Both CO and cold temperature can increase the load on the heart and thus increase the effect on cardiovascular morbidity (33). For other pollutants, including [SO.sub.2], the production of synergistic effects was biologically plausible, as both lower temperatures and high air pollutant concentrations were related to increased blood viscosity. Changes in blood theology may be caused by an inflammatory process in the lung induced by air pollution or by thermoregulatory adjustment to mild surface cooling in cold weather (34,35). This study in Hang Kong is the first to show that all of the oxidant pollutants under study increased effects (p [is less than] 0.07) on cardiovascular mortality in the cool season. It is also the first study to demonstrate the relationship between pollutant concentration and mortality stratified by cool and warm seasons, on the basis of statistical models with offset on expected counts from the same core model, thus ensuring comparability in the effect estimates between the two seasons. Overall, during the cool season there was a 5-10% (p [is less than] 0.038) increase in nonaccidental and cardiovascular mortality; this results in an increase from the 10th to the 90th percentile (from linear extrapolation) for each of the oxidant pollutants under study. The nonsignificant relative risk estimate for [SO.sub.2] on respiratory disease may be due to the small change in concentration from the warm season to the cool season. In a sensitivity analysis using the method with offset on expected counts, the estimated increases were consistent but lower, with increases of 2-7% (p [is less than] 0.046) (data not shown). The larger p-values may be due to lower power in stratified analysis when the sample size was halved.
Except for respiratory mortality, no strong effects of particulate pollutants were observed in Hong Kong, unlike in other places, although the levels were high. This should be investigated further. The difference may arise from the use of a time-series study, in that the magnitude of the effect estimates depends on the day-to-day covariation of the daily health outcomes and pollutant concentrations instead of the absolute levels of the pollutant concentrations. However, [PM.sub.10] was found to have a significant effect on respiratory mortality (RR = 1.05; p = 0.028) in both seasons combined, but the effects remain approximately the same (RR = 1.05-1.06) although nonsignificant (p [is greater than] 0.054) in the by-season estimates. The importance of [PM.sub.10] should not be diminished by this finding. [NO.sub.2] is important because of increasing volumes of vehicular traffic on the roads. [SO.sub.2] continues to have a strong effect, even though the concentrations have decreased and have been maintained at low levels. The formation of [O.sub.3] in the ambient air depends on a series of complicated photochemical reactions of oxygen, nitrogen oxides, and reactive hydrocarbons in the presence of sunlight. [O.sub.3] had been increasing until recently, and it is difficult and costly to control as a regional pollutant. HoWever, for most of the pollutants, the effects were nonsignificant after adjusting for a copollutant; this may arise from a problem of multicollinearity, except for [NO.sub.2] in cardiovascular mortality and [SO.sub.2] in both nonaccidental and cardiovascular mortality. These observations, together with strong positive exposure-response relationships for [NO.sub.2] and [SO.sub.2], suggested that [NO.sub.2] and [SO.sub.2] have independent effects and may be better indicators of effects on mortality in this subtropical city. For [O.sub.3] all of the RR estimates were not significantly greater than unity; the RR estimates were insensitive to adjustment either for autocorrelation or copollutant (Tables 2 and 3) and were insensitive to the use of maximum, minimum, or mean temperature in the model (data not shown).
In setting air pollution control policy from a public health viewpoint, it is important to identify the health effects of air pollutants from local data. Because of the lack of data, there are few studies based on daily hospital admissions and mortality in the Asian Pacific region. For hospital admissions, there has been only one study in Australia (36) and two in Hong Kong (30,37). For mortality studies, there have been one in Beijing, China (38) based on 1-year daily data, two in Australia (36,39), and two in Korea (40,41). Our report should contribute to the understanding of the effects of air pollutants in this region and may clarify the differences in effects and mechanisms between Western and Eastern populations.
Local data on health effects of air pollution are required for setting standards and objectives for air pollution controls. When local data are not available, foreign data may be helpful, but they may not be relevant or applicable because of a difference in climate or other conditions. Our findings in this study provide information to support a review of air quality objectives with consideration of their effects on health (10).
REFERENCES AND NOTES
(1.) Thurston GD. A critical review of [PM.sub.10] mortality time-series studies. J Expo Anal Environ Epidemiol 6:3-21 (1996).
(2.) Morgenstern H. Use of ecologic analysis in epidemiologic research. Am J Public Health 139:1336-1344 (1982).
(3.) Schimmel H, Murawski TJ. Evidence for possible acute health effects of ambient air pollution from time series analysis: methodological questions and some new results based on New York City daily mortality 1963-1976. Bull NY Acad Med 54:1052-1058 (1978).
(4.) Spix C. Daily time series of mortality counts: estimating the harvesting effect. Stat Mad (in press).
(5.) Zeger SL, Dominici F, Samet J. Harvesting-resistant estimates of air pollution effects on mortality. Epidemiology 10:171-175 (1999).
(6.) Rossi G, Zanobetti A, Marchi M. Time series analysis in environmental epidemiology: short-term effects of air pollution on mortality and morbidity. Epidemiol Prey 19:90-98 (1995).
(7.) Schwartz J, Spix C, Touloumi G, Bacharova L, Barumamdzadeh T, le Tertre A, Piekarksi T, Ponce de Leon A, Ponka A, Rossi G, et al. Methodological issues in studies of air pollution and daily counts of deaths or hospital admissions. J Epidemiol Community Health 50(suppl 1):S3-S11 (1996).
(8.) McMichael AJ, Anderson HR, Brunekreef B, Cohen AJ. Inappropriate use of daily mortality analyses to estimate long-term mortality effects of air pollution. Int J Epidemiol 27:450-453 (1998).
(9.) Ostro B, Sanchez JM, Aranda C, Eskeland GS. Air pollution and mortality: results from a study of Santiago, Chile. J Expo Anal Environ Epidamiol 6:97-114 (1996).
(10.) Hadley AJ, Wong CM, Wong TW, Tam CM, Yu ITS, Tam E. Health Effects of Air Pollution. Report from the SubWorking Group on the Review of Hong Kong's Air Quality Objectives. Hong Kong:Environmental Protection Department, Hong Kong Government, 1999.
(11.) Kinney PL, Ozkaynak H. Associations of daily mortality and air pollution in Los Angeles County. Environ Res 54:99-120 (1991).
(12.) Anderson HR, Ponce de Leon A, Bland JM, Bower JS, Strachan DP. Air pollution and daily mortality in London: 1987-92. Br Med J 312:665-669 (1996).
(13.) Sunyer J, Castellsague J, Saez M, Tobias A, Anto JM. Air pollution and mortality in Barcelona. J Epidemiol Community Health 50(suppl 1):S76-80 (1996).
(14.) Wietlisbach V, Pope CA, Ackermann-Liebrich U. Air pollution and daily mortality in three Swiss urban areas. Soz Praventivmed 41:107-115 (1996).
(15.) Touloumi G, Katsouyanni K, Zmirou D, Schwartz J, Spix C, Ponce de Leon A, Tobias A, Quennel P, Rabczenko D, Bacharova L, et al. Short-term effects of ambient oxidant exposure on mortality: a combined analysis within the APHEA project. Am J Epidemiol 146:177-185 (1997).
(16.) Schwartz J. Air pollution and daily mortality: a review and meta-analysis. Environ Rev 64:36-52 (1994).
(17.) Katsouyanni K, Touloumi G, Spix C, Schwartz J, Balducci F, Medina S, Rossi G, Wojtyniak B, Sunyer J, Bacharova L, et al. Short term effects of ambient sulphur dioxide and particulate matter on mortality in 12 European cities: results from time series data from the APHEA project. Air Pollution and Health: a European Approach. Br Meal J 314:1658-1663 (1997).
(18.) Dockery DW, Pope AC III, Xu X, Spengler JD, Ware JH, Fay ME, Ferris BG Jr, Speizer FE. An association between air pollution and mortality in six US cities. N Engl J Mad 329:1753-1759 (1993).
(19.) Thurston GD, Ito K, Lippmann M, Hayes C. Re-examination of London, England, mortality in relation to exposure to acidic aerosols during 1963-1972 winters. Environ Health Perspect 79:73-82 (1989).
(20.) Dockery DW, Schwartz J, Spengler JD. Air pollution and daily mortality: associations with particulates and acid aerosols. Environ Res 59:362-373 (1992).
(21.) Ito K, Thurston GD, Hayes C, Lippmann M. Associations of London, England, daily mortality with particulate matters, sulphur dioxide and acidic aerosol pollution. Arch Environ Health 48:21 3-220 (1993).
(22.) Borja-Aburto VH, Loomis DP, Bangdiwala SI, Shy CM, Rascon-Pacheco RA. Ozone, suspended particulates and daily mortality in Mexico City. Am J Epidemiol 145:258-268 (1997).
(23.) WHO. International Classification of Diseases, 1975 Revision. Geneva:World Health Organization, 1977.
(24.) Atkinson R. Personal communication.
(25.) Hastie TJ, Tibshirani RJ. Generalized Additive Models. London:Chapman and Hall, 1990.
(26.) Breslow NE, Day NR. Statistical Methods in Cancer Research, Vol II. The Design and Analysis of Cohort Studies. IARC Sci Pub182:151-152 (1987).
(27.) Planning, Environment and Lands Bureau. Clean Air for Hong Kong. Hong Kong:Hong Kong Special Administrative Region Government, 1999. Available: http://www.pelb.gov.hk/cleanair/index.htm [cited February 2001].
(28.) Kiely P, Yap D, DeBrou GD, Fraser D, Dong W, Brown A. A Comparative Study of Toronto's and Selected World Cities. Presented at the 90th Air and Waste Management Association, Toronto, Ontario, June 1997.
(29.) Bremner SA, Anderson HR, Atkinson RW, McMichael AJ, Strachan DP, Bland JM, Bower-JS. Short term associations between outdoor air pollution and mortality in London 1992-94. Occup Environ Med 56:237-244 (1999).
(30.) Wong CM, Ma S, Hedley AJ, Lam TH. Does ozone have any effect on daily hospital admissions for circulatory diseases? J Epidemiol Community Health 53:580-581 (1999).
(31.) Touloumi G, Samoli E, Katsouyanni K. Daily mortality and air pollution from particulate matter, sulphur dioxide, and carbon monoxide in Athens, Greece, 1987-91. A time-series analysis within the APHEA project. J Epidemiol Community Health 50 (suppl 1):S47-S51 (1996).
(32.) Zmirou D, Schwartz J, Saez M, Zanobetti A, Wojtyniak B, Touloumi G, Spix C, Ponce de Leon A, Le Moullec Y, Bacharova L, et al. Time-series analysis of air pollution and cause-specific mortality. Epidemiology 9:495-503 (1998).
(33.) Morris RD, Naumova EN. Carbon monoxide and hospital admissions for congestive heart failure: evidence of an increased effect at low temperatures. Environ Health Perspect 106:649-653 (1998).
(34.) Peters A, Doring A, Wichmann HE, Koenig W. Increased plasma viscosity during an air pollution episode: a link to mortality? Lancet 349:1582-1587 (1997).
(35.) Keatinge WR, Coleshaw SRK, Cotter F, Mattock M, Murphy M, Chelliah R. Increases in platelet and red cell counts, blood viscosity, and arterial pressure during mild surface cooling: factors in mortality from coronary and cerebral thrombosis in winter. Br Mad J 289:1405-1406 (1964).
(36.) Morgan G, Corbett S, Wlodarczyk J, Lewis P. Air pollution and daily mortality in Sydney, Australia, 1989 through 1993. Am J Public Health 88:759-764 (1998).
(37.) Wong TW, Lau TS, Yu TS, Neller A, Wong SL, Tam W, Pang SW. Air pollution and hospital admissions for respiratory and cardiovascular diseases in Hong Kong. Occup Environ Med 56:679-663 (1999).
(38.) Xu X, Gao J, Dockery DW, Chen Y. Air pollution and daily mortality in residential areas of Belling, China. Arch Environ Health 49:21 6-222 (1994).
(39.) Simpson RW, Williams G, Petroeschevsky A, Morgan G, Rutherford S. Associations between outdoor air pollution and daily mortality in Brisbane, Australia. Arch Environ Health 52:442-454 (1997).
(40.) Lee J-T, Shin D, Chung Y. Air pollution and daily mortality in Seoul and Ulsan, Korea. Environ Health Perspect 107:149-154 (1999).
(41.) Hong Y-C, Leem J-H, Ha E-H, Christiani DC. [PM.sub10] exposure, gaseous pollutants, and daily mortality in Inchon, South Korea. Environ Health Perspect 107:873-878 (1999).
Chit-Ming Wong, Stefan Ma, Anthony Johnson Hedley, and Tai-Hing Lam
Department of Community Medicine, The University of Hong Kong, Hong Kong, China
Address correspondence to A.J. Hedley, Department of Community Medicine, The University of Hong Kong, Patrick Manson Building South Wing, 7 Sassoon Road, Hong Kong, China. Telephone: (852) 2819 9282 / 2819 9280. Fax: (852) 2855 9528. E-mail: email@example.com
We thank R. Atkinson and R. Anderson for their advice on data analysis; M. Chi for secretarial support, and P. Chau for data analysis support.
Received 12 May 2000; accepted 24 October 2000.