key: cord-0807394-3n4ixehq authors: Gao, Hao-Yu; Liu, Xiao-Ling; Lu, Ya-Ke; Liu, Yu-Hong; Hu, Li-Kun; Li, Yan-Ling; Feng, Xiao-Dong; Yan, Yu-Xiang title: Short-term effects of gaseous air pollutants on outpatient visits for respiratory diseases: a case-crossover study in Baotou, China date: 2022-02-27 journal: Environ Sci Pollut Res Int DOI: 10.1007/s11356-022-19413-6 sha: 9f1fcc00985c8ce8243f1e3a33d25e7d1c4a7c5b doc_id: 807394 cord_uid: 3n4ixehq Air pollution is a major public health problem throughout the world. Although there have been several studies in this field, most of them have focused on particulate matter and only covered a few key cities. This study aimed to assess a potential association between exposure to gaseous air pollutants and outpatient visits for respiratory diseases in Baotou, China. Daily outpatient visits for respiratory diseases and daily averages of air pollutants and meteorological parameters from 2015 to 2020 were obtained. Time-stratified case-crossover design and restricted cubic splines were used to perform the analyses. Stratified analyses were performed in different hospital departments and districts. Significant association between the concentrations of air pollutants and outpatient visits for respiratory diseases was observed. The odds ratios of outpatient visits for respiratory diseases associated with per 10 μg/m(3) increases in concentrations of NO(2) and SO(2), and per 10 mg/m(3) increases in concentrations of CO were 1.033 (95% CI: 1.018 to 1.049), 0.965 (95% CI: 0.954 to 0.976), and 1.038 (95% CI: 1.006 to 1.071), respectively. Short-term exposure to NO(2), SO(2), and CO was positively associated with outpatient visits for respiratory diseases, with stronger effects among children. The relationship between O(3) and respiratory diseases varied at different concentrations. Ambient air pollution has been one of the most important public health concerns around the world (Cohen et al. 2017) . In October 2013, the International Agency for Research on Cancer (IARC) classified outdoor air pollution as carcinogenic to humans (Loomis et al. 2013 ). Contributing to 3.7 million deaths worldwide annually, air pollution has been identified as the worldwide largest single environmental health risks (Amann et al. 2020; Nwanaji-Enwerem et al. 2016 ). Ambient air pollutants include gaseous air pollutants and particulate matter. At present, the correlation between particulate matter and human health has been widely studied, which suggests that gaseous air pollutants are also closely associated with human health (Katsouyanni et al. 1997) . In recent years, as a country embarking on urbanization and industrialization, China has been struggling immensely with severe air pollution. "Brief Introduction of National Ecological Environment Quality in 2019" pointed out that among 337 prefecture-level cities in China, the proportion of air quality exceeding the standard is as high as 53.4% (China 2020). Respiratory diseases are diseases of the airways and other structures of the lung, among which chronic respiratory diseases are the leading causes of morbidity and mortality worldwide (Soriano et al. 2020 ). The respiratory system is exposed to the external environment directly, making it more susceptible to the effects and influences of the surrounding environment. Subsequently, air pollution has been widely acknowledged as a major influence and exacerbating factor in various respiratory diseases (Mokoena et al. 2019) . Despite the continuous effort and significant improvements in ambient air quality over the past decade, respiratory diseases in China remain increasingly serious (Bao et al. 2020) . Recently, research on the association between air pollution and respiratory diseases has received increased attention. But these studies were mainly focused on particulate matter, research on gaseous air pollutants has been relatively sparse Zhang et al. 2020 ). In addition, previous research mainly covered a few key cities or urban agglomerations (Chang et al. 2020; Phung et al. 2016; Zhang et al. 2019a) . As an important industrial city in Inner Mongolia, Baotou has been experiencing serious air pollution over the past few years. Nevertheless, limited studies have been conducted to evaluate the relationships between exposure to air pollutants and outpatient visits for respiratory diseases in this city. Given the above, more related studies are needed to be initiated there. The outpatient services are open to all diseases with various severities and are not restricted by bed availability. Thus, outpatient visits could better reflect the medical need of most people compared to emergency department visits and hospitalizations . The aim of this study was to assess the potential associations between exposure to gaseous air pollutants and outpatient visits for respiratory diseases in Baotou. Baotou is a city in the central and western parts of Inner Mongolia. As an important industrial base in China, Baotou's pillar industries are manufactured of steel and smelting of nonferrous metals. Baotou has a middle temperate semiarid continental monsoon climate, leading to lack of precipitation and seasonal temperature fluctuations. Compared with most of the cities in northern China, Baotou has a longer heating period and relies on coal-fired heating. Through the adoption of a series of control measures, the air quality in Baotou has gradually improved over recent years. However, the improvement in air quality is insignificant, and air pollution problem in Baotou is still fraught. Daily records of outpatient visits were collected from two general hospitals and two community hospitals in Baotou from January 1, 2015, to December 31, 2020. Four hospitals were selected from Donghe District and Qingshan District of Baotou. These two districts were defined as seriously polluted area and lightly polluted area according to the Ambient Air Quality Index, respectively. Total outpatient visits in all of these four hospitals had not reached saturation. Respiratory diseases were included in this study and were coded according to the International Classification of Diseases, Tenth Revision (ICD-10), chapter X: Diseases of respiratory system (J00-J99). Information about patients such as consultation date, hospital, and outpatient department was recorded. Nitrogen dioxide (NO 2 ), sulfur dioxide (SO 2 ), carbon monoxide (CO), ozone (O 3 ), inhalable particulate matter (PM 10 ), and fine particulate matter (PM 2.5 ) measurements in Baotou during January 1, 2015, to December 31, 2020, were retrieved from Baotou Environmental Monitoring Station. The Environmental Monitoring Station comprised eight monitoring sites in six different districts, and air pollutants data from two monitoring sites in Donghe District and two monitoring sites in Qingshan District was included in our study. For each district, we obtained daily mean concentrations of NO 2 , SO 2 , and CO and maximum daily 8-h mean concentrations of O 3 . Meteorological data, including mean temperature, relative humidity, atmospheric pressure, and wind speed, were obtained from Baotou Meteorological Bureau. A time-stratified case-crossover design was used to assess the acute effects of exposure to gaseous air pollutants on outpatient visits for respiratory diseases. With each case acting as its own control, case-crossover design is a variant of matched case-control study (Maclure 1991) and has been broadly adopted to investigate the associations of short-term environmental exposures with mortality and morbidity outcomes (Carracedo-Martinez et al. 2010) . The unique advantages of case-control study include controlling the effects of day of week, season, and time-invariant confounders such as individual characteristics. In this study, case period was defined as the day of outpatient visit. Control periods were identified by matching the day of the week within the same year and month in order to eliminate the potential confounding effects of long-term trend. Conditional logistic regression model was used to perform the analyses in our study. We estimated the odds ratios (ORs) and 95% confidence intervals (CIs) of outpatient visits for respiratory diseases associated with a 10.0 μg/ m 3 increase in daily concentrations of NO 2 , SO 2 , and O 3 , and a 1.0 mg/m 3 increase in daily concentrations of CO. Since the effects of environmental risk factors are sometimes not observed instantly but are delayed, a lag analysis was conducted to observe the delayed effects of gaseous air pollutants on outpatient visits. For the lag analysis, we applied both single-day lag structures (lag0 to lag3) and multi-day lag structures (lag01 to lag03). For example, lag1 corresponds to the air pollutant concentrations measured at 1 day prior to outpatient visit, while lag01 corresponds to the air pollutant concentrations averaged from current day's and lag1 day's measurements. To adjust for nonlinear confounding effects of meteorological factors, we applied a restricted cubic spline (RCS) function with 3 degrees of freedom (df) for mean temperature, relative humidity, atmospheric pressure, and wind speed. Public holidays were controlled through the use of a binary variable (coded as public holiday = 1 and no public holiday = 0) in the model. All results were from models including mean temperature, relative humidity, atmospheric pressure, wind speed, and public holidays. Because air pollutants possibly interact with each other and have a combined effect on health outcomes, we carried out analyses using multi-pollutant models. A principal component was introduced into the multi-pollutant models to obtain more reliable results without the impacts of collinearity between various air pollutants. Subgroup analyses stratified by types of hospitals and departments (department of internal medicine in general hospitals, department of pediatrics in general hospitals, general outpatient in community hospitals), and the district where the hospital is located (Donghe District and Qingshan District). Sensitivity analyses were conducted to check the robustness of the results. First, we used 4 df for all meteorological factors. Second, we excluded the data during the COVID-19 pandemic (data in 2020) for further analyses. All analyses were conducted using R software (version 4.1.0). The "survival" package and "splines2" package were used for conditional logistic regression analysis. And "rms" package was used to fit the restricted cubic splines. All statistical tests were two-sided, and the statistical significance was defined as p value less than 0.05. A total number of 168,201 patients visited the outpatient department for respiratory diseases in our participating hospitals between the years 2015 to 2020. The mean number of daily outpatient visits was 76.7, including 31.0 for the department of internal medicine in general hospitals, 32.1 for the department of pediatrics in general hospitals, and 13.7 for the general outpatient in community hospitals ( Table 1) . The daily average concentrations of air pollutants and meteorological variables in two districts between the years 2015 to 2020 are presented in Table 2 . The means (standard deviation, SD) of air pollutants were 39.2 (16.8) μg/m 3 for NO 2 , 26.7 (18.3) μg/m 3 for SO 2 , 1.2 (0.7) mg/m 3 for CO, and 93.8 (41.7) μg/m 3 for O 3 , respectively. Means (SD) of meteorological variables were 8.0 (12.7) °C for temperature, 56.0 (16.4)% for relative humidity, 902.3 (6.9) kPa for atmospheric pressure, and 2.9 (1.2) m/s for wind speed. In Donghe District, the means (SD) of air pollutants were 38.8 (16.8) μg/m 3 for NO 2 , 27.2 (18.2) μg/m 3 for SO 2 , 1.2 (0.7) mg/m 3 for CO, and 94.4 (40.9) μg/m 3 for O 3 , respectively. In Qingshan District, the means (SD) of air pollutants were 39.6 (16.9) μg/m 3 for NO 2 , 26.2 (18.4) μg/m 3 for SO 2 , 1.2 (0.7) mg/m 3 for CO, and 93.3 (42.6) μg/m 3 for O 3 , respectively. Spearman correlation coefficients for correlations among the air pollutants and meteorological parameters are presented in Table 3 . NO 2 , SO 2 , and CO showed positive correlations with each other (r = 0.79 to 0.83, p < 0.01) and particulate matter (r = 0.45 to 0.62, p < 0.01), and negative correlations with O 3 (r = − 0.52 to − 0.39, p < 0.01). Correlations were also observed between air pollutants and meteorological parameters. NO 2 , SO 2 , and CO were positively correlated with relative humidity (r = 0.11 to 0.36, p < 0.01) and atmospheric pressure (r = 0.24 to 0.36, p < 0.01), but negatively correlated with temperature (r = − 0.54 to − 0.38, p < 0.01) and wind speed (r = − 0.54 to − 0.28, p < 0.01). Conversely, O 3 was positively correlated with temperature (r = 0.84, Table 4 shows the ORs (95% CI) of outpatient visits for total respiratory diseases at various exposure days associated with per 10 μg/m 3 increase in exposure to NO 2 , SO 2 , and O 3 , and per 1 mg/m 3 increases in exposure to CO. Positive results were observed for NO 2 among all lag structures, for SO 2 at lag1-3 and lag02-03, and for CO at lag0-2 and lag01-03, where the highest was at lag0, lag03, and lag01 respectively. Negative results were observed for O 3 at lag1 and lag01-02. Figure 1 summarizes the subgroup of different departments' analyses results for lag days of 0-3 and 01-03. We found NO 2 was associated with outpatient visits for respiratory diseases in all subgroups and were associated most strongly among children (OR = 1.005 to 1.027). SO 2 was found associated with outpatient visits for respiratory diseases in department of internal medicine in general hospitals at lag1 and lag3, in department of pediatrics in general hospitals at lag1-2 and lag02-03, and in general outpatient in community hospitals at lag3. CO was found associated with outpatient visits for respiratory diseases in department of internal medicine in general hospitals at lag1-3 and lag01-03 and in department of pediatrics in general hospitals at lag1 and lag01-02. Unlike other gaseous pollutants, O 3 was found negatively associated with outpatient visits in department of internal medicine in general hospitals at lag1 and lag01-03. Figure 2 summarizes the subgroup of different districts' analyses results for lag days of 0-3 and 01-03. NO 2 was found more associated with outpatient visits for respiratory diseases in Qingshan District, which had larger ORs and narrower confidence intervals. SO 2 and CO were found more associated with outpatient visits for respiratory diseases in Donghe District. Also, O 3 were found negatively associated with outpatient visits in Qingshan District in most lag structures. Figure 3 shows the exposure-response relationships of outpatient visits for respiratory diseases with 4-day (lag0-3) moving average concentrations of NO 2 , SO 2 , CO, and O 3 in different pollutant models. The exposure-response curves suggested an approximately linear increase in outpatient visits risks with daily changes in NO 2 , SO 2 , and CO in both single-pollutant models and multi-pollutant models. For O 3 , there was a linear increase in outpatient visits risks at higher O 3 concentrations, with a threshold of around 89 μg/m 3 for a 8-h average exposure. Fig. 4 The association between air pollutants' concentrations and outpatient visits for respiratory diseases by adjusting the df (s) of meteorological factors The results of sensitivity analyses are shown in Fig. 4 and Fig. 5 . The main findings of our study did not change after adjusting the degrees of freedom of temperature, relative humidity, atmospheric pressure, and wind speed. The effects of SO 2 on outpatient visits for respiratory diseases decreased after excluding the data in 2020. In this time-stratified case-crossover study, exposures to NO 2 , SO 2 , and CO were found positively associated with outpatient visits for respiratory diseases. Besides, an approximately linear exposure-response relationship of outpatient visits for respiratory diseases at higher O 3 levels was observed. Short-term effects of gaseous air pollutants on outpatient visits for respiratory diseases varied slightly among different departments and districts. Globally, air pollution has become a serious public health concern (Chen et al. 2020), and various epidemiological studies report that short-term exposure to air pollutants may cause adverse health effects (Guan et al. 2016) . Our results suggested an approximately linear exposure-response relationship of respiratory diseases with short-term NO 2 , SO 2 , and CO exposures, and they were positively associated with outpatient visits for respiratory diseases. These results were consistent with the previous findings worldwide, in spite of heterogeneity in effect estimates across various regions of the world Tramuto et al. 2011; Yang et al. 2021) . Evidence of the effect of air pollutants on respiratory health has been accumulating (Dockery et al. 2013) . Putative biological mechanisms linking air pollution to respiratory diseases involve direct effect of air pollutants on the lung and indirect effects that mediated by oxidative stress and pulmonary inflammation (Zhang et al. 2019b) . Other studies have shown that meteorological factors have an effect on the concentrations of air pollutants, which is similar with some other studies Zhang et al. 2015) . The concentration of O 3 had a positive association with temperature and a negative association with relative humidity, because sunshine can promote the production process of O 3 , and relative humidity might have negative effects on sunshine duration (Jiang et al. 2020) . The concentrations of NO 2 , SO 2 , and CO were negatively correlated with wind speed and temperature, indicating that horizontal dispersion plays an important role in modulating air pollutants concentration (Zhang et al. 2015) . In China, low temperature might promote coal burning and heating, and thus, low temperature has a negative impact on air pollutants concentrations . In view of the fact that meteorological factors can influence air pollutants, thereby impacting health, we introduced meteorological factors into each models. In comparison to studies carried out in other cities, results of this study had similar trends. A study done in Fuzhou, China, associated exposure to NO 2 and outpatient visits for respiratory diseases with estimated OR of 1.051 (95% CI: 1.031 to 1.070) (Jiang et al. 2020) , while ours was 1.033 (95% CI: 1.018 to 1.049). Another study in Southeastern China associated exposure to SO 2 and outpatient visits for respiratory diseases with estimated relative risk of 1.035 (95% CI: 1.029 to 1.041) at lag1 (Mo et al. 2018) , while our OR was 1.008 (95% CI: 1.003 to 1.013). At the same time, our results showed that, although SO 2 was positively correlated with outpatient visits for respiratory diseases in several lag structures, they were negatively related to each other at lag0. Similar results to ours have been reported in other studies (Dong et al. 2021) , and one of the possible reasons is that irritating gaseous air pollutants including SO 2 can cause inflammation in the respiratory tract after several days (Chang et al. 2020) . Quite a few studies stated that the effect of CO on respiratory diseases was positive and has statistically significant effect (Chang et al. 2020) ; we also obtained the same result (OR = 1.038, 95% CI: 1.006 to 1.072). Unlike other gaseous air pollutants, we found O 3 was negatively correlated with outpatient visits for respiratory diseases at low concentration while positively correlated with that at high concentration. This might be due to the negative association between O 3 and other gaseous air pollutants, which were more strongly associated with respiratory diseases. Meanwhile, high concentration of O 3 led to bronchial inflammation and hyper-responsiveness, and this further led to respiratory diseases (Kim et al. 2020) . Single-day lag structures (lag0-3) and multi-day lag structures (lag01-03) were adopted in our study to investigate the lag relationship of gaseous air pollutants and outpatient visits for respiratory diseases. The effect estimates of multi-day lag structures were slightly larger and more stable than those of single-day lag structures for all four gaseous air pollutants. These findings are consistent with some other studies, which suggest that accumulated exposure of air pollutants increases the risk of respiratory diseases (Lin et al. 2017; Mokoena et al. 2019) . We also observed that the OR value of SO 2 was largest at lag2, while the OR values of NO 2 and CO were largest at lag0. The subgroup analyses showed that the positive associations of outpatient visits for respiratory diseases with NO 2 and O 3 appear to be stronger among children, and these were in accordance with previous studies (Mo et al. 2018) . These may suggest that children will be more sensitive than adults to air pollutants due to their narrower airways, developing lung functions, immature immune system, and longer staytime outdoors . Additionally, some multicity studies in China reported that stronger associations between air pollutants and respiratory diseases among areas with lowest air pollutants concentrations (Tian et al. 2019; Yang et al. 2021; Zhu et al. 2017 ). However, we did not find this trend in this study. As shown in Table 2 , the average concentrations of SO 2 and CO in Donghe District were higher than those in Qingshan District. Contrary to expectations, these two kinds of gaseous pollutants were found more associated with outpatient visits for respiratory diseases in Donghe District. These inconsistent results could be attributed to differences in both location and patient composition between hospitals in different districts. Results from sensitivity analyses by adjusting the degrees of freedom of meteorological parameters were essentially unchanged from primary analyses. After excluding the data in 2020, we found that effects of SO 2 on outpatient visits for respiratory diseases decreased. This may be attributed that more attention has been given to respiratory symptoms during the COVID-19 pandemic. And people are more willing to go to the hospital seeking treatment when they experience respiratory symptoms. There are several strengths of this study. The time-stratified case-crossover design using in our study could effectively control for confounding due to individual-level characteristics. Moreover, respiratory diseases in our study were diagnosed and classified based on ICD-10 code in order to eliminate the possibility misclassification of cases. We are aware of some limitations of our study. Air pollutant data in our study were based on district-wide averages rather than personal exposure measurements. In addition, the concentrations of air pollutants were collected on the day that patients visited the hospital not on the day of the onset of disease. Besides, we could not do subgroup analysis in different age groups and in subgroups of ICD-10, and analyze the differences in patient populations between four hospitals because these data are unavailable. These limitations might affect the results of our study to a certain extent. NO 2 , SO 2 , and CO were positively associated with outpatient visits for respiratory diseases. Stronger associations between them were among children. An approximately linear exposure-response relationship of outpatient visits for respiratory diseases with O 3 was observed at higher O 3 concentrations, with a threshold of around 89 μg/m 3 for an 8-h average exposure. Reducing global air pollution: the scope for further policy interventions Modeling the influence of nonclinic visits on the transmission of respiratory diseases Case-crossover analysis of air pollution health effects: a systematic review of methodology and application Ambient air pollution and daily hospital admissions for respiratory system-related diseases in a heavy polluted city in Northeast China impairs macrophage functions to exacerbate pneumococcus-induced pulmonary pathogenesis Acute respiratory response to individual particle exposure (PM(1.0), PM(2.5) and PM(10)) in the elderly with and without chronic respiratory diseases Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study Effect of air pollution control on mortality and hospital admissions in Ireland Revalue associations of short-term exposure to air pollution with respiratory hospital admissions in Lanzhou, China after the control and treatment of current pollution Impact of air pollution on the burden of chronic respiratory diseases in China: time for urgent action Temporal cross-correlations between air pollutants and outpatient visits for respiratory and circulatory system diseases in Fuzhou Short-term effects of ambient sulphur dioxide and particulate matter on mortality in 12 European cities: results from time series data from the APHEA project. Air Pollution and Health: a European Approach Health effects of ozone on respiratory diseases Air pollution characteristics in China during 2015-2016: spatiotemporal variations and key meteorological factors Hourly peak concentration measuring the PM 2.5 -mortality association: results from six cities in the Pearl River Delta study Short-term effects of ambient air pollution on pediatric outpatient visits for respiratory diseases in Yichang city The carcinogenicity of outdoor air pollution The case-crossover design: a method for studying transient effects on the risk of acute events Acute effects of air pollution on respiratory disease mortalities and outpatients in Southeastern China Ministry of Ecology and Environment of the People's Republic of China. Brief introduction of national ecological environment quality in 2019 Ambient air pollution and respiratory mortality in Xi'an, China: a time-series analysis Long-term ambient particle exposures and blood DNA methylation age: findings from the VA normative aging study Air pollution and risk of respiratory and cardiovascular hospitalizations in the most populous city in Vietnam Prevalence and attributable health burden of chronic respiratory diseases, 1990-2017: a systematic analysis for the Global Burden of Disease Study Fine particulate air pollution and adult hospital admissions in 200 Chinese cities: a time-series analysis Urban air pollution and emergency room admissions for respiratory symptoms: a case-crossover study in Palermo The impact of outdoor air pollutants on outpatient visits for respiratory diseases during 2012-2016 in Jinan What Are the Driving Forces of Urban CO2 Emissions in China? Int J Environ Res Public Health, A refined scale analysis between national and urban agglomeration levels Short term effects of air pollutants on hospital admissions for respiratory diseases among children: a multi-city time-series study in China Relationships between meteorological parameters and criteria air pollutants in three megacities in China Association of particulate matter air pollution and hospital visits for respiratory diseases: a time-series study from Impacts of event-specific air quality improvements on total hospital admissions and reduced systemic inflammation in COPD patients Shortterm effects of ambient PM1 and PM2.5 air pollution on hospital admission for respiratory diseases: Case-crossover evidence from Shenzhen, China Short-term effects of ambient air pollution and childhood lower respiratory diseases Acknowledgements We sincerely thank Baotou Center for Disease Control and Prevention for data collection and management.Author contribution Hao-Yu Gao performed data analysis and drafted the manuscript. Xiao-Ling Liu participated in the data collection. Ya-Ke Lu and Yu-Hong Liu participated in the data cleaning. Li-Kun Hu and Yan-Ling Li participated in data analysis. Xiao-Dong Feng and Yu-Xiang Yan developed the study design and edited the manuscript. All authors read and approved the final version of the manuscript. The data that support the findings of this study are available from the Baotou Center for Disease Control and Prevention, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of Baotou Center for Disease Control and Prevention. Ethics approval and consent to participate Not applicable. Competing interest The authors declare no conflict of interest.