Next Article in Journal
Assessing the Impact of Urbanization on Direct Runoff Using Improved Composite CN Method in a Large Urban Area
Previous Article in Journal
The Impact of Coworkers’ Safety Violations on an Individual Worker: A Social Contagion Effect within the Construction Crew
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multicity Analysis of the Short-Term Effects of Air Pollution on the Chronic Obstructive Pulmonary Disease Hospital Admissions in Shandong, China

1
Department of Epidemiology and Biostatistics, School of Public Health, Shandong University, 44, Wenhuaxi Street, Jinan 250012, China
2
Health and Family Planning Information Center of Shandong Province, 75, Yuhan Street, Jinan 250002, China
3
Shandong Academy of Environmental Science, 50, Lishan Street, Jinan 250013, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(4), 774; https://doi.org/10.3390/ijerph15040774
Submission received: 10 February 2018 / Revised: 8 April 2018 / Accepted: 12 April 2018 / Published: 17 April 2018
(This article belongs to the Section Environmental Health)

Abstract

:
Although there is growing evidence linking chronic obstructive pulmonary disease (COPD) hospital admissions to the exposure to ambient air pollution, the effect can vary depending on the local geography, pollution type, and pollution level. The number of large-scale multicity studies remains limited in China. This study aims to assess the short-term effects of ambient air pollution (PM2.5, PM10, SO2, NO2) on chronic obstructive pulmonary disease hospital admissions from 2015 to 2016, with a total of 216,159 records collected from 207 hospitals in 17 cities all over the Shandong province, east China. Generalized additive models and penalized splines were applied to study the data whilst controlling for confounding meteorological factors and long-term trends. The air pollution was analyzed with 0–6 day lag effects and the percentage change of hospital admissions was assessed for a 10-μg/m3 increase in the air pollution levels. We also examined the percentage changes for different age groups and gender, respectively. The results showed that air pollution was significantly associated with adverse health outcomes and stronger effects were observed for females. The air pollution health effects were also impacted by geographical factors such that the air pollution had weaker health effects in coastal cities.

1. Introduction

Many studies across the world have shown that poor air quality adversely impacts people’s health [1,2,3,4,5]. The Global Burden of Disease assessments indicated that ca. 7 million deaths per year and 3.1% of the global disease burden could be attributed to ambient particulate matter pollution, placing it among the top health risk factors globally [6,7]. Exposure to high levels of air pollutant concentration is particularly associated with respiratory and cardiovascular problems [8,9,10,11,12]. In respiratory diseases, Chronic Obstructive Pulmonary Disease (COPD) is presently the fourth leading cause of death and it is predicted to become the third leading cause by 2030, making this disease one of the major health challenges in the future [13,14]. COPD led to ca. 3 million deaths in 2004, 90% of which were from developing countries [15].
Due to the fast development in industrialization and urbanization with a marked incline in oil and coal burning, poor air quality has become a severe problem in China [16,17,18]. Heavy smog caused by a combination of industrial pollution sources and adverse weather conditions has been appearing more and more often in northern China. Considering the very large number of people who may be exposed to substantial amounts of air pollution, this is an important and urgent health issue [19,20,21,22].
Compared with extreme conditions such as mortality or morbidity, hospital admission is more directly related to air pollution exposure. Most of the existing studies were conducted in developed countries such as Europe and the United States [23,24,25,26,27,28,29,30,31]. The number of studies concentrating on the relationship between air pollution levels in multicities and the corresponding hospital admissions in China is limited. There are a few studies focused only on a single city [32,33,34,35,36].
COPD, as a type of respiratory disease, can be exacerbated by air pollution. In this study, we aimed to assess the short-term association between air pollution and chronic obstructive pulmonary disease (COPD) hospital admissions collected from all the tertiary-level and secondary-level hospitals in the 17 cities of the Shandong province. We then explored the influence of different geography, air pollution types, and pollution levels on human health; hence, providing public health evidence for preventative measures and health management for different ambient pollution environments. Shandong is located in east China; the eastern part of the province is by the sea and the western and middle parts are hilly areas. In addition, Shandong is one of the most prosperous and heavily polluted provinces in the country, with significant coal mining, oil refineries, and metallurgical and mechanical sectors.

2. Materials and Methods

2.1. Hospital Admissions

Hospital admission data were obtained from the Health and Family Planning Information Center of the Shandong Province. Daily counts of COPD hospital admissions between 1 January 2015 and 31 December 2016 were constructed from individual records of all tertiary-level and secondary-level hospitals in Shandong; there is at least one tertiary-level hospital in each city. Each record contains the hospital name, date of service, age, gender, date of birth, current address, and disease ICD code (International Classification of Diseases, Tenth Revision (ICD-10) codes). Chronic obstructive pulmonary disease outcomes were considered based on the ICD-10 codes J40–44.

2.2. Air Pollution and Meteorological Data

There are 93 air quality monitoring stations established all over Shandong by the China National Environmental Protection Administration. The number of air quality monitoring sites in each city are shown in Figure 1. Each air quality monitoring station provides hourly readings of PM10, PM2.5, SO2, and NO2 levels. Since the hospital admissions data were available daily, we aggregated the pollution data to daily counts by taking the mean value of all 24 hourly collected air pollution concentrations on each day. The meteorological data were obtained from the China Meteorological Administration Data Center (URL https://www.data.cma.cn/). These data were collected daily and included the location and altitude of the station, the daily mean, and the highest and lowest value of the monitored objects. In our study, the daily mean temperature, the relative humidity, and the wind speed were used to eliminate the confounding effects of weather.

2.3. Method of Analysis

We applied generalized additive Quasi-Poisson models to examine the associations between air pollution (PM2.5, PM10, SO2, NO2) and day-to-day variations in the COPD hospital admissions for the total patients, for males, for females, for patients aged over 65, and for patients under 65. The Quasi-Poisson distribution was adopted in order to overcome the overdispersion of hospital admissions data. Generalized additive models allow for highly flexible fitting as the outcome is assumed to be dependent on a sum of the smoothed and linear functions of the predictor variables. A smoothed function captures the nonlinear relationship between daily hospital admissions and the time-varying covariates such as the temperature and the calendar time [37,38]. The nonlinear term in generalized additive models can be estimated by using smoothing splines, which transform a possible nonlinear relationship into a linear form by separating the data into subintervals based on the basis functions. The smoothness of the splines depends on the number of basis functions: too many subintervals may lead to oversmoothness, whilst an inadequate number of subintervals leads to a rough model fitting. A penalized spline uses an overly large number of basis functions and penalizes excess curvature by using a penalty term. The trade-off between model fitting and model smoothness is controlled by the smoothing parameter which may be chosen by using a data-driven criterion such as generalized cross-validation (GCV) [39,40].
For each of the 17 cities in Shandong, the daily ambient air pollution (PM10, PM2.5, SO2, and NO2) concentrations for the 0–6 day lags were analyzed separately in association with daily hospitalizations. Briefly, we fitted the data into the following model:
L o g [ E ( Y i ) ] = α + β X i + p s ( t i m e ) + j p s ( Z i j ) + D O W
where E ( Y i ) is the expected number of daily COPD hospital admissions at day i , α is the intercept term, β is the regression coefficient, and X i is the pollutant measurements on day i . p s ( t i m e ) denotes the penalized splines of calendar time and p s ( Z i j ) denotes the penalized splines of the meteorological variables Z j such as the temperature, relative humidity, and wind speed, respectively. The dummy variable DOW is the indicator for the day of the week. All statistical analyses were performed with R V.3.4.3 (URL http://www.R-project.org) using the package mgcv (V.1.8–17). Then all results were presented as the percentage change in the relative risk (RR) of hospitalizations and its 95% confidence interval (CI) in association with a 10 μg/m3 increase in daily air pollutants.

3. Results

The locations of the 17 cities in Shandong Province are shown in Figure 1. There are 5 coastal cities: Dongying, Yantai, Weihai, Qingdao, and Rizhao. In this study, the total number of COPD hospital admissions during the study period was 216,159, the numbers of male and female records were 131,499 and 84,660, respectively. There were 162,601 records for patients over 65 years old, and 53,558 records for those under 65. The largest number of total hospital admissions during the two-year study time was 26,653 in the city of Weifang, while the smallest number was 4026 from the city of Dongying. The details of the hospital admissions in the 17 cities in Shandong can be found in Table 1.
Table A1, Table A2, Table A3, Table A4, Table A5, Table A6 and Table A7 in Appendix A present the minimum, maximum, 25% percentile, 50% percentile, 75% percentile, and mean values of the air pollution levels and meteorological factors from 1 January 2015 to 31 December 2016 in the 17 cities. Figure 2 shows the daily maximum air pollution levels in the 17 cities of the Shandong Province. Most of the daily maximum values of PM2.5, PM10, and NO2 in the 17 cities were 2–5 times that of China’s ambient air quality secondary standards (PM2.5: 75 μg/m3, PM10: 150 μg/m3, SO2: 150 μg/m3, and NO2: 80 μg/m3) [41], but the daily maximum concentrations of SO2 in 8 cities (Heze, Liaocheng, Linyi, Qingdao, Rizhao, Weihai, Ynatai, and Zaozhuang) were within the standard. There was an obvious spatial heterogeneity among the distribution of air pollution such that the air pollution concentrations in the eastern coastal cities like Qingdao, Yantai, Weihai, and Rizhao were generally lower than the measurements in the noncoastal cities.
Table 2 outlines the most statistically significant (p < 0.05) effect estimates and corresponding 95% confidence intervals (CIs) of the percentage changes of relative risk (RR) for the total COPD admissions per 10-μg/m3 increase in PM2.5, PM10, SO2, and NO2 concentrations for 0–6 day lags in Shandong. There were significant associations between PM2.5 and the total COPD admissions in 8 cities: Dezhou, Jinan, Jining, Laiwu, Linyi, Taian, Weifang, and Zaozhuang. The most severe effect was 0.978 (0.105, 1.858) by cause per 10-μg/m3 increase in PM2.5 in the city of Laiwu. The increase in PM10 concentrations led to a trend very similar to PM2.5, having a significant effect on the total admissions in 9 noncoastal cities, the largest one being 1.099 (0.458, 1.743), also in Laiwu. Furthermore, a 10-μg/m3 increase of SO2 was significantly associated with the increase of the total COPD hospital admissions in 12 cities, including two coastal cities: Dongying and Qingdao. Specifically, Qingdao had the most serious health effect: 3.164 (1.302, 5.059). Finally, the increase in NO2 levels had significant effects on hospitalizations in only 6 cities: Dezhou, Dongying, Liaocheng, Qingdao, Zaozhuang, and Zibo. The largest significant effect caused by NO2 in terms of the percentage change in the relative risk of hospitalizations for total patients was 3.515 (1.394, 5.590). This was in the city of Dezhou. These results for the total number of patients are also shown in the statistical map of Figure 3. The region colored in gray represents an insignificant effect in the area, while regions with a yellow to dark brown color demonstrate a significant percentage change in the relative risk of total COPD hospitalizations per 10-μg/m3 increase in that area. Figure 4 gives a forest plot of these results.
The health effects of air pollution were small in some cities with high pollution levels. From Table A1 and Table A2, we observed the highest levels of fine particulate matter (PM2.5 and PM10) occurred in the city of Liaocheng, which is located in the west of the province. However, the effects of PM2.5 and PM10 on the COPD hospital admissions in Liaocheng were 0.388 (−0.053, 0.830) and 0.262 (−0.052, 0.576), respectively; neither was significant. Similarly, during the study period, the most SO2- and NO2-polluted city was Zibo, located in the middle of the province. Although the effects of SO2 and NO2 were significant in Zibo (1.210 (0.375, 2.052) and 1.602 (0.130, 3.053), respectively), they were the weakest significant values among all the significant effects in the province. In contrast with this, SO2 and NO2 concentrations in Qingdao were relatively low, but the health effects from them were very strong: 3.164 (1.302, 5.059) and 2.824 (1.004, 4.677), respectively. In particular, the health effect caused by SO2 was the strongest among all the cities in Shandong.
Table A8, Table A9, Table A10 and Table A11 in Appendix A presents the results of the subgroup analysis. Regarding gender, stronger effects were observed for females. For example, the SO2 effect on male patients in Zaozhuang was insignificant: 1.279 (−0.346, 2.931), while the effect on female patients was significant: 2.570 (0.517, 4.666). Regarding age, although the effects on the aged patients were generally stronger, the effects on the younger patients were significant in some cities while the ones on older patients were not. For example, the PM2.5 effect on aged patients in Weifang was insignificant: 0.463 (−0.105, 1.035), but the effect on patients under 65 years of age was significant: 0.915 (0.226, 1.608).

4. Discussion

In this study, we examined the relationship between four air pollution (PM2.5, PM10, SO2, NO2) levels and the daily changes in COPD hospital admissions in 17 cities all over the Shandong province in China. We observed significant associations between ambient air pollution and COPD hospital admissions in Shandong. Epidemiological studies require reliable datasets. Therefore, for the health data, we adopted records from tertiary-level and secondary-level hospitals. All the patients went through a rigorous diagnosis and treatment during the admission time; hence, the accuracy of diagnosis can be guaranteed. Moreover, air pollution and meteorological data were collected from national monitoring stations. These datasets have the strength of analyzing the regional data using the national standard, avoiding the potential for publication bias in which positive findings are selectively reported.
From the results, we can see that the associations between air pollution and COPD hospitalizations were spatially heterogeneous across the province and that air pollution had obvious adverse health effects in the middle and west of the province. The cities in this area contained a large area of heavily polluting industrial sectors such as coal mining, the chemical industry, and the heavy industry. For example, Laiwu has the largest steel factory in the province, which made its PM10 and PM2.5 effects the strongest in the province. Geographical factors are also decisive. A large part of the middle and west of the Shandong province is hilly. The emitted ambient pollution is hard to blow away from these areas, consequently leading to a substantial exposure of the surrounding population to air pollution. Although the air pollution effects were much weaker in the coastal cities, SO2 and NO2 concentrations had significant influences on the hospitalizations in two coastal cities: Qingdao and Dongying. The most prominent source of NO2 and SO2 are fossil fuels burning from the industrial process with significant amounts of secondary aerosol, and the next most prominent sources are vehicle and domestic heating system exhaust emissions [42]. Qingdao is one of the most prosperous cities in China. It has the largest number of vehicles in Shandong: 2.3 million vehicles by the end of 2016, while there were only 1.8 million in the capital city of Jinan. Next, Dongying has the third largest oil field in China; the burning from the oil refinery process could be a major source of the SO2 and NO2 emissions which strongly affected the COPD hospital admissions.
We also found the effects of air pollution were small in some cities with high pollution levels. This could be caused by the active avoidance to exposure from the people living in heavily polluted cities, such as reducing outdoor activities or wearing masks. In addition, the body functions, especially the respiratory process, may have adjusted themselves to the long-term exposure to high-level pollution, hence, becoming less sensitive to toxic components. These results suggest that the studies focusing on a particular city may be unable to provide enough information for estimating the true effects, possibly even leading to biased significant results.
In the subgroup analysis, the relative risks of exposure to all of the pollutants for females were stronger for COPD. This is consistent with the results from previous studies that suggested that females were more vulnerable to ambient air pollution [33,43]. Furthermore, although the significant effects on the patients over 65 years of age were widespread in more cities, they were not overwhelmingly greater than the effects on younger patients. A possible explanation may be that the younger people would have more outdoor activities and are more likely to be exposed to air pollution if they were in the construction, manufacturing, or coal mining industries.
There are limitations in our analysis. Information on ambient concentrations in each city often comes from a number of monitoring sites, each of which may be subject to measurement errors and may contain periods of missing data. There is supposed to be a true underlying pollution surface which will form the basis of the exposures experienced by the population at risk. However, this surface is not directly observable and instead, measurements are taken at locations over space and time. The main limitation of this study is the lack of precise exposure estimates. Our exposure modeling approach is to simply take the average over all the air pollution data collected at the monitoring stations in the same area. This method may eliminate the spatial and temporal features of the ambient air pollution data; hence, the corresponding percentage change is not adequate. We need to apply a more sophisticated modeling method to the air pollution data in future studies.

5. Conclusions

This study assessed the associations between four types of air pollution and COPD hospital admissions using the daily time series data from 17 cities in the Shandong province, China. Our analysis showed that air pollution was significantly associated with adverse health outcomes. We found that the air pollution health effects on COPD were substantially impacted by the geographical factors and pollution type, such that air pollution had weaker impacts on COPD hospital admissions in coastal cities. Air pollution also had a stronger adverse effect on females. These results can help policymakers to determine the region-specific prevention measures for specific types of air pollution.

Supplementary Materials

Supplementary File 1

Acknowledgments

This study was supported by the Natural Science Foundation of Shandong Province (ZR2016HB49) and the Fundamental Research Funds of Shandong University (2016HW001).

Author Contributions

Yi Liu and Fuzhong Xue has full access to the data in this study and take responsibility for the integrity of the data and the accuracy of the data analysis. Yi Liu and Xue were responsible for the conception and study design. Jingjie Sun and Yannong Gou were responsible for the hospitalization data collection. Yi Liu, Xiubin Sun, Xiujun Li and Zhongshang Yuan conducted statistical analysis. Lizhi Kong was responsible for the air pollution data collection.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Descriptive statistics of the daily mean PM2.5 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table A1. Descriptive statistics of the daily mean PM2.5 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionMin25%50%75%MaxMean
Binzhou7.742.265.496.0336.775.1
Dezhou5.842.067.698.6307.179.4
Dongying7.031.052.079.3269.160.4
Heze9.248.770.1108.2427.485.8
Jinan13.345.063.595.3382.176.0
Jining11.842.960.893.1366.174.1
Laiwu16.145.768.9101.2323.680.8
Liaocheng9.750.974.8112.3529.889.7
Linyi11.234.456.289.3345.071.5
Qingdao7.723.637.257.9308.446.4
Rizhao4.824.040.466.0279.352.3
Taian12.137.655.181.9346.267.0
Weifang9.738.256.584.5336.067.4
Weihai2.319.027.946.2232.335.3
Yantai5.120.431.550.8245.139.5
Zaozhuang8.644.167.4100.4300.680.0
Zibo16.042.961.186.8335.870.7
Table A2. Descriptive statistics of the daily mean PM10 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table A2. Descriptive statistics of the daily mean PM10 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionMin25%50%75%MaxMean
Binzhou13.470.6107.7157.1486.9121.1
Dezhou9.081.7122.4177.2584.3138.0
Dongying15.961.694.0142.8546.6108.6
Heze13.788.7125.8184.6549.5143.6
Jinan21.191.3125.0176.1600.3141.9
Jining9.075.8105.9155.1589.2123.4
Laiwu12.183.3120.0159.0471.6129.6
Liaocheng10.497.2132.9186.4577.9151.0
Linyi19.273.2117.4170.9499.5134.2
Qingdao19.651.872.9107.0454.786.2
Rizhao8.545.673.8117.9414.789.5
Taian14.171.9105.1146.1486.1117.8
Weifang27.884.2113.9162.3460.8128.0
Weihai8.737.851.678.5464.564.3
Yantai13.640.257.183.9369.467.0
Zaozhuang23.190.3129.9178.8542.6142.4
Zibo18.780.4117.5154.3489.4125.6
Table A3. Descriptive statistics of the daily mean SO2 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table A3. Descriptive statistics of the daily mean SO2 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionMin25%50%75%MaxMean
Binzhou6.224.338.760.3248.546.2
Dezhou3.418.728.142.7181.235.0
Dongying3.024.137.356.1196.043.4
Heze9.723.132.846.7135.436.9
Jinan8.220.531.149.7185.239.3
Jining9.229.042.563.9166.449.1
Laiwu2.427.843.164.2217.551.1
Liaocheng3.720.031.147.0123.435.9
Linyi3.214.226.542.5139.032.3
Qingdao5.015.220.127.6130.324.0
Rizhao4.310.716.024.6110.020.2
Taian7.421.532.147.9149.937.7
Weifang6.718.829.150.0213.338.9
Weihai1.29.112.817.892.515.4
Yantai3.611.416.825.577.320.1
Zaozhuang10.425.840.857.6136.044.1
Zibo7.036.955.574.6208.260.8
Table A4. Descriptive statistics of the daily mean NO2 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table A4. Descriptive statistics of the daily mean NO2 levels (μg/m3) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionMin25%50%75%MaxMean
Binzhou7.824.335.249.1124.438.9
Dezhou5.723.632.243.7117.735.9
Dongying2.719.430.544.9103.134.0
Heze5.825.534.346.5102.437.6
Jinan10.931.342.557.7155.145.9
Jining8.029.240.653.6106.242.9
Laiwu8.231.944.458.6116.147.0
Liaocheng4.627.538.050.2115.841.4
Linyi8.529.641.957.2125.745.1
Qingdao7.121.630.841.1107.032.7
Rizhao3.919.530.644.198.233.1
Taian7.528.438.949.4107.540.3
Weifang8.722.631.944.196.535.0
Weihai3.813.519.727.992.722.1
Yantai8.020.729.139.884.631.8
Zaozhuang8.921.729.238.8101.431.4
Zibo13.738.748.858.3102.449.7
Table A5. Descriptive statistics of the daily mean temperature (°C) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table A5. Descriptive statistics of the daily mean temperature (°C) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionMin25%50%75%MaxMean
Binzhou−15.33.816.123.331.713.9
Dezhou−14.33.615.422.931.813.5
Dongying−13.23.916.524.032.614.4
Heze−9.25.816.423.331.314.7
Jinan−12.46.017.424.233.015.3
Jining−10.74.816.223.432.514.5
Laiwu−13.14.015.722.632.013.7
Liaocheng−10.95.216.423.531.514.5
Linyi−11.75.516.623.732.514.8
Qingdao−11.55.515.021.929.413.8
Rizhao−11.45.714.922.231.414.1
Taian−22.8−1.68.414.822.76.7
Weifang−13.04.516.323.633.214.3
Weihai−11.84.914.922.530.613.8
Yantai−11.93.715.022.631.413.3
Zaozhuang−11.05.816.423.432.214.7
Zibo−13.15.217.323.932.914.9
Table A6. Descriptive statistics of the daily mean relative humidity (%) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table A6. Descriptive statistics of the daily mean relative humidity (%) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionMin25%50%75%MaxMean
Binzhou255164779863.7
Dezhou2657708210069.3
Dongying234962739861.7
Heze306071839870.2
Jinan1942557210057.4
Jining3261738210070.9
Laiwu234963779762.3
Liaocheng245467799866.2
Linyi245368789965.5
Qingdao2157728310069.6
Rizhao225572859969.1
Taian7456787.510065.9
Weifang245164749662.6
Weihai185264769763.1
Yantai225466769664.2
Zaozhuang286375839872.9
Zibo204457729858
Table A7. Descriptive statistics of the daily mean wind speed (m/s) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table A7. Descriptive statistics of the daily mean wind speed (m/s) in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionMin25%50%75%MaxMean
Binzhou0.41.41.92.66.42.1
Dezhou0.31.52.23.17.22.4
Dongying0.51.62.02.65.42.2
Heze0.41.52.12.96.32.3
Jinan0.21.52.02.86.42.3
Jining0.21.01.41.94.21.5
Laiwu0.01.21.62.15.71.7
Liaocheng0.11.31.82.67.52.0
Linyi0.21.11.51.95.31.6
Qingdao1.02.32.93.78.43.2
Rizhao0.82.22.93.711.33.1
Taian1.04.15.97.720.56.2
Weifang0.21.41.92.55.92.0
Weihai0.52.13.04.49.63.4
Yantai0.72.02.73.88.73.0
Zaozhuang0.31.31.72.36.01.8
Zibo0.31.41.82.77.32.2
Table A8. Percentage changes of the relative risk of male COPD hospitalizations and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
Table A8. Percentage changes of the relative risk of male COPD hospitalizations and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
CaptionPM2.5PM10SO2NO2
Binzhou−0.512 (−1.24, 0.221)−0.301 (−0.892, 0.294)−1.149 (−2.420, 0.139)−1.540 (−3.690, 0.658)
Dezhou0.807 (0.051, 1.557) *0.612 (0.105, 1.117) *1.824 (−0.278, 3.969)3.467 (0.689, 6.167) *
Dongying−0.681 (−1.834, 0.485)−0.369 (−1.130, 0.398)2.142 (0.204, 4.118) *−1.422 (−4.223, 1.460)
Heze0.199 (−0.355, 0.755)0.315 (−0.067, 0.700)2.407 (0.834, 4.004) *−1.230 (−3.470, 1.062)
Jinan0.533 (−0.097, 1.167)0.373 (−0.017, 0.764)1.167 (−0.279, 2.633)1.939 (0.052, 3.924) *
Jining0.754 (0.182, 1.330) *0.624 (0.257, 0.993) *1.796 (0.666, 2.938) *1.652 (−0.781, 4.145)
Laiwu0.989 (−0.099, 2.090)0.650 (−0.143, 1.450)1.989 (−0.01, 4.029)3.134 (−0.101, 6.474)
Liaocheng0.492 (−0.030, 1.016)−0.247 (−0.619, 0.127)1.235 (−0.256, 2.747)−1.928 (−3.946, 0.133)
Linyi0.628 (0.141, 1.118) *0.404 (0.088, 0.722) *2.033 (0.817, 3.263) *0.752 (−0.535, 2.056)
Qingdao0.425 (−0.291, 1.146)0.391 (−0.076, 0.859)2.933 (0.609, 5.311) *2.715 (0.481, 4.998) *
Rizhao0.631 (−0.608, 1.886)0.574 (−0.329, 1.485)2.840 (−0.614, 6.413)2.175 (−1.222, 5.689)
Taian0.716 (−0.035, 1.473)0.271 (−0.165, 0.709)1.448 (−0.171, 3.094)0.851 (−1.346, 3.096)
Weifang0.414 (−0.097, 0.929)0.360 (0.016, 0.706) *0.612 (−0.299, 1.532)0.744 (−0.755, 2.266)
Weihai0.282 (−1.049, 1.631)−0.165 (−0.914, 0.590)−1.844 (−5.290, 1.727)−2.375 (−5.177, 0.510)
Yantai0.732 (−0.185, 1.657)0.538 (−0.092, 1.173)2.216 (−0.314, 4.809)1.011 (−0.993, 3.056)
Zaozhuang0.747 (0.041, 1.457) *0.501 (0.026, 0.979) *1.279 (−0.346, 2.931)2.034 (−0.571, 4.707)
Zibo0.627 (−0.084, 1.343)0.329 (−0.200, 0.862)1.050 (0.027, 2.082) *1.096 (−1.001, 3.238)
Table A9. Percentage changes of the relative risk of female COPD hospitalizations and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
Table A9. Percentage changes of the relative risk of female COPD hospitalizations and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
CaptionPM2.5PM10SO2NO2
Binzhou−0.796 (−1.723, 0.140)0.390 (−0.112, 0.895)1.616 (0.273, 2.978) *2.889 (0.542, 5.291) *
Dezhou1.033 (0.265, 1.795) *0.641 (0.139, 1.140) *2.196 (0.211, 4.220) *3.981 (1.204, 6.680) *
Dongying−1.328 (−2.664, 0.027)−0.321 (−1.186, 0.551)3.681 (1.000, 6.290) *4.568 (1.312, 7.716) *
Heze0.897 (0.120, 1.679) *0.753 (0.206, 1.303) *2.667 (0.609, 4.767) *3.086 (0.461, 5.779) *
Jinan1.132 (0.452, 1.818) *0.528 (0.087, 0.971) *1.840 (0.264, 3.440) *1.969 (0.099, 3.745) *
Jining0.773 (0.067, 1.485) *0.617 (0.161, 1.075) *2.025 (0.619, 3.450) *1.804 (0.105, 3.531) *
Laiwu1.826 (0.442, 3.229) *1.916 (0.893, 2.950) *2.819 (0.147, 5.564) *4.186 (−0.074, 8.628)
Liaocheng−0.325 (−0.845, 0.198)0.373 (0.004, 0.743) *1.228 (−0.551, 3.040)2.082 (0.361, 3.772) *
Linyi0.672 (0.202, 1.144) *0.426 (0.119, 0.734) *2.551 (1.291, 3.827) *1.669 (0.037, 3.326) *
Qingdao0.767 (−0.147, 1.689)0.561 (−0.079, 1.205)3.501 (0.930, 6.137) *2.798 (0.207, 5.456) *
Rizhao0.774 (−0.285, 1.844)0.535 (−0.234, 1.310)−3.286 (−6.891, 0.458)−1.864 (−4.632, 0.984)
Taian1.030 (0.197, 1.871) *0.575 (0.007, 1.147) *2.444 (0.500, 4.426) *2.703 (0.215, 5.129) *
Weifang0.713 (0.101, 1.329) *0.561 (0.147, 0.977) *1.071 (0.018, 2.135) *1.944 (0.231, 3.687) *
Weihai−1.775 (−3.721, 0.211)0.625 (−0.353, 1.613)4.056 (−0.496, 8.816)−2.327 (−6.040, 1.534)
Yantai1.719 (0.549, 2.903) *1.121 (0.283, 1.967) *4.304 (0.952, 7.767) *−1.751 (−4.356, 0.925)
Zaozhuang1.071 (0.213, 1.937) *0.712 (0.134, 1.292) *2.570 (0.517, 4.666) *4.861 (1.480, 8.354) *
Zibo−0.782 (−1.624, 0.066)−0.450 (−0.993, 0.096)1.398 (0.203, 2.607) *2.180 (−0.320, 4.742)
Table A10. Percentage changes of the relative risk of COPD hospitalizations for patients aged over 65 years and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
Table A10. Percentage changes of the relative risk of COPD hospitalizations for patients aged over 65 years and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
CaptionPM2.5PM10SO2NO2
Binzhou−0.691 (−1.447, 0.071)−0.409 (−0.914, 0.099)1.185 (0.100, 2.259) *1.256 (−0.624, 3.173)
Dezhou1.065 (0.425, 1.701) *0.569 (0.144, 0.992) *1.659 (0.019, 3.272) *4.125 (1.805, 6.390) *
Dongying−0.714 (−1.695, 0.278)0.400 (−0.276, 1.08)2.548 (0.545, 4.510) *2.723 (0.333, 5.056) *
Heze0.606 (0.119, 1.09) *0.419 (0.064, 0.776) *2.891 (1.439, 4.363) *2.389 (0.448, 4.367) *
Jinan0.728 (0.157, 1.295) *0.348 (0.001, 0.697) *1.347 (0.064, 2.646) *1.843 (0.120, 3.596) *
Jining0.620 (0.062, 1.182) *0.646 (0.303, 0.991) *1.711 (0.654, 2.779) *1.723 (−0.134, 3.614)
Laiwu0.927 (−0.082, 1.946)1.135 (0.394, 1.881) *2.030 (0.124, 3.972) *2.619 (−0.460, 5.793)
Liaocheng0.416 (−0.053, 0.887)−0.239 (−0.522, 0.046)1.250 (−0.103, 2.621)2.063 (0.502, 3.598) *
Linyi0.476 (0.017, 0.938) *0.351 (0.051, 0.652) *2.365 (1.222, 3.520) *0.930 (−0.536, 2.417)
Qingdao0.655 (−0.078, 1.393)0.504 (−0.008, 1.018)3.186 (1.065, 5.352) *3.218 (1.150, 5.330) *
Rizhao−0.532 (−1.490, 0.435)0.449 (−0.196, 1.098)2.163 (−0.845, 5.262)−1.107 (−3.428, 1.270)
Taian0.785 (0.148, 1.426) *0.431 (0.002, 0.865) *1.805 (0.325, 3.307) *1.498 (−0.448, 3.482)
Weifang0.463 (−0.105, 1.035)0.388 (0.043, 0.734) *0.559 (−0.325, 1.451)−1.100 (−2.478, 0.299)
Weihai−0.565 (−1.997, 0.887)0.218 (−0.493, 0.933)−2.280 (−6.180, 1.782)−1.765 (−4.492, 1.040)
Yantai0.573 (−0.304, 1.456)0.669 (0.050, 1.291) *1.995 (−0.389, 4.436)−1.428 (−3.290, 0.470)
Zaozhuang0.847 (0.148, 1.55) *0.484 (0.009, 0.961) *1.766 (0.181, 3.377) *3.761 (1.148, 6.442) *
Zibo0.583 (−0.097, 1.267)0.346 (−0.092, 0.787)1.044 (0.126, 1.971) *1.137 (−0.543, 2.845)
Table A11. Percentage changes of the relative risk of COPD hospitalizations for patients aged under 65 years and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
Table A11. Percentage changes of the relative risk of COPD hospitalizations for patients aged under 65 years and their 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
CaptionPM2.5PM10SO2NO2
Binzhou−0.988 (−2.092, 0.128)−0.635 (−1.369, 0.104)2.391 (0.395, 4.347) *3.063 (−0.921, 7.207)
Dezhou1.395 (0.320, 2.458) *1.103 (0.390, 1.810) *2.134 (−0.595, 4.938)5.770 (2.022, 9.374) *
Dongying−1.874 (−4.052, 0.352)−0.785 (−2.115, 0.562)−2.630 (−6.456, 1.352)−4.267 (−9.071, 0.791)
Heze0.465 (−0.496, 1.434)0.425 (−0.178, 1.032)1.804 (−0.725, 4.397)−3.105 (−6.244, 0.140)
Jinan0.671 (−0.216, 1.566)0.460 (−0.102, 1.025)1.400 (−0.683, 3.527)−1.546 (−3.882, 0.847)
Jining1.227 (0.419, 2.041) *0.641 (0.117, 1.167) *2.033 (0.429, 3.661) *−1.926 (−4.304, 0.512)
Laiwu1.492 (−0.193, 3.205)1.011 (−0.222, 2.259)3.206 (0.461, 6.027) *7.241 (1.718, 13.064) *
Liaocheng0.602 (−0.028, 1.236)0.400 (−0.105, 0.908)1.790 (−0.442, 4.071)−0.933 (−3.473, 1.674)
Linyi0.934 (0.340, 1.532) *0.476 (0.061, 0.893) *2.054 (0.468, 3.666) *1.806 (0.140, 3.501) *
Qingdao−0.924 (−1.994, 0.157)−0.316 (−1.030, 0.403)2.965 (−0.340, 6.380)2.849 (0.144, 5.480) *
Rizhao0.854 (−0.850, 2.587)0.754 (−0.216, 1.733)2.222 (−2.557, 7.235)2.463 (−1.233, 6.297)
Taian0.852 (−0.181, 1.896)−0.413 (−1.034, 0.212)1.710 (−0.516, 3.985)−1.951 (−4.626, 0.800)
Weifang0.915 (0.226, 1.608) *0.711 (0.244, 1.181) *−1.502 (−3.008, 0.027)2.676 (0.610, 4.785) *
Weihai−0.718 (−2.688, 1.291)0.703 (−0.297, 1.713)1.144 (−3.569, 6.087)−2.033 (−5.888, 1.980)
Yantai2.192 (0.843, 3.558) *1.584 (0.656, 2.521) *5.282 (1.384, 9.330) *2.677 (−0.473, 5.926)
Zaozhuang1.073 (0.006, 2.151) *0.750 (0.034, 1.471) *1.681 (−0.791, 4.213)2.202 (−1.837, 6.406)
Zibo−0.994 (−2.180, 0.207)−0.780 (−1.591, 0.037)1.606 (−0.108, 3.348)−2.753 (−5.729, 0.318)

References

  1. Katsouyanni, K. Ambient air pollution and health. Br. Med. Bull. 2003, 68, 143–156. [Google Scholar] [CrossRef] [PubMed]
  2. Brauer, M.; Freedman, G.; Frostad, J.; Van Donkelaar, A.; Martin, R.V.; Dentener, F.; Dingenen, R.V.; Estep, K.; Amini, H.; Apte, J.S.; et al. Ambient air pollution exposure estimation for the global burden of disease 2013. Environ. Sci. Technol. 2015, 50, 79–88. [Google Scholar] [CrossRef] [PubMed]
  3. Lelieveld, J.; Evans, J.S.; Fnais, M.; Giannadaki, D.; Pozzer, A. The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature 2015, 525, 367–371. [Google Scholar] [CrossRef] [PubMed]
  4. Raaschou-Nielsen, O.; Andersen, Z.J.; Beelen, R.; Samoli, E.; Stafoggia, M.; Weinmayr, G.; Hoffmann, B.; Fischer, P.; Nieuwenhuijsen, M.J.; Brunekreef, B.; et al. Air pollution and lung cancer incidence in 17 European cohorts: Prospective analyses from the European Study of Cohorts for Air Pollution Effects (ESCAPE). Lancet Oncol. 2013, 14, 813–822. [Google Scholar] [CrossRef]
  5. Brunekreef, B.; Holgate, S.T. Air pollution and health. Lancet 2002, 360, 1233–1242. [Google Scholar] [CrossRef]
  6. World Health Organization. Burden of Disease Attributable to Outdoor Air Pollution; World Health Organization: Geneva, Switzerland, 2011. [Google Scholar]
  7. Lim, S.S.; Vos, T.; Flaxman, A.D.; Danaei, G.; Shibuya, K.; Adair-Rohani, H.; AlMazroa, M.A.; Amann, M.; Anderson, H.R.; Andrews, K.G.; et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: A systematic analysis for the Global Burden of Disease Study 2010. Lancet 2013, 380, 2224–2260. [Google Scholar] [CrossRef]
  8. Stieb, D.M.; Judek, S.; Burnett, R.T. Meta-Analysis of Time-Series Studies of Air Pollution and Mortality: Effects of Gases and Particles and the Influence of Cause of Death, Age, and Season. J. Air Waste Manag. Assoc. 2002, 52, 470–484. [Google Scholar] [CrossRef] [PubMed]
  9. Sun, Q.; Hong, X.; Wold, L.E. Cardiovascular Effects of Ambient Particulate Air Pollution Exposure. Circulation 2010, 121, 2755–2765. [Google Scholar] [CrossRef] [PubMed]
  10. Dominici, F.; Peng, R.D.; Bell, M.L.; Pham, L.; McDermott, A.; Zeger, S.L.; Samet, J.M. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA 2006, 295, 1127–1134. [Google Scholar] [CrossRef] [PubMed]
  11. Brook, R.D.; Rajagopalan, S.; Pope, C.A.; Brook, J.R.; Bhatnagar, A.; Diez-Roux, A.V.; Holguin, F.; Hong, Y.; Luepker, R.V.; Mittleman, M.A.; et al. Particulate matter air pollution and cardiovascular disease. Circulation 2010, 121, 2331–2378. [Google Scholar] [CrossRef] [PubMed]
  12. Schikowski, T.; Sugiri, D.; Ranft, U.; Gehring, U.; Heinrich, J.; Wichmann, H.E.; Krämer, U. Long-term air pollution exposure and living close to busy roads are associated with COPD in women. Respir. Res. 2005, 6, 1–10. [Google Scholar] [CrossRef] [PubMed]
  13. Rabe, K.F.; Hurd, S.; Anzueto, A.; Jones, P.W.; Vogelmeier, C.; Anzueto, A.; Barnes, P.J.; Fabbri, L.M.; Martinez, F.J.; Nishimura, M.; et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am. J. Respir. Crit. Care Med. 2013, 187, 347–365. [Google Scholar] [CrossRef] [PubMed]
  14. World Health Organization. Chronic Obstructive Pulmonary Disease (COPD). Available online: http://www.who.int/respiratory/copd/en/ (accessed on 10 February 2018).
  15. World Health Organization. The Global Burden of Disease. 2004 Update; World Health Organization: Geneva, Switzerland, 2008. [Google Scholar]
  16. Chan, C.K.; Yao, X. Air pollution in mega cities in China. Atmos. Environ. 2008, 42, 1–42. [Google Scholar] [CrossRef]
  17. Rohde, R.A.; Muller, R.A. Air pollution in China: Mapping of concentrations and sources. PLoS ONE 2015, 10, e0135749. [Google Scholar] [CrossRef] [PubMed]
  18. Huang, R.J.; Zhang, Y.; Bozzetti, C.; Ho, K.F.; Cao, J.J.; Han, Y.; Daellenbach, K.R.; Slowik, J.G.; Platt, S.M.; Canonaco, F.; et al. High secondary aerosol contribution to particulate pollution during haze events in China. Nature 2014, 514, 218–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Cao, J.; Yang, C.; Li, J.; Chen, R.; Chen, B.; Gu, D.; Kan, H. Association between long-term exposure to outdoor air pollution and mortality in China: A cohort study. J. Hazard. Mater. 2011, 186, 1594–1600. [Google Scholar] [CrossRef] [PubMed]
  20. Chen, R.; Kan, H.; Chen, B.; Huang, W.; Bai, Z.; Song, G.; Pan, G.; CAPES Collaborative Group. Association of Particulate Air Pollution with Daily Mortality: The China Air Pollution and Health Effects Study. Am. J. Epidemiol. 2012, 175, 1173–1181. [Google Scholar] [CrossRef] [PubMed]
  21. Chen, Z.; Wang, J.N.; Ma, G.X.; Zhang, Y.S. China tackles the health effects of air pollution. Lancet 2013, 382, 1959–1960. [Google Scholar] [CrossRef]
  22. Zhou, M.; He, G.; Liu, Y.; Yin, P.; Li, Y.; Kan, H.; Fan, M.; Xue, A.; Fan, M. The associations between ambient air pollution and adult respiratory mortality in 32 major Chinese cities, 2006–2010. Environ. Res. 2015, 137, 278–286. [Google Scholar] [CrossRef] [PubMed]
  23. Schwartz, J.; Morris, R. Air pollution and hospital admissions for cardiovascular disease in Detroit, Michigan. Am. J. Epidemiol. 1995, 142, 23–35. [Google Scholar] [CrossRef] [PubMed]
  24. Schwartz, J. Air pollution and hospital admissions for heart disease in eight US counties. Epidemiology 1999, 10, 17–22. [Google Scholar] [CrossRef] [PubMed]
  25. Bell, M.L.; Ebisu, K.; Peng, R.D.; Samet, J.M.; Dominici, F. Hospital admissions and chemical composition of fine particle air pollution. Am. J. Respir. Crit. Care Med. 2009, 179, 1115–1120. [Google Scholar] [CrossRef] [PubMed]
  26. Peng, R.; Chang, H.; Bell, M.; McDermott, A.; Zeger, S.; Samet, J.; Dominici, F. Coarse particulate matter air pollution and hospital admissions for cardiovascular and respiratory diseases among Medicare patients. JAMA 2008, 299, 2172–2179. [Google Scholar] [CrossRef] [PubMed]
  27. Samoli, E.; Atkinson, R.W.; Analitis, A.; Fuller, G.W.; Green, D.C.; Mudway, I.; Anderson, H.R.; Kelly, F.J. Associations of short-term exposure to traffic-related air pollution with cardiovascular and respiratory hospital admissions in London, UK. Occup. Environ. Med. 2016, 73, 300–307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Atkinson, R.W.; Kang, S.; Anderson, H.R.; Mills, I.C.; Walton, H.A. Epidemiological time series studies of PM2.5 and daily mortality and hospital admissions: A systematic review and meta-analysis. Thorax 2014, 69, 660–665. [Google Scholar] [CrossRef] [PubMed]
  29. Dockery, D.W.; Rich, D.Q.; Goodman, P.G.; Clancy, L.; Ohman-Strickland, P.; George, P.; Kotlov, T. Effect of air pollution control on mortality and hospital admissions in Ireland. Res. Rep. Health Eff. Inst. 2013, 176, 3–109. [Google Scholar]
  30. Lanki, T.; Pekkanen, J.; Aalto, P.; Elosua, R.; Berqlind, N.; D’lppoliti, D.; Kulmala, M.; Nyberg, F.; Peters, A.; Picciotto, S.; et al. Associations of traffic related air pollution levels with hospitalization for first acute myocardial infarction: The HEAPSS study. Occup. Environ. Med. 2006, 63, 844–851. [Google Scholar] [CrossRef] [PubMed]
  31. Gan, W.Q.; FitzGerald, J.M.; Carlsten, C.; Sadatsafavi, M.; Brauer, M. Associations of ambient air pollution with chronic obstructive pulmonary disease hospitalization and mortality. Am. J. Respir. Crit. Care Med. 2013, 187, 721–727. [Google Scholar] [CrossRef] [PubMed]
  32. Xiong, Q.; Zhao, W.; Gong, Z.; Zhao, W.; Tang, T. Fine Particulate Matter Pollution and Hospital Admissions for Respiratory Diseases in Beijing, China. Int. J. Environ. Res. Public Health 2015, 12, 11880–11892. [Google Scholar] [CrossRef] [PubMed]
  33. Tao, Y.; Mi, S.; Zhou, S.; Wang, S.; Xie, X. Air pollution and hospital admissions for respiratory diseases in Lanzhou, China. Environ. Pollut. 2014, 185, 196–201. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, Z.; Wang, J.; Chen, L.; Chen, X.; Sun, G.; Zhong, N.; Kan, H.; Lu, W. Impact of haze and air pollution-related hazards on hospital admissions in Guangzhou, China. Environ. Sci. Pollut. Res. 2014, 21, 4236–4244. [Google Scholar] [CrossRef] [PubMed]
  35. Chen, R.; Chu, C.; Tan, J.; Cao, J.; Song, W.; Xu, X.; Jiang, C.; Ma, W.; Yang, C.; Chen, B.; et al. Ambient air pollution and hospital admission in Shanghai, China. J. Hazard. Mater. 2010, 181, 234–240. [Google Scholar] [CrossRef] [PubMed]
  36. Xiang, H.; Mertz, K.J.; Arena, V.C.; Brink, L.L.; Xu, X.; Bi, Y.; Talbott, E.O. Estimation of short-term effects of air pollution on stroke hospital admissions in Wuhan, China. PLoS ONE 2013, 8, e61168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hastie, T.J.; Tibshirani, R.J. Generalized additive models. In Statistical Models in S; CRC Press: Boca Raton, FL, USA, 1990. [Google Scholar]
  38. Wood, S. Generalized Additive Models: An Introduction with R; CRC Press: Boca Raton, FL, USA, 2006. [Google Scholar]
  39. Wahba, G. Spline Models for Observational Data; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1990. [Google Scholar]
  40. Gu, C. Smoothing Spline ANOVA Models; Springer: New York, NY, USA, 2013. [Google Scholar]
  41. Ministry of Environmental Protection of the People’s Republic of China. 2016. Available online: http://kjs.mep.gov.cn/hjbhbz/bzwb/dqhjbh/dqhjzlbz/201203/t20120302_224165.shtml (accessed on 10 February 2018).
  42. World Health Organization. Air Quality Guidelines: Global Update 2005: Particulate Matter, Ozone, Nitrogen Dioxide, and Sulfur Dioxide; World Health Organization: Geneva, Switzerland, 2006. [Google Scholar]
  43. Kan, H.; London, S.J.; Chen, G.; Zhang, Y.; Song, G.; Zhao, N.; Jiang, L.; Chen, B. Season, sex, age, and education as modifiers of the effects of outdoor air pollution on daily mortality in shanghai, china: The public health and air pollution in Asia (papa) study. Environ. Health Perspect. 2008, 116, 1183–1188. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The locations of the 17 cities in the Shandong Province; the coastal cities are shown in the green color. The numbers represent the count of air quality monitoring stations in each city.
Figure 1. The locations of the 17 cities in the Shandong Province; the coastal cities are shown in the green color. The numbers represent the count of air quality monitoring stations in each city.
Ijerph 15 00774 g001
Figure 2. The daily maximum air pollution levels in the 17 cities of the Shandong Province.
Figure 2. The daily maximum air pollution levels in the 17 cities of the Shandong Province.
Ijerph 15 00774 g002
Figure 3. The percentage change in the relative risk of the total COPD hospitalizations per 10-μg/m3 increase in air pollution levels in the 17 cities of the Shandong province.
Figure 3. The percentage change in the relative risk of the total COPD hospitalizations per 10-μg/m3 increase in air pollution levels in the 17 cities of the Shandong province.
Ijerph 15 00774 g003
Figure 4. A forest plot of the percentage change in the relative risk of the total COPD hospitalizations with a 95% confidence interval (CI) per 10-μg/m3 increase in the air pollution levels in the 17 cities of the Shandong province.
Figure 4. A forest plot of the percentage change in the relative risk of the total COPD hospitalizations with a 95% confidence interval (CI) per 10-μg/m3 increase in the air pollution levels in the 17 cities of the Shandong province.
Ijerph 15 00774 g004
Table 1. Descriptive statistics of the daily hospital admissions grouped by sex and age in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
Table 1. Descriptive statistics of the daily hospital admissions grouped by sex and age in the 17 cities of the Shandong province (from 1 January 2015 to 31 December 2016).
CaptionTotalMaleFemale≥65<65
Binzhou11,1035700540389722131
Dezhou62883134315447801508
Dongying4026229917273320706
Heze16,05511,007504812,5933462
Jinan13,7118383532810,7073004
Jining19,49312,884660914,6114882
Laiwu2823180010232227596
Liaocheng15,3689632573611,7833585
Linyi25,50813,83011,67818,4027106
Qingdao15,7979613618411,6954102
Rizhao56753118255741281547
Taian10,7527052370077093043
Weifang26,65315,83210,82119,4397214
Weihai72104612259849292281
Yantai11,8487738411087853063
Zaozhuang12,6237896472795373086
Zibo11,0256848417788222203
Total215,958131,37884,580162,43953,519
Table 2. The percentage change in the relative risk of the total COPD hospitalizations and its 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
Table 2. The percentage change in the relative risk of the total COPD hospitalizations and its 95% confidence interval (CI) per 10-μg/m3 increase in air pollutants in the 17 cities of the Shandong province, where * indicates statistically significant estimates (p < 0.05).
CityPM2.5LagPM10LagSO2LagNO2Lag
Binzhou−0.541 (−1.23, 0.153)00.172 (−0.202, 0.548)31.374 (0.383, 2.356) *51.168 (−0.555, 2.921)3
Dezhou0.895 (0.312, 1.474) *20.481 (0.095, 0.865) *21.427 (0.089, 2.970) *03.515 (1.394, 5.590) *2
Dongying−0.69 (−1.604, 0.232)3−0.305 (−0.902, 0.294)32.787 (0.926, 4.612) *03.316 (0.573, 5.984) *0
Heze0.279 (−0.203, 0.764)10.357 (0.022, 0.694) *32.558 (1.186, 3.948) *21.534 (−0.301, 3.402)3
Jinan0.501 (0.003, 1.018) *40.332 (0.008, 0.658) *41.405 (0.203, 2.621) *31.583 (−0.025, 3.218)0
Jining0.729 (0.232, 1.229) *20.605 (0.285, 0.925) *21.713 (0.734, 2.701) *41.164 (−0.309, 2.659)4
Laiwu0.978 (0.105, 1.858) *31.099 (0.458, 1.743) *31.530 (−0.132, 3.219)22.233 (−0.427, 4.965)2
Liaocheng0.388 (−0.053, 0.830)00.262 (−0.052, 0.576)01.217 (0.036, 2.485) *31.697 (0.247, 3.126) *5
Linyi0.501 (0.081, 0.922) *10.372 (0.099, 0.646) *12.276 (1.234, 3.328) *20.599 (−0.511, 1.721)2
Qingdao0.536 (−0.111, 1.188)00.388 (−0.064, 0.843)03.164 (1.302, 5.059) *02.824 (1.004, 4.677) *0
Rizhao0.635 (−0.179, 1.455)20.510 (−0.065, 1.088)31.916 (−0.772, 4.677)41.187 (−1.447, 3.890)0
Taian0.561 (0.004, 1.140) *20.314 (0.018, 0.678) *31.490 (0.167, 2.830) *2−1.501 (−3.191, 0.218)5
Weifang0.448 (0.007, 0.890) *20.406 (0.109, 0.704) *50.656 (−0.156, 1.475)21.245 (−0.052, 2.559)2
Weihai−0.667 (−1.907, 0.590)00.254 (−0.361, 0.873)2−1.717 (−5.102, 1.789)0−1.756 (−4.111, 0.657)2
Yantai0.658 (−0.120, 1.442)50.470 (−0.069, 1.012)21.768 (−0.357, 3.938)20.765 (−0.948, 2.508)2
Zaozhuang0.904 (0.256, 1.555) *10.542 (0.103, 0.983) *11.484 (0.015, 2.973) *53.013 (0.596, 5.488) *6
Zibo0.480 (−0.102, 1.066)3−0.189 (−0.567, 0.190)41.210 (0.375, 2.052) *31.602 (0.130, 3.053) *6

Share and Cite

MDPI and ACS Style

Liu, Y.; Sun, J.; Gou, Y.; Sun, X.; Li, X.; Yuan, Z.; Kong, L.; Xue, F. A Multicity Analysis of the Short-Term Effects of Air Pollution on the Chronic Obstructive Pulmonary Disease Hospital Admissions in Shandong, China. Int. J. Environ. Res. Public Health 2018, 15, 774. https://doi.org/10.3390/ijerph15040774

AMA Style

Liu Y, Sun J, Gou Y, Sun X, Li X, Yuan Z, Kong L, Xue F. A Multicity Analysis of the Short-Term Effects of Air Pollution on the Chronic Obstructive Pulmonary Disease Hospital Admissions in Shandong, China. International Journal of Environmental Research and Public Health. 2018; 15(4):774. https://doi.org/10.3390/ijerph15040774

Chicago/Turabian Style

Liu, Yi, Jingjie Sun, Yannong Gou, Xiubin Sun, Xiujun Li, Zhongshang Yuan, Lizhi Kong, and Fuzhong Xue. 2018. "A Multicity Analysis of the Short-Term Effects of Air Pollution on the Chronic Obstructive Pulmonary Disease Hospital Admissions in Shandong, China" International Journal of Environmental Research and Public Health 15, no. 4: 774. https://doi.org/10.3390/ijerph15040774

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop