Part 2. Association of daily mortality with ambient air pollution, and effect modification by extremely high temperature in Wuhan, China.
Fewer studies have been published on the association between daily mortality and ambient air pollution in Asia than in the United States and Europe. This study was undertaken in Wuhan, China, to investigate the acute effects of air pollution on mortality with an emphasis on particulate matter (PM*). There were three primary aims: (1) to examine the associations of daily mortality due to all natural causes and daily cause-specific mortality (cardiovascular [CVD], stroke, cardiac [CARD], respiratory [RD], cardiopulmonary [CP], and non-cardiopulmonary [non-CP] causes) with daily mean concentrations (microg/m3) of PM with an aerodynamic diameter--10 pm (PM10), sulfur dioxide (SO2), nitrogen dioxide (NO2), or ozone (O3); (2) to investigate the effect modification of extremely high temperature on the association between air pollution and daily mortality due to all natural causes and daily cause-specific mortality; and (3) to assess the uncertainty of effect estimates caused by the change in International Classification of Disease (ICD) coding of mortality data from Revision 9 (ICD-9) to Revision 10 (ICD-10) code. Wuhan is called an "oven city" in China because of its extremely hot summers (the average daily temperature in July is 37.2 degrees C and maximum daily temperature often exceeds 40 degrees C). Approximately 4.5 million residents live in the core city area of 201 km2, where air pollution levels are higher and ranges are wider than the levels in most cities studied in the published literature. We obtained daily mean levels of PM10, SO2, and NO2 concentrations from five fixed-site air monitoring stations operated by the Wuhan Environmental Monitoring Center (WEMC). O3 data were obtained from two stations, and 8-hour averages, from 10:00 to 18:00, were used. Daily mortality data were obtained from the Wuhan Centres for Disease Prevention and Control (WCDC) during the study period of July 1, 2000, to June 30, 2004. To achieve the first aim, we used a regression of the logarithm of daily counts of mortality due to all natural causes and cause-specific mortality on the daily mean concentrations of the four pollutants while controlling for weather, temporal factors, and other important covariates with generalized additive models (GAMs). We derived pollutant effect estimations for 0-day, 1-day, 2-day, 3-day, and 4-day lagged exposure levels, and the averages of 0-day and 1-day lags (lag 0-1 day) and of 0-day, 1-day, 2-day, and 3-day lags (lag 0-3 days) before the event of death. In addition, we used individual-level data (e.g., age and sex) to classify subgroups in stratified analyses. Furthermore, we explored the nonlinear shapes ("thresholds") of the exposure-response relations. To achieve the second aim, we tested the hypothesis that extremely high temperature modifies the associations between air pollution and daily mortality. We developed three corresponding weather indicators: "extremely hot," "extremely cold," and "normal temperatures." The estimates were obtained from the models for the main effects and for the pollutant-temperature interaction for each pollutant and each cause of mortality. To achieve the third aim, we conducted an additional analysis. We examined the concordance rates and kappa statistics between the ICD-9-coded mortality data and the ICD-10-coded mortality data for the year 2002. We also compared the magnitudes of the estimated effects resulting from the use of the two types of ICD-coded mortality data. In general, the largest pollutant effects were observed at lag 0-1 day. Therefore, for this report, we focused on the results obtained from the lag 0-1 models. We observed consistent associations between PM10 and mortality: every 10-microg/m3 increase in PM10 daily concentration at lag 0-1 day produced a statistically significant association with an increase in mortality due to all natural causes (0.43%; 95% confidence interval [CI], 0.24 to 0.62), CVD (0.57%; 95% CI, 0.31 to 0.84), stroke (0.57%; 95% CI, 0.25 to 0.88), CARD (0.49%; 95% CI, 0.04 to 0.94), RD (0.87%; 95% CI, 0.34 to 1.41), CP (0.52%; 95% CI, 0.27 to 0.77), and non-CP (0.30%; 95% CI, 0.05 to 0.54). In general, these effects were stronger in females than in males and were also stronger among the elderly (> or = 65 years) than among the young. The results of sensitivity testing over the range of exposures from 24.8 to 477.8 microg/m3 also suggest the appropriateness of assuming a linear relation between daily mortality and PM10. Among the gaseous pollutants, we also observed statistically significant associations of mortality with NO, and SO2, and that the estimated effects of these two pollutants were stronger than the PM10 effects. The patterns of NO2 and SO2 associations were similar to those of PM10 in terms of sex, age, and linearity. O3 was not associated with mortality. In the analysis of the effect modification of extremely high temperature on the association between air pollution and daily mortality, only the interaction of PM10 with temperature was statistically significant. Specifically, the interaction terms were statistically significant for mortality due to all natural (P = 0.014), CVD (P = 0.007), and CP (P = 0.014) causes. Across the three temperature groups, the strongest PM10 effects occurred mainly on days with extremely high temperatures for mortality due to all natural (2.20%; 95% CI, 0.74 to 3.68), CVD (3.28%; 95% CI, 1.24 to 5.37), and CP (3.02%; 95% CI, 1.03 to 5.04) causes. The weakest effects occurred at normal temperature days, with the effects on days with low temperatures in the middle. To assess the uncertainty of the effect estimates caused by the change from ICD-9-coded mortality data to ICD-10-coded mortality data, we compared the two sets of data and found high concordance rates (> 99.3%) and kappa statistics close to 1.0 (> 0.98). All effect estimates showed very little change. All statistically significant levels of the estimated effects remained unchanged. In conclusion, the findings for the aims from the current study are consistent with those in most previous studies of air pollution and mortality. The small differences between mortality effects for deaths coded using ICD-9 and ICD-10 show that the change in coding had a minimal impact on our study. Few published papers have reported synergistic effects of extremely high temperatures and air pollution on mortality, and further studies are needed. Establishing causal links between heat, PM10, and mortality will require further toxicologic and cohort studies.
Qian Z
,He Q
,Lin HM
,Kong L
,Zhou D
,Liang S
,Zhu Z
,Liao D
,Liu W
,Bentley CM
,Dan J
,Wang B
,Yang N
,Xu S
,Gong J
,Wei H
,Sun H
,Qin Z
,HEI Health Review Committee
... -
《-》
Multicity study of air pollution and mortality in Latin America (the ESCALA study).
The ESCALA* project (Estudio de Salud y Contaminación del Aire en Latinoamérica) is an HEI-funded study that aims to examine the association between exposure to outdoor air pollution and mortality in nine Latin American cities, using a common analytic framework to obtain comparable and updated information on the effects of air pollution on several causes of death in different age groups. This report summarizes the work conducted between 2006 and 2009, describes the methodologic issues addressed during project development, and presents city-specific results of meta-analyses and meta-regression analyses.
The ESCALA project involved three teams of investigators responsible for collection and analysis of city-specific air pollution and mortality data from three different countries. The teams designed five different protocols to standardize the methods of data collection and analysis that would be used to evaluate the effects of air pollution on mortality (see Appendices B-F). By following the same protocols, the investigators could directly compare the results among cities. The analysis was conducted in two stages. The first stage included analyses of all-natural-cause and cause-specific mortality related to particulate matter < or = 10 pm in aerodynamic diameter (PM10) and to ozone (O3) in cities of Brazil, Chile, and México. Analyses for PM10 and O3 were also stratified by age group and O3 analyses were stratified by season. Generalized linear models (GLM) in Poisson regression were used to fit the time-series data. Time trends and seasonality were modeled using natural splines with 3, 6, 9, or 12 degrees of freedom (df) per year. Temperature and humidity were also modeled using natural splines, initially with 3 or 6 df, and then with degrees of freedom chosen on the basis of residual diagnostics (i.e., partial autocorrelation function [PACF], periodograms, and a Q-Q plot) (Appendix H, available on the HEI Web site). Indicator variables for day-of-week and holidays were used to account for short-term cyclic fluctuations. To assess the association between exposure to air pollution and risk of death, the PM10 and O3 data were fit using distributed lag models (DLMs). These models are based on findings indicating that the health effects associated with air pollutant concentrations on a given day may accumulate over several subsequent days. Each DLM measured the cumulative effect of a pollutant concentration on a given day (day 0) and that day's contribution to the effect of that pollutant on multiple subsequent (lagged) days. For this study, exposure lags of up to 3, 5, and 10 days were explored. However, only the results of the DLMs using a 3-day lag (DLM 0-3) are presented in this report because we found a decreasing association with mortality in various age-cause groups for increasing lag effects from 3 to 5 days for both PM10 and O3. The potential modifying effect of socioeconomic status (SES) on the association of PM10 or O3 concentration and mortality was also explored in four cities: Mexico City, Rio de Janeiro, São Paulo, and Santiago. The methodology for developing a common SES index is presented in the report. The second stage included meta-analyses and metaregression. During this stage, the associations between mortality and air pollution were compared among cities to evaluate the presence of heterogeneity and to explore city-level variables that might explain this heterogeneity. Meta-analyses were conducted to combine mortality effect estimates across cities and to evaluate the presence of heterogeneity among city results, whereas meta-regression models were used to explore variables that might explain the heterogeneity among cities in mortality risks associated with exposures to PM10 (but not to O3).
The results of the mortality analyses are presented as risk percent changes (RPC) with a 95% confidence interval (CI). RPC is the increase in mortality risk associated with an increase of 10 microg/m3 in the 24-hour average concentration of PM10 or in the daily maximum 8-hour moving average concentration of O3. Most of the results for PM10 were positive and statistically significant, showing an increased risk of mortality with increased ambient concentrations. Results for O3 also showed a statistically significant increase in mortality in the cities with available data. With the distributed lag model, DLM 0-3, PM10 ambient concentrations were associated with an increased risk of mortality in all cities except Concepci6n and Temuco. In Mexico City and Santiago the RPC and 95% CIs were 1.02% (0.87 to 1.17) and 0.48% (0.35 to 0.61), respectively. PM10 was also significantly associated with increased mortality from cardiopulmonary, respiratory, cardiovascular, cerebrovascular-stroke, and chronic obstructive lung diseases (COPD) in most cities. The few nonsignificant effects generally were observed in the smallest cities (Concepción, Temuco, and Toluca). The percentage increases in mortality associated with ambient O3 concentrations were smaller than for those associated with PM10. All-natural-cause mortality was significantly related to O3 in Mexico City, Monterrey, São Paulo and Rio de Janeiro. Increased mortality risks for some specific causes were also observed in these cities and in Santiago. In the analyses stratified by season, different patterns in mortality and O3 were observed for cold and warm seasons. Risk estimates for the warm season were larger and significant for several causes of death in São Paulo and Rio de Janeiro. Risk estimates for the cold season were larger and significant for some causes of death in Mexico City, Monterrey, and Toluca. In an analysis stratified by SES, the all-natural-cause mortality risk in Mexico City was larger for people with a medium SES; however we observed that the risk of mortality related to respiratory causes was larger among people with a low SES, while the risk of mortality related to cardiovascular and cerebrovascular-stroke causes was larger among people with medium or high SES. In São Paulo, the all-natural-cause mortality risk was larger in people with a high SES, while in Rio de Janeiro the all-natural-cause mortality risk was larger in people with a low SES. In both Brazilian cities, the risks of mortality were larger for respiratory causes, especially for the low- and high-SES groups. In Santiago, all-natural-cause mortality risk did not vary with level of SES; however, people with a low SES had a higher respiratory mortality risk, particularly for COPD. People with a medium SES had larger risks of mortality from cardiovascular and cerebrovascular-stroke disease. The effect of ambient PM10 concentrations on infant and child mortality from respiratory causes and lower respiratory infection (LRI) was studied only for Mexico City, Santiago, and São Paulo. Significant increased mortality risk from these causes was observed in both Santiago (in infants and older children) and Mexico City (only in infants). For O3, an increased mortality risk was observed in Mexico City (in infants and older children) and in São Paulo (only in infants during the warm season). The results of the meta-analyses confirmed the positive and statistically significant association between PM10 and all-natural-cause mortality (RPC = 0.77% [95% CI: 0.60 to 1.00]) using the random-effects model. For mortality from specific causes, the percentage increase in mortality ranged from 0.72% (0.54 to 0.89) for cardiovascular disease to 2.44% (1.36 to 3.59) for COPD, also using the random-effects model. For O3, significant positive associations were observed using the random-effects model for some causes, but not for all natural causes or for respiratory diseases in people 65 years or older (> or = 65 years), and not for COPD and cerebrovascular-stroke in the all-age and the > or = 65 age groups. The percentage increase in all-natural-cause mortality was 0.16% (-0.02 to 0.33). In the meta-regression analyses, variables that best explained heterogeneity in mortality risks among cities were the mean average of temperature in the warm season, population percentage of infants (< 1 year), population percentage of children at least 1 year old but < 5 years (i.e., 1-4 years), population percentage of people > or = 65 years, geographic density of PM10 monitors, annual average concentrations of PM10, and mortality rates for lung cancer.
The ESCALA project was undertaken to obtain information for assessing the effects of air pollutants on mortality in Latin America, where large populations are exposed to relatively high levels of ambient air pollution. An important goal was to provide evidence that could inform policies for controlling air pollution in Latin America. This project included the development of standardized protocols for data collection and for statistical analyses as well as statistical analytic programs (routines developed in R by the ESCALA team) to insure comparability of results. The analytic approach and statistical programming developed within this project should be of value for researchers carrying out single-city analyses and should facilitate the inclusion of additional Latin American cities within the ESCALA multicity project. Our analyses confirm what has been observed in other parts of the world regarding the effects of ambient PM10 and 03 concentrations on daily mortality. They also suggest that SES plays a role in the susceptibility of a population to air pollution; people with a lower SES appeared to have an increased risk of death from respiratory causes, particularly COPD. Compared with the general population, infants and young children appeared to be more susceptible to both PM10 and O3, although an increased risk of mortality was not observed in these age groups in all cities. (ABSTRACT TRUNCATED)
Romieu I
,Gouveia N
,Cifuentes LA
,de Leon AP
,Junger W
,Vera J
,Strappa V
,Hurtado-Díaz M
,Miranda-Soberanis V
,Rojas-Bracho L
,Carbajal-Arroyo L
,Tzintzun-Cervantes G
,HEI Health Review Committee
... -
《-》