Copyright ©ERS Journals Ltd 2006 Short-term effects of nitrogen dioxide on mortality: an analysis within the APHEA project Fikfak111 Dept of Hygiene and Epidemiology, University of Athens, Athens, Greece. 2 Harvard School of Public Health, Boston, MA, USA. 3 Dept of Public Health and Clinical Medicine, Umeå University, Umeå, Sweden. 4 Regional Observatory of Health, Ile-de-France, France. 5 Dept of Public Health and Clinical Nutrition, University of Kuopio, and 6 National Public Health Institute, Unit of Environmental Epidemiology, Kuopio, Finland. 7 National Institute of Hygiene, Dept of Medical Statistics, Population Studies Laboratory, Warsaw, Poland. 8 Institute of Social and Preventive Medicine, University of Basel, Basel, Switzerland. 9 Dept of Environmental Health, Institute of Public Health Bucharest, Bucharest, Romania. 10 Center for Quality & Safety, Hadassah University Hospital, Jerusalem, Israel. 11 Institute of Public Health, Ljubljana, Slovenia. CORRESPONDENCE: E. Samoli, Dept of Hygiene and Epidemiology, University of Athens Medical School, 75 Mikras Asias Street, 115 27 Athens, Greece, Fax: 30 2107462205. E-mail: esamoli{at}med.uoa.gr Keywords: Air pollution, heterogeneity, modelling, mortality, nitrogen dioxide
Received: December 7, 2005
The short-term effects of nitrogen dioxide (NO2) on total, cardiovascular and respiratory mortality in 30 European cities participating in the Air Pollution on Health: a European Approach (APHEA)-2 project were investigated. The association was examined using hierarchical models implemented in two stages. In the first stage, data from each city were analysed separately, whereas in the second stage, the city-specific air pollution estimates were regressed on city-specific covariates to obtain overall estimates and to explore sources of possible heterogeneity. A significant association of NO2 with total, cardiovascular and respiratory mortality was found, with stronger effects on cause-specific mortality. There was evidence of confounding in respiratory mortality with black smoke and sulphur dioxide. The effect of NO2 on total and cardiovascular mortality was observed mainly in western and southern European cities, and was larger when smoking prevalence was lower and household gas consumption was higher. The effect of NO2 on respiratory mortality was higher in cities with a larger proportion of elderly persons in the population and higher levels of particulate matter with a 50% cut-off aerodynamic diameter of 10 µm. The results of this large study are consistent with an independent effect of nitrogen dioxide on mortality, but the role of nitrogen dioxide as a surrogate of other unmeasured pollutants cannot be completely ruled out. Many epidemiological studies have documented adverse short-term effects of different types of air pollution on health outcomes in recent years 15. The pollution indicators used were mainly ambient particles 4, 6, 7, but gaseous pollutants, such as nitrogen dioxide (NO2), ozone (O3) and carbon monoxide (CO), have also been shown to have adverse effects on mortality and morbidity 3, 810. The results from epidemiological studies have led to revisions of air quality guidelines and standards and scheduled dates for regular revisions in the future 1113. Nitrogen oxides (NOx), primarily nitric oxide (NO), are produced from high-temperature combustion, such as when fuel is burned in motor vehicles and power plants. Once emitted, NO relatively rapidly reacts with O3 or oxygen to form NO2. This usually dominating part is known as secondary NO2, but the transformation mainly takes place close to the source. Primary NO2 emissions are particularly important from diesel vehicles, and this fraction seems to be increasing in Europe. Due to nonlinearities in the oxidation reactions and variations in the background O3 levels, NOx has a stronger temporal correlation than NO2 has with other combustion products (exhaust particles, CO, sulphur dioxide (SO2)) emitted at the same time from cars and other sources. Indoor air can also be contaminated with high levels of NO2, since unventilated heaters and gas stoves also emit substantial amounts of NO2.
In spite of laboratory, clinical and epidemiological research, the health effects of NO2 exposure on humans are not well understood. The toxicological evidence suggests that increased susceptibility to infection, functional deficits from effects on airways, and deterioration of the status of persons with chronic respiratory conditions, including asthmatics, are of potential concern. NO2 is a highly reactive, nitrogen-centred free radical, poorly water-soluble gas deposited peripherally in the lungs. It is absorbed along the entire respiratory tract, but exposure studies indicate that the major target site for the action of NO2 is the terminal bronchioles. The main mechanism of NO2 toxicity has been suggested to involve lipid peroxidation in cell membranes and various actions of free radicals on structural and functional molecules. Antioxidants such as ascorbic acid and Most epidemiological studies on the health effects of NO2 have focused on morbidity rather than mortality. The short-term effects of NO2 on mortality within the Air Pollution and Health: A European Approach (APHEA)-2 project, which uses an extensive European database from 30 European cities, were investigated. Special attention was paid to efforts to distinguish the effects of NO2 per se from confounding or modifying effects of other pollutants, such as particles.
APHEA-2 is a multicentre project including 30 cities across Europe and associated regions, which studies short-term health effects of air pollution. Data were collected on daily counts of all-cause mortality, excluding deaths from external causes (International Classification of Disease (ICD)-9: >800), cardiovascular mortality (ICD-9: 390459) and respiratory mortality (ICD-9: 460519). The data covered at 3 consecutive years for each city within the years 19901997. In all, 2,893,430 deaths occurred in the cities studied during that period. Details of the data have been published elsewhere 5. Daily air pollution measurements were provided by the monitoring networks established in each town participating in the APHEA-2 project 5. Time series data on daily temperature (°C, daily mean) and relative humidity (%) were used to control for the potential confounding effects of weather. External information on influenza epidemics or other unusual events (heat waves, strikes, etc.) was also collected, if available 5.
Table 1
For the four Polish cities included in the analysis (i.e. Cracow, Poznan, Lodz, Wroclaw), only NO2 cumulative 24-h measurements were available. In order to include them in the analysis, the maximum hourly NO2 concentration of each day (NO2 1-h values) was estimated as 1.64 times the 24-h values, where 1.64 is the average of the ratio between the two measurements in the cities that provided both. There was a substantial variability among all cities in the study of the levels of all pollutants, as well as in the mean daily temperature and humidity.
Within the APHEA-2 project the maximum hourly NO2 concentration of each day (NO2 1-h) was analysed rather than the average concentration of NO2 over 24 h (NO2 24-h), since more cities provided measurements for the former. The maximum daily 1- and 24-h NO2 concentrations are highly correlated (the correlation coefficient ranges 0.800.94, with a median of 0.90). In the present study, it was decided a priori to use the average of lags 0 and 1 for NO2 for all cities, since there is evidence that the average of 2 days pollution correlates better with mortality than a single days exposure 19. Furthermore, this approach avoids potential bias that could result from selectively reporting the most significant lags. In this analysis, a linear doseresponse relationship between NO2 1-h and mortality was assumed. This assumption was based on previously published results from the APHEA-2 project indicating that the doseresponse relationship could be adequately approximated by a linear association (E. Samoli; personal communication) 20. To investigate the effect of NO2 over a larger number of days and examine the shape of the association with each analysed health outcome, polynomial distributed lag models for the NO2 were fitted, using lags 0 to 5 19. For the analysis, a hierarchical modelling approach was used. First, regression models in each city were fitted separately to allow specific control for seasonal effects, weather and other potential confounders. The individual city results were then used in a second-stage analysis to obtain overall estimates and to investigate potential effect modifiers.
Individual city analysis The general methodological guidelines developed within the framework of the APHEA-2 project were followed, and are described in detail elsewhere 21. One additional feature is the use of penalised regression splines instead of the nonparametric function locally weighted nonparametric smoothers (loess) as smoothing functions to control for possible confounding. The smooth functions of time serve as a proxy for any time-dependent outcome predictors or confounders with long-term trends and seasonal patterns not explicitly included in the model. Hence long-term trends and seasonal patterns were removed from the data to guard against this confounding by omitted variables. Weather variables, which are potential confounders, were included explicitly. In particular, same-day temperature and humidity, and lagged values of these meteorological variables were included in the models. Thin-plate regression splines were used as basis functions for the penalised regression splines 23. Based on experiences from the previous analyses of the APHEA-2 data, the number of basis functions was chosen to be 40 for the time variable and 10 for the weather variables. The smoothing parameters that minimised the absolute value of the sum of partial autocorrelations (PACs) of the residuals from lags 3 to 30 were then chosen. To account for serial correlation in the residuals where it remained in the final model, autoregressive terms were added into the model as appropriate 24. In the special case of the small cities (and especially in cause-specific mortality), a minimum of one degree of freedom per year was required. The APHEA-2 method was used for influenza control, including a dummy variable taking the value of one when the 7-day moving average of the respiratory mortality was greater than the 90th percentile of its city-specific distribution. Since the influenza control, as described, was based on the distribution of respiratory mortality, the influenza dummy variable was included only when the total and cardiovascular mortality were analysed 21. Based on previously published results 25, 26, there is no indication that omitting control for influenza when respiratory mortality was analysed would influence the association between air pollution and mortality. It is unclear why the specific time within a winter that an epidemic occurs in a particular city should have much to do with air pollution levels and hence confound the relation under investigation. Nevertheless, to further explore the potential confounding effects of influenza on respiratory mortality, respiratory mortality was also analysed for days below the 90th percentile of its distribution. To evaluate how sensitive the present results are to the choice of the degree of smoothing, models with ±25% of the degrees of freedom for the time smoothing that were chosen based on the PAC criterion were also applied. To investigate potential confounding effects of the daily levels of other pollutants, two pollutant models were also applied, i.e. NO2 and alternatively PM10 (24-h mean), SO2 (24-h mean) or O3 (maximum of 8-h means) were included in the model, and were also applied. Some studies have demonstrated that the effects of air pollution on mortality are spread over >2 days 19, 27. To examine this question, a separate distributed lag model was fitted for each city. A cubic polynomial distributed lag with lags to 5 days before the deaths was used, which has proven adequate in past studies 19. Therefore, the coefficients for the cubic polynomial that define the shape of the distributed lag were estimated in each city.
Second-stage analysis
Table 2
Figure 1
Figure 2
Figure 3
Figure 4
Table 3
The observed heterogeneity in the effect estimates of NO2 was examined using potential-effect modifiers in second-stage regression models. Potential-effect modifiers used in the APHEA-2 analysis included variables describing the air pollution level and mix in each city, the health status of the population, the geographical area and the climatic conditions 5. Table 4
When investigating the source of heterogeneity in the association between NO2 and total and cardiovascular mortality, the most important effect modifier was the geographical area (defined as western, southern and central eastern European cities), followed by the prevalence of smoking in the city. More specifically, in cities with a lower prevalence of smoking the effect of NO2 on total and cardiovascular mortality was greater. The highest effect of NO2 on total and cardiovascular disease mortality was in north-western cities, followed closely by the effect in southern European cities, while there was only a small and nonsignificant effect in eastern cities. Concerning cardiovascular mortality, there was additional evidence that household consumption of natural gas acted as an effect modifier, with higher NO2 effect where the consumption was higher. Finally, in the association with respiratory mortality the most important effect modifiers were the median levels of PM10 followed by the proportion of the elderly (i.e. >65 yrs of age) in the population of the city. In cities with high median PM10 levels and a high proportion of elderly people, there was a stronger effect of NO2 on respiratory mortality.
In studies published during the recent decades, NO2 has been associated with a decrease in lung function, an increase in respiratory symptoms and an increase in asthma and COPD hospital admissions. Most of the epidemiological studies of short-term effects of NO2 on health were focused on symptoms reported in diaries, on hospitalisation for respiratory diseases, and on decrease of pulmonary function 3033. A few studies have been conducted on the effects of photochemical air pollution on mortality 3437. The most extensive European database available today to investigate the potential effects of NO2 on mortality was used. Significant adverse health effects of NO2 on total, cardiovascular and respiratory mortality, with stronger effects on cause-specific mortality have been found. These findings complement those previously reported from APHEA-1, which was the first part of the APHEA project and included a smaller number of cites (i.e. only six compared with 30 cities included in the APHEA-2 analysis). As part of that analysis, Touloumi et al. 3 have reported significant positive associations between NO2 and daily total deaths. In that first part of the APHEA analysis, an increase per 10 µg·m-3 in NO2 was associated with 0.26% increase in total mortality compared with 0.30% increase reported in the present results. Kinney and Ozkaynak 34 also found a significant association of NO2 with total and cardiovascular mortality in Los Angeles County, USA. However, no significant association was observed for respiratory deaths. The authors suggested that the small number of deaths from respiratory causes may have limited the power to detect small pollution effects. In the present analysis, where power is gained by use of multiple locations, there is a stronger effect on respiratory mortality compared with total mortality. However, the associated standard error is larger and the effect decreases and becomes nonsignificant when controlling for BS and SO2. Finally, in the analysis of 20 US cities within the National Morbidity, Mortality, and Air Pollution Study (NMMAPS) project, no consistent pattern of association between total mortality and NO2 was found 38. The investigators found a positive, but not statistically significant, effect of NO2 at lag 0 and at lag 1, and a highly statistically significant result at lag 2. The difference between the NMMAPS and APHEA findings may be attributed to the varying air pollution sources and mixture in Europe and the US. In this study, the cumulative effect over 6 days was larger by about 22% for total and cardiovascular mortality and by 45% for respiratory mortality compared with the average exposure over 2 days. This indicates that previous estimates of NO2 effects may in fact represent an underestimation if they take into account only very short-term health effects. Moreover, when the shape of the association between mortality and NO2 is considered, two different patterns can be distinguished: total and cardiovascular mortality is clearly more affected by NO2 levels on the 2 previous days, whereas for respiratory mortality the effects are more evenly distributed over the 6 previous days. This difference may be explained by differences in the biological mechanisms underlying the health effects of NO2. Although the above results indicate adverse effects of NO2 on mortality, the independence of its effect from those of other pollutants is still unclear. A single pollutant could act as a marker of a pollution mixture. Hence, NO2 could be a marker of other pollutants generated by vehicle exhausts such as particles. An attempt has been made to estimate the independent effects of pollutants using two pollutant models, for those pollutants for which data were available. No evidence of confounding for total and cardiovascular mortality has been found. There was evidence of confounding by BS and SO2 in the effect on respiratory mortality. A possible explanation is that the BS- and SO2-adjusted effect estimates of NO2 on respiratory mortality may reflect, to a larger extent, the effects from sources other than traffic. BS is more specific for traffic-related particles than PM10 and provides a means of addressing the question of particle composition. It was impossible to control for indices for which no measurements were available, for example, the number of particles. Further study focusing on exposure to mixtures including NO2, BS and SO2 is needed to further understanding of the aetiological mechanism through which the pollutants affect mortality and, in particular, respiratory mortality. In order to further investigate the independence of the NO2 effects and contribute to the ongoing discussion on the possible confounding effect of smaller particle sizes, the PM2.5 median level for nine of the analysed cities as potential effect modifier has been used 39. The PM2.5 median level did not act as significant effect modifier and reduced the heterogeneity by <10% in each of the studied outcomes. Besides, if NO2 truly reflects PM2.5, then in cities with a higher PM2.5/PM10 ratio, a higher NO2/PM10 level would be expected. This scenario was not verified in four cities with relevant data (Athens, Birmingham, Helsinki, Amsterdam; the Relationship between Ultrafine and fine Particulate matter in Indoor and Outdoor air and respiratory Health (RUPIOH) project, unpublished data). Finally, the correlation between PM10 and PM2.5 in these four cities ranged 0.580.95. These high correlations indicate that controlling for PM10, as in the two-pollutants analysis, is largely equivalent to controlling for PM2.5. However, the data used in this approach are limited to only a few cities. Seaton et al. 40 suggest that NO2 is a surrogate for particle numbers. Their hypothesis is illustrated using data only from Aberdeen where the correlation of NO2 with PM10 is 0.45, with PM2.5, 0.55, and with particle numbers, 0.89. However, this correlation is fairly high. In a recent three-city study 41, the correlation between ultrafine particles and NO2 were 0.49, 0.72 and 0.82. Therefore, in many cases it is possible to separate the effects of NO2 and ultrafine particles. Using information from four European cities (the RUPIOH project, unpublished data) that provided relevant data for the present analysis, the correlation between NO2 and PM10 ranged 0.380.54, between NO2 and PM2.5, 0.430.63, between NO2 and particle numbers, 0.530.72, and between PM10 and particle numbers, 0.190.54. Hence, in the present data, the associations between NO2 and particle numbers are not so high as to make the distinction of the effects impossible, as in the case of Seaton et al. 40. It should be noted that among these four cities, the one with the highest NO2 particle numbers correlation (Helsinki) had the lowest estimated NO2 effect on total mortality. The present authors investigated heterogeneity in the effect estimates between cities by taking into account several city characteristics, and environmental and social factors as potential effect modifiers. It was found that when smoking prevalence is lower, the effect of NO2 on total and cardiovascular mortality is larger. A possible explanation may be that smoking acts as a competing risk, harvesting the population of susceptible individuals, but this issue needs further investigation. Furthermore, the effect of NO2 on total and cardiovascular mortality was observed mainly in western and southern European cities. In more than half of the eastern European cities involved in the analysis, the NO2 1-h levels were in the lower range of the observed distribution. Domestic gas consumption as an indicator of domestic NO2 exposure acted as potential effect modifier in cardiovascular mortality, indicating that the effect of outdoor NO2 exposure increased with higher indoor exposure. The risk contrast is not so pronounced and it is difficult to believe that indoor pollution levels act as an effect modifier for the association of outdoor NO2 levels and mortality, unless the slope becomes steeper at higher cumulative exposures. Another likely explanation may lie in the association of higher gas consumption in the more developed cities with characteristic air pollution sources (e.g. more traffic). Finally, the effect of NO2 on respiratory mortality was higher in cities with a larger proportion of elderly persons in the population, for whom there is evidence of increased susceptibility to other pollutants. For cities with higher PM10 levels, larger NO2 effects were estimated, which may mean that the two pollutants act synergistically. One limitation of the present studys second-stage analysis is that there is a noticeable difference in the effects between Eastern European and other cities. Hence, any variables, such as the mortality rate, etc. that strongly differ between the two regions are likely to appear as significant explanatory factors of the heterogeneity in the effect estimates. Or alternatively, if the effect modifiers identified are indeed true modifiers, their distribution across Europe may result in these marked geographical differences. In conclusion, the present results confirm those previously reported concerning the adverse effects of nitrogen dioxide on mortality and complement them by investigating potential confounding by other pollutants and possible effect modification. The results of this large study are consistent with an independent effect of nitrogen dioxide on mortality, but the role of nitrogen dioxide as a surrogate of other unmeasured pollutants cannot be completely ruled out. Since the short-term effects of nitrogen dioxide on respiratory mortality may be confounded by other vehicle-derived pollutants, further investigation is needed to enhance understanding of the underlying biological mechanism.
This article has been cited by other articles:
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||