Space-time analysis of the effect of air pollution on children’s health

Emerson Pedreira Matos Valdério Anselmo Reisen Faradiba Sarquis Serpa Paulo Roberto Prezotti Filho Maria de Fátima Silva Leite About the authors

Abstract

The study aimed to investigate the short-term association between air pollution and emergency treatments for respiratory diseases in children 0 to 6 years of age. This was an ecological space-time study in Greater Metropolitan Vitória, Espírito Santo State, Brazil. A Poisson regression general additive model (GAM) used the number of daily treatments for respiratory diseases as the dependent variable, and the independent variables were daily concentrations of air pollutants (PM10, SO2, NO2, O3, and CO), temperature, humidity, and precipitation. Average daily concentrations were used to make estimates for the entire metropolitan area and in loco analyses considering children residing in a 2km radius around 8 air quality monitoring stations. An increase of 10μg/m3 in the concentration of air pollutants increased the risk of emergency treatment for respiratory disease. In the overall area, for PM10, the increase was 2.43%, 2.73%, and 3.29% in the cumulative values at 5, 6, and 7 days, respectively. For SO2, the increase was 4.47% on the day of exposure, 5.26% two days later, and 6.47%, 8.8%, 8.76%, and 7.09% for the cumulative values at days 2, 3, 4, and 5 days, respectively. CO showed a significant association for residents around two stations, and O3 for only one. Even within the limits set by the World Health Organization, the pollutants PM10, SO2, NO2, and O3 are associated with increased risk of treatment for respiratory diseases in children 0 to 6 years of age, and some effects were only identified when disaggregating by neighborhood, i.e., in loco, which allows capturing greater variation in the data.

Keywords:
Air Pollution; Respiratory Tract Diseases; Child


Introduction

Respiratory diseases are the principal cause of morbidity and mortality in children and can be caused or aggravated by exposure to air pollutants emitted on a large scale throughout the planet 11. World Health Organization. Inheriting a sustainable world? Atlas on children's health and the environment. Geneva: World Health Organization; 2017.,22. Holgate ST. Every breath we take: the lifelong impact of air pollution - a call for action. Clin Med (Lond) 2017; 17:8-12.. Emissions of particulate matter (PM), nitrogen oxides (NOx), volatile organic compounds (VOCs), sulfur dioxide (SO2), and photochemical pollutants such as ozone (O3) have increased in recent decades due to growth of the motor vehicle fleet and increasing industrialization 22. Holgate ST. Every breath we take: the lifelong impact of air pollution - a call for action. Clin Med (Lond) 2017; 17:8-12.. Studies in large urban areas prove a significant association between levels of these pollutants and the number of emergency treatments and hospitalizations from respiratory causes 33. Cifuentes LA, Vega J, Köpfer K, Lave LB. Effect of the fine fraction of particulate matter versus the coarse mass and other pollutants on daily mortality in Santiago, Chile. J Air Waste Manag Assoc 2000; 50:1287-98.,44. Gouveia N, Mendonça GAS, Ponce de Leon A, Correia JEM, Junger WL, Freitas CU, et al. Poluição do ar e efeitos na saúde nas populações de duas grandes metrópoles brasileiras. Epidemiol Serv Saúde 2003; 12:29-40.,55. Daumas RP, Mendonça GAS, Ponce de León A. Poluição do ar e mortalidade em idosos no município do Rio de Janeiro: análise de série temporal. Cad Saúde Pública 2004; 20:311-9.,66. Bakonyi SM, Danni-Oliveira IM, Martins LC, Braga AL. Poluição atmosférica e doenças respiratórias em crianças na Cidade de Curitiba - PR. Rev Saúde Pública 2004; 38:695-700.,77. Nascimento LF, Pereira LA, Braga AL, Módolo MCC, Carvalho Jr. JA. Effects of air pollution on children's health in a city in southeastern Brazil. Rev Saúde Pública 2006; 40:77-82.,88. Braga ALF, Pereira LAA, Procópio M, André PA, Saldiva PHN. Associação entre poluição atmosférica e doenças respiratórias e cardiovasculares na cidade de Itabira, Minas Gerais, Brasil. Cad Saúde Pública 2007; 23 Suppl 4:S570-8.,99. Chen R, Chu C, Tan J, Cao J, Song W, Xu X, et al. Ambient air pollution and hospital admission in Shanghai, China. J Hazard Mater 2010; 181:234-40.,1010. Samoli E, Nastos PT, Paliatsos AG, Katsouyanni K, Priftis KN. Acute effects of air pollution on pediatric asthma exacerbation: evidence of association and effect modification. Environ Res 2011; 111:418-24.,1111. Almeida SP, Casimiro E, Calheiros J. Short-term association between exposure to ozone and mortality in Oporto, Portugal. Environ Res 2011; 111:406-10.,1212. Souza JB, Reisen VA, Santos JM, Franco GC. Principal components and generalized linear modeling in the correlation between hospital admissions and air pollution. Rev Saúde Pública 2014; 48:451-8.,1313. Souza JB, Reisen VA, Franco GC, Ispány M, Bondon P, Santos JM. Generalized additive models with principal component analysis: an application to time series of respiratory disease and air pollution data. J R Stat Soc Series B Stat Methodol 2018; 67:453-80.,1414. Nascimento AP, Santos JM, Mill JG, Souza JB, Reis Júnior NC, Reisen VA. Associação entre concentração de partículas finas na atmosfera e doenças respiratórias agudas em crianças. Rev Saúde Pública 2017; 51:3..

A recent systematic review on Brazilian children’s health and the environment identified only 17 time-series studies that correlated air pollution and respiratory causes, all conducted in Southeast Brazil. These studies showed a positive association between levels of pollutants and emergency treatments, hospitalizations, and decreased lung function 1515. Froes Asmus CI, Camara VM, Landrigan PJ, Claudio L. A systematic review of children's environmental health in Brazil. Ann Glob Health 2016; 82:132-48.. Two of these studies were conducted in Greater Metropolitan Vitória, Espírito Santo State, Brazil 1212. Souza JB, Reisen VA, Santos JM, Franco GC. Principal components and generalized linear modeling in the correlation between hospital admissions and air pollution. Rev Saúde Pública 2014; 48:451-8.,1616. Castro HA, Hacon S, Argento R, Junger WL, Mello CF, Castiglioni Júnior N. Air pollution and respiratory diseases in the Municipality of Vitória, Espírito Santo State, Brazil. Cad Saúde Pública 2007; 23 Suppl 4:S630-42..

Greater Metropolitan Vitória is predominantly urban and highly industrialized and is undergoing urban sprawl. Air quality is affected by the emission of pollutants from motor vehicles and large-scale industrial activities, in addition to the logistics sector, with the existence of a large port and airport complex 1717. Instituto Estadual de Meio Ambiente e Recursos Hídricos. Relatório da qualidade do ar: Grande Vitória 2013. Vitória: Instituto Estadual de Meio Ambiente e Recursos Hídricos; 2013.. Emissions from these mobile and fixed sources directly influence the area’s air quality 1717. Instituto Estadual de Meio Ambiente e Recursos Hídricos. Relatório da qualidade do ar: Grande Vitória 2013. Vitória: Instituto Estadual de Meio Ambiente e Recursos Hídricos; 2013..

Detection of the association between concentrations of pollutants and respiratory conditions requires statistical techniques that allow isolating the effects of air pollution, since there may be various confounding factors such as temperature and humidity. Time-series models have allowed more precise analyses of such associations.

Epidemiological time-series studies generally use a single fixed monitoring station or the average of stations to represent the entire population’s exposure to pollutants. However, this approach may not reflect people’s true exposure 1818. Namdeo A, Tiwary A, Farrow E. Estimation of age-related vulnerability to air pollution: assessment of respiratory health at local scale. Environ Int 2011; 37:829-37.. Thus, a more realistic analysis might consider a small geographic area close to the monitoring station 1919. Zeger SL, Thomas D, Dominici F, Samet JM, Schwartz J, Dockery D, et al. Exposure measurement error in time-series studies of air pollution: concepts and consequences. Environ Health Perspect 2000; 108:419-26.. Still, a limited number of epidemiological studies use this type of assessment, although there is some evidence that exposure to classification error in time-series analysis tends towards a downward bias in estimates and in this sense does not limit the results’ importance for public health 1919. Zeger SL, Thomas D, Dominici F, Samet JM, Schwartz J, Dockery D, et al. Exposure measurement error in time-series studies of air pollution: concepts and consequences. Environ Health Perspect 2000; 108:419-26..

Thus, methodologies in which the areas covered by measurement are better discriminated spatially are important for obtaining more precise inferences. In this context, the current study aimed to assess, via GAM, the short-term relationship between the number of emergency treatments for respiratory problems in children under six years of age and levels of air pollutants in Greater Metropolitan Vitória, considering the temporal variables measured on site (area around monitoring stations) and the average between stations (overall average for Greater Metropolitan Vitória).

Methods

Study area

This was an ecological space-time study conducted in Greater Metropolitan Vitória from January 1, 2005, to December 31, 2010. Greater Metropolitan Vitória consists of 7 municipalities (Vitória, Vila Velha, Cariacica, Serra, Viana, Guarapari, and Fundão), with a territory of 2,318.917km2 and approximately 1.7 million inhabitants, as one of the main urban and industrial development hubs in the State of Espírito Santo (Instituto Brasileiro de Geografia e Estatística. Censo populacional 2010. https://ww2.ibge.gov.br/home/estatistica/populacao/censo2010, accessed on 19/Oct/2018). The principal polluting activities in Greater Metropolitan Vitória include roadways, various industries (steel, pelletizing, mining, cement mills), ports, airports, and residential and commercial emissions 1717. Instituto Estadual de Meio Ambiente e Recursos Hídricos. Relatório da qualidade do ar: Grande Vitória 2013. Vitória: Instituto Estadual de Meio Ambiente e Recursos Hídricos; 2013..

Air quality monitoring in Greater Metropolitan Vitória is done by nine stations which jointly comprise the Automatic Air Quality Monitoring Network (RAMQAr, in Portuguese) managed by the Espírito Santo State Institute for the Environment and Water Resources 1717. Instituto Estadual de Meio Ambiente e Recursos Hídricos. Relatório da qualidade do ar: Grande Vitória 2013. Vitória: Instituto Estadual de Meio Ambiente e Recursos Hídricos; 2013.. The spatial analysis considers each station’s area of influence, as proposed by Santolim 2020. Santolim L. Rede otimizada para o monitoramento da qualidade do ar da Grande Vitória, Espírito Santo [Dissertação de Mestrado]. Vitória: Universidade Federal do Espírito Santo; 1991.. The area corresponds to a circle around the station with a radius of approximately 2km. At the time of the study, eight of the nine RAMQAr stations (Figure 1) were functioning.

Figure 1
Neighborhoods with air quality monitoring stations.

Health outcome

Data on emergency treatments for respiratory conditions in children under 6 years of age were obtained from the emergency departments of two hospitals, one public and the other private, the Nossa Senhora da Glória Children’s Hospital and the Unimed Vitória Integrated Healthcare Center, respectively. Respiratory diseases were coded according to the 10th revision of the International Classification of Diseases (ICD-10: J00-J99). We selected the treatments of children that lived in the neighborhoods located in the area of influence of each of the RAMQAr stations.

Environmental pollutants and meteorological variables

The concentrations of particulate matter with aerodynamic volume up to 10 microns (PM10), SO2, nitrogen dioxide (NO2), O3, and carbon monoxide (CO) for the target period were provided by the Espírito Santo State Institute for the Environment and Water Resources (IEMA), which is responsible for the nine RAMQAr stations: Laranjeiras, Carapina, Cidade Continental, Jardim Camburi, Enseada do Suá, Vitória-Centro, Ibes, Vila Velha-Centro, and Cariacica. These automatic stations continuously collect and analyze air samples and process the data as hourly averages, on-site and in real time. The measurement method for PM10 is Tapered Element Oscillating Microbalance (TEOM), which measures the mass concentration continuously. Measurement of SO2 uses the principle of ultraviolet ray fluorescence. Measurement of NO2 uses a combination of dual crossflow modulation type, also used to measure CO and O3, with the principle of chemiluminescence and the differential calculus method 1717. Instituto Estadual de Meio Ambiente e Recursos Hídricos. Relatório da qualidade do ar: Grande Vitória 2013. Vitória: Instituto Estadual de Meio Ambiente e Recursos Hídricos; 2013.. We thus calculated the daily averages for PM10 and SO2 and used maximum hourly concentrations of NO2 and maximum daily concentrations of the eight-hour moving averages for CO and O3. All pollutants in all the stations were measured in micrograms per cubic meter (µg/m3).

The target meteorological variables were: temperature in degrees Celsius (ºC); relative humidity in percentages (%), recorded by the Carapina and Cariacica stations; and precipitation in millimeters (mm), recorded by the Carapina monitor, which are the stations that record these measurements. We used these data to calculate the arithmetic means of the recorded measurements of humidity and temperature (maximum, average, and minimum) and precipitation represented by daily rainfall. We also included for each day the temperature and humidity values from the current medical care day (lag0), the previous day (lag1), two days previously (lag2), and moving averages of the two and three previous days (mm01, mm02).

Missing data

Failures in monitoring of pollutants in the RAMQAr network during the period, both on single days and on consecutive days, caused gaps in the records of concentrations and were corrected by imputation, according to the methodology described by Junger 2121. Junger WL. Análise, imputação de dados e interfaces computacionais em estudos de séries temporais epidemiológicas [Tese de Doutorado]. Rio de Janeiro: Universidade do Estado do Rio de Janeiro; 2008.. The estimates obtained with this method are explained by the spatial correlation between the different levels of the same pollutant in the different monitors and by the autocorrelation in the pollutant’s levels in the same monitor, over time. The temporal patterns are modeled with a Gaussian autoregressive integrated moving average (ARIMA) model. This algorithm is especially adapted to weather data with missing measurements in some monitors across a given region. The upper and lower limits in the imputation of missing data were the maximum and minimum concentrations, respectively, observed in the historical series for each pollutant’s concentrations.

Statistical analysis

The modeling strategy consisted of defining a central model with all the known information (trend, seasonality, days of the week, holidays, and meteorological conditions) in order to explain the variability in the number of treatments for respiratory diseases, except for the concentration of pollutants.

The choice of variables and covariables for the model was based on tests and diagnoses in each stage of the modeling process. The diagnoses were based on residual analysis and the Akaike criterion (AIC) 2222. Hastie TJ, Tibshirani RJ. Monographs on statistics and applied probability 43: generalizes additive models. New York: Chapman & Hall; 1990..

Generalized additive model

The daily number of medical treatments represents a counting process, and the generalized additive model (GAM) with Poisson marginal distribution was the statistical tool used to estimate the shape of the curve in the relationship between the health outcome and air pollution 55. Daumas RP, Mendonça GAS, Ponce de León A. Poluição do ar e mortalidade em idosos no município do Rio de Janeiro: análise de série temporal. Cad Saúde Pública 2004; 20:311-9.,1313. Souza JB, Reisen VA, Franco GC, Ispány M, Bondon P, Santos JM. Generalized additive models with principal component analysis: an application to time series of respiratory disease and air pollution data. J R Stat Soc Series B Stat Methodol 2018; 67:453-80.,1414. Nascimento AP, Santos JM, Mill JG, Souza JB, Reis Júnior NC, Reisen VA. Associação entre concentração de partículas finas na atmosfera e doenças respiratórias agudas em crianças. Rev Saúde Pública 2017; 51:3..

Let {Υt} ≡ {Υt}t∈Z be a counting series, that is, yt ∈ {0,1,...}. The conditional distribution Υt, given the past Ft-1 that contains all the available information thus far t - 1, is denoted by:

p(yt;μt|Ft-1)=e-μtμtytyt!(1)

in which μt represents the expected value (average) of Yt. Thus, given a sample Y,...,Yn, consisting of “n” mutually conditional independent random variables belonging to “Yt”, the conditional log-likelihood function is given by:

l(μ)=t=1nlnp(yt;μt|Ft-1)t=1nYtlnμt-μt(2)

in which the vector μ = (μ1,...,μn) depends on the parameters and the process {Yt}. Let Xt = [X1t,..., Xpt]T be the vector of covariables in the dimension p in time t, in which T denotes the transpose, which can include past values of Yt and other ancillary information such as pollutants and confounding variables (trend, seasonality, and meteorological variables, among others). In this study, the sequence X1t,..., Xqt denotes the concentrations of pollutants PM10, SO2, NO2, O3, and CO, thus, q = 5, and X(q+1)t,..., Xpt indicates the confounding variables in time t, (p > q).

The relationship between the vector is given by:

lnμt=qj=0βjXjt+pj=q+1fjXjt(3)

in which (β0,β), with β = (β1,..., βq)T is the vector of the coefficients to be estimated (βj it is the j-th covariable coefficient) and fj is the smoothing function for the j-th confounding variable. In addition, β0 indicates the curve’s intercept and is associated with X0t= 1 for every t. The entire modeling process was performed in the R software (http://www.r-project.org) with the ARES package 2323. Junger W, Leon AP. Ares: environment air pollution epidemiology. A library for time series analysis. Vienna: The Comprehensive R Archive Network; 2011. (CRAN Archives)..

The relative risk (RR) of a pollutant covariable Xi,j = 1,...,q is given as the relative variation in the expected count of respiratory disease events by the unit variation ξ in the covariable while maintaining the other covariables constant. According to Baxter et al. 2424. Baxter LA, Finch SJ, Lipfert FW, Yu Q. Comparing estimates of the effects of air pollution on human mortality obtained using different regression methodologies. Risk Anal 1997; 17:273-8., formula (8), the RR is given by:

RRXj(ξ)=E(Y|Xj=ξ,Xi,ij)E(Y|Xj=0,Xi,ij)(4)

For the Poisson regression, the RR does not depend on the xi, i ≠ j values of the other covariables and can be expressed as:

RRX1ξ=exp(βj^ξ)(5)

For the GAM with Poisson marginal distribution, the RR and its approximate confidence interval (CI, or IC in Portuguese), at level of significance α of a covariable, Xj, j = 1...q, is estimated as follows:

RRXj(ξ)Xj=exp(βj^ξ)(6)

IC(RRXj(ξ))=exp(β^jξ±zα/2β^jξ)(7)

β^j is the estimated coefficient associated with pollutant Xj in a study with standard error sβ^j and Zα/2 is the quantile α/2 of the standard normal distribution. At level of significance α, the hypothesis to be tested is defined as H0: RRxj = 1 against H0: RRxj > 1, in which RRxj = RRxj(1), that is, RR of the unit variation in Xj. The rejection of H0 implies statistically that the respective pollutant has a significant adverse effect on health.

In this study, calculations of the RR(xj) values correspond to an increase of 1,000μg/m3 in the CO levels and of 10μg/m3 for the other pollutants. The results are presented as percentage increases in the number of medical treatments and are calculated by the expression:

%RRx=RRx-1x100(8)

Lag

The biological manifestation of the effects of pollution on human health apparently display a behavior that shows a lag in relation to the individual’s exposure to the pollutants, that is, events that occurred on a given day may be associated with the pollution levels on that day and/or on previous days. Thus, according to previous studies 2525. Zhang H, Niu Y, Yao Y, Chen R, Zhou X, Kan H. The impact of ambient air pollution on daily hospital visits for various respiratory diseases and the relevant medical expenditures in Shanghai, China. Int J Environ Res Public Health 2018; 15:E425.,2626. Chen R, Yin P, Meng X, Liu C, Wang L, Xu X. Fine particulate air pollution and daily mortality. a nationwide analysis in 272 Chinese cities. Am J Respir Crit Care Med 2017; 196:73-81., we decided to investigate the association between the number of treatments for respiratory diseases and the pollution levels on the day of the emergency treatment (lag0) and on the previous days (lag1, lag2, lag3). The cumulative effect was assessed with the moving averages from two to eight days (MA01, MA02, MA03, MA04, MA05, MA06, MA07) 2727. Katsouyanni K, Schwartz J, Spix C, Touloumi G, Zmirou D, Zanobetti A, et al. Short term effects of air pollution on health: a european approach using epidemiologic time series data: the APHEA protocol. J Epidemiol Community Health 1996; 50 Suppl 1:S12-8.,2828. Qu Y, Pan Y, Niu H, He Y, Li M, Li L, et al. Short-term effects of fine particulate matter on non-accidental and circulatory diseases mortality: a time series study among the elder in Changchun. PLoS One 2018; 13:e0209793..

The project was approved by the Institutional Review Board of the Health Sciences Center of the Federal University of Espírito Santo, under number 04/11, on May 14, 2011.

Results

During the study period there were 46,421 emergency treatments recorded for respiratory diseases in children 0 to 6 years of age living in the area of influence of the eight RAMQAr monitoring stations. The average daily number of treatments in Greater Metropolitan Vitória was 21.19 (SD = 9.90), ranging from 1.72 to 4.84 per day, with higher monthly averages from March to June, corresponding to autumn and early winter. The highest number of treatments was in children living in the Enseada do Suá neighborhood.

Average temperature during the period varied from 20.85 to 29.36°C, and rainfall varied from 0mm to 117.80mm (average = 3.78mm) and relative humidity varied from 61.79% to 97.27% (average = 77.47%).

Concentrations of pollutants did not display a uniform behavior between the different RAMQAr stations. The highest average concentrations recorded in the period were: PM10 (43.06μg/m3) in Cariacica, SO2 (16.32μg/m3) in Enseada do Suá, O3 (38.66μg/m3) in Ibes, and NO2 (44.10μg/m3) and CO (1,730.91μg/m3) in Vitória-Centro. Table 1 shows the pollutants’ minimum, average, and maximum concentrations and 25%, 50%, and 75% quartiles and standard deviations (SD).

Table 1
Descriptive statistics of treatments for respiratory diseases, by daily averages of PM10 (μg/m3), SO2 (μg/m3), and NO2 (μg/m3) and 8-hour moving averages of O3 (μg/m3) and CO (μg/m3) in each Automatic Air Quality Monitoring Network (RAMQAr, in Portuguese), 2005-2010.

The series of daily counts of treatments in the overall region was smoothed with a “spline” with 12 degrees of freedom, defined by the Akaike modeling criterion and residual analysis. Non-parametric analysis evidenced seasonality and a downward trend over time, confounding factors that were included in the modeling process. The period from autumn (March) to winter (June) displayed a seasonality with an increase in the number of treatments for respiratory diseases.

Figure 2 shows the goodness-of-fit analysis for the model described previously via the results from Poisson regression for estimation of the effect of PM10 on the 6-day moving average for Greater Metropolitan Vitória.

Figure 2
Diagnosis of the central model for Greater Metropolitan Vitória, Espírito Santo State, Brazil.

The estimated coefficient for PM10 was 0.0032, with a standard error of 0.0015 and p-value of 0.0301, which we represent as [MA6, Greater Metropolitan Vitória] (0.0032; 0.00145; 0.0301).

The model’s fit to the data for the areas in loco resulted in [MA6, Vila Velha-Centro] (0.0086; 0.0036; 0.0173), [MA6, Jardim Camburi] (0.0114; 0.0045; 0.0112), [MA6, Cariacica] (0.0044, 0.0021, 0.0331), and [MA6, Carapina] (0.0116; 0.0043; 0.0074). The results for the other sites were not statistically significant.

With the values estimated by the models, we calculated the relative risks corresponding to the estimates and proceeded to the analyses as presented.

We observed no empirical evidence of the model’s poor fit, that is, the residuals are not correlated and are approximately normal. The periodogram also proved that the residuals presented characteristics of white noise. These residual graphic analyses provide the necessary support for the model’s good fit and thus for performing inferences. Thus, the quality of the fitted model was guaranteed by the empirical properties shown by the residuals.

We calculated the relative risks (RR) and confidence intervals (CI) for each model fitted to the station’s neighborhood and for the overall region for each pollutant. Table 2 shows the estimated relative risks for an increase of 10 in the concentration of each of the air pollutants for the overall region. When we analyzed the patterns shown in this table, the gradients clearly evidenced the significant effects of PM10 on the treatments, for cumulative lags. For SO2, in the simple lag, the RR values were significant for lag0 and lag2. NO2, CO, and O3 did not show statistical significance in the RR for Greater Metropolitan Vitória.

Table 2
Relative risk (RR) of treatments for respiratory diseases in children under 6 years of age for 10µg/m3 increases in PM10 (μg/m3), SO2 (μg/m3), NO2 (μg/m3), O3 (μg/m3), and CO (μg/m3) in Greater Metropolitan Vitória, Espírito Santo State, Brazil, 2005 to 2010.

Table 3 presents a summary of the RR values for the effects of each pollutant whose calculations showed statistical significance (the results in bold print are those with the highest magnitudes). For the pollutants PM10 and SO2, the relative risks were higher when estimated for each station separately, compared to the RR for the entire Greater Metropolitan Vitória, that is, the adjusted model in which the covariables correspond to the average concentrations. As seen previously, the effects of NO2, CO, and O3 for the aggregate area (Greater Metropolitan Vitória) did not show statistical significance in the RR, which can be considered a spurious result, since the effects of these pollutants for the disaggregated sites led to highly significant relative risks.

Table 3
Percentage increase and 95% confidence interval (95%CI) in emergency pediatric treatments for respiratory problems. Greater Metropolitan Vitória, Espírito Santo State, Brazil, 2005-2010.

The entire analysis was based on the comparison of RR for the different pollutants and spatial points and compared to an average

Discussion

This study analyzed the effect of air pollution on children’s health in Greater Metropolitan Vitória by estimating relative risk in a GAM. Concentrations of PM10, SO2, NO2, and O3, even complying with Brazil’s prevailing legislation and WHO guidelines, showed a significant association with increases in emergency treatments for respiratory diseases in children 0 to 6 years of age, proving that there is no limit on these pollutants that is safe for human health 11. World Health Organization. Inheriting a sustainable world? Atlas on children's health and the environment. Geneva: World Health Organization; 2017..

The number treatments increased from March to June, corresponding to autumn to early winter. This expected increase is due to different factors such as low temperatures that predispose to aggravation of preexisting respiratory diseases, higher incidence of viral respiratory diseases, and increased concentration of primary pollutants determined by scarce rainfall and thermal inversion.

The neighborhood with the most children treated at the emergency departments was Enseada do Suá. This neighborhood is heavily impacted by mobile polluting winds, since it is close to roadways with intense traffic, including a bridge connecting the city of Vitória to Vila Velha, with up to 70,000 motor vehicles circulating per day, as well as fixed sources, mainly an industrial complex inside the urban grid. A recent study showed that the resident population in this neighborhood of Vitória reported the most discomfort and aggravation from sedimented dust 2929. Melo MM, Santos JM, Frere S, Reisen VA, Reis Jr. NC, Leite MDFS. Annoyance caused by air pollution: a comparative study of two industrialized regions. International Journal of Environmental and Ecological Engineering 2015; 9:182-7..

During the study period, the concentrations of pollutants varied between the different RAMQAr stations, which can be explained by each neighborhood’s specific polluting activities. The average pollutant levels in these stations are based on on-site data and are usually quite close and significantly different from the overall average. This empirical evidence shows that the local averages, E(μL t), (L = 1,...,8) differ from the Greater Metropolitan Vitória average E(μR t). This result is expected, since the process is not stationary in time, that is, E(μL t) ≠ E(μR t), for L=1,…,8, and justifies the proposed study in the sense of comparing the health outcome to the pollutants in spatially discriminated areas.

PM10 and SO2 were the pollutants that showed the most consistent association with the increase in treatments for respiratory problems. PM10 showed a substantial effect on treatments in all the lags and in nearly all the neighborhoods, with the highest magnitude in Jardim Camburi and Carapina. This effect pattern is consistent with other studies. Samoli et al. 1010. Samoli E, Nastos PT, Paliatsos AG, Katsouyanni K, Priftis KN. Acute effects of air pollution on pediatric asthma exacerbation: evidence of association and effect modification. Environ Res 2011; 111:418-24., in Greece, found that an increment of 10μg/m3 in PM10 and SO2 was associated with increases of 2.54% and 5.98% in the number of hospitalizations from respiratory diseases, respectively. In Itabira, Minas Gerais, Brazil, an increase of 10μg/m3 in the concentration of PM10 was associated with a 4% increase in emergency department treatments for respiratory diseases in children for lag0 and lag18. In Greater Metropolitan Vitória, Souza et al. 1313. Souza JB, Reisen VA, Franco GC, Ispány M, Bondon P, Santos JM. Generalized additive models with principal component analysis: an application to time series of respiratory disease and air pollution data. J R Stat Soc Series B Stat Methodol 2018; 67:453-80. observed a significant relationship between the concentrations of all the pollutants monitored by the RAMQAr (PM10, SO2, NO2, O3, and CO) and the number of emergency treatments for respiratory problems in children up to 6 years of age, through hybrid modeling, vector autoregression (VAR), principal components analysis (PCA), and GAM.

We observed high concentrations of PM10 in the Cariacica station, which can be explained by the station’s location inside the Espírito Santo Farm Produce Clearinghouse (CEASA/ES), which has intense internal truck traffic and is located close to two federal highways with heavy truck traffic (BR-262 and BR-101).

The higher effect of PM10 and SO2 in the Jardim Camburi, Carapina, and Laranjeiras neighborhoods can be explained by their proximity to heavily travelled roadways and the largest industrial complex in the state of Espírito Santo, both of which are potential sources of PM10 and SO2.

We found the highest RR for SO2 around the Jardim Camburi and Laranjeiras stations. In Enseada do Suá, the pollutant with the highest average was SO2, but due to the lack of available data, the RR was only calculated for PM10 in this neighborhood, representing a shortcoming in the study.

The adjustments to the model for each of the neighborhoods around the RAMQAr monitoring stations allowed better observation of the effect of all the pollutants on the number of emergency treatments for respiratory diseases. Evidence has shown that in time-series studies, measurements from a central station, compared to statistical or dispersion models or satellite-derived estimates, can be an adequate exposure metric for pollutants distributed homogeneously over the area 3030. Dionisio KL, Baxter LK, Burke LK, Özkaynak H. The importance of the exposure metric in air pollution epidemiology studies: when does it matter, and why? Air Qual Atmos Health 2016; 9:495-502., which does not apply to the current study.

Comparing the estimates when we simulated for the entire Greater Metropolitan Vitória to the analyses around the RAMQAr stations, some effects were only perceptible when the analyses were done in the sites disaggregated by neighborhood, signaling a greater effect compared to the estimate for the entire Greater Metropolitan Vitória. The explanation is that the analysis by the average for all the stations tends to smooth the data and thus decrease their variability, disguising some effects and underestimating air pollution’s effects on health. Wilson et al. 3131. Wilson JG, Kingham S, Pearce J, Sturman AP. A review of intraurban variations in particulate air pollution: implications for epidemiological research. Atmos Environ 2005; 39:6444-62. also concluded that consideration of intraurban variations in concentrations in epidemiological studies allows minimizing exposure errors and uncertainties in the relative risk, but they found that this only applies to long-term studies. However, using this approach we observed short-term effects of NO2, CO, and O3 on the number of treatments for respiratory diseases in young children, especially in the Jardim Camburi, Vila Velha-Centro, and Laranjeiras neighborhoods.

Despite the association between the various air pollutants and the risk of emergency treatment for respiratory diseases in children in Greater Metropolitan Vitória, we observed a downward trend in the levels of pollutants over time, which can be explained by better local air pollution control.

For future studies, other susceptible groups could be investigated in the same neighborhood in order to allow elaborating a complete picture of the acute effects of air pollution on the population’s health. As an alternative methodology, bootstrap techniques could be used to obtain intervals with the same precision but with a narrower sampling range. Another potential methodology is estimation of variance with heteroscedastic models. The GLARMA 3232. Davis RA, Dunsmuir WTM, Streett SB. Maximum likelihood estimation for an observation driven model for poisson counts. Methodol Comput Appl Probab 2005; 7:149-59. and PINAR 3333. Monteiro M, Scotto MG, Pereira I. Integer-valued autoregressive processes with periodic structure. J Stat Plan Inference 2010; 140:1529-41. models, with greater structural complexity, are statistical tools that could be addressed in these studies.

The consistency of associations and the magnitude of the observed effects in the neighborhoods, even in an environment with pollution levels that comply with limits set by regulatory agencies, are extremely relevant to public health. The findings provide backing for measures to minimize health risks, besides contributing to environmental and urban health planning and improvement of public policies.

Acknowledgments

The authors wish to acknowledge the support of various research funding agencies: Brazilian National Research Council (CNPq), Graduate Studies Coordinating Board (Capes), Espírito Santo State Research Foundation (FAPES), Vitoria Science and Technology Support Fund (FACITEC), Espírito Santo State Institute for the Environment (IEMA), Nossa Senhora da Glória Children’s Hospital (HINSG), and the Unimed Vitória Integrated Healthcare Center (CIAS).

References

  • 1
    World Health Organization. Inheriting a sustainable world? Atlas on children's health and the environment. Geneva: World Health Organization; 2017.
  • 2
    Holgate ST. Every breath we take: the lifelong impact of air pollution - a call for action. Clin Med (Lond) 2017; 17:8-12.
  • 3
    Cifuentes LA, Vega J, Köpfer K, Lave LB. Effect of the fine fraction of particulate matter versus the coarse mass and other pollutants on daily mortality in Santiago, Chile. J Air Waste Manag Assoc 2000; 50:1287-98.
  • 4
    Gouveia N, Mendonça GAS, Ponce de Leon A, Correia JEM, Junger WL, Freitas CU, et al. Poluição do ar e efeitos na saúde nas populações de duas grandes metrópoles brasileiras. Epidemiol Serv Saúde 2003; 12:29-40.
  • 5
    Daumas RP, Mendonça GAS, Ponce de León A. Poluição do ar e mortalidade em idosos no município do Rio de Janeiro: análise de série temporal. Cad Saúde Pública 2004; 20:311-9.
  • 6
    Bakonyi SM, Danni-Oliveira IM, Martins LC, Braga AL. Poluição atmosférica e doenças respiratórias em crianças na Cidade de Curitiba - PR. Rev Saúde Pública 2004; 38:695-700.
  • 7
    Nascimento LF, Pereira LA, Braga AL, Módolo MCC, Carvalho Jr. JA. Effects of air pollution on children's health in a city in southeastern Brazil. Rev Saúde Pública 2006; 40:77-82.
  • 8
    Braga ALF, Pereira LAA, Procópio M, André PA, Saldiva PHN. Associação entre poluição atmosférica e doenças respiratórias e cardiovasculares na cidade de Itabira, Minas Gerais, Brasil. Cad Saúde Pública 2007; 23 Suppl 4:S570-8.
  • 9
    Chen R, Chu C, Tan J, Cao J, Song W, Xu X, et al. Ambient air pollution and hospital admission in Shanghai, China. J Hazard Mater 2010; 181:234-40.
  • 10
    Samoli E, Nastos PT, Paliatsos AG, Katsouyanni K, Priftis KN. Acute effects of air pollution on pediatric asthma exacerbation: evidence of association and effect modification. Environ Res 2011; 111:418-24.
  • 11
    Almeida SP, Casimiro E, Calheiros J. Short-term association between exposure to ozone and mortality in Oporto, Portugal. Environ Res 2011; 111:406-10.
  • 12
    Souza JB, Reisen VA, Santos JM, Franco GC. Principal components and generalized linear modeling in the correlation between hospital admissions and air pollution. Rev Saúde Pública 2014; 48:451-8.
  • 13
    Souza JB, Reisen VA, Franco GC, Ispány M, Bondon P, Santos JM. Generalized additive models with principal component analysis: an application to time series of respiratory disease and air pollution data. J R Stat Soc Series B Stat Methodol 2018; 67:453-80.
  • 14
    Nascimento AP, Santos JM, Mill JG, Souza JB, Reis Júnior NC, Reisen VA. Associação entre concentração de partículas finas na atmosfera e doenças respiratórias agudas em crianças. Rev Saúde Pública 2017; 51:3.
  • 15
    Froes Asmus CI, Camara VM, Landrigan PJ, Claudio L. A systematic review of children's environmental health in Brazil. Ann Glob Health 2016; 82:132-48.
  • 16
    Castro HA, Hacon S, Argento R, Junger WL, Mello CF, Castiglioni Júnior N. Air pollution and respiratory diseases in the Municipality of Vitória, Espírito Santo State, Brazil. Cad Saúde Pública 2007; 23 Suppl 4:S630-42.
  • 17
    Instituto Estadual de Meio Ambiente e Recursos Hídricos. Relatório da qualidade do ar: Grande Vitória 2013. Vitória: Instituto Estadual de Meio Ambiente e Recursos Hídricos; 2013.
  • 18
    Namdeo A, Tiwary A, Farrow E. Estimation of age-related vulnerability to air pollution: assessment of respiratory health at local scale. Environ Int 2011; 37:829-37.
  • 19
    Zeger SL, Thomas D, Dominici F, Samet JM, Schwartz J, Dockery D, et al. Exposure measurement error in time-series studies of air pollution: concepts and consequences. Environ Health Perspect 2000; 108:419-26.
  • 20
    Santolim L. Rede otimizada para o monitoramento da qualidade do ar da Grande Vitória, Espírito Santo [Dissertação de Mestrado]. Vitória: Universidade Federal do Espírito Santo; 1991.
  • 21
    Junger WL. Análise, imputação de dados e interfaces computacionais em estudos de séries temporais epidemiológicas [Tese de Doutorado]. Rio de Janeiro: Universidade do Estado do Rio de Janeiro; 2008.
  • 22
    Hastie TJ, Tibshirani RJ. Monographs on statistics and applied probability 43: generalizes additive models. New York: Chapman & Hall; 1990.
  • 23
    Junger W, Leon AP. Ares: environment air pollution epidemiology. A library for time series analysis. Vienna: The Comprehensive R Archive Network; 2011. (CRAN Archives).
  • 24
    Baxter LA, Finch SJ, Lipfert FW, Yu Q. Comparing estimates of the effects of air pollution on human mortality obtained using different regression methodologies. Risk Anal 1997; 17:273-8.
  • 25
    Zhang H, Niu Y, Yao Y, Chen R, Zhou X, Kan H. The impact of ambient air pollution on daily hospital visits for various respiratory diseases and the relevant medical expenditures in Shanghai, China. Int J Environ Res Public Health 2018; 15:E425.
  • 26
    Chen R, Yin P, Meng X, Liu C, Wang L, Xu X. Fine particulate air pollution and daily mortality. a nationwide analysis in 272 Chinese cities. Am J Respir Crit Care Med 2017; 196:73-81.
  • 27
    Katsouyanni K, Schwartz J, Spix C, Touloumi G, Zmirou D, Zanobetti A, et al. Short term effects of air pollution on health: a european approach using epidemiologic time series data: the APHEA protocol. J Epidemiol Community Health 1996; 50 Suppl 1:S12-8.
  • 28
    Qu Y, Pan Y, Niu H, He Y, Li M, Li L, et al. Short-term effects of fine particulate matter on non-accidental and circulatory diseases mortality: a time series study among the elder in Changchun. PLoS One 2018; 13:e0209793.
  • 29
    Melo MM, Santos JM, Frere S, Reisen VA, Reis Jr. NC, Leite MDFS. Annoyance caused by air pollution: a comparative study of two industrialized regions. International Journal of Environmental and Ecological Engineering 2015; 9:182-7.
  • 30
    Dionisio KL, Baxter LK, Burke LK, Özkaynak H. The importance of the exposure metric in air pollution epidemiology studies: when does it matter, and why? Air Qual Atmos Health 2016; 9:495-502.
  • 31
    Wilson JG, Kingham S, Pearce J, Sturman AP. A review of intraurban variations in particulate air pollution: implications for epidemiological research. Atmos Environ 2005; 39:6444-62.
  • 32
    Davis RA, Dunsmuir WTM, Streett SB. Maximum likelihood estimation for an observation driven model for poisson counts. Methodol Comput Appl Probab 2005; 7:149-59.
  • 33
    Monteiro M, Scotto MG, Pereira I. Integer-valued autoregressive processes with periodic structure. J Stat Plan Inference 2010; 140:1529-41.

Publication Dates

  • Publication in this collection
    07 Oct 2019
  • Date of issue
    2019

History

  • Received
    24 July 2018
  • Reviewed
    20 Jan 2019
  • Accepted
    05 Apr 2019
Escola Nacional de Saúde Pública Sergio Arouca, Fundação Oswaldo Cruz Rio de Janeiro - RJ - Brazil
E-mail: cadernos@ensp.fiocruz.br