ABSTRACT

OBJECTIVE

To evaluate the impact of air pollution on hospitalizations for respiratory and cardiovascular diseases in the largest Brazilian metropolis.

METHODS

This study was carried out at the Metropolitan Region of São Paulo, Brazil. Environmental data were obtained from the network of monitoring stations of nine municipalities. Air pollution exposure was measured by daily means of PM₁₀ (particles with a nominal mean aerodynamic diameter ≤ 10 μm) per municipality, while daily counts of hospitalizations for respiratory and cardiovascular diseases within the Brazilian Unified Health System were the outcome. For each municipality a time series analysis was carried out in which a semiparametric Poisson regression model was the framework to explain the daily fluctuations on counts of hospitalizations over time. The results were combined in a meta-analysis to estimate the overall risk of PM₁₀ in hospitalizations for respiratory and cardiovascular diseases at the Metropolitan Region of São Paulo.

RESULTS

Regarding hospitalizations for respiratory diseases, the effect estimates were statistically significant (p < 0.05) for all municipalities, except Santo André and Taboão da Serra. The RR (Relative Risk) of this outcome for an increase of 10 µg/m³ in the levels of PM₁₀ ranged from 1.011 (95%CI 1.009–1.013) for São Paulo to 1.032 (95%CI 1.024–1.040) in São Bernardo do Campo. The RR of hospitalization for respiratory diseases in children for an increase of 10 µg/m3 of PM₁₀ ranged from 1.009 (95%CI 1.001–1.017) in Santo André to 1.077 (95%CI 1.056–1.098) in Mauá. Only São Paulo and São Bernardo do Campo presented positive and statistically significant results for hospitalizations for cardiovascular diseases.

CONCLUSIONS

This is the first study to estimate the risk of illness from air pollution in the set of municipalities of the Metropolitan Region of São Paulo, Brazil. Global estimates of the effect of exposure to pollution in the region indicated associations only with respiratory diseases. Only São Paulo and São Bernardo do Campo showed an association between the levels of PM₁₀ and hospitalizations for cardiovascular diseases.

Air Pollution, adverse effects; Respiratory Diseases, epidemiology; Cardiovascular Diseases, epidemiology; Meta-analysis

INTRODUCTION

The RMSP (Metropolitan Region of São Paulo), also known as Greater São Paulo, accounts for thirty-nine municipalities of the state of São Paulo, comprising an area of 7,946 km², with a population of approximately 20 million inhabitants in 2010 (Figure 1)^aa Fundação Sistema Estadual de Análise de Dados. Portal de Estatísticas do Estado de São Paulo: tabelas da Região Metropolitana do Estado de São Paulo. São Paulo: SEADE; 2010 [cited 2016 Jul 1]. Available from: http://www.imp.seade.gov.br/frontend/#/tabelas . The region had an estimated circulating fleet of 6.9 million vehicles in 2012, which corresponds to 49% of the fleet of the state, and it has the largest concentration of industries in Brazil, which put it in the spotlight in relation to air pollution and the study of its impacts from mobile and fixed sources^bb Companhia de Tecnologia e Saneamento Ambiental de São Paulo. Relatório da Qualidade do Ar 2013. São Paulo: CETESB; 2014 [cited 2016 Jul 5]. Available from: http://www.cetesb.sp.gov.br/ar/qualidade-do-ar/31-publicacoes-e-relatorios .

These sources of pollution accounted for the emission of approximately 132,000 metric tons per year of CO (carbon monoxide), 42,000 metric tons per year of HC (hydrocarbons), 77,000 metric tons per year of NO_x (nitrogen oxides), 4,500 metric tons per year of PM (particulate matter), and 11,000 metric tons per year of SO_x (sulfur oxides); in addition, vehicles account for 97% of CO emissions, 81% of HC emissions, 80% of NO_x emissions, 48% of SO_x emissions, and 40% of PM emissions^cc Companhia de Tecnologia e Saneamento Ambiental de São Paulo. Emissões veiculares no Estado de São Paulo 2013. São Paulo: CETESB; 2014 [cited 2016 Jul 5]. Available from: http://veicular.cetesb.sp.gov.br/wp-content/uploads/sites/35/2013/12/relatorio-emissoes-veiculares-2013.pdf . We highlight that the RMSP is responsible for approximately 40% of the CO, NMHC (non-methane hydrocarbon), and RCHO (aldehyde) emissions and 30% of the NO_x, PM, and SO₂ emissions in the state of São Paulo^cc Companhia de Tecnologia e Saneamento Ambiental de São Paulo. Emissões veiculares no Estado de São Paulo 2013. São Paulo: CETESB; 2014 [cited 2016 Jul 5]. Available from: http://veicular.cetesb.sp.gov.br/wp-content/uploads/sites/35/2013/12/relatorio-emissoes-veiculares-2013.pdf .

Thus, in order to follow the evolution of air quality and to guide environmental control measures, CETESB (Environmental Company of the State of São Paulo) manages an air quality-monitoring network in the RMSP, with twenty-six fixed stations, in ten municipalities (Carapicuíba, Diadema, Guarulhos, Mauá, Osasco, Santo André, São Bernardo do Campo, São Caetano do Sul, São Paulo, and Taboão da Serra). Twenty-three of these stations measure particles with a nominal mean aerodynamic diameter ≤ 10 μm (PM₁₀)^bb Companhia de Tecnologia e Saneamento Ambiental de São Paulo. Relatório da Qualidade do Ar 2013. São Paulo: CETESB; 2014 [cited 2016 Jul 5]. Available from: http://www.cetesb.sp.gov.br/ar/qualidade-do-ar/31-publicacoes-e-relatorios .

Numerous studies have shown the impact of exposure to PM₁₀ on health, with emphasis on respiratory and cardiovascular diseases. Among the respiratory outcomes, exposure to PM₁₀ is associated with increased hospitalizations for asthma, pneumonia, and chronic obstructive pulmonary disease¹⁰10. Gouveia N, Freitas CU, Martins LC, Marcilio IO. Hospitalizações por causas respiratórias e cardiovasculares associadas à contaminação atmosférica no Município de São Paulo, Brasil. Cad Saude Publica. 2006;22(12):2669-77. https://doi.org/10.1590/S0102-311X2006001200016
https://doi.org/10.1590/S0102-311X200600... , as well as exacerbation of symptoms associated with respiratory allergic diseases¹⁷17. Nicolussi FH, Santos APM, André SCS, Veiga TB, Takayanagui AMM. Poluição do ar e doenças respiratórias alérgicas em escolares. Rev Saude Publica. 2014;48(2):326-30. http://doi.org/10.1590/S0034-8910.2014048004940
http://doi.org/10.1590/S0034-8910.201404... . Regarding cardiovascular diseases, exposure to PM₁₀ increases hospitalizations and death from cerebrovascular accident¹1. Amancio CT, Nascimento LF. Environmental pollution and deaths due to stroke in a city with low levels of air pollution: ecological time series study. Sao Paulo Med J. 2014;132(6):353-8. https://doi.org/10.1590/1516-3180.2014.1326733
https://doi.org/10.1590/1516-3180.2014.1... ^,1515. Nascimento LFC, Francisco JB, Patto MBR, Antunes AM. Environmental pollutants and stroke-related hospital admissions. Cad Saude Publica. 2012;28(7):1319-24. https://doi.org/10.1590/S0102-311X2012000700010
https://doi.org/10.1590/S0102-311X201200... and congestive heart failure⁷7. Evo CPR, Ulrych BK, Takegawa B, Soares G, Nogueira G, Oliveira LO, et al. Poluição do ar e internação por insuficiência cardíaca congestiva em idosos no município de Santo André. Arq Bras Cienc Saude. 2011 [cited 2017 Jul 10];36(1):6-9. Available from: http://files.bvs.br/upload/S/1983-2451/2011/v36n1/a1917.pdf
http://files.bvs.br/upload/S/1983-2451/2... , hospitalizations and emergency services for hypertension¹¹11. Guo Y, Tong S, Zhang Y, Barnett AG, Jia Y, Pan X. The relationship between particulate air pollution and emergency hospital visits for hypertension in Beijing, China. Sci Total Environ. 2010;408(20):4446-50. https://doi.org/10.1016/j.scitotenv.2010.06.042
https://doi.org/10.1016/j.scitotenv.2010... and arrhythmia²⁰20. Santos UP, Terra-Filho M, Lin CA, Pereira LA, Vieira TC, Saldiva PH, et al. Cardiac arrhythmia emergency room visits and environmental air pollution in Sao Paulo, Brazil. J Epidemiol Community Health. 2008;62(3):267-72. https://doi.org/10.1136/jech.2006.058123
https://doi.org/10.1136/jech.2006.058123... , and hospitalizations for acute myocardial infarction¹⁰10. Gouveia N, Freitas CU, Martins LC, Marcilio IO. Hospitalizações por causas respiratórias e cardiovasculares associadas à contaminação atmosférica no Município de São Paulo, Brasil. Cad Saude Publica. 2006;22(12):2669-77. https://doi.org/10.1590/S0102-311X2006001200016
https://doi.org/10.1590/S0102-311X200600... .

In addition to the health effects described above, an adverse association has been reported between PM₁₀ and cancer incidence (skin, lung, thyroid, larynx, and bladder), deaths from lung cancer²³23. Yanagi Y, Assunção JV, Barrozo LV. The impact of atmospheric particulate matter on cancer incidence and mortality in the city of São Paulo, Brazil. Cad Saude Publica. 2012;28(9):1737-48. https://doi.org/10.1590/S0102-311X2012000900012
https://doi.org/10.1590/S0102-311X201200... , incidence of deep venous thrombosis²2. Baccarelli A, Martinelli I, Zanobetti A, Grillo P, Hou LF, Bertazzi PA, et al. Exposure to particulate air pollution and risk of deep vein thrombosis. Arch Intern Med. 2008;168(9):920-7. https://doi.org/10.1001/archinte.168.9.920
https://doi.org/10.1001/archinte.168.9.9... , premature births¹²12. Lima TAC, Nascimento LFC, Medeiros APP, Santos VP. Association between maternal exposure to particulate matter and premature birth. Rev Ambient Agua. 2014;9(1):27-36. https://doi.org/10.4136/ambi-agua.1262
https://doi.org/10.4136/ambi-agua.1262... , low birth weight⁹9. Gouveia N, Bremner SA, Novaes HMD. Association between ambient air pollution and birth weight in São Paulo, Brazil. J Epidemiol Community Health. 2004;58(1):11-7. https://doi.org/10.1136/jech.58.1.11
https://doi.org/10.1136/jech.58.1.11... ^,1414. Medeiros A, Gouveia N. Relação entre baixo peso ao nascer e a poluição do ar no Município de São Paulo. Rev Saude Publica. 2005;39(6):965-72. https://doi.org/10.1590/S0034-89102005000600015
https://doi.org/10.1590/S0034-8910200500... , increased neonatal mortality¹³13. Lin CA, Pereira LAA, Nishioka DC, Conceição GMS, Braga ALF, Saldiva PHN. Air pollution and neonatal deaths in São Paulo, Brazil. Braz J Med Biol Res. 2004;37(5):765-70. https://doi.org/10.1590/S0100-879X2004000500019
https://doi.org/10.1590/S0100-879X200400... . Furthermore, PM₁₀ is an important risk factor in the early development of atherosclerosis¹⁸18. Poursafa P, Kelishadi R, Lahijanzadeh A, Modaresi M, Javanmard SH, Assari R, et al. The relationship of air pollution and surrogate markers of endothelial dysfunction in a population-based sample of children. BMC Public Health. 2011;11:115. https://doi.org/10.1186/1471-2458-11-115
https://doi.org/10.1186/1471-2458-11-115... and in the onset of painful episodes in patients with sickle cell anemia³3. Barbosa SMM, Farhat SCL, Martins LC, Pereira LAA, Saldiva PHN, Zanobetti A, et al. Air pollution and children’s health: sickle cell disease. Cad Saude Publica. 2015;31(2):265-75. https://doi.org/10.1590/0102-311X00013214
https://doi.org/10.1590/0102-311X0001321... .

Among the municipalities of the RMSP that have air quality monitoring, several studies on the deleterious effects of atmospheric contamination were carried out in the municipality of São Paulo. More recently, studies in the municipality of Santo André have shown an association between pollution and hospitalizations for heart failure⁷7. Evo CPR, Ulrych BK, Takegawa B, Soares G, Nogueira G, Oliveira LO, et al. Poluição do ar e internação por insuficiência cardíaca congestiva em idosos no município de Santo André. Arq Bras Cienc Saude. 2011 [cited 2017 Jul 10];36(1):6-9. Available from: http://files.bvs.br/upload/S/1983-2451/2011/v36n1/a1917.pdf
http://files.bvs.br/upload/S/1983-2451/2... ^,1616. Negrete BR, Rosa CC, Ikeuti DH, Delena PJ, Borba TM, Braga ALF. Poluição atmosférica e internações por insuficiência cardíaca congestiva em adultos e idosos em Santo André (SP). Arq Bras Cienc Saude. 2010 [cited 2017 Jul 10];35(3):208-12. Available from: http://files.bvs.br/upload/S/1983-2451/2010/v35n3/a1688
http://files.bvs.br/upload/S/1983-2451/2... , low birth weight¹⁹19. Romão R, Pereira LAA, Saldiva PHN, Pinheiro PM, Braga ALF, Martins LC. The relationship between low birth weight and exposure to inhalable particulate matter. Cad Saude Publica. 2013;29(6):1101-8. https://doi.org/10.1590/S0102-311X2013000600007
https://doi.org/10.1590/S0102-311X201300... , and increased blood pressure in traffic controllers⁵5. Chiarelli PS, Pereira LAA, Saldiva PHN, Ferreira Filho C, Garcia MLB, Braga ALF, et al. The association between air pollution and blood pressure in traffic controllers in Santo André, São Paulo, Brazil. Environ Res. 2011;111(5):650-5. https://doi.org/10.1016/j.envres.2011.04.007
https://doi.org/10.1016/j.envres.2011.04... .

Thus, in order to evaluate the impact of air pollution in the RMSP, a study was carried out with data from these monitoring stations, and PM₁₀ was the indicator of pollution and hospitalizations for respiratory and cardiovascular diseases were the indicators of effect. The results were then combined in a meta-analysis to estimate the overall risk of PM₁₀ in the region. Estimates of the global impact of air pollution on the population of the RMSP are important in order to subsidize discussions about interventions aimed at controlling this exposure.

METHODS

The data sets were extracted from the website of information whose construction was financed by the Ministry of Health for the Project named Projeto de Avaliação de Impacto da Poluição do Ar nas Cidades Brasileiras (http://www.observandosaopaulo.com.br). In this website, data on AIH (Hospitalization Authorizations) were grouped per day according to municipality of residence, from 2000 to 2008. Information on AIH is systematized by DATASUS (Department of Informatics of the Brazilian Universal Health System). In addition, this website also stores data on air pollutants collected by the monitoring stations of CETESB, as well as meteorological data from IAG (Instituto de Astronomia e Geofísica) of the Universidade de São Paulo and CETESB. The dependent variables were the daily counts of hospitalizations for: all respiratory diseases, i.e., ICD-10 J00-J99, in all ages, hereafter TRD, all respiratory diseases in children younger than five years old, hereafter CRD, and cardiovascular diseases, i.e., ICD-10 I00-I99, in adults aged 39 years old or more, hereafter CVD. Some municipalities presented great variability in the volume of hospitalizations, with a very small daily average and with many days without events, which makes the analysis rather complex. Thus, as a general rule for the inclusion of a municipality in the Projeto de Avaliação de Impacto da Poluição do Ar nas Cidades Brasileiras, it was established that:

• municipalities with daily average of three or more events were analyzed;

• for municipalities with a daily average of less than three events but more than one, the model fit was regarded more carefully. If there was no adjustment, the corresponding result would be examined with caution; and

• municipalities with daily average of less than one event were not analyzed.

In the RMSP, not all municipalities monitor every regulated pollutant, but they all monitor PM₁₀. Therefore, the focus of this study was to analyze exposure to PM₁₀ in all municipalities as well as to summarize its effect estimates in the RMSP. Hourly PM₁₀ measures were aggregated to daily averages. The municipality of São Paulo has thirteen air quality monitoring stations that measure PM₁₀ and the municipalities of Santo André and Guarulhos have two stations each. In the other municipalities, only one air quality monitoring station of the CETESB measures PM₁₀.

For each municipality, overall exposure was the daily average over the set of available monitoring stations in the municipality. Time series of a pollutant measured in a monitoring station may have gaps, often with large periods without information. Therefore, all databases were evaluated for the completeness of the data on PM₁₀. A minimum period of analysis of three consecutive years was required, and the total sum of gaps could not account for more than 15% of the days. Consequently, there were variable periods of analysis of the municipalities to the detriment of what we initially proposed, which was to evaluate the impact of air pollution on the health of the RMSP population for the period from 2000 to 2008. Furthermore, the municipality of Carapicuíba was not included in the analysis, as it did not fulfill the above mentioned requirements in the period from 2000 to 2008.

Temperature and relative humidity of the air were used as controls in every time series analysis; however, most of the municipalities do not have meteorological stations. Therefore, the temperature and relative humidity indicators used for the municipalities were the same as those used for São Paulo. They were the average of the daily averages, the average of the daily minima, and the average of the daily maxima of the CETESB and IAG stations, all located in São Paulo. Even though there are meteorological monitoring stations in Guarulhos and São Caetano do Sul, their datasets presented large gaps, which hindered the use of these data.

The time series analyses were carried out using the ARES library, developed for the R platform^dd Junger WL. Análise, imputação de dados e interfaces computacionais em estudos de séries temporais epidemiológicas [thesis]. Rio de Janeiro: Instituto de Medicina Social da Universidade Estadual do Rio de Janeiro; 2008 [cited 2017 Jul 10]. Available from: http://www.bdtd.uerj.br/tde_arquivos/44/TDE-2013-05-16T162258Z-3184/Publico/Washington%20Leite%20Junger-tese.pdf . A common methodology was applied to choose the models for the sites that joined the “Projeto de Avaliação de Impacto da Poluição do Ar nas Cidades Brasileiras”.

We constructed explanatory models for daily counts of hospitalizations for each cause over time for each municipality participating in the study. The proposed models belong to the class of GAM (Generalized Additive Models), with the option of Poisson regression, according to the equation:

where Y_t and X_1t are the numbers of morbid events (outcome) and the level of a given pollutant (exposure) in day t, respectively, X_it are the predictor variables, which include time, and S_i are smoothing functions, here natural splines. We also added indicator variables (dichotomous) for the days of the week and national or local holidays, testing their significance. Holidays with significance of up to 0.09 were grouped according to the direction of their effect: positive or negative. In the time series modeling process, we sought to minimize the AIC (Akaike Information Criterion) and optimize the PACF (partial autocorrelation function).

A working model was developed by municipality (Core Model). Under each model, all control variables were included as described above, followed by a thorough examination of the properties of the residuals, after which the model was regarded as appropriate, or not. If the residual diagnostics revealed any anomaly, the Core model was modified accordingly. The exposure time series – daily PM₁₀ – was introduced in lags of up to five days (single lag) and we also verified the cumulative effect in that period, using a polynomial distributed lag model. This model, besides considering the latency of the effect of the pollutants, minimizes the instability in the estimation process, typical of analyses that use multiple lags²¹21. Schwartz J. The distributed lag between air pollution and daily deaths. Epidemiology. 2000;11(3):320-6. https://doi.org/10.1097/00001648-200005000-00016
https://doi.org/10.1097/00001648-2000050... .

The risk estimates were calculated after introducing the exposure variable in the Core model via a linear term (see above equation). From the statistical model, we can obtain the percentage relative risk for each increase of 10µg/m³ in the levels of PM_10. In all analyses, we considered a significance level of 5%.

The following step was to perform the meta-analysis of the individual findings (municipalities), under the assumption of random effects (DerSimonian and Laird method⁶6. DerSimonian R, Laird N. Meta-analysis in clinical trials. Control Clin Trials. 1986;7(3):177-88. https://doi.org/10.1016/0197-2456(86)90046-2
https://doi.org/10.1016/0197-2456(86)900... ) and the estimation of heterogeneity as proposed by Mantel-Haenszel, since the hypothesis of homogeneity was rejected in all cases. The meta-analysis was performed in Stata 12, using the macro metan, which allows the meta-analysis of fixed or random effects.

RESULTS

Table 1 shows the average number of hospitalizations for respiratory and cardiovascular diseases, as well as the descriptive analysis of the levels of PM₁₀, for each of the nine municipalities analyzed. With the exception of São Paulo, the other municipalities presented relatively small numbers of daily hospitalizations; however, only São Caetano do Sul had an average of less than one hospitalization per day for respiratory diseases in children. Guarulhos and Osasco presented the highest average values of PM₁₀. The percentage of missing PM₁₀ values was rather low for all municipalities.

Given that the core model was developed for each municipality, there was a great variation in the number of degrees of freedom used to control time trends and seasonality, mainly because of the given period of analysis and the model choice criteria such as parsimony (smaller AIC) and optimization of the PACF. As an example, for total respiratory diseases in São Paulo, with a longer time period of analysis, 36 degrees of freedom were regarded to account for seasonality and trend, while for Taboão da Serra only 12 were needed. The daily value of each indicator, the average of two days, or the simple lag up to two days was used to control the temperature and humidity conditions in each municipality. For these indicators, few degrees of freedom were necessary, reaching a maximum of five in Mauá.

The relative risk estimates for the cumulative effects of zero to five days using the polynomial distributed lag model are presented in Table 2. All municipalities presented statistically significant estimates for hospitalizations for respiratory diseases with the exception of TRD in Santo André and CRD in Taboão da Serra. The relative risk of hospitalization for TRD associated to an increase of 10 µg/m³ in the levels of PM₁₀ ranged from 1.011 (95%CI 1.009–1.013) for São Paulo to 1.032 (95%CI 1.024–1.040) in São Bernardo do Campo. The risk of hospitalization for CRD ranged from 1.009 (95%CI 1.001–1.017) in Santo André to 1.077 (95%CI 1.056–1.098) in Mauá.

On the other hand, the estimates for hospitalizations for CVD were statistically significant only for São Paulo, São Bernardo do Campo, and Santo André. In Santo André, we observed a relative risk of less than one. The other estimates presented great variation in addition to not being statistically significant (Table 2).

The meta-analysis of the three outcomes studied indicated heterogeneity among the individual estimates of each municipality; thus, random effects models were estimated to obtain summarized estimates. The results of the meta-analysis are presented in Figures 2, 3, and 4.

Note that the overall estimates of the effect of PM₁₀ for an increase of 10 μg/m³ were statistically significant and of similar magnitude for hospitalizations for respiratory diseases and all ages (RR = 1.014, 95%CI 1.009–1.019) and for children younger than five years old (RR = 1.019, 95%CI 1.012–1.025). On the other hand, the summary estimate for hospitalizations for CVD was approximately one and not statistically significant.

DISCUSSION

This is the first study to estimate the risk of illness from air pollution in a large set of municipalities of the RMSP. We observed positive and statistically significant associations between exposure to PM₁₀ and hospitalizations for total respiratory diseases and in children younger than five years old in almost all nine municipalities that have air quality monitors. On the other hand, only São Paulo and São Bernardo do Campo showed an association between the levels of PM₁₀ and hospitalizations for CVD. Similarly, global estimates of the effect of exposure to pollution in the region indicated associations only with respiratory diseases. In summary, there was an increase of 1.4% in hospitalizations for total respiratory diseases for each increase of 10 μg/m³ in the levels of PM₁₀. For children younger than five years old, the effect is slightly higher, with a 1.9% increase in hospitalizations.

The magnitude of the effects, both for the global estimates and those related to each municipality analyzed, was consistent with similar studies, both national⁸8. Freitas CU, Junger W, Ponce de Leon AC, Grimaldi R, Silva MAFR, Gouveia N. Poluição do ar em cidades brasileiras: selecionando indicadores de impacto na saúde para fins de vigilância. Epidemiol Serv Saude. 2013;22(3):445-54. https://doi.org/10.5123/S1679-49742013000300009
https://doi.org/10.5123/S1679-4974201300... ^,1010. Gouveia N, Freitas CU, Martins LC, Marcilio IO. Hospitalizações por causas respiratórias e cardiovasculares associadas à contaminação atmosférica no Município de São Paulo, Brasil. Cad Saude Publica. 2006;22(12):2669-77. https://doi.org/10.1590/S0102-311X2006001200016
https://doi.org/10.1590/S0102-311X200600... and international⁴4. Carugno M, Consonni D, Randi G, Catelan D, Grisotto L, Bertazzi PA, et al. Air pollution exposure, cause-specific deaths and hospitalizations in a highly polluted Italian region. Environ Res. 2016;147:415-24. https://doi.org/10.1016/j.envres.2016.03.003
https://doi.org/10.1016/j.envres.2016.03... ^,1111. Guo Y, Tong S, Zhang Y, Barnett AG, Jia Y, Pan X. The relationship between particulate air pollution and emergency hospital visits for hypertension in Beijing, China. Sci Total Environ. 2010;408(20):4446-50. https://doi.org/10.1016/j.scitotenv.2010.06.042
https://doi.org/10.1016/j.scitotenv.2010... . For example, a study carried out in the municipality of São Paulo has found an increase of 2.4% in hospitalizations for respiratory diseases in children younger than five years old associated with an increase of 10 μg/m³ in the levels of PM₁₀¹⁰10. Gouveia N, Freitas CU, Martins LC, Marcilio IO. Hospitalizações por causas respiratórias e cardiovasculares associadas à contaminação atmosférica no Município de São Paulo, Brasil. Cad Saude Publica. 2006;22(12):2669-77. https://doi.org/10.1590/S0102-311X2006001200016
https://doi.org/10.1590/S0102-311X200600... . This same study has found a statistically significant association for hospitalizations for CVD and with a magnitude very close to what we observed in this analysis for São Paulo.

We can also observe that respiratory diseases seem to be the best indicator of the health effects associated with air pollution. Freitas et al.⁸8. Freitas CU, Junger W, Ponce de Leon AC, Grimaldi R, Silva MAFR, Gouveia N. Poluição do ar em cidades brasileiras: selecionando indicadores de impacto na saúde para fins de vigilância. Epidemiol Serv Saude. 2013;22(3):445-54. https://doi.org/10.5123/S1679-49742013000300009
https://doi.org/10.5123/S1679-4974201300... , in a study carried out in several Brazilian municipalities, have identified hospitalizations for respiratory diseases in children younger than five years old, followed by hospitalizations for total respiratory diseases, as the best indicators of the health effects of air pollution for surveillance purposes.

We can also note the high degree of atmospheric contamination to which the population of the RMSP is exposed. Two municipalities, Guarulhos and Osasco, presented an annual average of PM₁₀ higher than the air quality standard established by the National Environment Council (CONAMA), equivalent to 50 μg/m³. However, we highlight that the State Environment Council (CONSEMA) has ruled that the state of São Paulo should adopt stricter standards for air quality²²22. Valentim LSO. Consema aprova novos padrões de qualidade do ar para São Paulo. BEPA Bol Epidemiol Paulista. 2011 [cited 2017 Jul 10];8(91):27-9. Available from: http://www.cvs.saude.sp.gov.br/up/Consema%20aprova%20novos%20%20%20padr%C3%B5es%20de%20qualidade%20do%20ar%20para%20S%C3%A3o%20Paulo.pdf
http://www.cvs.saude.sp.gov.br/up/Consem... . In May 2011, new air quality standards were adopted for the state of São Paulo, aiming to reach the limits recommended by the World Health Organization, of 20 μg/m³, as the maximum annual average for PM₁₀. We highlight that all the municipalities of the RMSP would be outside this value according to this new parameter not yet fully in force in the state.

As any epidemiological study, this analysis presents some limitations, such as the fact that only 10 of the 39 municipalities of the RMSP monitor air quality. In addition, many of the municipalities have only one monitor, which may not adequately reflect the exposure of the whole population living there. Another limitation is the fact that many residents of these municipalities tend to spend most of the day in other locations, especially in the more central regions of São Paulo, because of work or study.

On the other hand, the estimates obtained by the meta-analysis should be reflecting the true impact of air pollution on the residents of the RMSP. It is important to know and quantify the estimates of the global impact of air pollution on the population living there to help design interventions that aim to reduce exposure to air pollution, as well as to estimate the economic costs of the health effects of this exposure. These data can also support the monitoring of variations in the health effects resulting from these possible interventions. As air quality gains more and more prominence as an important determinant of population health, cities must seek to ensure a better quality of air and, consequently, a better quality of life for their inhabitants.

REFERENCES

Publication Dates

History