Introduction
Air pollution is an important environmental issue due to its serious impacts on human health. Air pollutants such as CO, NO
2, SO
2, O
3, and PM are mainly emitted from fossil fuels' combustion
1. According to the World Health Organization (WHO), approximately 3.5 million deaths occur due to air pollution annually. Many air pollutants such as CO, NO
2, SO
2, O
3, particulate matter (PM
10), heavy metals, radon, etc. have negative health effects
2. This rate increased by 60% from the year 2000 to 2010 and 65% of the deaths were attributed to the Asian countries
3. In the recent two decades, epidemiological studies showed that air pollution caused cardiovascular and respiratory problems, lung function defect, chronic bronchitis, and death
4. Moreover, lack of birth weight, pulmonary disorders, acute heart attacks, cardiovascular diseases, and other illnesses were related to urban air pollution
5. Among these diseases, outbreak of pulmonary diseases and daily hospital admissions due to respiratory disease were among the most significant indicator diseases associated with air pollution. In this regard, epidemiological and laboratorial observations proved specific roles of the suspended particles, O
3, NO
2, and SO
2 on human health
6. These factors have a significant role in stimulation of the respiratory system, oxidation, reduction in the lung capacity, and lung cell immunity against biological factors. As a result, a variety of diseases raise in the respiratory tract, which may lead to death
7. The consequent increased rate of hospital admission can lead to an additional burden on the health system. Similar studies in other cities of Iran showed that environmental pollutants affected the mortality rate and pulmonary diseases with an increasing trend in urban areas
8. Some studies investigated the relationship between air pollution and the number of patients admitted to hospitals, but this relationship has not been determined in Yazd city. Epidemiological and human studies confirmed the relationship between air pollutants and types of diseases
7,9,10. As a result, the authorities are recommended to try to reduce the disease burden related to air pollution. Yazd is located in the central desert regions of Iran, where the amount of suspended particles has increased with increase of desertification trend and dust streams. So, the quality of air is highly affected in this region by natural phenomena such as dust stream
11. The current study aimed to investigate the correlation between air pollution and respiratory disease by using time-series analysis in Yazd city, but according to similar studies did not investigate the confounder factors
10,12-14. However, statistical models are among the effective methods in health studies caused by air pollution. Nevertheless, statistical models such as regression can adjust the confounders; time series models can determine the trend, behavior, and time variations of the parameters. Other benefits of time series include modeling and forecasting of variation trend based on the past data. Considering the effect of air pollution, especially suspended particles, on outbreak of the respiratory diseases and the gap in the relationship between these two factors
13, this study was conducted. Moreover, variation trend of the air pollution and daily admission hospital rate were determined.
Materials and Method
Location study
Iran is a high plateau located at latitudes within the range of 25- 40° in an arid zone in the northern hemisphere of the Earth. The world’s most oppressive deserts are located in center of the large and closed cradle of Iran's plateau, at the center of Yazd province, which is surrounded by large deserts, such as Loot and Dashte Kavir
15.
Air pollution and health data
Daily average concentrations of the five significant pollutants (SO
2NO
2O
3COPM
10), were taken from two monitoring air pollutants stations of Yazd Environmental Protection Agency in 2017 (March 25 to September 22, 2017). All five pollutants was precisely monitored at both stations during one year and full information of the indicator pollutants was collected and analyzed by a time series method. Information about the daily hospital admission due to pulmonary diseases was collected from Shahid Sadoughi Hospital and Shahid Rahnemoon Hospital of Yazd. , the total number of pulmonary diseases were classified monthly.
Statistical analyses
Descriptive statistics (including mean, frequency, minimum, maximum, percent, and Standard) and chi-square test were run to analyze the epidemiological data. The data trend of time variation (the daily admission of patients versus the concentration of pollutants) was investigated and modeled within the studied time by time series using autoregression moving average (ARMA). Autocorrelation Function (ACF) and partial autocorrelation functions (PACF) are two suitable tools for determining the primary q (the order of the MA part) and p (order of AR part), respectively. The Root mean square error (RMSE) criterion derived from fitting of the model and the Box-Pierce test were used to determine the best model. Accordingly, AR and ARMA models were checked with different parameters and the suitable models were selected. Among the selected models, the model with lower RMSE was selected as the best model. The RMSE ranged from 3.82 to 10.81 (3.72, 3.8, 10.24, and 10.81); therefore, ARMA (1,1) was selected as best model. The lags between the daily admitted patients due to pulmonary disease and the concentration of air pollutants in =12, 24, 36, 48 hr. were evaluated using this model. In this study, the Poisson regression model was used to describe and predict the average daily admission of patients with pulmonary diseases. The Pearson correlation statistic (r) was also applied to determine the significant relationship between air pollutants and the daily admission of patients due to pulmonary disease. Moreover, regression analysis was used to describe the average parameters. Time series modeling using Minitab
17.1, Poisson regression analysis using Stata
14.2, and correlation test were run using SPSS
20. The significance level was less than 0.05.
Ethical issues
This study was conducted with approval of Shahid Sadoughi University of Medical Sciences and Health Services, Medical Ethics Committee. Code: IR. SSU. SPH.REC.1397.164.
Results
The average daily admission of patients due to pulmonary diseases was 9.7%, with an average age of 34.5. This study aimed to investigate the relationship between air pollution and daily hospital admissions due to respiratory disease using the time series model.
In this study, the average Air Quality Index (AQI)
16 was in a good state, only PM
10 was at a moderate level. The suspended particle concentration was in this range (in a good state, according to AQI index) during 96% of days. Moreover, the suspended particles' concentrations were unhealthy for sensitive, unhealthy, and at-risk individuals in 10%, 2%, and 1.6% of the total days, respectively.
The relationship between other pollutants and daily admission patients was insignificant (P > 0.05).
On the contrary, exposure to PM
10 was associated with death and cardio-respiratory diseases. However, lags (3 days) were observed between the PM concentration and its health effect.
Regression analysis
Since a disease outbreak depends on the interaction between several mechanisms and factors, regression analysis was used to determine the relationship between variables and their effects on each other
14. Regarding the specific accountability of the responses, Multiple Poison Regression was used (Table 1). In this study, the daily admission of the patients to hospital was considered as a dependent variable and the pollutants' concentration as an independent variable.
Age is one of the confounders adjusted in our regression models. According to the regression model, the relationship between O3 and the daily admission of patients to hospital due to pulmonary diseases was statistically significant (p < 0.049). In other words, 1 unit of increase in O3 concentration increases the risk of pulmonary diseases about 0.986. Variations of patients' daily admission to the hospital versus O3 concentrations are shown in Figure 1.
According to Table 2, the ARMA (1, 1) model was the best model for explaining the effect of air pollution on daily admission of patients to hospital due to respiratory disease, because its mean was less than the other models. Furthermore, these models are suitable for shock-induced systems such as the daily admission patient to hospital due to pulmonary disease. AR (3) and AR (4) models were statistically significant (p < 0.001), but their coefficient was not statistically significant.
In this study, time series' models were also evaluated in lag = 12, 24, 36, and 48. Results showed that exposure to air pollutants during 48 hours (lag = 4) was the most effective agent on the daily admission of patients to hospital due to pulmonary disease in Yazd. However, the statistical difference was significant in the modified Box-pierce for AR (3) and AR (4) models. According to Table 2, with increase of lags, chi-square and P increased in ARMA models (P < 0.02) considering the delay in admission of patients to hospital due to respiratory diseases. With regard to the results, ARMA (1, 1) was selected as the best model and trend of the predicted variation based on this model was shown by the residual plots (Figure 2).
As shown in figure 2, the fitting was perfectly suitable and aptness of the proposed model was confirmed. Adequacy data also showed that autocorrelation and partial autocorrelation were near to zero, which showed insignificancy of the residuals in all lags. In the following, the plot of observed versus fitted data during the study period is shown in figure 3.
According to Figure 3 and 4 difference between the fitted and observed values is low. However, due to the study short duration, the frequency of variation was not determined in time series. Considering the important effect of air pollutants on outbreak of the respiratory diseases and the depended-variations' trend of admission patient to hospital with air pollution, further investigations are required at higher periods. This relationship can be investigated regarding other diseases such as cardiovascular diseases.
Discussion
A total of 1770 patients were admitted in the studied hospitals due to respiratory diseases in Yazd. Among the patients, 808 were women and were 962 men. Only the annual average of PM
10 exceeded from the WHO standard
19, while the other factors were within the acceptable range. High concentration of PM was observed in many Asian cities due to the specific properties of this geographic reign
20. According to the meteorological data, Yazd City is located in the warm and dry area of Iran in a desert region, where many local and regional dust steams occur. The correlation between CO and SO
2 was more than the other factors (r
2 = 0.145, p < 0.05), because they are the primary contaminants emitted to air from combustion of fuel. Therefore, their correlation can be due to similarity of the emission source and the effective environmental conditions that involve these pollutants in the atmosphere
21. Recently, many studies have been conducted on potential mechanisms of the exposure to air pollutants that lead to respiratory disease. For example, air pollutants can cause damage to epithelial cells and macrophages, which may reduce the body's defense mechanisms against pulmonary infections or increase the reactivity of respiratory system cells. The severity of cell damages varied according to the concentration of contaminants, exposure duration, and type of pollutant. However, our results showed a significant and positive correlation between air pollution and prevalence of respiratory diseases
22. So, the outbreak of pulmonary illnesses, mortality rate, and admission of patient to hospital increased noticeably
18. In the current study, was tried to determine the correlation between air pollution in Yazd city and daily admission of the patient to hospital due to pulmonary diseases using the Pearson test. The findings indicated that health effects of air pollution on human health can be observed as the rate of patients' daily admission to hospital.
The correlation between CO and SO
2 was more than other factors (r
2 = 0.145, p < 0.05). Results showed that the correlation between daily admission of patients to hospital and the concentration of NO
2 was higher than other pollutants (r
2 = 0.4) (p < 0.001), because NO
2 stimulated the pulmonary system in higher depth. Moreover, the secondary particulate matters, such as PANs are generated due to the influence of NO
2 on allergic reactions and acute disease
23,24. In other studies, the 4% NO
2 increased acute respiratory infections
22. Moreover, toxicological evidence showed that NO
2 had an indirect effect on the pulmonary illnesses and infections, because it may increase the sensitivity of respiratory cell to bacterial and viral infections
25. Moreover, this correlation between pulmonary illnesses and NO
2 may be due to the photolysis process of the No
x cycle. This is due to the fact that NO is emitted from the moving resources as primary pollutant and is involved in atmospheric reactions related to light. This reaction was improved in hot and dry climate such as Yazd. On the other hand, presence of a high level of humidity can limit it
26. Moreover, the correlation between patients' admission and CO and PM
10 was negative (p < 0.05) since the health effect and toxicity of CO are related to indoor space and climatic pattern that generated COHb and reduced the transportation capacity of oxygen in the blood
27. Therefore, CO is not related to the mortality rate and hospital admission due to the pulmonary disease. PM
10 can penetrate into the human pulmonary system and create inflammatory and allergenic responses. Moreover, various chemical and microbial contaminants aggregated on its surface, which cause synergistic effects
25.
Considering that these pollutants do not have protective effects on the health, application of models able to determine the co-effects and the variation trend during time is necessary
13,14.
About 92% of the diseases are multifactorial and depend on many factors, such as socio-economic, genetic, and other factors. Moreover, modeling and determining the trend of changes during the period of study were used by time-series models with regard to the synergistic effects of contaminants on creating health complications.
A significant difference was found between the observed and fitted values, which indicates inappropriateness of used models in this study. On the other hand, P in ARMA grade 1 and 3 was statistically significant, indicating that the model was suitable.
With regard to the delayed effects of air pollution on health, identification of the health risks is significantly necessary. However, the findings are controversial regarding the effect of lags on the relationship between air pollutants and health effects
10.
Conclusion
In this study, the relationship between air pollution and the number of pulmonary patients in Yazd was studied.
- The air quality index for Yazd was good; because the average concentration of air pollutants was less than the EPA and WHO standards (except PM10 that exceeded from this limit during 4% of the study days).
- The daily admission of patients to hospital due to pulmonary disease was highly correlated with NO2 and O3, while its correlation with CO, PM10, and SO2 was insignificant. This was due to the fact that NO2 and O3 were identified as oxidation factors stimulating the respiratory system.
- According to the regression model and after adjusting the age as a confounding factor, O3 was the most effective pollutant on daily admission of patients to hospital due to pulmonary disease.
- The variation trend of air pollution and daily admission of patient to hospital due to pulmonary disease were determined using time series analyses. According to this model, ARMA (1,1) was the best model showing better results in lag = 48 hours since the air pollutants, especially PMs delayed the effects on human health. Moreover, they are able to weaken the immunity system and develop other diseases.
Acknowledgments
Thanks are owed to the Department of Environmental Health Engineering, Abarkouh paramedical school, Yazd Shahid Sadoughi University of Medical Sciences for their help in conducting this study.
Funding
This study was funded by Yazd, Shahid Sadoughi University of Medical Sciences, Iran.
Conflict of interest
The authors declare no conflict of interest regarding the publication of this research.
This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work for commercial use.
References
1. Abdolahnejad A, Jafari N, Mohammadi A, et al. Mortality and morbidity due to exposure to ambient NO2, SO2, and O3 in Isfahan in 2013–2014. Int J Prev Med. 2018;9:11.
2. Farrokhzadeh H, Jafari N, Sadeghi M, et al. Estimation of spatial distribution of PM10, lead, and radon concentrations in Sepahanshahr, Iran using Geographic Information System (GIS). Journal of Mazandaran University of Medical Sciences. 2018;27(159):84-96.
3. Brunekreef B, Holgate ST. Air pollution and health. The lancet. 2002;360(9341):1233-42.
4. Fattore E, Paiano V, Borgini A, et al. Human health risk in relation to air quality in two municipalities in an industrialized area of Northern Italy. Environ Res. 2011; 111(8):1321-7.
5. Kermani M, Aghaei M. Estimation of cardiovascular death, myocardial infarction and chronic obstructive pulmonary disease (COPD) attributed to PM and SO2 in the air of Tehran metropolis. Journal of Research in Environmental Health. 2016;2(2):116-26.
6. Bahrami Asl F, Kermani M, Aghaei M, et al. Estimation of diseases and mortality attributed to NO2 pollutant in five metropolises of Iran using AirQ model in 2011-2012. Journal of Mazandaran University of Medical Sciences. 2015;24(121): 239-49.
7. Szyszkowicz M, Kousha T, Castner J, et al. Air pollution and emergency department visits for respiratory diseases: A multi-city case crossover study. Environ Res. 2018;163:263-9.
8. Khaniabadi YO, Sicard P, Takdastan A, et al. Mortality and morbidity due to ambient air pollution in Iran. Clin Epidemiol Glob Health. 2019;7(2):222-7.
9. Geravandi S, Goudarzi G, Soltani F, et al. Sulfur dioxide pollutant and its effects on disease incidence and death among the citizens of Bushehr city. Iranian South Medical Journal. 2016;19(4):598-607.
10. Maturo A, Moretti V. Sociological theories on air pollution: between environmental justice and the risk society approach. Clinical Handbook of Air Pollution-Related Diseases: Springer; 2018: 603-20.
11. Modarres R, Sadeghi S. Spatial and temporal trends of dust storms across desert regions of Iran. Nat Hazards. 2018;90(1):101-14.
12. Zhang Y, Ni H, Bai L, et al. The short-term association between air pollution and childhood asthma hospital admissions in urban areas of Hefei City in China: A time-series study. Environ Res. 2019;169:510-6.
13. Lin H, Tao J, Qian ZM, et al. Shipping pollution emission associated with increased cardiovascular mortality: A time series study in Guangzhou, China. Environ Pollut. 2018;241:862-8.
14. Fernández-Niño JA, Astudillo-García CI, Rodríguez-Villamizar LA, et al. Association between air pollution and suicide: a time series analysis in four Colombian cities. Environmental Health. 2018;17(1):47.
15. Keshtkaran P. Harmonization between climate and architecture in vernacular heritage: A case study in Yazd, Iran. Procedia Eng. 2011;21:428-38.
16. James KA, Strand M, Hamer MK, et al. Health services utilization in asthma exacerbations and PM10 levels in rural Colorado. Ann Am Thorac Soc. 2018;15(8):947-54.
17. Swartzlander MD, Blakney AK, Amer LD, et al. Immunomodulation by mesenchymal stem cells combats the foreign body response to cell-laden synthetic hydrogels. Biomaterials. 2015;41:79-88.
18. Karimi B, Samadi S. Mortality and hospitalizations due to cardiovascular and respiratory diseases associated with air pollution in Iran: A systematic review and meta-analysis. Atmos Environ. 2019;198:438-47.
19. Zeng Q, Ni Y, Li G, et al. The quantitative assessment of the public excess disease burden advanced by inhalable particulate matter under different air quality standard targets in Tianjin, China. Environ Sci Pollut Res. 2019;26(7):6931-8.
20. Sousa S, Pires J, Martins E, et al. Short-term effects of air pollution on respiratory morbidity at Rio de Janeiro-Part II: health assessment. Environment international. 2012;43:1-5.
21. Galanopoulos C, Yan J, Li H, et al. Impacts of acidic gas components on combustion of contaminated biomass fuels. Biomass Bioenergy. 2018;111:263-77.
22. Shakerkhatibi M, Dianat I, Jafarabadi MA, et al. Air pollution and hospital admissions for cardiorespiratory diseases in Iran: artificial neural network versus conditional logistic regression. Int J Environ Sci Technol. 2015;12(11):3433-42.
23. Karoui A, Crochemore C, Mekki M, et al. Comparison of cardiovascular responses to nitrogen dioxide and diesel exhausts. An experimental study with controlled exposures to air pollutants. Archives of Cardiovascular Diseases Supplements. 2015;7(2):142.
24. Brugha R, Grigg J. Urban air pollution and respiratory infections. Paediatr Respir Rev. 2014;15(2):194-9.
25. Gadsby N, McHugh M, Russell C, et al. Development of two real-time multiplex PCR assays for the detection and quantification of eight key bacterial pathogens in lower respiratory tract infections. Clin Microbiol Infect. 2015;21(8):788.
26. Zarei ME, Ashkezari SFM, Yari M. The investigation of the function of the central courtyard in moderating the harsh environmental conditions of a hot and dry climate (Case study: City of Yazd, Iran). Spatium. 2017(38):1-9.
27. Sobhakumari A, Poppenga RH, Pesavento JB, et al. Pathology of carbon monoxide poisoning in two cats. BMC Vet Res. 2018;14(1):67.