Next Article in Journal
Dog and Cat Ownership Predicts Adolescents’ Mental Well-Being: A Population-Based Longitudinal Study
Previous Article in Journal
The Impact of Online Health Information on Patient Health Behaviours and Making Decisions Concerning Health
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Short-Term Effects of Air Pollution on Cardiovascular Disease Using Bayesian Spatio-Temporal Models

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 250014, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(3), 879; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17030879
Submission received: 16 January 2020 / Revised: 25 January 2020 / Accepted: 26 January 2020 / Published: 31 January 2020
(This article belongs to the Section Environmental Health)

Abstract

:
There has been an increasing number of clinical and epidemiologic research projects providing supporting evidence that short-term exposure to ambient air pollution contributes to the exacerbation of cardiovascular disease. However, few studies consider measurement error and spatial effects in the estimate of underlying air pollution levels, and less is known about the influence of baseline air pollution levels on cardiovascular disease. We used hospital admissions data for cardiovascular diseases (CVD) collected from an inland, heavily polluted city and a coastal city in Shandong Province, China. Bayesian spatio-temporal models were applied to obtain the underlying pollution level in each city, then generalized additive models were adopted to assess the health effects. The total cardiovascular disease hospitalizations were significantly increased in the inland city by 0.401% (0.029, 0.775), 0.316% (0.086, 0.547), 0.903% (0.252, 1.559), and 2.647% (1.607, 3.697) per 10 μg/m3 increase in PM2.5, PM10, SO2, and NO2, respectively. The total cardiovascular diseases hospitalizations were increased by 6.568% (3.636, 9.584) per 10μg/m3 increase in the level of NO2. Although the air pollution overall had a more significant adverse impact on cardiovascular disease hospital admissions in the heavily polluted inland city, the short-term increases in air pollution levels in the less polluted coastal areas led to excessive exacerbations of cardiovascular disease.

1. Introduction

Cardiovascular diseases (CVD) cause 17.7 million deaths every year (31% of all global deaths), and are the most common cause of death and disability around the world [1,2,3,4,5]. Many studies across the world have shown that ambient air pollution contributes to the increase of mortality, morbidity, and hospitalizations of CVD [6,7,8,9,10,11,12,13]. Without a precise control of the sources of pollution and the application of strict policies, air pollution has become one of the major public health concerns in China [14,15,16,17,18,19].
Researchers have previously conducted a series of analyses on the effect of ambient air pollution on CVD in China [20,21,22,23,24,25,26,27]; however, few studies have considered measurement error and spatial effects in the estimate of the underlying air pollution levels. Health data are generally only available in an aggregated form for an entire administrative area, but daily measurements of air pollution consist of concentrations from a set of different monitoring sites, meaning that a corresponding single daily measure of pollution is required as a predictor in the regression model. The most commonly used approach is to simply average the measurements over spatial locations, and to use the resulting summary in a health model. However, this simple summary ignores any measurement error or spatial variability in the underlying pollution field and, as such, subsequent estimates of the health effects may be biased [28]. Different baselines of air pollution concentration may lead to different health outcomes when there is a short-term increase in the air pollution level, a topic few existing studies have focused on.
The inland city of Jinan is located in a hilly area of Shandong province, East China. It is an industrial city that suffers from serious air pollution, largely due to the oil and coal burning from industrial processes and vehicle emissions. In contrast, the coastal city of Weihai is one of the least polluted regions in China. It has been named the “National Hygienic City” and the “National Garden City”. In 2003, the Weihai Municipal Government was awarded the “UN-Habitat Scroll of Honor”, for outstanding improvements in shelter and the urban environment [29].
In this study, we proposed a Bayesian spatio-temporal model to accurately estimate the underlying air pollution level in the study regions, then assessed the short-term effect of air pollutants, including PM2.5, PM10, SO2, and NO2, on the exacerbation of CVD by analyzing the data of hospital admissions collected from all the tertiary-level and secondary-level hospitals in Jinan and Weihai, both in Shandong Province. Based on the analysis results, we further looked into the baseline influence on the health effect of the air pollution. This study thus provides public health evidence for developing policies and health management strategies for different ambient pollution environments.

2. Materials and methods

2.1. Hospital Admissions

Hospital admission records for CVD between 1st January, 2014 and 31st December, 2016 were obtained from all tertiary-level and secondary-level hospitals in Jinan and Weihai. Each record contained 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). Cardiovascular disease outcomes were considered, based on the ICD-10 codes: I00–99.

2.2. Air Pollution and Meteorological Data

The monitoring stations established by both the China National Environmental Protection Administration and Shandong local government were considered in this study. There were two sites excluded from the study in the city of Jinan, both of which were new monitoring stations established in May, 2015. Hence, in total, 15 and 9 stations in Jinan and Weihai, respectively, were included in the analysis, as shown in Figure 1. The maximum distance between monitors was 74.64 km and 87.64 km in Jinan and Weihai, respectively. Each air quality monitoring station provided hourly readings of PM2.5, PM10, SO2, and NO2 levels, between 1st January, 2014 and 31st December, 2016. 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 daily mean, and the highest and lowest values of the monitored objects. We considered the daily mean temperature, relative humidity, wind speed, and pressure as confounding factors of weather in the analysis.

2.3. Method of Analysis

In this paper, we proposed a two-stage approach: firstly, we modeled the underlying pollution concentrations using a Bayesian spatio-temporal model; secondly, we assessed their relationship with the health data, using a frequentist generalized additive model. We implemented the first stage Bayesian spatio-temporal model using a Markov chain Monte Carlo (MCMC) simulation. There were three parts of this hierarchical model:
ln ( W s t ) N ( Z t + m s , σ ε 2 )
(1)
The measurement error level: We modeled the air pollution data on the log scale because the pollution concentrations were nonnegative and often skewed to the right. The combination of true underlying process Z t and spatial effect m s is linked with the observed pollution concentrations using a measurement error model, where the measurement errors are independent with zero mean and constant variance of σ ε 2
Z t N ( μ + ρ ( Z t 1 μ ) , σ z 2 ) m s M V N ( 0 , σ m 2 Σ ( ϕ ) )
(2)
The underlying process level: The true levels were modeled by a first order autoregressive process which had variance σ z 2 and lag one correlation coefficient ρ . The spatial structure was represented by a set of zero-mean Gaussian random effects m s , which had variance σ m 2 and correlation matrix Σ ( ϕ ) . The latter was constructed using the Matern class of functions, with a smoothness parameter fixed at 0.5. This gave an exponential correlation structure C o r r [ s , s ] = exp ( s s ϕ ) , where . denotes Euclidean distance.
σ ε 2 I n v e r s e G a m m a ( 0.01 , 0.01 ) σ z 2 I n v e r s e G a m m a ( 0.01 , 0.01 ) σ m 2 I n v e r s e G a m m a ( 0.01 , 0.01 ) ϕ D i s c r e t e U n i f o r m ( a ϕ , b ϕ )
(3)
The parameter level: The spatial range parameter ϕ was assigned a discrete uniform prior with possible values ranging between the minimum and maximum distance between monitoring locations [30]. The variance parameters σ ε 2 , σ z 2 , and σ m 2 were assigned weekly informative priors [31].
In order to draw inference using this Bayesian model, the MCMC algorithm was run for 40,000 iterations, discarding the first 10,000 as “burn in” samples. Finally, the posterior mean values were summarized as estimates of the true underlying air pollution levels and plugged into the second stage health models.
In the second stage, we applied generalized additive Poisson models in frequentist settings to examine the associations between air pollution (PM2.5, PM10, SO2, and NO2) for 0–6 day lags and day-to-day variations in the CVD hospital admissions for the total number of patients, stratified by sex and age. We also fitted multipollutant models to assess the stability of air pollution health effects, where daily air pollution exposures were included jointly, at the same lag as in the models. A quasi-Poisson distribution was adopted, in order to overcome the overdispersion of hospital admissions data. The smoothed functions captured the nonlinear relationship between daily hospital admissions and the time-varying covariates [32,33], while the degrees of freedom of natural splines were determined by the Akaike’s information criterion [34]; therefore, calendar time (7 df per year) and the meteorological variables, including temperature (4 df), relative humidity (4 df), wind speed (3 df), and pressure (5 df), were adopted. In addition, the day of the week (DOW) was also included as a dummy variable in the models. The MGVC package (v.1.8–17) [33] in R v.3.4.3 [35] was applied, to construct the health models. Then, all the 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

In this study, the number of CVD hospital admissions during the study period was 306,963 and 91,232 in Jinan and Weihai, respectively. In total, there were 219,121 male patients and 179,074 female patients. Moreover, there were 209,554 patients that were over 65 years old, and 188,641 records were for patients aged under 65. The details of the hospital admissions in the two cities can be found in Table 1. Figure 2 and Figure 3 present the boxplot of the daily mean air pollution levels and the meteorological factors during the study period in Jinan and Weihai. The daily maximum values of the air pollution concentrations in Jinan were above 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) [36], while the daily maximum concentrations of SO2 and NO2 in Weihai were below the standard. There was an obvious spatial heterogeneity among the distribution of ambient air pollution levels, such that the pollutant measurements in the coastal city of Weihai were much lower than the ones in the inland city of Jinan.
Figure 4 presents the results for lags over 0–6 days from the single-pollutant models, including the effect estimates and corresponding 95% confidence intervals (CIs) of the percentage changes of relative risk (RR) for the total CVD admissions per 10μg/m3 increase in PM2.5, PM10, SO2, and NO2 concentrations, in Jinan and Weihai. There were significant associations between all types of pollutants and the CVD admissions in Jinan. The total CVD hospitalizations significantly increased by 0.401% (0.029, 0.775), 0.316% (0.086, 0.547), 0.903% (0.252, 1.559), and 2.647% (1.607, 3.697) per 10μg/m3 increase in PM2.5, PM10, SO2, and NO2, respectively. The health effects of air pollution were less dramatic in the coastal city of Weihai, although the total CVD hospitalizations increased by 6.568% (3.636, 9.584) per 10 μg/m3 increase in NO2 at lag 0.
The results of the multipollutant models are shown in Table 2. For total CVD hospitalizations in Jinan, the effect of PM2.5, PM10, and SO2 become insignificant when adjusted for NO2. Then, the effect of PM2.5 and PM10 becomes insignificant when adjusted for each other. Moreover, the association between an increase in the levels of SO2 and NO2 and the total hospital admissions was consistent when adjusted for PM2.5 and PM10. For total CVD hospital admissions in Weihai, the effect of NO2 was consistent when adjusted for other pollutants, such that the total hospital admissions were significantly increased by 9.606% (5.887, 13.456), 8.453% (5.048, 11.97), and 8.415% (4.448, 12.532), per 10 μg/m3 increase in NO2, when adjusted for PM2.5, PM10, and SO2, respectively.
Table 3 shows the results of the subgroup analysis. Regarding gender, stronger effects were observed for females in Jinan, except for PM2.5, such that the significant effects of PM2.5 were 0.404% (0.009, 0.801) and 0.398% (0.013, 0.783) on male and female patients, respectively, both at lag 0. However, we found the opposite results in Weihai, where the air pollution had a stronger effect on male patients. For example, the SO2 effect on male patients was significant (2.859% (0.340, 5.441) at lag 0), while the effect on female patients was insignificant. In addition, the significant NO2 effects were 7.419% (4.031, 10.918) and 5.535% (2.070, 9.119) for male and female patients, respectively, at lag 0. Regarding age, the results varied between the two cities. The air pollution effects on the younger patients were stronger in Jinan, for example, the PM2.5 effect on aged patients in Jinan was insignificant, but the effect on patients under 65 years of age was significant: 0.485% (0.074, 0.899) at lag 0. In Weihai, the air pollution effects on the aged patients were more severe. The younger patients’ hospitalizations were significantly affected only by the increase in NO2 levels, but PM2.5, SO2, and NO2 effects on aged patients were significant. The NO2 effects were 5.236% (1.591, 9.011) and 7.612% (4.321, 11.006) on younger and aged patients, respectively.

4. Discussion

In this study, we assessed the effects of short-term exposure to ambient air pollution on the exacerbation of CVD, using hospital admissions data from a coastal and an inland city in Shandong province, China. Since the health data were collected from tertiary-level and secondary-level hospitals, all the patients went through a reasonable diagnosis and treatment during the admission time; therefore, the accuracy of diagnosis and the reliability of the datasets can be guaranteed.
The public health burden of the present exposure level is substantial. As air pollution is experienced by everybody, even a small increase in risk can be associated with very large numbers of people being ill [37,38,39]. It is therefore extremely important that the risks to health are estimated accurately, and for that to happen there needs to be accurate recording of the levels of exposure that might be experienced by populations at risk. In this analysis, a Bayesian spatio-temporal model was applied to obtain the underlying pollution level. This method effectively captures the spatial and temporal features of the ambient air pollution data, and also considers the measurement error term in the data. As a result, the corresponding estimated percentage change from the health model is more adequate than the ones obtained from the model when simply taking the average over all of the air pollution data collected at the monitoring stations in the same area.
There are likely to be computational considerations associated with jointly fitting the health and exposure models in Bayesian settings, especially if the latter uses large amounts of data over space and time. When the exposure model is complicated, or when one is interested in running multiple candidate epidemiological models with different sets of covariates, a single model is not going to provide an efficient method of investigation. We therefore adopted two-stage model settings in order to ease the computational burden of running a combined model (this has been adopted in a number of cases [28,40,41,42,43]). The two-stage approach has the advantage that the exposure model, which is likely to be the most computationally demanding, does not require a refit when running multiple health effect analyses. Moreover, we implemented the health models in a frequentist setting, as generalized additive models are the most commonly used models in epidemiological studies. The frequentist modeling obtained the confidence intervals of the estimates, instead of credible intervals deducted from a Bayesian method, which benefited us when comparing our results with the ones from other epidemiological studies.
We found that the associations between air pollution and CVD hospitalizations were spatially heterogeneous, and that air pollution contributed, overall, more prominent adverse health effects in the inland city. Jinan contains significant oil refineries and metallurgical and mechanical sectors and, moreover, vehicle emissions are one of the major causes of the heavy air pollution in the city. Furthermore, geographical factors are also influential, as Jinan is situated in a typical valley basin and is surrounded by the largest mountain group in East China, making it difficult for the emitted ambient pollution to drift away from the city and consequently leading to frequent cases of heavy smog.
In some cases, we found that the adverse health effect is smaller in a highly polluted area when the air pollution level goes up dramatically in a short period of time—this is consistent with the results from previous studies [24,27,44,45]. We found that only the effect of an increase of 10 μg/m3 of NO2 (lag 0) on CVD hospital admissions was significant in the coastal city of Weihai. Comparing the effects of NO2 in the two cities, we noticed that the effect in Weihai was much stronger than that in Jinan, being 6.568% (3.636, 9.584) in Weihai and 2.647% (1.607, 3.697) in Jinan. However, Figure 2 shows that the median NO2 concentration in Weihai during the study period was 16, and for Jinan was 44. This might be because in cities with higher air pollution such as Jinan, the public are less sensitive to air pollution, as their body functions, especially the underlying biochemical processes, may have become adjusted to exposure to a high level of the toxic component. Moreover, the estimates for NO2 remained significant after adjustment for other pollutants, suggesting that NO2 is important for the air pollution mixture in these two cities.
The subgroup analysis also shows the spatial heterogeneity. For example, the adverse effects of air pollution on the hospitalizations of younger patients were stronger than the aged patients in Jinan, but the results were opposite in Weihai. This could be caused by aged people actively adopting defensive measures when living in heavily polluted cities, such as purchasing air filters, reducing outdoor activities, or wearing masks. On the contrary, young people have more outdoor activities and are more likely to be exposed to air pollution, especially if they are in the construction or manufactural industry. In contrast, people that live in less polluted regions, such as Weihai, may have a very weak sense of air pollution self-protection and precautions, which could expalin why the air pollution effect on aged patients was stronger.
This study has a limitation. We excluded two monitoring sites from the study in the city of Jinan, which were both newly established in May, 2015. Since our study period was from 2014 to 2016, the missing rates of these two sites were about 44%. Given that these missing values appeared sequentially since the beginning of the study period, the accuracy of the predicted values using the Bayesian method may be an issue. In our study, we firstly integrated the spatial and temporal information to obtain the underlying air pollution levels, then plugged such values into the health models to estimate the health effects. Undoubtedly, the error introduced through predictions in the first stage can lead to bias in the estimate of the health effects; therefore, we need to apply a more efficient modeling method to predict this type of missing data in future studies.

5. Conclusions

Our analysis showed that air pollution was significantly associated with CVD hospital admissions in the inland city, but the short-term increases in air pollution levels, especially NO2, in the less polluted coastal city can also lead to the excessive exacerbation of CVD. Therefore, we suggest that not only the public living in heavily polluted cities need defensive activities but that the people inhabiting less polluted cities should also adopt self-protection measures when there is an obvious incline in the air pollution in a short time period.

Author Contributions

Data curation, J.S. and Y.G.; Methodology, X.S.; Supervision, F.X.; Visualization, D.Z.; Writing—original draft, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of Shandong Province: ZR2016HB49 and Fundamental Research Fund of Shandong University: 2016HW001.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Cardiovascular Diseases. Available online: http://www.who.int/cardiovascular_diseases/en/ (accessed on 13 January 2020).
  2. Lozano, R.; Naghavi, M.; Foreman, K.; Lim, S.; Shibuya, K.; Aboyans, V.; Abraham, J.; Adair, T.; Aggarwal, R.; Ahn, S.Y.; et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: A systematic analysis for the Global Burden of Disease Study 2010. Lancet 2012, 380, 2095–2128. [Google Scholar] [CrossRef]
  3. 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 2012, 380, 2224–2260. [Google Scholar] [CrossRef] [Green Version]
  4. Liu, J.; Liu, Y.; Wang, L.; Yin, P.; Liu, S.; You, J.; Zeng, X.; Zhou, M. The disease burden of cardiovascular and circulatory diseases in China, 1990 and 2010. Chin. J. Prev. Med. 2015, 49, 315–320. [Google Scholar]
  5. Mathers, C. The Global Burden of Disease: 2004 Update; World Health Organization: Geneva, Switzerland, 2008. [Google Scholar]
  6. Xie, W.; Li, G.; Zhao, D.; Xie, X.; Wei, Z.; Wang, W.; Wang, M.; Li, G.; Liu, W.; Sun, J.; et al. Relationship between fine particulate air pollution and ischaemic heart disease morbidity and mortality. Heart 2015, 101, 257–263. [Google Scholar] [CrossRef] [PubMed]
  7. Brook, R.D.; Rajagopalan, S.; Pope, C.A., III; 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: An update to the scientific statement from the American Heart Association. Circulation 2010, 121, 2331–2378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Franchini, M.; Mannucci, P.M. Air pollution and cardiovascular disease. Thrombosis Res. 2012, 129, 230–234. [Google Scholar] [CrossRef]
  9. 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] [Green Version]
  10. Sun, Q.; Hong, X.; Wold, L.E. Cardiovascular effects of ambient particulate air pollution exposure. Circulation 2010, 121, 2755–2765. [Google Scholar] [CrossRef]
  11. 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. J. Am. Med. Assoc. 2006, 295, 1127–1134. [Google Scholar] [CrossRef] [Green Version]
  12. Pope, C.A., III; Burnett, R.T.; Turner, M.C.; Cohen, A.; Krewski, D.; Jerrett, M.; Gapstur, S.M.; Thun, M.J. Lung cancer and cardiovascular disease mortality associated with ambient air pollution and cigarette smoke: Shape of the exposure–response relationships. Environ. Health Perspect. 2011, 119, 1616–1621. [Google Scholar] [CrossRef] [Green Version]
  13. Peters, A.; Liu, E.; Verrier, R.L.; Schwartz, J.; Gold, D.R.; Mittleman, M.; Baliff, J.; Oh, J.A.; Allen, G.; Monahan, K.; et al. Air pollution and incidence of cardiac arrhythmia. Epidemiology 2000, 11, 11–17. [Google Scholar] [CrossRef] [PubMed]
  14. Chan, C.K.; Yao, X. Air pollution in mega cities in China. Atmos. Environ. 2008, 42, 1–42. [Google Scholar] [CrossRef]
  15. Rohde, R.A.; Muller, R.A. Air pollution in China: Mapping of concentrations and sources. PLoS ONE 2015, 10, e0135749. [Google Scholar] [CrossRef] [PubMed]
  16. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Kan, H.; Chen, R.; Tong, S. Ambient air pollution, climate change, and population health in China. Environ. Int. 2012, 42, 10–19. [Google Scholar] [CrossRef] [PubMed]
  18. Matus, K.; Nam, K.M.; Selin, N.E.; Lamsal, L.N.; Reilly, J.M.; Paltsev, S. Health damages from air pollution in China. Glob. Environ.Chang. 2012, 22, 55–66. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, M.; Huang, Y.; Ma, Z.; Jin, Z.; Liu, X.; Wang, H.; Liu, Y.; Wang, J.; Jantunen, M.; Bi, J.; et al. Spatial and temporal trends in the mortality burden of air pollution in China: 2004–2012. Environ. Int. 2017, 98, 75–81. [Google Scholar] [CrossRef] [Green Version]
  20. Guo, Y.; Jia, Y.; Pan, X.; Liu, L.; Wichmann, H.E. The association between fine particulate air pollution and hospital emergency room visits for cardiovascular diseases in Beijing, China. Sci. Total Environ. 2009, 407, 4826–4830. [Google Scholar] [CrossRef]
  21. Kan, H.; Huang, W.; Chen, B.; Zhao, N. Impact of outdoor air pollution on cardiovascular health in Mainland China. CVD Prev. Control. 2009, 4, 71–78. [Google Scholar] [CrossRef]
  22. Zhang, P.; Dong, G.; Sun, B.; Zhang, L.; Chen, X.; Ma, N.; Yu, F.; Guo, H.; Huang, H.; Lee, Y.L.; et al. Long-term exposure to ambient air pollution and mortality due to cardiovascular disease and cerebrovascular disease in Shenyang, China. PLoS ONE 2011, 6, e20827. [Google Scholar] [CrossRef]
  23. Wong, T.W.; Tam, W.S.; Yu, T.S.; Wong, A.H. Associations between daily mortalities from respiratory and cardiovascular diseases and air pollution in Hong Kong, China. Occ. Environ. Med. 2002, 59, 30–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Li, H.; Chen, R.; Meng, X.; Zhao, Z.; Cai, J.; Wang, C.; Yang, C.; Kan, H. Short-term exposure to ambient air pollution and coronary heart disease mortality in 8 Chinese cities. Int. J. Cardiol. 2015, 197, 265–270. [Google Scholar] [CrossRef] [PubMed]
  25. Ye, X.; Peng, L.; Kan, H.; Wang, W.; Geng, F.; Mu, Z.; Zhou, J.; Yang, D. Acute effects of particulate air pollution on the incidence of coronary heart disease in Shanghai, China. PLoS ONE 2016, 11, e0151119. [Google Scholar] [CrossRef] [PubMed]
  26. Xu, M.; Guo, Y.; Zhang, Y.; Westerdahl, D.; Mo, Y.; Liang, F.; Pan, X. Spatiotemporal analysis of particulate air pollution and ischemic heart disease mortality in Beijing, China. Environ. Health 2014, 13, 109. [Google Scholar] [CrossRef] [Green Version]
  27. Yin, P.; He, G.; Fan, M.; Chiu, K.Y.; Fan, M.; Liu, C.; Xue, A.; Liu, T.; Pan, Y.; Mu, Q.; et al. Particulate air pollution and mortality in 38 of China’s largest cities: Time series analysis. Br. Med. J. 2017, 356, j667. [Google Scholar] [CrossRef] [Green Version]
  28. Lee, D.; Shaddick, G. Spatial modeling of air pollution in studies of its short-term health effects. Biometrics 2010, 66, 1238–1246. [Google Scholar] [CrossRef]
  29. Yan, H. Heritage Tourism in China: Modernity, Identity and Sustainability; Channel View Publications: Bristol, UK, 2017. [Google Scholar]
  30. Diggle, P.J.; Tawn, J.A.; Moyeed, R.A. Model-based geostatistics. J. R. Stat. Soc. Ser. C App. Stat. 1998, 47, 299–350. [Google Scholar] [CrossRef]
  31. Gelman, A. Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper). Bayesian Anal. 2006, 1, 515–534. [Google Scholar] [CrossRef]
  32. Hastie, T.J. Generalized additive models. In Statistical Models in S; Routledge: London, UK, 2017; pp. 249–307. [Google Scholar]
  33. Wood, S.N. Generalized Additive Models: An Introduction with R; CRC Press: London, UK, 2017. [Google Scholar]
  34. Akaike, H. Factor analysis and AIC. In Selected Papers of Hirotugu Akaike; Springer: New York, NY, USA, 1987; pp. 371–386. [Google Scholar]
  35. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018; Available online: https://www.R-project.org/ (accessed on 13 January 2020).
  36. Ministry of Environmental Protection of the People’s Republic of China. China’s Ambient air Quality Secondary Standards. Available online: http://kjs.mep.gov.cn/hjbhbz/bzwb/dqhjbh/dqhjzlbz/201203/t20120302_224165.shtml (accessed on 13 January 2020).
  37. Cohen, A.J.; Ross Anderson, H.; Ostro, B.; Pandey, K.D.; Krzyzanowski, M.; Künzli, N.; Gutschmidt, K.; Pope, A.; Romieu, I.; Samet, J.M.; et al. The global burden of disease due to outdoor air pollution. J. Toxicol. Environ. Health Part A 2005, 68, 1301–1307. [Google Scholar] [CrossRef]
  38. 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]
  39. Pope, C.A., III; Dockery, D.W. Health effects of fine particulate air pollution: Lines that connect. J. Air Waste Manag. Assoc. 2006, 56, 709–742. [Google Scholar] [CrossRef] [PubMed]
  40. Carlin, B.P.; Xia, H.; Devine, O.; Tolbert, P.; Mulholland, J. Spatio-temporal hierarchical models for analyzing Atlanta pediatric asthma ER visit rates. In Case Studies in Bayesian Statistics; Springer: New York, NY, USA, 1999; pp. 303–320. [Google Scholar]
  41. Zhu, L.; Carlin, B.P.; Gelfand, A.E. Hierarchical regression with misaligned spatial data: Relating ambient ozone and pediatric asthma ER visits in Atlanta. Environ. Off. J. Int. Environ. Soc. 2003, 14, 537–557. [Google Scholar] [CrossRef]
  42. Chang, H.H.; Peng, R.D.; Dominici, F. Estimating the acute health effects of coarse particulate matter accounting for exposure measurement error. Biostatistics 2011, 12, 637–652. [Google Scholar] [CrossRef] [PubMed]
  43. Peng, R.D.; Bell, M.L. Spatial misalignment in time series studies of air pollution and health data. Biostatistics 2010, 11, 720–740. [Google Scholar] [CrossRef]
  44. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Chen, R.; Kan, H.; Chen, B.; Huang, W.; Bai, Z.; Song, G.; Pan, G. 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]
Figure 1. The location of Jinan and Weihai in the Shandong Province; the two cities are shown in blue. The numbers represent the number of air quality monitoring stations in each city.
Figure 1. The location of Jinan and Weihai in the Shandong Province; the two cities are shown in blue. The numbers represent the number of air quality monitoring stations in each city.
Ijerph 17 00879 g001
Figure 2. A boxplot of daily mean air pollution levels in Jinan and Weihai.
Figure 2. A boxplot of daily mean air pollution levels in Jinan and Weihai.
Ijerph 17 00879 g002
Figure 3. A boxplot of meteorological factors in Jinan and Weihai.
Figure 3. A boxplot of meteorological factors in Jinan and Weihai.
Ijerph 17 00879 g003
Figure 4. A plot of the percentage change in the relative risk of the total cardiovascular diseases hospitalizations with a 95% confidence interval (CI) per 10 μg/m3 increase in the air pollution levels, in the cities of Jinan and Weihai.
Figure 4. A plot of the percentage change in the relative risk of the total cardiovascular diseases hospitalizations with a 95% confidence interval (CI) per 10 μg/m3 increase in the air pollution levels, in the cities of Jinan and Weihai.
Ijerph 17 00879 g004
Table 1. The daily hospital admissions grouped by sex, age, and sub-disease in the cities of Jinan and Weihai (from 1st January, 2014 to 31st December, 2016).
Table 1. The daily hospital admissions grouped by sex, age, and sub-disease in the cities of Jinan and Weihai (from 1st January, 2014 to 31st December, 2016).
CityHospitalization (Cases/Per Day)MinP25P50P75Max
Total8195291346580
Male5106158191308
JinanFemale389131157274
≥653103148180361
<65592138169296
Total1678096200
Male0374554113
WeihaiFemale129364389
≥650374554109
<651283544104
Table 2. The percentage change in the relative risk of the total cardiovascular disease hospitalizations with a 95% confidence interval (CI) per 10 μg/m3 increase in the air pollution levels estimated from multipollutant models, where * indicates statistically significant estimates (p < 0.05).
Table 2. The percentage change in the relative risk of the total cardiovascular disease hospitalizations with a 95% confidence interval (CI) per 10 μg/m3 increase in the air pollution levels estimated from multipollutant models, where * indicates statistically significant estimates (p < 0.05).
Jinan
PM2.5PM10SO2NO2
Single pollutant0.401 (0.029, 0.775) *Single pollutant0.316 (0.086, 0.547) *Single pollutant0.903 (0.252, 1.559) *Single pollutant2.647 (1.607, 3.697) *
+PM10−0.279 (−1.106, 0.555)+PM2.50.472 (−0.045, 0.992)+PM2.50.889 (0.156, 1.628) *+ PM2.53.164 (1.835, 4.509) *
+SO20.360 (−0.057, 0.778)+SO20.307 (0.050, 0.565) *+PM100.859 (0.115, 1.609) *+ PM103.018 (1.631, 4.423) *
+NO2−0.295 (−0.764, 0.177)+NO2−0.124 (−0.429, 0.183)+NO20.298 (−0.584, 1.188)+ SO23.879 (2.483, 5.295) *
+All−0.342 (−1.181, 0.504)+All0.071 (−0.475, 0.620)+All0.316 (−0.568, 1.208)+All4.229 (2.564, 5.922) *
Weihai
PM2.5PM10SO2NO2
Single pollutant 0.161 (−0.356, 0.680)Single pollutant0.143 (−0.210, 0.497)Single pollutant2.149 (−0.031, 4.377)Single pollutant6.568 (3.636, 9.584) *
+PM10−0.001 (−0.808, 0.812)+PM2.50.144 (−0.409, 0.699)+PM2.53.580 (0.472, 6.785) *+PM2.59.606 (5.887, 13.456) *
+SO2−0.239 (−0.949, 0.477)+SO2−0.021 (−0.447, 0.407)+ PM102.875 (0.264, 5.554) *+PM108.453 (5.048, 11.970) *
+NO2−0.536 (−1.204, 0.136)+NO2−0.167 (−0.582, 0.249)+NO2−1.980 (−4.767, 0.887) *+SO28.415 (4.448, 12.532) *
+All−0.654 (−1.597, 0.298)+All0.098 (−0.457, 0.656)+All0.035 (−3.277, 3.460)+All9.648 (5.518, 13.940) *
Table 3. The changes in the relative risk of cardiovascular disease hospitalizations for patients by gender and age, and their 95% confidence interval (CI) per 10 μg/m3 increase in air pollutants in Jinan and Weihai, where * indicates statistically significant estimates (p < 0.05).
Table 3. The changes in the relative risk of cardiovascular disease hospitalizations for patients by gender and age, and their 95% confidence interval (CI) per 10 μg/m3 increase in air pollutants in Jinan and Weihai, where * indicates statistically significant estimates (p < 0.05).
JinanWeihai
PollutantClassEstimateLagEstimateLag
PM2.5Male0.404 (0.009, 0.801) *0−0.419 (−0.992, 0.157)2
Female0.398 (0.013, 0.783) *0−0.420 (−1.010, 0.174)3
≥650.323 (−0.056, 0.704)00.660 (0.109, 1.214) *5
<650.485 (0.074, 0.899) *0−0.577 (−1.198, 0.048)3
PM10Male0.308 (0.063, 0.553) *00.119 (−0.260, 0.499)3
Female0.326 (0.088, 0.565) *00.283 (−0.138, 0.706)1
≥650.279 (0.045, 0.514) *00.306 (−0.068, 0.682)5
<650.356 (0.100, 0.612) *0−0.257 (−0.684, 0.172)4
SO2Male0.749 (0.058, 1.445) *12.859 (0.340, 5.441) *0
Female1.089 (0.415, 1.768) *11.281 (−1.292, 3.921)0
≥650.709 (0.041, 1.381) *12.438 (0.344, 4.576) *5
<651.116 (0.398, 1.839) *11.924 (−0.807, 4.730)0
NO2Male2.516 (1.414, 3.630) *07.419 (4.031, 10.918) *0
Female2.803 (1.725, 3.892) *05.535 (2.070, 9.119) *0
≥652.284 (1.225, 3.355) *07.612 (4.321, 11.006) *0
<653.033 (1.880, 4.199) *05.236 (1.591, 9.011) *0

Share and Cite

MDPI and ACS Style

Liu, Y.; Sun, J.; Gou, Y.; Sun, X.; Zhang, D.; Xue, F. Analysis of Short-Term Effects of Air Pollution on Cardiovascular Disease Using Bayesian Spatio-Temporal Models. Int. J. Environ. Res. Public Health 2020, 17, 879. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17030879

AMA Style

Liu Y, Sun J, Gou Y, Sun X, Zhang D, Xue F. Analysis of Short-Term Effects of Air Pollution on Cardiovascular Disease Using Bayesian Spatio-Temporal Models. International Journal of Environmental Research and Public Health. 2020; 17(3):879. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17030879

Chicago/Turabian Style

Liu, Yi, Jingjie Sun, Yannong Gou, Xiubin Sun, Dandan Zhang, and Fuzhong Xue. 2020. "Analysis of Short-Term Effects of Air Pollution on Cardiovascular Disease Using Bayesian Spatio-Temporal Models" International Journal of Environmental Research and Public Health 17, no. 3: 879. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17030879

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