Association between long-term exposure to air pollutants and cardiopulmonary mortality rates in South Korea

DOI: https://doi.org/10.21203/rs.2.23984/v4

Abstract

Background: The association between long-term exposure to air pollutants, including nitrogen dioxide (NO2), carbon monoxide (CO), sulfur dioxide (SO2), ozone (O3), and particulate matter 10 μm or less in diameter (PM10), and mortality by ischemic heart disease (IHD), cerebrovascular disease (CVD), pneumonia (PN), and chronic lower respiratory disease (CLRD) is unclear. We investigated whether living in an administrative district with heavy air pollution is associated with an increased risk of mortality by the diseases through an ecological study using South Korean administrative data over 19 years.

Methods: A total of 249 Si-Gun-Gus, unit of administrative districts in South Korea were studied. In each district, the daily concentrations of CO, SO2, NO2, O3, and PM10 were averaged over 19 years (2001–2018). Age-adjusted mortality rates by IHD, CVD, PN and CLRD for each district were averaged for the same study period. Multivariate beta-regression analysis was performed to estimate the associations between air pollutant concentrations and mortality rates, after adjusting for confounding factors including altitude, population density, higher education rate, smoking rate, obesity rate, and gross regional domestic product per capita. Associations were also estimated for two subgrouping schema: Capital and non-Capital areas (77:172 districts) and urban and rural areas (168:81 districts).

Results: For IHD, higher SO2 concentrations were significantly associated with a higher mortality rate, whereas other air pollutants had null associations. For CVD, SO2 and PM10 concentrations were significantly associated with a higher mortality rate. For PN, O3 concentrations had significant positive associations with a higher mortality rate, while SO2, NO2, and PM10 concentrations had significant negative associations. For CLRD, O3 concentrations were associated with an increased mortality rate, while CO, NO2, and PM10 concentrations had negative associations. In the subgroup analysis, positive associations between SO2 concentrations and IHD mortality were consistently observed in all subgroups, while other pollutant-disease pairs showed null, or mixed associations.

Conclusion: Long-term exposure to high SO2 concentration was significantly and consistently associated with a high mortality rate nationwide and in Capital and non-Capital areas, and in urban and rural areas. Associations between other air pollutants and disease-related mortalities need to be investigated in further studies. 

Background

There is increasing evidence on the harmful associations between air pollution and cardiopulmonary mortality [1-9]. Many short-term studies have reported compelling evidence on such associations [2, 3, 5, 6, 9]; however, relatively limited number of long-term studies were performed. This may be partly because collecting and analyzing long-term air pollution and cardiopulmonary mortality data together are relatively difficult than collecting and analyzing short-term data together.

Although a meta-analysis by Vodonos et al. [8], and a recent large representative cohort study by Pope et al. [4] provide compelling evidence on long-term associations between air pollution and cardiopulmonary mortality in a cohort design, these studies have only focused on exposure to fine particulate matter.

However, long-term association studies on cardiopulmonary mortality performed in South Korea [10-12] have focused on particulate matter 10 μm or less in diameter (PM10). Kim et al. [11] used the National Health Insurance Service sample cohort representing the general population in South Korea and estimated the individual exposure to PM10 as a 5-year average (2002–2006); they found positive but insignificant associations between PM10 exposure and cardiopulmonary diseases. Tran et al. [12] found associations between pneumonia mortality and PM10 concentrations (2005–2015), and Kim et al. [10] reported the cardiopulmonary mortality benefits of PM10 reduction. Both studies were conducted in 25 districts in Seoul, the capital of South Korea.

However, there is a knowledge gap regarding whether exposure to higher concentrations of air pollutants, including carbon monoxide (CO), sulfur dioxide (SO2), nitrogen dioxide (NO2), ozone (O3), and PM10 in a residential district in South Korea over a long term, such as 19 years, would be associated with higher cardiopulmonary mortality. We investigated a total of 249 districts in South Korea from 2001 to 2018 to evaluate the associations between air pollutants including CO, SO2, NO2, O3, and PM10, and age-adjusted mortality rates related to ischemic heart disease (IHD), cerebrovascular disease (CVD), pneumonia (PN), and chronic lower respiratory disease (CLRD) nationwide after adjusting for altitude, population density, higher education rate, smoking rate, obesity rate, and gross regional domestic product per capita (GRDP). Because there may be uncaptured socioeconomic or cultural differences between the capital and non-capital areas, and urban and rural areas, we also investigated whether the associations found in the nationwide setting remained qualitatively similar in subgroups.

Methods

Study design and ethics

The study used an ecological design. Ethical approval was not required because the study used only publicly accessible, national statistics database.

Air pollution

CO, SO2, NO2, O3, PM10, and PM2.5 concentrations measured by the National Ambient Air Quality Monitoring Information System, are publicly accessible via the AirKorea website. In South Korea, there are 332 measurement stations. Due to a shortage in measurement stations before 2015, PM2.5 concentrations were not assessed in the current study. The average concentrations of each pollutant per day were collected for each station. The air pollution measurement station system was not directly matched to the Si-Gun-Gu district system, on which populations and mortality statistics dataset were based. Longitudes and latitudes of all air pollution measurement stations and all districts’ administrative authorities offices were obtained, and then the average air pollutant concentrations throughout the study period for each administrative office were estimated by linearly interpolating air pollutant measurements from the nearest three stations. Python programming language version 2.7 (Python Software Foundation, Beaverton, Oregon, United States) was used in this procedure. For each air pollutant and district, the average air pollutant concentrations throughout the study period (2001–2018) was computed. The above mentioned method is largely similar to that used in our previous study [13].

Mortality statistics

According to the 10th revision of the International Classification of Diseases, age-adjusted mortality rates were obtained from death certificates and population census data were obtained from the Korean Statistical Information Service (KOSIS) during the study period (2001–2018). In detail, mortality rates of IHD (I20–I25), CVD (I60–I69), PN (J12–J18), and CLRD (J40–J47) were obtained. As of 2018, there were 250 Si-Gun-Gus in South Korea as of 2018. Si-Gun-Gu is a level in the Korean administrative area system, which is comparable to counties in the United States. All Si-Gun-Gus in South Korea were included in this study, except Sejong-Si, which was newly designated in 2012. The mortality rates per 100,000 were age-adjusted by using the standard population as of July 1, 2010 in South Korea. The age-adjusted mortality rates were extracted from the KOSIS database and calculated as follows:

Confounding factors

The annual average of confounding factors including altitude, smoking rate (rate of current smokers adjusted for the age of the national standard population), higher education rate (rate of >15-year-old persons with college education or more in the district), obesity rate (rate of persons with body mass index > 25 kg/m2), population density based on the 2010 Census, and gross regional domestic production (GRDP) per capita as of 2011 were accessed for all districts using the KOSIS.

Statistical analysis

Data are presented as median, interquartile range, and 95% confidence interval (95% CI) where applicable. Per interquartile increase of air pollutant concentrations, multivariable beta regression [14, 15] models were built, and the odds ratio of each air pollutant to the mortality rates were estimated while adjusting for the confounding factors. A basic bootstrap method was utilized to estimate the 95% CIs for the odds ratios. Statistical analyses were performed by using R statistics software version 3.6.3 (R Foundation for Statistical Computing, Vienna, Austria).

Subgrouping schema

Two subgrouping schemas were applied: capital and non-capital areas and urban and rural areas. Among the 249 districts, the capital area included 77 districts in Seoul, the capital city of South Korea, Incheon, and Gyeonggi-do. These 77 districts are geographically in the vicinity of the capital and linked to each other by public transportation such as the subway system. The capital area contains 49% of the total South Korean population. The non-Capital area consists of the remaining 172 districts. The urban subgroup contained 168 districts identified as Gu or Si, whereas the rural subgroup contained 81 districts identified as Gun.

Results

Table 1 shows the medians and interquartile ranges of the mortality rates of the four diseases, concentrations of the five air pollutants, and confounding factors averaged from 2001 to 2018. The population of South Korea, as of 2010 (in the middle of the study period) was 50,515,666 persons. Throughout the study period, a total of 4,558,640 all-cause mortalities were recorded. Among them, 242,711 deaths were attributed to IHD, 509,740 deaths to CVD, 160,174 deaths to PN, and 138,271 deaths to CLRD.

Table 1. Characteristics of the study area 

Characteristics

Median (first–third quartile range)

Number of districts analyzed

249

IHDa: Age-adjusted mortality rate (per 100,000)

21.6 (19.4–25.3)

CVDb: Age-adjusted mortality rate (per 100,000)

50.7 (46.5–54.5)

PNc: Age-adjusted mortality rate (per 100,000)

12.2 (10.9–13.5)

CLRDd: Age-adjusted mortality rate (per 100,000)

13.8 (11.7–16.0)

Carbon monoxide (ppb) e

554 (476–624)

Nitrogen dioxide (ppb)

20.1 (14.6–26.8)

Sulfate dioxide (ppb)

4.89 (4.15–5.56)

Ozone (ppb)

24.3 (21.6–27.6)

PM10 (µg/m3)

50.2 (46.1–55.5)

Altitude (m)

124 (59.9–220)

Population density (per km2)

361 (109–6042)

Higher-education rate f (%)

34.4 (22.9–43.0)

Smoking rate g (%)

25.1 (23.5–26.8)

Obesity rate h (%)

22.4 (20.9–24.4)

GRDP (million won)

22.9 (16.7–29.0)

From 2001 to 2018, the population of all 249 South Korean administrative districts, Si-Gun-Gus, were studied. The South Korean population, as of 2010 (in the middle of the study period) was 50,515,666 persons. Throughout the study period, a total of 4,558,640 all-cause deaths were recorded. Among them, 242,711 deaths were attributed to IHD, 509,740 deaths to CVD, 160,174 deaths to PN, and 138,271 deaths to CLRD. 

Ischemic Heart Disease (IHD): Total number of deaths by IHD in 2001–2018 was 242,711. The annual raw mortality rates throughout the study period were adjusted for each district’s age distribution to the standard South Korean population in 2010.

Cerebrovascular Disease (CVD): Total number of deaths by CVD in 2001–2018 was 509,740. 

Pneumonia (PN): Total number of deaths by PN in 2001–2018 was 160,174. 

Chronic Lower Respiratory Disease (CLRD): Total number of deaths by CLRD in 2001–2018 was 138,271. 

Air pollution data from 2001 to 2018, in daily mean concentrations according to the positions of monitoring stations were accessed using the AirKorea database. An interpolation model based on a geographical information system was applied to yield the average air pollutant concentration throughout the study period of the corresponding districts. 

Rate of >15-year-old persons with college education or more in 2010

Rate of current smokers adjusted for the age of the national standard population in 2010

Rate of population with body mass index > 25 kg/m2 adjusted for the age of the national standard population in 2010

Gross Regional Domestic Product per capita (GRDP) in 2011


For IHD, an increased SO2 concentrations were significantly associated with a higher mortality rate (odds ratio per interquartile range [OR] 1.09; 95% CI, 1.05–1.12), whereas other air pollutants had null associations. For CVD, SO2 (OR 1.03; 95% CI 1.01–1.05) and PM10 (OR 1.04; 95% CI 1.02–1.07) concentrations had significant associations with a higher mortality rate. For PN, O3 (OR 1.06; 95% CI 1.02–1.09) concentrations had significant positive associations with a higher mortality rate, while SO2 (OR 0.968; 95% CI 0.943–0.994), NO2 (OR 0.893; 95% CI 0.861–0.923), and PM10 (OR 0.947; 95% CI 0.919–0.980) concentrations had significant negative associations. For CLRD, O3 concentrations were associated with an increased mortality rate (OR 1.08; 95% CI 1.01–1.13), while CO (OR 0.891; 95% CI 0.856–0.935), NO2 (OR 0.822; 95% CI 0.780–0.865), and PM10 (OR 0.934; 95% CI 0.902–0.977) concentrations had negative associations. Figure 1 shows the ORs and 95% CIs of the estimated associations.

In the subgroup analysis that divided the 249 districts into capital or non-capital areas (77:172 districts) and into urban or rural areas (168:81 districts), positive associations between SO2 concentrations and IHD mortality were consistently observed in all subgroups, while other pollutant-disease pairs showed null or mixed associations(Figure 2 and Figure 3). Table 2 summarizes the qualitative associations between disease mortality and air pollutant concentrations in the corresponding subgroup schema. A ‘+’ denote a significant positive association, a ‘–’ to negative, or blank to insignificant. For example, associations between CVD mortality and NO2 concentrations exhibited a paradoxical pattern in the subgroup analysis because a significant negative association was found in capital districts but positively associated in non-capital areas. However, null associations were found nationwide and in urban and rural areas. In contrast, significant negative associations were found between NO2 concentrations and CLRD mortality nationwide and in capital, non-capital, and urban areas; however, positive associations were found in rural areas.

Table 2. Significant associations between mortality rates and air pollutant concentrations nationwide and in capital, non-capital, urban, and rural areas. 

Ischemic heart disease mortality rate

Air pollutant

Nationwide

Capital

Non-capital

Urban

Rural

CO

 

 

 

 

SO2

+

+

+

+

+

NO2

 

 

+

 

 

O3

 

 

 

+

 

PM10

 

 

 

 

 

Cerebrovascular disease mortality rate

Air pollutant

Nationwide

Capital

Non-capital

Urban

Rural

CO

 

 

 

 

+

SO2

+

+

+

+

 

NO2

 

+

 

 

O3

 

 

 

 

PM10

+

 

+

 

+

Pneumonia mortality rate

Air pollutant

Nationwide

Capital

Non-capital

Urban

Rural

CO

 

 

 

 

 

SO2

 

 

 

 

NO2

 

 

O3

+

 

 

+

 

PM10

 

 

 

Chronic Lower Respiratory Disease mortality rate

Air pollutant

Nationwide

Capital

Non-capital

Urban

Rural

CO

 

 

SO2

 

 

 

 

+

NO2

+

O3

+

+

 

+

 

PM10

 

 

 

 

“+” denotes a positive association, which means increased air pollutant concentrations are associated with higher mortality rates in the study districts.  

“–” denotes a negative association, which that means the opposite of “+”. 

Blank denotes null associations. 

The significance of the associations was determined using 95% confidence inte


Discussion

In the nationwide analysis, we found significant positive associations between SO2 concentrations and IHD and CVD mortality, PM10 concentrations and CVD mortality, and O3 concentrations and PN and CLRD mortality, which were consistent with those reported in previous studies [1, 5, 9, 16]; however, direct comparisons of the effect sizes are not appropriate because of differences in the study design, area, and period (Figure 1). However, significant negative associations between SO2 concentrations and PN mortality, NO2 concentrations and PN and CLRD mortality, and PM10 concentrations and PN and CLRD mortality have not been reported before and are hard to explain intuitively. In subgroup analysis, we consistently found positive associations between SO2 concentrations and IHD mortality regardless of the subgrouping schema; hence, we can confidently state that long-term (19 years) exposure to increased SO2 concentrations is associated with increased IHD mortality. However, for other disease mortality-air pollutant pairs, it is precarious to conclude that there is a positive or negative association because the beta-regression results differed among the subgroups. Associations between SO2 concentrations and CVD mortality were significantly positive nationwide and in capital, non-capital, and urban areas and marginally positive in rural areas; hence, there probably is a positive correlation. For CLRD mortality-NO2 pair, the associations are negative in some subgroups and positive in others, which is paradoxical. This paradoxical association patterns among subgroups may imply that there is an important but unidentified confounding factor that was not incorporated in the regression model. For example, in both urban and rural areas, NO2 concentrations may be harmful for respiratory health and there may have been an increase medical service usage such as emergency department visits or hospitalization related to CLRD, but accessibility or the quality of medical service may be different in urban and rural areas. In our previous study on the association between air pollution and the incidence and mortality rates of breast cancer, air pollution was positively associated with the incidence rates but not with the mortality rates [13]. This conjecture could be resolved with data on the incidence rates of CLRD per district, which is not presently available.

There are suggested pathways linking long-term exposure to air pollution and cardiopulmonary disease mortality. In a study by Hoek et al., PM10 concentrations were associated with a significant increase in blood pressure and induced infection and inflammation in circulatory and respiratory diseases [17]. Hiraiwa et al. have suggested that excess cytokines such as interleukin (IL)-1, IL-6, and tumor necrosis factor (TNF) can induce vascular events in patients with chronic obstructive pulmonary disease via systemic oxidative stress and inflammation in the lung to promote endothelial dysfunction and atherosclerotic plaque rupture, possibly leading to acute cardiac events or stroke [18]. According to Mukae et al, human alveolar macrophages, when exposed to high PM10 concentrations, can phagocytose these particles and produce an array of cytokines such as TNFa and IL-1b, which are part of the innate immune response [19]. Hence, long-term exposure to PM10 may aggravate premature mortality from CVD and CLRD.

The consistent positive associations between SO2 and IHD and CVD mortality that we found agree with previous publications those reported similar associations in short-term. In a systematic review on air pollution and stroke by Shah et al [20] reported significant positive associations between SO2 and mortality and hospital admissions by stroke. Hong et al [21] found a significant positive association between ischemic stroke mortality and SO2 and total suspended particulates in the short term (0–3 lagged days) in Seoul, South Korea. Qian et al [22] reported significant associations of cardiovascular disease mortality with PM10 and SO2 in Wuhan, China. Moolgavkar et al [23] also found out similar associations in Los Angeles, United States. Wichmann et al [24] reported associations between cardiovascular and cerebrovascular mortality with SO2 in Cape Town, South Africa, too. Amancio et al [25] also found positive associations between short term SO2 exposure and circulatory disease and stroke mortality in Brazil with ecological study design. Moreover, Chung et al [26] found significant positive associations between PM10 and SO2 and cardioembolic stroke incidence on the basis of the Clinical Research Center for Stroke 5th division centers registry data in South Korea. There are studies suggesting etiological links between SO2 exposure and IHD and CVD. Routledge et al [27] found out that SO2 exposure reduce cardiac vagal control, a response that would be expected to increase susceptibility to ventricular arrhythmia. Szyszkowicz et al [28] suggested that SO2-derived acidic compounds may penetrate the brain barrier to mediate abnormal brain neural activity or brain ischemia.  

The present study is an ecological analysis rather than an individual-level cohort study nor a case-control study. From the ecological nature of our study arises a limitation. Because the unit of analysis was a district and not an individual, we did not obtain patient-specific information, such as comorbidities, medication, occupational history, and patient-specific exposure to air pollution, that are of considerable importance in cardiopulmonary mortality. Another limitation is the lack of migration history data. We contend that migrations would not significantly impact the current study because only approximately 10% of the population moved between different districts in South Korea between 2003 and 2013 [11, 29].  In addition, there were differences in the density of distributed stations per km2(6 times denser in the capital area than in non-capital areas), which could cause potential bias of the daily measurements of the air pollutant concentrations.

Conclusion

Long-term exposure (19 years) to high SO2 concentrations was consistently and significantly associated with a high mortality rate nationwide and in capital and non-capital areas and in urban and rural areas. Associations between other air pollutants and disease-related mortalities need to be investigated in further studies.

Abbreviations

IHD: ischemic heart disease

CVD: cerebrovascular disease

PN: pneumonia

CLRD: chronic lower respiratory disease

CO: carbon monoxide

SO2: sulfur dioxide

NO2: nitrogen dioxide

O3: ozone

PM10: particulate matter 10 μm or less in diameter

GRDP: gross regional domestic product per capita

KOSIS: Korean Statistical Information Service

IQR: interquartile range

OR: odds Ratio

IL: interleukin

TNF: tumor necrosis factor

Declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Availability of data and materials

All datasets used during the current study are publicly available via AirKorea (https://www.airkorea.or.kr/) and Korean Statistical Information Service (https://kosis.kr/).

Competing interests

The authors declare that they have no competing interests.

Funding

None

Author’s contributions

JK, JH, and NK conceived the study design. JK, JH, and HJB have made contributions to the acquisition of data. JK, JH, HY, HJB and NK did the statistical analysis and data interpretation. MJ and NK checked the integrity of the manuscript including consistency of the analysis results and data interpretation. JH, JK and NK wrote the first draft. JH, HY, HJB, and MJ critically revised the draft. All authors have reviewed and approved the final draft.

Acknowledgements

None

References

  1. Cromar KR, Gladson LA, Ghazipura M, Ewart G: Estimated Excess Morbidity and Mortality Associated with Air Pollution above American Thoracic Society-recommended Standards, 2013-2015 American Thoracic Society and Marron Institute Report. Ann Am Thorac Soc 2018, 15(5):542-551.
  2. Di Q, Dai LZ, Wang Y, Zanobetti A, Choirat C, Schwartz JD, Dominici F: Association of Short-term Exposure to Air Pollution With Mortality in Older Adults. Jama-J Am Med Assoc 2017, 318(24):2446-2456.
  3. Lu F, Zhou L, Xu Y, Zheng TZ, Guo YM, Wellenius GA, Bassig BA, Chen XD, Wang HC, Zheng XY: Short-term effects of air pollution on daily mortality and years of life lost in Nanjing, China. Sci Total Environ 2015, 536:123-129.
  4. Pope CA, Lefler JS, Ezzati M, Higbee JD, Marshall JD, Kim SY, Bechle M, Gilliat KS, Vernon SE, Robinson AL et al: Mortality Risk and Fine Particulate Air Pollution in a Large, Representative Cohort of U.S. Adults (vol 127, 077007, 2019). Environ Health Persp 2019, 127(9).
  5. Ren M, Li N, Wang Z, Liu YS, Chen X, Chu YY, Li XY, Zhu ZM, Tian LQ, Xiang H: The short-term effects of air pollutants on respiratory disease mortality in Wuhan, China: comparison of time-series and case-crossover analyses. Sci Rep-Uk 2017, 7.
  6. Sun SZ, Tian LW, Cao WN, Lai PC, Wong PPY, Lee RSY, Mason TG, Kramer A, Wong CM: Urban climate modified short-term association of air pollution with pneumonia mortality in Hong Kong. Sci Total Environ 2019, 646:618-624.
  7. Thurston GD, Kipen H, Annesi-Maesano I, Balmes J, Brook RD, Cromar K, De Matteis S, Forastiere F, Forsberg B, Frampton MW et al: A joint ERS/ATS policy statement: what constitutes an adverse health effect of air pollution? An analytical framework. Eur Respir J 2017, 49(1).
  8. Vodonos A, Abu Awad Y, Schwartz J: The concentration-response between long-term PM2.5 exposure and mortality; A meta-regression approach. Environ Res 2018, 166:677-689.
  9. Vodonos A, Kloog I, Boehm L, Novack V: The impact of exposure to particulate air pollution from non-anthropogenic sources on hospital admissions due to pneumonia. Eur Respir J 2016, 48(6):1791-1794.
  10. Kim D, Kim J, Jeong J, Choi M: Estimation of health benefits from air quality improvement using the MODIS AOD dataset in Seoul, Korea. Environ Res 2019, 173:452-461.
  11. Kim OJ, Kim SY, Kim H: Association between Long-Term Exposure to Particulate Matter Air Pollution and Mortality in a South Korean National Cohort: Comparison across Different Exposure Assessment Approaches. Int J Env Res Pub He 2017, 14(10).
  12. Tran H, Kim J, Kim D, Choi M, Choi M: Impact of air pollution on cause-specificmortality in Korea: Results from Bayesian Model Averaging and Principle Component Regression approaches. Sci Total Environ 2018, 636:1020-1031.
  13. Hwang J, Bae H, Choi S, Yi H, Ko B, Kim N: Impact of air pollution on breast cancer incidence and mortality: a nationwide analysis in South Korea. Sci Rep 2020, 10(1):5392.
  14. Cribari-Neto F, Zeileis A: Beta Regression in R. J Stat Softw 2010, 34(2):1-24.
  15. Ferrari SLP, Cribari-Neto F: Beta regression for modelling rates and proportions. J Appl Stat 2004, 31(7):799-815.
  16. Cromar KR, Gladson LA, Ewart G: Trends in Excess Morbidity and Mortality Associated with Air Pollution above American Thoracic Society-Recommended Standards, 2008-2017. Ann Am Thorac Soc 2019, 16(7):836-845.
  17. Hoek G, Krishnan RM, Beelen R, Peters A, Ostro B, Brunekreef B, Kaufman JD: Long-term air pollution exposure and cardio- respiratory mortality: a review. Environ Health 2013, 12(1):43.
  18. Hiraiwa K, van Eeden SF: Nature and consequences of the systemic inflammatory response induced by lung inflammation. Lung Inflammation 2014:79.
  19. Mukae H, Vincent R, Quinlan K, English D, Hards J, Hogg JC, van Eeden SF: The effect of repeated exposure to particulate air pollution (PM10) on the bone marrow. Am J Respir Crit Care Med 2001, 163(1):201-209.
  20. Shah ASV, 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-British Medical Journal 2015, 350.
  21. Hong YC, Lee JT, Kim H, Kwon HJ: Air pollution - A new risk factor in ischemic stroke mortality. Stroke 2002, 33(9):2165-2169.
  22. Qian Z, He Q, Lin HM, Kong L, Zhou D, Liang S, Zhu Z, Liao D, Liu W, Bentley CM et al: Part 2. Association of daily mortality with ambient air pollution, and effect modification by extremely high temperature in Wuhan, China. Res Rep Health Eff Inst 2010(154):91-217.
  23. Moolgavkar SH: Air pollution and daily mortality in three U.S. counties. Environ Health Perspect 2000, 108(8):777-784.
  24. Wichmann J, Voyi K: Ambient air pollution exposure and respiratory, cardiovascular and cerebrovascular mortality in Cape Town, South Africa: 2001-2006. Int J Environ Res Public Health 2012, 9(11):3978-4016.
  25. Amancio CT, Nascimento LF: Association of sulfur dioxide exposure with circulatory system deaths in a medium-sized city in Brazil. Braz J Med Biol Res 2012, 45(11):1080-1085.
  26. Chung JW, Bang OY, Ahn K, Park SS, Park TH, Kim JG, Ko Y, Lee S, Lee KB, Lee J et al: Air Pollution Is Associated With Ischemic Stroke via Cardiogenic Embolism. Stroke 2017, 48(1):17-23.
  27. 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-227.
  28. Szyszkowicz M, Porada E, Tremblay N, Grafstein E: Sulfur dioxide and emergency department visits for stroke and seizure. Stroke Res Treat 2012, 2012:824724.
  29. Statistics Korea News [http://kostat.go.kr/portal/eng/news/3/index.board]