1 Introduction

The Deaths caused by cerebrovascular diseases (CEVD) account for a large proportion of the population mortality []. Recently, the associations between ambient air pollution exposure and CEVD mortality have gained more attention. Some current evidences suggest that short-term exposure to ambient air pollution can increase daily CEVD mortality on the same day or several days after the exposure [, , ]. But some studies also indicated no associations between air pollution and CEVD mortality [, ]. These divergences have not been adequately reviewed, but spatial variations were reported to be a key factor []. Spatial variations may include environmental factors, social-economic factors and other domestic characteristics of the study population. During the past few years, the gender-based differences in the vulnerability to ambient air pollution have been given more importance. There are already some evidences implying that the vulnerability to ambient air pollution between the two genders might be considerable different, and women could be more susceptible under most circumstances [, ]. So far, the age-based differences in susceptibility to air pollution have been frequently reported [, , ]. However, few studies about the gender-based differences on deaths caused by CEVD have been published. Moreover, the current environment protection and healthcare policies also pay more attention on the infants, the elderly or people with specific diseases, but seldom take the gender-based vulnerability into consideration. To provide more evidences for sifting susceptible population, further exploration of gender-based vulnerability to ambient air pollution are urgently needed.

On the other hand, low social-economic status are reported to aggravate the vulnerability to air pollution []. Therefore, although some relevant studies have been published, those who were carried out under worse air quality and among lower social-economic status population (for example, in developing countries) could yield different results. Therefore, to further explore the effects of short-term ambient air pollution exposure on CEVD mortality, the relevant studies in developing countries with worse air quality are still necessary.

In this study, a time-series analysis was performed. We aimed to: a. explore the associations between short-term ambient PM2.5, PM10 and O3 and daily CEVD mortality; b. detect the potential differences of the gender-based vulnerability to ambient air pollution.

2 Materials and Methods

2.1 Study area and population

The population of Shenyang, which is the largest city in Northeastern China, were selected. According to the yearbook of Liaoning Province of 2017, the number of Shenyang’s residents was 7.3 million. As a city dedicating to heavy industry, Shenyang was once named as the top 10 heavy polluted cities around the world by the World Bank in 2001. Although the environment has been significantly improved in the past decades, the levels of some air pollutants still kept above the national ambient air quality standard of China (NAAQS). According to the Shenyang Environmental Bulletin, for example, there are 174, 158 and 117 polluted days in Shenyang during 2014, 2015 and 2016, respectively [].

2.2 Data Sources

The daily CEVD (ICD-10, code: I60-I69.9) mortality data during 2014–2017 were obtained from Center of Disease Control (CDC) of Shenyang. The access to the ambient air pollution monitoring data and the meteorological data of Shenyang were presented in detail in a previous study []. In brief, daily mean concentration of Particulate matter with diameter <2.5μm (PM2.5), Particulate matter with diameter <10μm (PM10), Ozone (O3) were retrieved from 11 state-controlled environmental air quality automatic monitoring stations through the website of Shenyang Bureau of Ecology and Environment. The overall mean concentration of the 11 stations for each type of pollutant were calculated and used as the general exposure level of the Shenyang population [, , , ]. During 2015 April 4th to 9th, the air quality monitoring data were missing. Therefore, a total of 1455 days’ air quality data (four years total days of 1461 minus 6 days) were available. The source of the meteorological information was from National Oceanic and Atmospheric Administration (NOAA), and was retrieved via the R package ‘worldmet’. Factors such as daily average temperature and daily average relative humidity were extracted.

2.3 Statistical Methods

Similar statistical approaches in our previous study were adopted []. Briefly, we used a generalized additive model (GAM) with a time series approach, in which daily CEVD deaths were reckoned following a quasi-Poisson distribution. The relationships between daily CEVD deaths and air pollution were considered linear, and both single-pollutant and two-pollutant models were calculated [, ]. Meteorological factors including relative humidity and daily average temperature were also included.

The models were listed as below (Model 1):

Mi~quasiPossion(μi)Log[μi]=α+β1Pi+ns(RH, df1)+ns(T, df2)+ns(daysi, df3)+Factor(DOWi)

Where Mi denotes the expected daily CEVD death (including total, men or women) counts of the day i; α is the intercept; βj refers to the corresponding coefficient of ambient air pollutants listed in the models; meteorological factors including daily average temperature (T) and average relative humidity (RH); ns(X, x) is the natural spline function; Factor is a categorical datum; DOWi is the day of the week on dayi; dfs is the degree of freedom of the natural spline function, which is selected from the model with the smallest unbiased risk estimator score; and dayi is the consecutively numbered day of the study, in which January 1, 2014 is day1 and January 1, 2015 is day366.

Excess Relative Risk (ERR10), which means the excess risks caused by a 10 unit elevation of pollutants, as well as Relative Risk (RR) was used to represent the risks of the ambient air pollutants on the CEVD population. The lag effects during lag0-2, lag0-3, lag0-5, and lag0-7 were explored. To further detect the gender-based differences in the vulnerability to air pollution, we used the Distributed lag nonlinear models (DLNMs) to explore the seasonal exposure-response relationships (during lag0-7) []. Considering the exposure-response relationships were not always linear or monotonic, we selected a natural spline function whose dfs were chosen according to the minimum unbiased risk estimator (UBRE) value [, , ].

Cross-validations were performed with a leave-1-year-out cross-validation strategy to check the stability of the findings. All statistical processes were performed with the help of the R statistic (version 3.5.1, the comprehensive R Archive Network, http://cran.r-project.org/).

3 Results

3.1 CEVD Mortality

There were 26,781 CEVD deaths in Shenyang City during 2014–2017, in which the daily deaths of men (57.79%) were significantly more than women (42.21%, P < 0.01). More deaths occurred in cold period (November to February) than in other periods. As to the deaths in each month, CEVD deaths were at peak in January with a 4 year average monthly CEVD deaths of 664, and declined to the minimum of 517 deaths in August. Details about the CEVD deaths in Shenyang during 2014-2017 were presented in Figure 1 and Table 1.

Figure 1 

Monthly CEVD deaths during 2014–2017. Generally CEVD deaths of women were more than those of men, and CEVD deaths were more in cold season than in warm season in both men and women.

Table 1

Descriptives of daily CEVD deaths during 2014–2017.


CountsDaily MeanSDCountsDaily MeanSDCountsDaily MeanSD

All Time2678118.414.891547710.643.55113047.772.93

3.2 Associations between pollutants and CEVD deaths

In the two-pollutant models, we found short-term O3 exposure, instead of PM2.5 or PM10, was positively associated with daily CEVD deaths, implying O3 might be the most suspicious air pollutant that lead to the excess CEVD mortality. For example, with the ERR10 of total CEVD deaths during lag0-1, lag0-2 and lag0-3 were 1.02% (95%CI 0.63%–1.4%), 0.95% (95%CI 0.62%–1.27%) and 0.74% (95%CI 0.46%–1.03%), respectively. Women tended to be more vulnerable to O3 that the ERRs of women were generally larger than those of men. In addition, the ERRs of women decreased with the extension of the lag period, while the ERRs of men increased at lag0-3 and then decreased with time. No positive associations were found between PM2.5 or PM10 and CEVD deaths. The results of the single-pollutant models varied slightly from those in the two-pollutant models. The association between short-term ambient air pollution exposure and CEVD mortality were shown in Table 2 and Figure 2 (from the two-pollutant model) and Figure S1 (from the single-pollutant model) of supplementary files.

Table 2

Cumulative ERR10 of each pollutant at different lags (Unit: %).

LagsPollutantsSingle-pollutant modelTwo-pollutant model



lag0-1PM2.5 –0.07–0.03–0.01–0.07–0.08–0.04
PM10 0.05––0.090.01
O3 0.091.42*1.01*0.11.43*1.02*
lag0-2PM2.5 –0.13–0.04–0.06–0.13–0.09–0.1
PM10 –0.04–0.1–0.05–0.04–0.12–0.07
O3 0.081.32*0.94*0.091.33*0.95*
lag0-3PM2.5 –0.07–0.05–0.06–0.08–0.08–0.09
PM10 0.02–0.12–0.050.01–0.13–0.06
O3 0.49*1.07*0.74*0.5*1.08*0.74*
lag0-5PM2.5 –0.08–0.02–0.06–0.09–0.04–0.08
PM10 –0.02–0.05–0.04–0.03–0.06–0.05
O3 0.39*0.75*0.54*0.4*0.76*0.55*
lag0-7PM2.5 –0.13–0.03–0.09–0.14–0.05–0.1
PM10 –0.05–0.05–0.05–0.06–0.06–0.06
O3 0.31*0.62*0.43*0.32*0.62*0.44*
Figure 2 

The associations between ambient air pollution and daily CEVD mortality. O3 seems to be the most suspicious pollutant that leads to excess CEVD mortality. Generally women with CEVD were more vulnerable to ambient air pollutants.

3.3 Seasonal analyses

The Seasonal analyses by two-pollutant model indicated that the CEVD population were more vulnerable to air pollution in cold seasons (winter and spring). In cold seasons, for example, the ERRs of O3 among both men and women were mostly significantly above 0. However, we can hardly find significant positive associations in summer or autumn. Some negative associations were found between PMs and CEVD deaths. They might suggest non-linear relationships between short-term PMs exposure and CEVD deaths, since air pollution exposure cannot be a protective factor. Even in the seasonal analyses, the ERRs of women were generally larger than those of men, and this was more obvious in ERRs of O3 during winter. The details of seasonal analyses are shown in Table 3.

Table 3

Seasonal analyses of the associations between short-term ambient air pollution exposure and CEVD death.



Womenlag0-1PM2.5 –0.1–
PM10 –0.28–0.350.060.2
O3 1.39*0.071.51*4.12*
lag0-2PM2.5 0.19–0.110.03–0.24
PM10 –0.16–0.390.04–0.09
O3 1.65*0.320.93.2*
lag0-3PM2.5 0.06–0.540.05–0.26
PM10 –0.18–0.80.05–0.15
O3 1.25*0.340.452.78*
lag0-5PM2.5 0.2–0.660.13–0.41*
PM10 0.11–0.750.13–0.18
O3 1.36*0.29–0.22.23*
lag0-7PM2.5 0.13–0.720.08–0.25
PM10 0.02–0.9*0.09–0.04
O3 1.37*0.33–0.312.44*
Menlag0-1PM2.5 –0.29–1.680.04–0.12
PM10 0.67*–0.89–0.020.22
O3 1.29*–0.13–0.991.57
lag0-2PM2.5 –0.29–1.48–0.06–0.27
PM10 0.25–0.9–0.10.08
O3 1.24*–0.11–1.041.41
lag0-3PM2.5 –0.69–1.41*–0.1–0.12
PM10 –0.03–0.93*–0.110.17
O3 1.47*–0.05–1.14*1.04
lag0-5PM2.5 –0.61–2.01*–0.170.03
PM10 –0.02–1.47*–0.170.18
O3 0.99*–0.03–0.81.34*
lag0-7PM2.5 –0.6–1.16*–0.28*0.12
PM10 –0.02–1.33*–0.27*0.23*
O3 0.54–0.14–0.631.57*

The seasonal exposure-response relationships of PM2.5, PM10 and O3 by two-pollutant were shown in Figures 3, 4 and 5, respectively. Among these curves, some unignorable non-linear curves confirmed our hypothesis that the exposure-response relationships were not always linear. For example, the curves of O3 in autumn (Figure 5) of both genders are obviously not linear, which usually increased first and then decreased. Some curves also increased and then decreased, but they showed a generally decreasing trend, such as the PM10 curve of men in spring and winter. This can possibly explain some of the negative associations between PMs and CEVD mortality in the linear models. Another notable result is that the trends of curves sometimes differed significantly between men and women. For example, the O3 exposure-response curve of men in autumn was obviously bell-shaped, but the curve of women in autumn showed a slightly bell-shaped trend. And the PM2.5 or PM10 curves of men in winter generally increased, while those curves of women were generally decreased.

Figure 3 

Comparison of Seasonal exposure-response curves of PM2.5-CEVD mortality between men and women. Most curves of both men and women tended to be non-linear, and the RRs of women were generally larger than those of men. However, men CEVD mortality were shown to be more vulnerable to PM2.5 in winter.

Figure 4 

Comparison of Seasonal exposure-response curves of PM10-CEVD mortality between men and women. Some curves showed both non-linear and generally increasing or decreasing trend, which may provide a possible explanation for the insignificant associations between PMs and daily CEVD deaths. Generally women were more vulnerable to PM10 in spring and summer, but men were more vulnerable in autumn and winter.

Figure 5 

Comparison of Seasonal exposure-response curves of O3-CEVD mortality between men and women. Significant differences of gender-based vulnerability were observed. In spring and autumn, men were more vulnerable, but in the other seasons women were more vulnerable.

No significant variations were found when using the eliminated data of one-year each time, suggesting good robustness of the results.

4 Discussion

In this study, we found limited evidences that PMs could induce increased CEVD mortality when adopting linear relationships. However, O3 was shown to associate with increased daily CEVD mortality in both men and women, and women tended to be more vulnerable. These findings implied that O3 could be the primary risk pollutant of the excess CEVD deaths, and women were more vulnerable to ambient air pollutants than men.

We found the ERRs of women were larger than those of men in general. This accorded with the previous outcomes among Shanghai, Hong Kong and Japanese population [, , , ]. They are also supported by the mechanism studies indicating that women were more susceptible to oxidative stress (OS) []. Bell et al. suggested geographical variation in the magnitude of gender difference, which were larger in the Northeast than in other regions of the US []. This is also similar to our study that was performed in the North, since some studies that were performed in middle or southern China indicated no gender differences or men were more vulnerable. However, a recent comprehensive assessment by EPA noted there were not enough evidences to draw a conclusion on this issue []. Therefore more relevant studies are still in need.

In addition, we found men might be more vulnerable to PM2.5 in winter. Although more papers reported women as the susceptible gender to ambient air pollution, there are also studies indicated that the mortality of men was more sensitive to PM2.5 [, , ]. These variations of the gender-based vulnerability to PM2.5 across the studies may be due to different components of PMs as well as their proportions, since the components vary significantly across regions [, ]. The haze is also a key risk for circulatory mortality, and most of the haze days in Shenyang occurred in winter []. This may also aggravate the adverse effects of short-term air pollution on the CEVD population. As to short-term O3 exposure, it has been frequently reported to associate with elevated risk of the increased circulatory mortality. In this study, we further confirmed ambient O3 as a fatal risk for CEVD population, as well as the gender-based vulnerability differences. The mechanisms that O3 induces excess CEVD mortality are not yet clear, but the influences of O3 on circulatory system have been validated by previous studies. One possible path is by inducing hypertension, when ET-1 were generated after short term exposure to O3 and then further elevated the blood pressure []. And vasoconstriction caused by O3 exposure not only induces hypertension, but can also slow the vascular flow, which was also fatal for the CEVD population [].

Seasonal analyses suggested that both men and women CEVD were more vulnerable to air pollution in cold seasons (winter and spring), and It was more evident in the ERRs of O3. Some previous papers reported that O3 could elevate the cardiovascular mortality only in warm seasons, but more recent studies (especially those carried out in east Asia) reported more robust associations in cold seasons [, , , ]. In this study, the seasonal variations could largely attribute to the cold climate in Shenyang City. The mean of the daily lowest temperature was –13.5°C in Shenyang’ winter during 2014 to 2017. As the triggers of CEVD deaths, the increased blood pressure and blood viscosity can be induced by the exposure to low ambient temperatures []. Several studies also indicated that low ambient temperature can exacerbate the fatal effects of air pollution []. Therefore, even if average daily O3 concentration were significantly lower than in cold seasons than in warm seasons, the CEVD population of Shenyang could also be more sensitive.

It is interesting that the result of CEVD differed from those of CHD in Liu et al., that CHD mortality was associated with PM2.5 while CEVD mortality was associated with O3. This may be attributed to the negative correlation between O3 and PM2.5, and the different triggers of acute onset between CHD and CEVD. First, in this study the ERRs of O3 is higher than PMs, indicating that O3 induces more risks in the CEVD population of Shenyang. Second, the paths of air pollution-induced acute CHD and acute CEVD may be different. Many studies proved that short-term exposure to either PM2.5 or O3 would result in acute increase of blood pressure, which is possibly how short-term air pollution exposure increased CEVD mortality [, , ]. However, the role of acute increasing blood pressure in inducing the onset of acute CHD was not as strong as in inducing acute CEVD (such as hemorrhagic stroke). Third, the concentration of PM2.5 and O3 were negatively correlated in Shenyang (shown in supplementary files). In summary, it is possible that in the linear models O3 associated with CEVD mortality while PM2.5 not.

It is not rare that no associations between PMs and daily CEVD deaths were found in men or women. Several studies failed to prove the short-term PMs exposure as a risk factor that could elevate daily cardiovascular hospital admissions or mortality [, ]. One main reason for these inconsistencies is that the exposure-response relationships can be non-linear [, , ]. Moreover our previous study also indicated that this relationship between ASHD mortality and air pollution seemed to be non-linear []. Just as shown in Figures 3, 4 and 5, some exposure-response curves of PMs were bell-shaped, so if they were treated as linear the generally tendency can be decreasing. On the other hand, the inconsistencies may also be attributed to spatial variation; however the potential influences of detailed individual factors remained to be explored.

This study provided the further evidence that women CEVD population were more vulnerable to ambient air pollution, especially more vulnerable to O3. It also offered evidence for policy making. The public attention should not only focus on the total air quality, but also on each type of pollutants as well as their susceptible population. In addition, the concentrations of O3 are usually higher on sunny days with lower PMs concentrations, so people are more likely to misjudge the air quality and ignore its potential risks. And currently the most efforts to improve air quality focus on controlling the concentrations of PMs. Therefore the ambient O3 also deserved great attention in both environment protection and in public education.

There are also some deficiencies in this study. First, although over 26,000 CEVD death cases were included, the outcomes of the multi-center studies can be more convincible. Second, the findings in the exposure-response curves cannot be completely explained by the current evidences. Therefore, more studies on the gender-based differences in the vulnerability to ambient air pollution are needed.

5 Conclusions

Women were generally more vulnerable to short-term air pollution exposure than men in Shenyang. Ozone can induce excess daily CEVD mortality, which deserves more attention in public health education.

Additional File

The additional file for this article can be found as follows:

Supplementary file

Figure S1. DOI: https://doi.org/10.5334/gh.849.s1