Association of short-term exposure to sulfur dioxide and hospitalization for ischemic and hemorrhagic stroke in Guangzhou, China
BMC Public Health volume 20, Article number: 263 (2020)
In developing countries, ambient sulfur dioxide (SO2) is a serious air pollutant concern, but there is no enough and consistent epidemiological evidence about its health effects on stroke hospitalization.
We collected the daily air pollution data, meteorological data and number of daily hospital admissions for ischemic and hemorrhagic stroke, in Guangzhou from January 1st 2009 to December 31st 2014. Then we applied generalized additive model with a quasi-Poisson link to assess the relationship between short-term SO2 exposure and the total number of hospital admissions for ischemic and hemorrhagic stroke. In addition, we evaluated the effect of ambient SO2 by age (< 65 years and ≥ 65 years).
During the study period, a 24-h mean concentration of ambient SO2 of 27.82 μg/m3, a total of 58,473 ischemic stroke and 9167 hemorrhagic stroke hospital admissions hospital were recorded. Ambient SO2 was found to increase the risk for both ischemic and hemorrhagic stroke hospital admission in single pollutant model. The maximum value of percentage changes for ischemic and hemorrhagic stroke occurred in lag 0 day and lag 1 day, per 10 μg/m3 increase in SO2 concentrations was corresponded to a 1.27% (95% confidence interval (CI), 0.42–2.12%) and 1.55% (95%CI, 0.02–3.11%) increased risk, respectively. The association between SO2 and ischemic stroke hospitalization was robust to two pollutant model, but for hemorrhagic stroke it’s partially weakened after adjusting for co-pollutants. The effect of ambient SO2 on ischemic stroke appeared to be greater for people < 65 years old, but null effect on hemorrhagic stroke was identified for both age groups.
We found short-term exposure to ambient SO2 may significantly increase the risks of hospitalization for ischemic stroke. The findings may contribute to a better understanding of the health effects of low-levels of SO2.
Stroke is the second leading cause of death in the world . According to Global Burden of Disease (GBD) Study 2016, the highest age-standardized incidences of stroke were observed in east Asia, especially China (354 [95% UI 331–378] per 100,000 person-years) . Ischemic stroke and hemorrhagic stroke account for approximately 70 and 28% of all stroke cases, respectively . Moreover, the incidence of stroke continues to rise in China in recent years . Therefore, identification of modifiable risk factors for stroke has significant public health implications.
Hypertension, diabetes mellitus and dyslipidemia are major risk factors for cardiovascular diseases including stroke [4, 5]. The high prevalence of these risk factors is a major burden to most economies around the world. Ambient air pollution is also a well-known risk factor for stroke morbidity , but the majority of evidence are focused on particulate matter (PM) [7, 8]. Sulfur dioxide (SO2) is a non-flammable, colorless gas with a strong pungent odor at room temperature, which is primarily released from the combustion of sulfur containing fossil fuels at power plants (73%) and other industry facilities (20%) . As a strong respiratory irritant and bronchoconstrictor, SO2 could induce systemic inflammation and oxidative stress , which might take a role in stroke. Short-term exposure to ambient SO2 was associated with increased risk of stroke morbidity in a meta-analysis mainly based on developed countries . However, in developing countries like China, the evidence is limited and very few studies  distinguished between ischemic and hemorrhagic stroke. Most studies only focused on ischemic stroke but the results have not been conclusive yet, with some finding evidence of a positive association  but not in other studies [13,14,15,16].
Guangzhou is one of the most densely populated cities in China. The industrial production and traffic burden make SO2 one of the main atmospheric pollution sources in Guangzhou , which is severer than those in developed countries . Moreover, the stroke mortality in Guangzhou was relatively higher compared to other cities of China .
Therefore, we conducted a time-series study to examine the short-term association between ambient SO2 and hospital admissions of ischemic and hemorrhage stroke in Guangzhou, respectively.
Guangzhou is located in the South China, which has a typical subtropical humid-monsoon climate with an average annual temperature of 22 °C and an average rainfall of 1500–2000 mm . There are 12.7 million people lived in Guangzhou in 2010 . The city is one of the top-ranking commercial and manufacturing regions in China. As the economy develops in Guangzhou, it also introduces severe air pollution problems which might lead to risks to human health .
We obtained the air quality monitoring data from the public sharing system of Guangzhou Environmental Monitoring Center (http://126.96.36.199:8023/gzaqi_new/RealTimeDate.html), which includes data of all the 11 national air quality monitor stations in Guangzhou. The data of SO2, nitrogen dioxide (NO2) and suspended particulates smaller than 10 μm in aerodynamic diameter (PM10) covered the period from January 1, 2009 to December 31, 2014. But for suspended particulates smaller than 2.5 μm in aerodynamic diameter (PM2.5) and ozone (O3), the data was limited to the period of March 8, 2012 to December 31, 2014, because Guangzhou Environmental Monitoring Center started the surveillance of PM2.5 and O3 from March 8, 2012. All data from the monitoring stations were averaged in the present study. We used daily average concentration for SO2, NO2, PM10, PM2.5, and daily maximum 8-h average (O3-8h max) for O3.
Daily meteorological data (daily mean temperature and relative humidity) from January 1, 2009 to December 31, 2014 were obtained from the China Meteorological Data Sharing Service System. All meteorological stations were averaged in the analysis.
The data on air pollution and meteorology follow the quality control programs mandated by the Chinese government. These types of data have been widely used to explore their effects on morbidity and mortality in Guangzhou [22,23,24].
Daily hospital admissions data were obtained from eight hospitals in Guangzhou during January 1, 2009 to December 31, 2014. If there were recurrent strokes occurring within 28 days of a previous admission, then it was treated as a single stroke event . All the cases were coded under the World Health Organization’s International Classification of Diseases, the 10th version (ICD-10). Ischemic and hemorrhagic stroke hospital admissions were identified as I63 and I61, respectively.
This study does not involve experimental animals or individual information on human subjects. Ethical approval was not required for this study. Research projects that do not involve human participants, their data or tissue do not require ethical review by the ethics committee in Southern Medical University.
To assess the relationship between short-term exposure of SO2 and stroke hospital admissions, we applied a generalized additive model (GAM) with Poisson link in this study . In accordance with previous studies, we used 7 degrees of freedom (df) per year for the smooth function of calendar time to control for fluctuations in hospitalization on long time trend and seasonality , 6 df for the current day’s temperature, and 3 df for relative humidity . Day of the week (DOW) was also included as a dummy variable.
We first used the single pollutant model to examine the independent association between SO2 and ischemic/hemorrhagic stroke risk. Then two pollutant model including another pollutant was built to test if there are any potential cofounding effects of other pollutant (e.g. NO2) on the relationship between SO2 and stroke.
These models can be generically represented as:
Single pollutant model:
Log [E (Yt)] = α + βiXi + s (temperature, df = 6) + s (humidity, df = 3) + DOW + s (time, df = 7/year).
Two pollutants model:
Log [E (Yt)] = α + ∑βiXi + s (temperature, df = 6) + s (humidity, df = 3) + DOW + s (time, df = 7/year).
where E (Yt) stands for the expected number of stroke hospital admissions on day t. α represent the intercept. X is the concentrations of pollutant (SO2, PM10, PM2.5, NO2 and O3). βi stands for the coefficient for Xi. s() represented a smoother based on penalized smoothing splines, which captures the nonlinear relationships of the covariates of time trend, temperature and humidity with daily hospitalization. df represents the degree of freedom. Daily mean temperature and relative humidity were used in all models to control for confounding. Time is the time to adjust for long-term trend and seasonality.
To explore the potential modification, we examined outcomes by stratification of age (≥65 years and < 65 years). We also conducted sensitivity analyses to check the robustness of results by changing the degrees of freedom in the smooth function of time trend (7–9 per year). P value < 0.05 was considered as statistically significant. All the analyses were performed with R software (R Development Core Team, 2015).
Table 1 presented the stroke, meteorological and air pollution data in Guangzhou during January 1, 2009 to December 31, 2014. Overall, 58,473 ischemic stroke and 9167 hemorrhagic stroke hospital admissions were recorded, the distribution of age is presented in Supplementary Table 1. The daily average temperature was 21.99 °C while daily relative humidity was 76.30%. The average concentrations of air pollutants including SO2, PM10, PM2.5, NO2 and O3 were 27.82 μg/m3, 72.93 μg/m3, 46.07 μg/m3, 54.33 μg/m3 and 68.69 μg/m3, respectively.
The correlation coefficients between SO2 and other pollutants varied from 0.31 to 0.48. And SO2 is negatively correlated with humidity, which is opposite to temperature. (Supplementary Table 2).
For ischemic stroke, we observed the maximum effect of SO2 at lag0 day (Table 2). Every 10 μg/m3 increase of SO2 was associated with 1.27% (95%CI, 0.42–2.12%) increase of ischemic stroke (Table 2). In all the two-pollutant models, the effects of SO2 remained in line with that of single pollutant model (Table 4). When looking into the different age groups, a 10 μg/m3 increase of SO2 concentrations on lag 0 day represented 1.39% (95%CI, 0.35–2.45%) and 1.21% (95%CI, 0.02–2.42%) increment for ischemic stroke hospital admissions in < 65 and ≥ 65 years old group, respectively (Table 3). And we found SO2 still increased the risk of hemorrhagic stroke hospital admission with the adjustment of NO2 and O3 for < 65 years (Table 4). But null association existed in two pollutant model for ≥65 years (Table 4).
For hemorrhagic stroke, the effect of SO2 reached the maximum value on lag1 day (Table 2), with percentage change of 1.55% (95%CI, 0.02–3.11%) per 10 μg/m3. But neither of the effect is significant in < 65 and ≥ 65 years old group (Table 3). The effect of SO2 lost statistical significance after including NO2, PM10 and PM2.5 but it’s still significant when O3 was included (Table 4).
Sensitivity analyses were performed to check the robustness of results in terms of the degrees of freedom in the smooth function of time trend, which was similar to that in the main model (Supplementary Table 3).
This is one of the few studies specified not only association between SO2 and ischemic stroke but also hemorrhagic stroke hospital admissions in developing countries so far.
Our study suggested that ambient SO2 was significantly associated with increased hospital admissions for ischemic stroke. A 10 μg/m3 increase of the same day/lag 1 day SO2 was associated with increments of 1.27% (0.42–2.12%) and 0.93% (0.09–1.78%) ischemic stroke hospital admissions risk, respectively. The associations were robust to the adjustment of all the other air pollutants. For example, the most informative pollutant of traffic emissions NO2. The association were comparable to previously studies [6, 11, 12, 16].
However, the associations of SO2 and hemorrhagic stroke were more variable and imprecise. We observed significant associations between lag 1/lag 2 days SO2 exposures and hemorrhagic stroke hospitalization risk, with percentage change of 1.55% (0.02–3.11%) and 1.52% (0.02–3.04%), respectively. But when adjusted for other air pollutants in our models (except O3), the magnitude of risk due to SO2 exposure turned out to be statistically insignificant. Previous evidences on SO2-hemorrhagic stroke association are scarce and inconsistent [11, 29]. The discrepant results on hemorrhagic stroke may be explained by differences in average pollutant levels, meteorological patterns and population characteristics between studies . For example, in the study of Liu et al. , SO2 concentration level was 39.6 ± 41.2 μg/m3 with remarkably higher variation than that in Guangzhou (21.99 ± 6.37 μg/m3) in our study. The exposure levels and exposure patterns might partly explain the difference.
The identification of potentially vulnerable subpopulation has significant implications for public health , therefore we stratified the stroke cases into < 65 and ≥ 65 years groups. For ischemic stroke, we found that the effect estimates of SO2 were higher for the < 65 than that in the ≥65 age group, which is the opposite of Tian et al. . This may be due to the varied activity patterns for different age groups in different cities. The elderly in Guangzhou tend to stay in door or in the park nearby, while people under 65 years old usually take vehicles or public transportation to work and have more exposure on SO2, because the vehicle exhaust emission is one of the main sources of SO2. But for hemorrhagic stroke hospitalization, there was no evidence of effect in neither age subgroup in any lag structure. It might be related to the smaller number of daily hospital admissions of hemorrhagic stroke (4 ± 3 cases/day) in this study. The lower incidence of hemorrhagic stroke may lead to larger imprecision in the estimates . In addition, we noted that 70.14% hemorrhagic stroke cases were < 65 years old, while 59.75% ischemic stroke cases were ≥ 65 years old, which was similar to previous observations . But the different age distribution of ischemic and hemorrhagic stroke should not be ignored. It might also contribute to the different relationship between SO2 and ischemic and hemorrhagic stroke.
The mechanism behind the adverse stroke effects associated with SO2 is biological plausible. The inhalation of SO2 could alter heart rate variability, elevate oxidative and enhanced blood coagulation and thrombosis formation . Animal study also showed that SO2 elevated the levels of cyclooxygenase-2, interleukin-1b, tumor necrosis factor-a, intercellular adhesion molecule-1 mRNA and protein, which might also play a role in the development of stroke . These pathophysiologic changes may be related to the development and progression of stroke. Nevertheless, the exact biological pathways underlying the health impact of SO2 need further investigations.
Although this study found that the magnitude of elevated risk of stroke due to SO2 exposures in Guangzhou was relatively small, it should not be ignored. Coal related air pollution, characterized by high level of SO2, remains to be an important environmental problem in China, including Guangzhou . The increased risk of ischemic stroke and hemorrhagic stroke due to the exposures of SO2 estimated in this work highlights the need for continued vehicle emission control, desulfurization, and promotion of hydropower and nuclear power from the government  and vigilance for the health risks of air pollution. There are several limitations of this study. First, averaging air pollution exposures over Guangzhou may introduce exposure measurement error. Secondly, as a retrospective ecological study, personal exposure confounders could not be excluded in the model, and we have no access to individual air pollution exposure data. Further studies based on personal measurements or modeled residential concentrations are warranted to evaluate the health effects of SO2, especially at low levels or in the real world. Moreover, the potential misclassifications for the stroke diagnosis should also be taken into consideration. However, all the hospitals included in this study are top-ranked public hospitals with high quality of electronic diseases diagnosis records in Guangzhou. Another limitation of the present study is the lack of consideration of more age groups, as the distribution of ischemic and hemorrhagic stroke might vary across the age (e.g. 65–74, and ≥ 75 years), we recommend future studies to examine the association in subgroups as 65–74, and ≥ 75 years. Due to the limited number of hemorrhagic stroke hospitalizations, studies based on enlarged sample size or in multi-cites should be performed to examine the association in the future. In addition, hypertension, diabetes mellitus and dyslipidemia are important risk factors for stroke . Further studies taking the comorbidity of stroke and hypertension, diabetes mellitus and dyslipidemia into consideration could be helpful to identify effective ways to prevent stroke.
In summary, we found short-term exposure to ambient SO2 may significantly increase the risks of hospitalization for ischemic stroke. The findings may contribute to a better understanding of the health effects of SO2.
Availability of supporting data
The air pollution data that support the findings of this study are available on request from the corresponding author, Qing Chen. The stroke data are not publicly available due to the privacy restriction. Competing interests.
The authors declare they have no competing interests.
Day of the week
Generalized additive model
Global Burden of Disease
International Classification of Diseases, the 10th version
- O3 :
- SO2 :
Collaborators GBDLRoS, Feigin VL, Nguyen G, Cercy K, Johnson CO, Alam T, Parmar PG, Abajobir AA, Abate KH, Abd-Allah F, et al. global, regional, and country-specific lifetime risks of stroke, 1990 and 2016. N Engl J Med. 2018;379(25):2429–37.
Collaborators GBDS. Global, regional, and national burden of stroke, 1990-2016: a systematic analysis for the global burden of disease study 2016. Lancet Neurol. 2019;18(5):439–58.
Wu S, Wu B, Liu M, Chen Z, Wang W, Anderson CS, Sandercock P, Wang Y, Huang Y, Cui L, et al. Stroke in China: advances and challenges in epidemiology, prevention, and management. Lancet Neurol. 2019;18(4):394–405.
Beaney TBL, Castillo RR, Charchar FJ, Cro S, Damasceno A, Kruger R, Nilsson PM, Prabhakaran D, Ramirez AJ, Schlaich MP, Schutte AE, Tomaszewski M, Touyz R, Wang JG, Weber MA, Poulter NR. May Measurement Month 2018: A pragmatic global screening campaign to raise awareness of blood pressure by the International Society of Hypertension. Eur Heart J. 2019(40):2006–17.
Chia YCCS, Chew BN, Devaraj NK, Siti Suhaila MY, Tay CL, Kang PS, Verna Lee KM, Kong SZ, Teoh SW, Nurjasmine AJ. May measurement month 2017 blood pressure screening: findings from Malaysia—South-East Asia and Australasia. Eur Heart J Suppl. 2019;21(April):D77–9.
Shah AS, Lee KK, McAllister DA, Hunter A, Nair H, Whiteley W, Langrish JP, Newby DE, Mills NL. Short term exposure to air pollution and stroke: systematic review and meta-analysis. BMJ. 2015;350:h1295.
Yitshak Sade M, Novack V, Ifergane G, Horev A, Kloog I. Air pollution and ischemic stroke among young adults. Stroke. 2015;46(12):3348–53.
Zeng W, Zhang Y, Wang L, Wei Y, Lu R, Xia J, Chai B, Liang X. Ambient fine particulate pollution and daily morbidity of stroke in Chengdu, China. PloS One. 2018;13(11):e0206836.
Ge E, Fan M, Qiu H, Hu H, Tian L, Wang X, Xu G, Wei X. Ambient sulfur dioxide levels associated with reduced risk of initial outpatient visits for tuberculosis: a population based time series analysis. Environ Pollut. 2017;228:408–15.
Wang LJ, Liu C, Meng X, Niu Y, Lin ZJ, Liu YN, Liu JM, Qi JL, You JL, Tse LA, et al. Associations between short-term exposure to ambient sulfur dioxide and increased cause-specific mortality in 272 Chinese cities. Environ Int. 2018;117:33–9.
Liu H, Tian YH, Xu Y, Huang Z, Huang C, Hu YH, Zhang J. Association between ambient air pollution and hospitalization for ischemic and hemorrhagic stroke in China: a multicity case-crossover study. Environ Pollut. 2017;230:234–41.
Tian Y, Liu H, Zhao Z, Xiang X, Li M, Juan J, Song J, Cao Y, Wang X, Chen L, et al. Association between ambient air pollution and daily hospital admissions for ischemic stroke: a nationwide time-series analysis. PLoS Med. 2018;15(10):e1002668.
Dong H, Yu Y, Yao S, Lu Y, Chen Z, Li G, Yao Y, Yao X, Wang SL, Zhang Z. Acute effects of air pollution on ischaemic stroke onset and deaths: a time-series study in Changzhou, China. BMJ Open. 2018;8(7):e020425.
Maheswaran R, Pearson T, Beevers SD, Campbell MJ, Wolfe CD. Air pollution and subtypes, severity and vulnerability to ischemic stroke-a population based case-crossover study. PLoS One. 2016;11(6):e0158556.
Xiang H, Mertz KJ, Arena VC, Brink LL, Xu X, Bi Y, Talbott EO. Estimation of short-term effects of air pollution on stroke hospital admissions in Wuhan, China. PloS One. 2013;8(4):e61168.
Guo P, Wang Y, Feng W, Wu J, Fu C, Deng H, Huang J, Wang L, Zheng M, Liu H. Ambient Air Pollution and Risk for Ischemic Stroke: A Short-Term Exposure Assessment in South China. Int J Environ Res Public Health. 2017;14(9).
Zhou K, Ye YH, Liu Q, Liu AJ, Peng SL. Evaluation of ambient air quality in Guangzhou, China. J Environ Sci (China). 2007;19(4):432–7.
Yang J, Yin P, Zhou M, Ou CQ, Li M, Li J, Liu X, Gao J, Liu Y, Qin R, et al. The burden of stroke mortality attributable to cold and hot ambient temperatures: epidemiological evidence from China. Environ Int. 2016;92-93:232–8.
Lin H, Tao J, Du Y, Liu T, Qian Z, Tian L, Di Q, Rutherford S, Guo L, Zeng W, et al. Particle size and chemical constituents of ambient particulate pollution associated with cardiovascular mortality in Guangzhou, China. Environ Pollut. 2016;208(Pt B):758–66.
Yu ITS, Zhang Y, San Tam WW, Yan QH, Yj X, Xj X, Wu W, Ma WJ, Tian LW, Tse LA, et al. Effect of ambient air pollution on daily mortality rates in Guangzhou, China. Atmos Environ. 2012;46:528–35.
Huang D, Xu J, Zhang S. Valuing the health risks of particulate air pollution in the Pearl River Delta, China. Environ Sci Policy. 2012;15(1):38–47.
Liu T, Zhang YH, Xu YJ, Lin HL, Xu XJ, Luo Y, Xiao J, Zeng WL, Zhang WF, Chu C. The effects of dust–haze on mortality are modified by seasons and individual characteristics in Guangzhou, China. Environ Pollut. 2014;187:116–23.
Li X, Xiao J, Lin H, Liu T, Qian Z, Zeng W, Guo L, Ma W. The construction and validity analysis of AQHI based on mortality risk: A case study in Guangzhou, China. Environ Pollut. 2017;220(Pt A):487–94.
He JR, Liu Y, Xia XY, Ma WJ, Lin HL, Kan HD, Lu JH, Feng Q, Mo WJ, Wang P, et al. Ambient temperature and the risk of preterm birth in Guangzhou, China (2001-2011). Environ Health Perspect. 2016;124(7):1100–6.
Guo P, Zheng M, Wang Y, Feng W, Wu J, Deng C, Luo G, Wang L, Pan B, Liu H. Effects of ambient temperature on stroke hospital admissions: results from a time-series analysis of 104,432 strokes in Guangzhou, China. Sci Total Environ. 2017;580:307–15.
Peng RD, Dominici F, Louis TA. Model choice in time series studies of air pollution and mortality. J Royal Stat Soc. 2006;169(2):179–203.
Peng RD, Dominici F, Pastor-Barriuso R, Zeger SL, Samet JM. Seasonal analyses of air pollution and mortality in 100 US cities. Am J Epidemiol. 2005;161(6):585–94.
Wong CM, Vichit-Vadakan N, Kan H, Qian Z. Public health and air pollution in Asia (PAPA): a multicity study of short-term effects of air pollution on mortality. Environ Health Perspect. 2008;116(9):1195–202.
Wellenius GA, Schwartz J, Mittleman MA. Air pollution and hospital admissions for ischemic and hemorrhagic stroke among medicare beneficiaries. Stroke. 2005;36(12):2549–53.
Chen R, Yin P, Meng X, Wang L, Liu C, Niu Y, Lin Z, Liu Y, Liu J, Qi J, et al. Associations between ambient nitrogen dioxide and daily cause-specific mortality: evidence from 272 Chinese cities. Epidemiol. 2018;29(4):482–9.
Routledge HC, Manney S, Harrison RM, Ayres JG, Townend JN. Effect of inhaled Sulphur dioxide and carbon particles on heart rate variability and markers of inflammation and coagulation in human subjects. Heart. 2006;92(2):220–7.
Sang N, Yun Y, Li H, Hou L, Han M, Li G. SO2 inhalation contributes to the development and progression of ischemic stroke in the brain. Toxicol Sci. 2010;114(2):226–36.
AK-y L. Limits to vehicle emission control: a case of Guangzhou. Procedia Environ Sci. 2016;36:106–13.
Ethics approval and consent to participate
The protocol of study and accessing the hospital admission data was approved by the ethics committee of School of Public Health of Southern Medical University, and was performed in accordance with the ethical standards noted in the 1964 Declaration of Helsinki and its later amendments.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The distribution of ischemic and hemorrhagic stroke events by age. Supplementary Table 2. Correlations between air pollutants and meteorological factors. Supplementary Table 3. The percent change of stroke risk associated with 10 μg/m3 and per IQR increase of SO2 with different df (sensitivity analysis).
About this article
Cite this article
Shen, S., Li, X., Yuan, C. et al. Association of short-term exposure to sulfur dioxide and hospitalization for ischemic and hemorrhagic stroke in Guangzhou, China. BMC Public Health 20, 263 (2020). https://doi.org/10.1186/s12889-020-8354-0