Short-term effects of desert and non-desert PM(10) on mortality in Sicily, Italy.
Increased PM10 concentrations are commonly observed during Saharan dust advections. Limited epidemiological evidence suggests that PM10 from anthropogenic and desert sources increase mortality. We aimed to evaluate the association between source-specific PM10 (non-desert and desert) and cause-specific mortality in Sicily during 2006-2012 period.
Daily PM10 concentrations at 1-km2 were estimated in Sicily using satellite-based data, fixed monitors and land use variables. We identified Saharan dust episodes using meteorological models, back-trajectories, aerosol maps, and satellite images. For each dust day, we estimated desert and non-desert PM10 concentrations. We applied a time-series approach on 390 municipalities of Sicily to estimate the association between PM10 (non-desert and desert) and daily cause-specific mortality.
33% of all days were affected by Saharan dust advections. PM10 concentrations were 8 μg/m3 higher during dust days compared to other days. We found positive associations of both non-desert and desert PM10 with cause-specific mortality. We estimated percent increases of risk (IR%) of non-accidental mortality equal to 2.3% (95% Confidence Interval [CI]: 1.4, 3.1) and 3.8% (3.2, 4.4), per 10 μg/m3 increases in non-desert and desert PM10 at lag 0-5, respectively. We also observed significant associations with cardiovascular (2.4% [1.3, 3.4] and 4.5% [3.8, 5.3]) and respiratory mortality (8.1% [6.8, 9.5], and 6.3% [5.4, 7.2]). We estimated higher effects during April-September, with IR% = 4.4% (3.2, 5.7) and 6.3% (5.4, 7.2) for non-desert and desert PM10, respectively.
Our results confirm previous evidence of harmful effects of desert PM10 on non-accidental and cardio-respiratory mortality, especially during the warm season.
Renzi M
,Forastiere F
,Calzolari R
,Cernigliaro A
,Madonia G
,Michelozzi P
,Davoli M
,Scondotto S
,Stafoggia M
... -
《-》
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
... -
《-》
[Meta-analysis of the Italian studies on short-term effects of air pollution].
In recent years, much attention has been given to review reports on the early effects of air pollution on health, measured through daily series of deaths and/or hospital admissions. A number of large planned meta-analyses (in which methods for data retrieval and processing are commonly planned a priori for all participating centers) are on going both in the US and in Europe. The National Mortality, Morbidity and Air Pollution Study included data from 90 US cities, whereas APHEA (Air Pollution and Health, a European Approach) considers data from about 30 european cities. The present paper summarizes methods and findings of MISA, a meta-analysis of data from 8 Italian cities. It belongs to an ad hoc supplement of Epidemiologia & Prevenzione (Epidemiol Prev 2001; 25 (2) Suppl: 1-72), the official Journal of the Italian Association of Epidemiology, which contains a full description of the study. MISA was launched on March 2000, within the project "Statistics, Environment and Health" (GRASPA), funded by the Italian Ministry of Education. Additional support was given by the Authorities of the 8 participating cities (from North to South: Turin, Milan, Verona, Ravenna, Bologna, Florence, Rome and Palermo). DAILY HEALTH DATA: Deaths certificate and hospital admission data have been collected respectively from the Local Health Authority and regional files. The same programme for retrieval of data on selected hospital admissions for acute conditions was used in the 8 cities. Main data are summarized in Table 1. DAILY CONCENTRATION OF POLLUTANTS: Most data were obtained from Regional Environmental Protection Agencies, which are responsible for environmental monitoring since 1993. Verona, Palermo and Milan (1990-94) data were obtained from local sources. Monitors with more than 25% of missing data were excluded. Meteorological data were collected by the same monitors and completed with data from monitors situated in the suburbs or (in Milan and Bologna) in the airport. The monitors were selected by a group of experts to ensure comparability. For SO2 and NO2 daily averages of hourly measurements were used, whereas concentrations of ozone and CO were estimated as the maximum 8 hours moving average. Total suspended particulate or PM10 were measured as 24 hours deposition. All analyses used the whole range of observed values (Table 2). Daily data were considered as missing when more than 25% of hourly data were not available. Missing data in one monitor were imputed as average of data from the remaining monitors weighted by the ratio between the specific monitor's year average and the general year average of all the selected city monitors. Missing data in one day were imputed as average of four days (preceding and following day, the same day of the previous and following weeks). In the city of Florence and Palermo PM10 concentrations were available. For the other cities we applied a conversion factor from PTS to PM10 (0.6 for Turin and 0.8 for all the others) estimated through validation studies. Ozone concentrations were used only where background monitors were available (Turin, Verona, Bologna and Florence) and limited to the warm season (May through September).
A common protocol for the city-specific analyses was defined on the basis of a structured exploratory analysis. The adopted basic model was a Generalized Additive Model for Poisson data. Effect estimates were age-adjusted (0-64, 65-74, 75+) and formal tests of interaction pollutant-age were conducted. In the first two age groups, indicator variables for seasonality were specified, and cubic splines with fixed number of degree of freedom were specified for the last age group and for all age groups for the morbidity data. Model adequacy was checked by residual analysis and inspection of the partial autocorrelation function. In a sensitivity analysis non linear pollutant effects were considered and overdispersed [table: see text] transitional models were fitted; the analysis was conducted for all lags 0-3 and some distributed lags (0-1, 1-2, 0-3); no multipollutant models were fitted. The same model was fitted to the city data. No model selection was done: Table 3 describes the steps in model building. In the meta-analysis, for each outcome, the estimates for each pollutant and for each city were combined using fixed and random effects models. Heterogeneity of effects was tested according to DerSimonian and Laird. Results were checked using a hierarchical bayesian model, which was used to investigate heterogeneity across cities in a meta-regression phase. Non informative priors were used. Posterior distributions of parameters of interest have been obtained with WinBUGS. 10,000 iterations (excluding [table: see text] the first 2000) were retained, while for the meta-regression 100,000 iterations (excluding the first 4000) were stored. To approximate the marginal posteriors only one sample out of five were used. Achieved convergence was assessed using the Gelman and Rubin approach. In the meta-regression the models specified were the following: [formula: see text] i denotes city, j calendar period (1990-1994; 1995-1999). The first model includes only period as effect modifier, while the second model other potential variables. The ui terms (which do not vary with j) represent city specific random effects.
For each pollutant, the meta-analysis detected a statistically significant association with mortality for natural causes. But for ozone, positive associations were commonly found for death and hospital admissions for both cardiovascular and respiratory diseases. Indeed, the only estimates whose lower 95% confidence limit bore a negative sign regarded the association between PM10 and mortality from respiratory diseases. Ozone in the warm season was positively and significantly associated with daily mortality and mortality for cardiovascular diseases whereas other estimates did not reach statistical significance and some were negative (only lag 0-1 for external comparability are reported in Table 4). Risks were highest (up to 4%) for respiratory conditions (Table 4). They were more pronounced at lag 1-2 for mortality, and at lag 0-3 for hospital admissions. Age was an effect modifier for mortality, the elderly being more susceptible. In the random effect meta-analysis, at lag 1-2, excess risks for unit increase of the pollutants at age 75+ and at age 0-64 were respectively: 4.9% and -0.4% for SO2, 1.7% and 0.6% for NO2; 2.3% and 0.2% for CO. Corresponding figures for PM10 at lag 0-1 were 1.1% and 0.2%. The effect of PM10 on mortality [table: see text] was greater during the warm season (2.8% vs 0.8%). A complete analysis is reported in the Italian text. Here we provide some details on the effects of PM10, about which the residual heterogeneity across cities was highest (Table 4). In addition, the epidemiological evidence on the hazards from this fraction of particulate matter is more controversial. Table 5 reports the excess risk estimated through the meta-analysis in 1995-99 for a 10 micrograms/m3 increase of PM10 for some outcomes. Proper prior distributions (overdispersed normal and inverse gamma) were adopted in the final bayesian analyses. The sensitivity of results to the choice of the priors were investigated (we defined proper and improper uniform, student's t), obtaining comparable results. Total natural mortality was significantly heterogeneous across cities (Q = 18.96, 5 df, p < 0.001). City-specific estimates are represented graphically in Fig. 1. As expected, the confidence (credibility) intervals are widest [table: see text] for bayesian estimates, intermediate for those obtained under a random effects model, and narrowest for those found under a fixed effects model. Nevertheless, differences in point estimates are negligible. A North-South gradient in risk is obvious. Table 6 shows, for the cities for which mortality data were available, the improvement in precision and the shrinkage of effect estimates toward the overall mean introduced by the bayesian modelling. In the meta-regression, total mortality and a deprivation score were associated with greater effects. The excess risks on hospital admission were modified by the deprivation score and by the NO2/PM10 ratio. Overall, the risk estimates were greater in the calendar period 1995-99 and there was a North-South gradient, with larger effects in cities located in Central and Southern Italy (Florence, Rome, Palermo).
The meta-analysis of the Italian studies on short-term effects of air pollution in 8 cities, MISA, exhibits the following features: With the exception of Naples, all greatest Italian cities were included; overall a population of 7 million was enrolled. The study protocol was accurate with regard to the selection of hospital admissions for acute conditions. Monitored data of concentration of pollutant were carefully evaluated before their inclusion in the meta-analysis. City specific analyses were carried out according to a common protocol controlling for seasonality, influenza epidemics, age and meterological variables; [table: see text] the protocol derived from a structured exploratory analysis. The meta-analysis was done using fixed and random effects models; a hierarchical bayesian model was fitted in a sensitivity analysis. The heterogeneity of effects across cities was investigated using a hierarchical bayesian model for meta-regression. While mortality data are of good quality, hospital admission data are more problematic. Since the filing criteria for the latter changed around 1995, comparability of results before and after such date is limited. Moreover, hospital admissions rely on availability of beds, the offer of which may be restricted during the warm season. Comparability of pollutant concentration estimates among cities may have been influenced by differences in monitor characteristics. (ABSTRACT TRUNCATED)
Biggeri A
,Bellini P
,Terracini B
,Italian MISA Group
... -
《Epidemiologia & Prevenzione》