- Open Access
The effects of traffic-related air pollutants on chronic obstructive pulmonary disease in the community-based general population
Respiratory Research volume 22, Article number: 217 (2021)
Previous studies have shown inconsistent results regarding the impact of traffic pollution on the prevalence of chronic obstructive pulmonary disease (COPD). Therefore, using frequency matching and propensity scores, we explored the association between traffic pollution and COPD in a cohort of 8284 residents in a major agricultural county in Taiwan.
All subjects completed a structured questionnaire interview and health checkups. Subjects with COPD were identified using Taiwan National Health Insurance Research Databases. A hybrid kriging/LUR model was used to identify levels of traffic-related air pollutants (PM2.5 and O3). Multiple logistic regression models were used to calculate the prevalence ratios (PRs) of COPD and evaluate the role played by traffic-related indices between air pollutants and COPD. The distributed lag nonlinear model was applied in the analysis; we excluded current or ever smokers to perform the sensitivity analysis.
Increased PRs of COPD per SD increment of PM2.5 were 1.10 (95% CI 1.05–1.15) and 1.25 (95% CI 1.13–1.40) in the population with age and sex matching as well as propensity-score matching, respectively. The results of the sensitivity analysis were similar between the single and two pollutant models. PM2.5 concentrations were significantly associated with traffic flow including sedans, buses, and trucks (p < 0.01). The higher road area and the higher PM2.5 concentrations near the subject’s residence correlated with a greater risk of developing COPD (p for interaction < 0.01).
Our results suggest that long-term exposure to traffic-related air pollution may be positively associated with the prevalence of COPD.
Chiayi County is located in the southwest part of Taiwan, and the Tropic of Cancer runs through it. It has the third largest agricultural land area in Taiwan. Therefore, the population density and industrial area of Chiayi County are not very dense. However, according to the environmental resource database of the Taiwan Environmental Protection Administration (TWEPA, 2020), the population of Chiayi County decreased from 563,365 in 2001 to 507,068 in 2018. Nevertheless, the number of vehicles (including cars and motorcycles) increased from 448,824 in 2001 to 523,851 in 2018. In other words, the number of vehicles owned by each person increased from 0.8 in 2001 to 1.03 in 2018; indicating that everyone in Chiayi County has one vehicle. Under such circumstances, the air pollution caused by traffic and its effects on human health is a subject worth exploring. According to data from the Taiwan Ministry of Health and Welfare, the age-standardized mortality rate of chronic lower respiratory tract in Chiayi County in 2009 was 14.5 per 100,000 persons, and by 2018, it rose to 16.3 per 100,000 persons. In contrast, Taiwan’s national data showed that in 2009, it was 14.9 per 100,000 persons, and by 2018, it had dropped to 12.8 per 100,000 persons (Taiwan Ministry of Health and Welfare, 2020). These data show that the mortality of respiratory disease in Chiayi County is gradually worsening compared with that of the whole country.
According to data from the Taiwan Air Quality Monitoring Network operated by Taiwan Environmental Protection Administration (2020), the most important pollutant in Chiayi County in the past six years (2014–2019) was PM2.5. During these periods, there were approximately 28% of daily mean PM2.5 concentrations in Chiayi County that exceeded Taiwan’s air quality standard (35 μg/m3). It is suggested that traffic is one of the most important contributors of PM2.5 in Chiayi County.
Particulate matter (PM) is a complicated mixture of solid and liquid particles consisting of consists of organic chemicals, metals, sulfate, nitrate, and ammonium. In addition, ozone (O3) is a well-known strong oxidizing agent and a secondary pollutant produced by nitrogen oxides and volatile organic compounds. Since humans are exposed to air pollutants primarily by inhalation, the respiratory system is one of the most important target organs of the harmful effects of air pollutants. In vitro studies have indicated exposure to variety of air pollutants may cause damage to lung, trachea, or bronchus, especially for fine particulate matter [1, 2]. Studies have shown that these pollutants cause lung inflammation, alveolar epithelial damage, and impaired mitochondrial function of the bronchial epithelial cells [2, 3]. They are potential candidates that induce COPD in the residents of Chiayi County.
Previous epidemiological studies on the relationship between exposure to traffic pollution sources and COPD have shown inconsistent conclusions. Andersen et al.  used a Cox proportional hazards model to study the association between COPD and exposure to traffic-related air pollution in Denmark. Their results showed that long-term exposure to traffic-related air pollution was associated with development of COPD. A cross-sectional study was conducted in Germany to investigate the effect of long-term exposure to PM10 from traffic on COPD. The results indicated that chronic exposure to PM10, NO2 and living near a major road might increase the risk of developing COPD and can have a detrimental effect on lung function . Another cross-sectional study conducted in Southern Sweden also demonstrated that living close to traffic was associated with prevalence of COPD in adults . However, a British nationwide cross-sectional study showed that close residential proximity to main roads did not increase the health risks of asthma, COPD, or allergic disease . Likewise, another city-based British study showed that there was no evidence to suggest that living in close proximity to traffic is a major determinant of asthma, allergic disease, or COPD in adults . Possible reasons for the different findings may be related to the assessment method of exposure concentration, culture differences as well as individual susceptibility to air pollution. Recently, studies have suggested that location-based air pollution data frequently have a far lower resolution than location-based health data. It is very likely to cause bias in assessing the relationship between exposure and response [9, 10].
One of the indicators for the development of COPD is decreasing lung function. A large-scale longitudinal cohort study conducted in Taiwan showed that for every 5 μg/m3 increase in PM2.5 concentration, forced vital capacity (FVC) decreased by 1.18%, and forced expiratory volume in 1 s (FEV1) decreased by 1.46%. Compared with the participants exposed to PM2.5 in the first quartile group, hazard ratio of developing COPD in the fourth quartile group was 1.23 (95% CI 1.09–1.39) . Another population-based UK Biobank study found that exposure to PM2.5 is associated with decreased lung function and increased COPD prevalence, indicating that exposure to PM2.5 is a risk factor for COPD . Furthermore, a population-based cohort study in the United States demonstrated that even at relatively low levels, long-term exposure to traffic and PM2.5 decrease FEV1 and FVC. In addition, it accelerated the rate of lung function decline . Therefore, further study the relationship between exposure to traffic pollution and COPD is suggested.
Findings from recent studies indicate that outdoor air pollution and some health outcomes have a nonlinear exposure–response relationship [14,15,16]. Therefore, in this study, we used the hybrid kriging/land-use regression (LUR) model to obtain the levels of traffic-related pollutants (PM2.5 and O3) and further explored the effects of these pollutants on COPD risk in a large-scale community-based population under different study designs and statistical analysis. In addition, we adopted the distributed lag nonlinear model (DLNM) in the analysis of non-linear relationship and search for potential index for PM2.5 as well as O3 increments. Also, we evaluated the roles that traffic-related factors played in the association between air pollutants and COPD risk, such as road area as well as the traffic load of different kinds of vehicles. Considering that smoking is a significant risk factor for COPD and increased exposure to PM2.5; we executed a sensitivity analysis and approve the above results in the non-smokers.
Materials and methods
Study area and participants
The study population was derived from a Chiayi County Complex Health Screening (CCHS) program launched from 2012 to 2013 in Chiayi County, which comprises 18 townships and has the highest percentage of elderly people (18.6%) in Taiwan. This program aimed to explore the effects of air pollution on long-term health risk for a community-based general population ≥ 40 years. Community recruitment was conducted during 2012 and 2013 and all residents aged ≥ 40 years lived in Chiayi County were invited to participate in this health program by mailing leaflet. A total of 8284 community residents voluntarily participated in the study. As shown in Fig. 1, we excluded participants younger than 40 years of age (N = 462), those who did not live in Chiayi County (N = 975), and records with missing sex data (N = 55). The final analytic cohort consisted of 6792 participants. According to the previous Epidemiology and Impact of COPD (EPIC) Asia survey, prevalence of COPD in Taiwan was approximately 10% . About 3,416 of sample size are needed for all population aged ≥ 40 years old in Chiayi (N = 280,000) to have a confidence level of 95% that the real value is within ± 1% of the surveyed value. Therefore, our sample size is large enough to represent the population in Chiayi County. Written informed consent was obtained from all participants. This study was approved by the Research Ethics Committee of China Medical University Hospital, Taichung, Taiwan (DMR101-IRB061).
Collection of questionnaires and health examinations
Well-trained personnel carried out standardized personal interviews based on a structured questionnaire, which contained demographic data and lifestyle variables contained cigarette smoking and quantity of areca nut chewing, consumption of alcohol and other beverages, participation in sport, consumption of three regular meals per day, and personal and familial history of cancer or other related diseases. In addition, all study population received health checkups, including baseline anthropometric and general biochemistry examinations such as blood pressure, plasma levels of triglycerides, total cholesterol, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), blood glucose, and blood creatinine after an 8-h fasting period. The biochemical data and self-reported data in the questionnaires were combined. Diabetes was defined as fasting glucose ≥ 126 mg/dL or use of insulin or oral hypoglycemic medications; hypertension was defined as systolic blood pressure (SBP) ≥ 140 mmHg, diastolic blood pressure (DBP) ≥ 90 mmHg, or antihypertensive medication use. Hyperlipidemia was defined as total cholesterol ≥ 200 mg/dL or triglyceride ≥ 130 mg/dL or confirmed disease status. In addition, baseline information about demographic characteristics and lifestyle variables was acquired from the questionnaires. Quantified (yes/ no) or frequency (numbers/ week) data of cigarette smoking and areca nut chewing, consumption of alcohol and other beverages, sport and three regular meal habits, personal and familial history of cancer or other related diseases were also collected. The prevalent cases of heart diseases, arthritis, asthma, chronic kidney diseases, and cancer were defined if they self-reported “Yes” to “Have you ever had heart diseases?” in the medical history portion of the questionnaire.
Outcome Assessment and sensitivity analysis
COPD diagnoses were identified with ICD-9-CM codes (491, 492, and 496) through linkage of Taiwan National Health Insurance Research Databases (TNHIRD); almost all residents (> 99%) in Taiwan have been joined this TNHIRD program since 1995. The accuracy of COPD diagnosis recorded in the NHIRD has been validated , which indicating the accuracy of diagnoses was excellent. There were 668 COPD cases in the original population (N = 6792). We then adopted 1:4 frequency-matching with age and gender to the non-COPD group; there were 654 COPD cases and 2616 non-COPD controls. In addition, we constructed a propensity-score-matched population of 451 COPD cases and 1804 controls by matching age, gender, categorized BMI, ethnicity, level of education, marital status, and COPD-related comorbidity, including hypertension, diabetes, hyperlipidemia, heart diseases, arthritis, asthma, CKD, and cancer. Given that smoking habit is an important risk factor for COPD and increased exposure to PM2.5, we excluded current or ever smokers (N = 1272) and performed a sensitivity analysis through the similar methods of 1: 4 frequency-matching and propensity-score-matching (Fig. 1).
A hybrid kriging/land-use regression (LUR) model for PM 2.5 , and O 3 estimation
We adopted the air pollutant data (PM2.5 and O3) collected from 71 Taiwan Environmental Protection Administration air quality monitoring stations between 2006 and 2013 for PM2.5 and 2000–2013 for O3 to calibrate our previously developed hybrid kriging/LUR model [10, 19]. The model takes into consideration the land use patterns, such as green areas, major traffic roads, temple incense burning, and industrial areas to improve the accuracy of predicting the exposure concentration of pollutants. Besides the above factors, other important variables such as temperature, relative humidity, wind-related factors, and meteorological as well as topography factors were also considered in the sequent stepwise model selection. Variance inflation factor < 3 was selected for the collinearity test in the model and also the statistical criterion of stepwise variable selection used in this study was 0.1. Furthermore, to improve the accuracy of PM2.5 and O3 variation predictions, the LUR model included the predicted concentration level from the kriging interpolation as a variable. Therefore, the hybrid approach further included the kriging-based concentration estimations as a predictor variable to improve the prediction performance of LUR. For data validity, the cross-validated R2 and RMSE were 0.87 and 5.02 μg/m3 for PM2.5 and 0.70 and 0.04 μg/m3 for O3, respectively. Finally, we calculated the overall average values of PM2.5 and O3 of all residents from the earliest start year to their corresponding year of recruitment (2012 or 2013).
Baseline characteristics of the COPD and non-COPD groups were compared using chi-square test for categorical variables. Multiple logistic regression models were used to evaluate the prevalence ratios (PRs) and 95% confidence interval (CI) for the associations between PM2.5, and O3 (exposure variables) with respect to the PRs of COPD (outcome variable). Statistically significant variables of COPD shown in Tables 1 and 2 were considered as covariate variables adjusted for by including them in final multiple logistic regression models. Also, two-pollutant models were used to assess the associations of PM2.5 with COPD; the two-pollutant models included O3. Because the effect of O3 on COPD risk disappeared after adjustment for PM2.5, we only explored the role of PM2.5 in the following analysis. Either b-spline or natural cubic was used for fitting the non-linear dose–response relationship of PM2.5 and COPD under different degrees of freedom. Selection of the degrees of freedoms in the final model was determined using the minimum Akaike information criterion (DLNM package in the R program). To determine the important factors for PM2.5 increments, we contained all LUR-related indices with different circular buffers (500 m, 1000 m, 1500 m, 2000 m, 2500 m, 3000 m, 4000 m, and 5000 m) in the model selection. These variables included normalized difference vegetation index (NDVI), area of industrial land, number of temples, distance from residence to municipal waste incinerator, area of traffic road. Among these variables, the overall value of R2 of the four indices including within a 5000 m circular buffer of NDVI, area of local road area, area of industrial land within a 4000 m circular buffer, and amounts of temples within a 500 m circular buffer, was approximately 0.75 for PM2.5. We further evaluated the interactions of LUR-related indices and levels of PM2.5 on PRs of COPD in the propensity-score-matched population using multiple logistic regression models. The cut-off values were defined as median of LUR-related indices and PM2.5 in the non-COPD groups. Furthermore, we analyzed the associations between the daily traffic load of different types of vehicles and the annual PM2.5 level at traffic load monitoring stations in Chiayi County during the study period by repeated measurements analysis. Different correlation structures of repeated data were all used in the generalized estimating equation approach with calculating Pan’s quasilikelihood under the independence model information criterion (QIC). The correlation structure with the lowest QIC score was generally considered to be the best . All data were analyzed using the SAS statistical package (SAS, version 9.4, Cary, NC) and R software version 3.6.3. A two-sided p-value < 0.05 was considered statistically significant.
Basic characteristics and lifestyle-related variables in COPD and non COPD population
Descriptive characteristics between study participants with and without COPD are shown in Table 1. There was approximately a 1:1 sex ratio of males to females, with average age of 66 years. About half of the entire study population had a BMI that was categorized as overweight or obesity, 62% had an elementary level education or below, and 85% were married.
For comorbidities, the study population with COPD had a higher prevalence of diabetes, heart diseases, arthritis, asthma, and cancer (all p-values < 0.05). In addition, we constructed a propensity-score analysis, matching the study population with similar distributions of COPD-related variables described above between the COPD and non-COPD groups.
The associations between lifestyles- and dietary-related variables and prevalent rate ratios of COPD are shown in Table 2. Approximately 10%-30% of the study population had habits of cigarette smoking, alcohol consumption, tea and coffee drinking, and betel consumption. Most of the subjects reported consuming less than 3 bottles of sugary drinks per week. Half of study population consumed one to three bowls of vegetables per day; however, about 60% had less than one bowl of fruit per day. Cigarette smoking and betel consumption were significantly associated with COPD irrespective of age and sex matching as well as propensity-score matching. In addition, through age and gender matching, there was a significant difference in the distribution of alcohol drinking between the COPD group and the control group. However, by way of propensity score matching, this significant difference disappeared.
Linear and non-linear relationships of PM2.5, O3, and prevalence risk of COPD
We further explored the associations between exposure to air pollution and prevalence ratios of COPD by adjusting for risk factors shown in Tables 1 and 2. The results are presented in Fig. 2. The increased prevalence ratios of COPD per SD increment of PM2.5 were 1.10 (95% CI 1.05–1.15), and 1.25-fold (95% CI 1.13–1.40) in the population with age and sex matching as well as propensity-score matching, respectively (Fig. 2A). The significant results were also found in two-pollutant models, which included the levels of O3. In addition, for O3, there was about 11–12% significantly increased prevalence ratios of COPD per SD increment of O3 in two kinds of matching population; however, after adjusting for PM2.5 in two-pollutant models, the effects of O3 on COPD were not found. The detailed data of the above prevalence ratios are shown in Additional file 1: Table S1.
In the sensitivity analysis, we excluded the ever smokers from the original cohort of 6792 residents and performed a similar 1:4 matching to repeatedly explore the association between exposure to air pollution (PM2.5 and O3) and COPD. The results indicated that there were significant positive associations between levels of PM2.5 and COPD, irrespective of whether a single or two-pollutant mode was used (Fig. 2B). For two-pollutant models, there were significantly 1.25-fold and 1.19-fold risks of COPD prevalence (both p < 0.05).
We further executed the DLNM analysis, which was set as the simple b-splines function with five degrees of freedom and with the minimum AIC value (2198.85) (Additional file 1: Table S2). Also, we selected the lowest 5% of PM2.5 levels as a reference value and the results showed a non-linear dose–response relationship of significant COPD prevalence and PM2.5 above 35 μg/m3 among propensity-score matching population. The spline results are displayed graphically in Fig. 3.
The role of traffic-related variables for PM2.5 and prevalence ratios of COPD
Since the novelty effect of PM2.5 was observed in the above results compared to O3, we only explored the role of PM2.5 in the following analysis. For the four important contributors of PM2.5 increment from model selection (NDVI, area of industrial land, number of temples, and area of traffic road) we evaluated the interactions of individual contributor and PM2.5 levels on prevalence ratios of COPD in the propensity score matched population in Table 3. The results suggest statistically significant dose–response relationships between increasing levels of PM2.5 combined with increased amount of individual contributor for positive associations of COPD (all p < 0.01). Among the four aforementioned contributors, high PM2.5 levels combined with high road area or with high area of industrial land significantly interact on increased prevalence ratios of COPD (both p < 0.01) after adjusting for smoking and betel consumption. Furthermore, we calculated the nearest distance from residential address to the surrounding bus station and analyzed the associations between the levels of PM2.5 and distance to bus station through generalized linear regression model. For residents that lived within 1 km of the bus station, the results suggested residents who live farther away from the bus station have a low PM2.5 concentration after adjustment for townships and smoking habit (p = 0.004, data not shown).
For the role of area of traffic road on PM2.5 increment, we further analyzed the associations between the daily traffic load of different types of vehicles and annual PM2.5 at traffic flow monitoring stations in Chiayi County during 2009–2014 using repeated measurements analysis. The results showed an unstructured correlation structure had the lowest QIC score. It indicated a statistically significant positive association between traffic load and increased levels of PM2.5, especially for passenger cars, buses, and trucks (p < 0.01, p < 0.01, and p < 0.01, respectively) (Table 4). It demonstrated that the PM2.5 concentrations in Chiayi County were significantly related to emissions from sedans, buses, and trucks.
This present study adopted the hybrid kriging/LUR model included the kriging-based concentration estimations as a predictor variable to improve the prediction performance of LUR and provided the individual level of exposure to ambient air pollutants. The results indicated a positive correlation between exposure to PM2.5 and PRs of COPD in a community-based population even excluding all smokers from the original cohort. Besides, traffic-related variables including road area, traffic load of sedans, buses, and trucks, as well as living near a bus station were relative to PM2.5 increment.
Several epidemiological studies have shown that exposure to PM2.5 is related to COPD. It is estimated that long-term exposure to ambient PM2.5 contributes between 10.7% and 15.3% to COPD in Iran . A meta-analysis study demonstrated that exposure to PM2.5 is significantly associated with prevalence of COPD (OR: 2.32, 95% CI 1.91–2.82) . A Korean study showed that concentration of PM2.5 was associated with increased COPD-related hospital visits in Chuncheon . The results of this investigation are consistent with those of the above studies. In addition, we used the DLNM model to analyze the dose–response relationship between exposure to PM2.5 and prevalence of COPD. Our results demonstrated that there was a non-linear relationship between PM2.5 and prevalence of COPD. As the concentration of PM2.5 rose above 35 μg/m3, the prevalence of COPD increased significantly. This finding is very similar to another Taiwanese population-based study, which found that exposure to PM2.5 at concentrations greater than 38.98 μg/m3 increased susceptibility to COPD among nonsmokers .
The present study applied the hybrid kriging/LUR model to predict the exposure concentration of air pollutants at residential addresses. This model has been used extensively to assess the effects of exposure to PM2.5 on human health [25,26,27]. In the model, we consider the NDVI, the areas of the industrial land, the areas of the traffic road, and the number of temples, etc. to predict the exposure concentration of PM2.5 and ozone. We found that increased levels of PM2.5 were significantly associated with NDVI, areas of industrial land, traffic roads, and number of temples. Among them, the areas of industrial parks and traffic roads near the subject’s residency and the exposure concentration of PM2.5 have a multiplicative effect on the PRs of COPD. The results demonstrated that compared with the low traffic road area and low PM2.5 concentration, the population with high traffic road area and a high PM2.5 concentration near their homes had a 1.66 times greater risk (95% CI 1.21–2.28) of developing COPD. In addition, we verified this association by calculating the nearest distance from the residential address of subjects to surrounding bus stations. A statistically significant negative association between the distance and levels of PM2.5 was observed (data not shown). The number of vehicles in Chiayi County increasing yearly is considered to be correlated to gradually increasing the risk of respiratory diseases in Chiayi County; this is even higher than the situation in Taiwan (Taiwan Ministry of Health and Welfare, 2020). Our data show that traffic is one of the main factors contributing to the increasing development of COPD in Chiayi County.
We used the data from the traffic flow monitoring stations (N = 74) in Chiayi County to analyze the relationship between the daily traffic flow of various vehicles and the concentration of PM2.5 to explore the impact of different types of vehicles on the concentrations of PM2.5. Our results show that buses had the greatest impact on the variability of PM2.5 concentration, followed by sedans, and then trucks. The relationship between the traffic load of these three types of vehicles and PM2.5 showed a statistically significant correlation. Studies have suggested that diesel-fueled vehicles have a higher PM2.5 emission factor than that of gasoline-fueled vehicles. For example, a study in Taiwan showed that a sedan's emission factor is 1.25 mg/km, while those of diesel engine trucks are as high as 185 mg/km . Another Chinese study also showed that the emission factor for diesel-fueled vehicles is about 257 mg/km and for non-diesel-fueled vehicles is about 17 mg/km . Both buses and trucks are diesel-fueled vehicles. Therefore, the findings of our study showed that the traffic flow of the buses and trucks had a positive correlation with the concentration of PM2.5. Since buses need to allow passengers to get on and off within a certain distance, the operations of their engines vary greatly in such conditions, which may be attributed to buses having the greatest impact on PM2.5 emissions.
Previous large-scale community-based studies suggest there are high frequencies of additional comorbidities among elderly people with COPD [30, 31], which is consistent with our findings. A survey from the US Center for Disease Control and Prevention assistance found that the most common co-morbid chronic conditions among people with COPD included arthritis, asthma, cancer, and coronary artery disease . Oxidative stress and systemic inflammation were the common mechanisms in the development and progression of COPD as well as other comorbidities [33, 34]. COPD is a state of systemic inflammation  and high levels of inflammatory markers are associated with the severity of airflow obstruction and cardiac injury . Hypoxia induced expression of proinflammatory transcription factors lead to endothelial dysfunction and then to atherosclerosis . Lung infections were very important precipitating factors for acute exacerbation of COPD, and some of them accelerated atherosclerosis and precipitated acute coronary syndrome by causing plaque instability . COPD also had a detrimental effect on CAD through hypoxia, decreased respiratory muscle strength, and use of bronchodilators . Inflammatory mediators, such as IL-17 and anti–citrullinated protein antibodies, which play a role in arthritis, are also involved in the pathogenesis of COPD . COPD, arthritis, and coronary artery disease share many of the same risk factors, such as sex, age, tobacco use, obesity, and sedentary lifestyle . The abovementioned mechanisms must be considered to explain the association between COPD, CAD, and arthritis.
Our study has several limitations. We could not exclude the temporality of the association between air pollution and COPD. This may overlook pre-existing lung injury attributed to occupational exposure or lifestyle factors such as smoking and increased sensitivity to exposure to air pollution. In addition, empirical duration in outdoor or indoor was not acquired in the analysis. However, smoking is an important source of indoor exposure to air pollution. Therefore, we constructed a sensitivity analysis to excluded the possible effect of smoking and still found a positive association between levels of PM2.5 and prevalence ratios of COPD. Another limitation is that the study population was volunteers participating in a community-based health screening; therefore, most participants tended to be healthy and this could underestimate the prevalence ratios of COPD. In addition, study participants without sever COPD did not seek medical help and could not be acquired in our study, which may bias the present results. Moreover, about 70% of study population was elder with average aged 65 years old (in Table 1) and they had no work. The data of outdoor activities duration were not well evaluated for elder population. Therefore, these variables were not further considered in our study. We did not collect chemical composition data of PM2.5 such as polycyclic aromatic hydrocarbons or heavy metals in this study. However, it is suggested that PAHs or heavy metals in PM2.5 may have a positive association with development of COPD [42, 43]. We recommended that follow-up studies could further analyze the impact of PM2.5 components on COPD to understand the mechanisms of PM2.5 on COPD.
The positive association between levels of PM2.5 and PRs of COPD was observed in a community-based population even excluding all smokers from the original cohort. The results showed a non-linear relationship of PM2.5 above 35 μg/m3 and prevalence of COPD in the DLNM analysis. Areas of industrial land as well as roads individually interacted with levels of PM2.5 on increased prevalence ratios of COPD. In the analysis of traffic flow, we found that sedans, buses, and trucks all had a significant positive association with the variations in PM2.5 concentration. It is suggested that traffic emitted PM2.5 is an important factor in the development of COPD in Chiayi County. In the future, empirical measurements of indoor and outdoor air pollution should be further expanded and explored with regard to the role of air pollution on COPD incidence.
Chronic obstructive pulmonary disease
- PM2.5 :
- NO2 :
- O3 :
Low-density lipoprotein cholesterol
High-density lipoprotein cholesterol
Systolic blood pressure
Diastolic blood pressure
Taiwan National Health Insurance Research Databases
Hybrid kriging/land-use regression
Normalized difference vegetation index
Wang G, Zhao J, Jiang R, Song W. Rat lung response to ozone and fine particulate matter (PM2.5) exposures. Environ Toxicol. 2015;30(3):343–56.
Leclercq B, Kluza J, Antherieu S, Sotty J, Alleman LY, Perdrix E, Loyens A, Coddeville P, Lo Guidice JM, Marchetti P, et al. Air pollution-derived PM2.5 impairs mitochondrial function in healthy and chronic obstructive pulmonary diseased human bronchial epithelial cells. Environ Pollut. 2018;243(Pt B):1434–49.
Hansel NN, McCormack MC, Kim V. The effects of air pollution and temperature on COPD. COPD. 2016;13(3):372–9.
Andersen ZJ, Hvidberg M, Jensen SS, Ketzel M, Loft S, Sorensen M, Tjonneland A, Overvad K, Raaschou-Nielsen O. Chronic obstructive pulmonary disease and long-term exposure to traffic-related air pollution: a cohort study. Am J Respir Crit Care Med. 2011;183(4):455–61.
Schikowski T, Sugiri D, Ranft U, Gehring U, Heinrich J, Wichmann HE, Kramer U. Long-term air pollution exposure and living close to busy roads are associated with COPD in women. Respir Res. 2005;6:152.
Lindgren A, Stroh E, Montnemery P, Nihlen U, Jakobsson K, Axmon A. Traffic-related air pollution associated with prevalence of asthma and COPD/chronic bronchitis. A cross-sectional study in Southern Sweden. Int J Health Geogr. 2009;8:2.
Pujades-Rodriguez M, Lewis S, McKeever T, Britton J, Venn A. Effect of living close to a main road on asthma, allergy, lung function and chronic obstructive pulmonary disease. Occup Environ Med. 2009;66(10):679–84.
Pujades-Rodriguez M, McKeever T, Lewis S, Whyatt D, Britton J, Venn A. Effect of traffic pollution on respiratory and allergic disease in adults: cross-sectional and longitudinal analyses. BMC Pulm Med. 2009;9:42.
Alexeeff SE, Schwartz J, Kloog I, Chudnovsky A, Koutrakis P, Coull BA. Consequences of kriging and land use regression for PM2.5 predictions in epidemiologic analyses: insights into spatial variability using high-resolution satellite data. J Expo Sci Environ Epidemiol. 2015;25(2):138–44.
Wu CD, Zeng YT, Lung SC. A hybrid kriging/land-use regression model to assess PM2.5 spatial-temporal variability. Sci Total Environ. 2018;645:1456–64.
Doiron D, de Hoogh K, Probst-Hensch N, Fortier I, Cai Y, De Matteis S, Hansell AL. Air pollution, lung function and COPD: results from the population-based UK Biobank study. Eur Respir J. 2019;54(1):1802140.
Guo C, Zhang Z, Lau AKH, Lin CQ, Chuang YC, Chan J, Jiang WK, Tam T, Yeoh EK, Chan TC, et al. Effect of long-term exposure to fine particulate matter on lung function decline and risk of chronic obstructive pulmonary disease in Taiwan: a longitudinal, cohort study. Lancet Planet Health. 2018;2(3):e114–25.
Rice MB, Ljungman PL, Wilker EH, Dorans KS, Gold DR, Schwartz J, Koutrakis P, Washko GR, O’Connor GT, Mittleman MA. Long-term exposure to traffic emissions and fine particulate matter and lung function decline in the Framingham heart study. Am J Respir Crit Care Med. 2015;191(6):656–64.
Wang C, Feng L, Chen K. The impact of ambient particulate matter on hospital outpatient visits for respiratory and circulatory system disease in an urban Chinese population. Sci Total Environ. 2019;666:672–9.
Jung CR, Chung WT, Chen WT, Lee RY, Hwang BF. Long-term exposure to traffic-related air pollution and systemic lupus erythematosus in Taiwan: a cohort study. Sci Total Environ. 2019;668:342–9.
Li Z, Mao X, Liu Q, Song H, Ji Y, Xu D, Qiu B, Tian D, Wang J. Long-term effect of exposure to ambient air pollution on the risk of active tuberculosis. Int J Infect Dis. 2019;87:177–84.
Lim S, Lam DC, Muttalif AR, Yunus F, Wongtim S, le Lan TT, Shetty V, Chu R, Zheng J, Perng DW, et al. Impact of chronic obstructive pulmonary disease (COPD) in the Asia-Pacific region: the EPIC Asia population-based survey. Asia Pac Fam Med. 2015;14(1):4.
Su VY, Yang KY, Yang YH, Tsai YH, Perng DW, Su WJ, Chou KT, Su KC, Yen YF, Chen PC. Use of ICS/LABA combinations or LAMA Is associated with a lower risk of acute exacerbation in patients with coexistent COPD and asthma. J Allergy Clin Immunol Pract. 2018;6(6):1927-35 e3.
Hsu CY, Wu JY, Chen YC, Chen NT, Chen MJ, Pan WC, Lung SC, Guo YL, Wu CD. Asian culturally specific predictors in a large-scale land use regression model to predict spatial-temporal variability of ozone concentration. Int J Environ Res Public Health. 2019;16(7):1300.
Ballinger GA. Using generalized estimating equations for longitudinal data analysis. Organ Res Methods. 2004;7(2):127–50.
Faridi S, Shamsipour M, Krzyzanowski M, Kunzli N, Amini H, Azimi F, Malkawi M, Momeniha F, Gholampour A, Hassanvand MS, et al. Long-term trends and health impact of PM2.5 and O3 in Tehran, Iran, 2006–2015. Environ Int. 2018;114:37–49.
Han F, Yang X, Xu D, Wang Q, Xu D. Association between outdoor PM2.5 and prevalence of COPD: a systematic review and meta-analysis. Postgrad Med J. 2019;95(1129):612–8.
Jo YS, Lim MN, Han YJ, Kim WJ. Epidemiological study of PM2.5 and risk of COPD-related hospital visits in association with particle constituents in Chuncheon, Korea. Int J Chron Obstruct Pulmon Dis. 2018;13:299–307.
Huang HC, Lin FC, Wu MF, Nfor ON, Hsu SY, Lung CC, Ho CC, Chen CY, Liaw YP. Association between chronic obstructive pulmonary disease and PM2.5 in Taiwanese nonsmokers. Int J Hyg Environ Health. 2019;222(5):884–8.
Li CY, Wu CD, Pan WC, Chen YC, Su HJ. Association between long-term exposure to PM25 and incidence of type 2 diabetes in taiwan: a national retrospective cohort study. Epidemiology. 2019;30(Suppl 1):S67–75.
Wu CD, Chern YR, Pan WC, Lung SC, Yao TC, Tsai HJ, Spengler JD. Effects of surrounding environment on incidence of end stage renal disease. Sci Total Environ. 2020;723:137915.
Chen XY, Feng PH, Han CL, Jheng YT, Wu CD, Chou HC, Chen YY, Wu SM, Lee KY, Kuo HP, et al. Alveolar epithelial inter-alpha-trypsin inhibitor heavy chain 4 deficiency associated with senescence-regulated apoptosis by air pollution. Environ Pollut. 2021;278:116863.
Yang HH, Dhital NB, Wang LC, Hsieh YS, Lee KT, Hsu YT, Huang SC. Chemical characterization of fine particulate matter in gasoline and diesel vehicle exhaust. Aerosol Air Quality Res. 2019;19(6):1439–49.
Cheng Y, Lee SC, Ho KF, Chow JC, Watson JG, Louie PK, Cao JJ, Hai X. Chemically-speciated on-road PM(2.5) motor vehicle emission factors in Hong Kong. Sci Total Environ. 2010;408(7):1621–7.
Cazzola M, Bettoncelli G, Sessa E, Cricelli C, Biscione G. Prevalence of comorbidities in patients with chronic obstructive pulmonary disease. Respiration. 2010;80(2):112–9.
Schnell K, Weiss CO, Lee T, Krishnan JA, Leff B, Wolff JL, Boyd C. The prevalence of clinically-relevant comorbid conditions in patients with physician-diagnosed COPD: a cross-sectional study using data from NHANES 1999–2008. BMC Pulm Med. 2012;12:26.
Cunningham TJ, Ford ES, Rolle IV, Wheaton AG, Croft JB. Associations of self-reported cigarette smoking with chronic obstructive pulmonary disease and co-morbid chronic conditions in the United States. COPD. 2015;12(3):276–86.
Liu Y, Wheaton AG, Murphy LB, Xu F, Croft JB, Greenlund KJ. Chronic obstructive pulmonary disease and arthritis among US adults, 2016. Prev Chronic Dis. 2019. https://doi.org/10.5888/pcd16.190035.
Pastor MD, Nogal A, Molina-Pinelo S, Melendez R, Romero-Romero B, Mediano MD, Lopez-Campos JL, Garcia-Carbonero R, Sanchez-Gastaldo A, Carnero A, et al. Identification of oxidative stress related proteins as biomarkers for lung cancer and chronic obstructive pulmonary disease in bronchoalveolar lavage. Int J Mol Sci. 2013;14(2):3440–55.
Sin DD, Man SF. Why are patients with chronic obstructive pulmonary disease at increased risk of cardiovascular diseases? The potential role of systemic inflammation in chronic obstructive pulmonary disease. Circulation. 2003;107(11):1514–9.
Perotin JM, Adam D, Vella-Boucaud J, Delepine G, Sandu S, Jonvel AC, Prevost A, Berthiot G, Pison C, Lebargy F, et al. Delay of airway epithelial wound repair in COPD is associated with airflow obstruction severity. Respir Res. 2014;15:151.
Ryan S, Taylor CT, McNicholas WT. Systemic inflammation: a key factor in the pathogenesis of cardiovascular complications in obstructive sleep apnoea syndrome? Thorax. 2009;64(7):631–6.
Roversi S, Fabbri LM, Sin DD, Hawkins NM, Agusti A. Chronic obstructive pulmonary disease and cardiac diseases. An urgent need for integrated care. Am J Respir Crit Care Med. 2016;194(11):1319–36.
Singer J, Yelin EH, Katz PP, Sanchez G, Iribarren C, Eisner MD, Blanc PD. Respiratory and skeletal muscle strength in chronic obstructive pulmonary disease: impact on exercise capacity and lower extremity function. J Cardiopulm Rehabil Prev. 2011;31(2):111–9.
McGuire K, Avina-Zubieta JA, Esdaile JM, Sadatsafavi M, Sayre EC, Abrahamowicz M, Lacaille D. Risk of incident chronic obstructive pulmonary disease in rheumatoid arthritis: a population-based cohort study. Arthritis Care Res (Hoboken). 2019;71(5):602–10.
Liu Y, Pleasants RA, Croft JB, Lugogo N, Ohar J, Heidari K, Strange C, Wheaton AG, Mannino DM, Kraft M. Body mass index, respiratory conditions, asthma, and chronic obstructive pulmonary disease. Respir Med. 2015;109(7):851–9.
Huang BF, Chang YC, Han AL, Hsu HT. Metal composition of ambient PM2.5 influences the pulmonary function of schoolchildren: a case study of school located nearby of an electric arc furnace factory. Toxicol Ind Health. 2018;34(4):253–61.
Lag M, Ovrevik J, Refsnes M, Holme JA. Potential role of polycyclic aromatic hydrocarbons in air pollution-induced non-malignant respiratory diseases. Respir Res. 2020;21(1):299.
This study was supported by grants from the Ministry of Science and Technology of Taiwan under Grant No. MOST 105-2314-B-039-022 and Ministry of Education of Taiwan under Grant No. 11046H02. The funders did not participate in the study design, data collection and analysis, publication decision, or manuscript preparation.
Ethics approval and consent to participate
All participants were given detailed oral and written information on the study and gave written informed consent. This study was approved by the Research Ethics Committee of China Medical University Hospital, Taichung, Taiwan.
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1
: Table S1. Associations between indices of air pollutants and COPD risk: single- and two-pollutant models. Table S2. AIC values under different spline functions and degree of freedom.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Hsu, HT., Wu, CD., Chung, MC. et al. The effects of traffic-related air pollutants on chronic obstructive pulmonary disease in the community-based general population. Respir Res 22, 217 (2021). https://doi.org/10.1186/s12931-021-01812-x
- Chronic obstructive pulmonary disease
- Land-use regression model
- Air pollution