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

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)
GRDPi (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

aIschemic 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

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

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

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

eAir 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

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

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

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

iGross 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.

Fig. 1

Associations between air pollutant concentrations (□: CO, ○: SO2, △: NO2, ◇: O3, ⊠: PM10) and a ischemic heart disease (IHD), b cerebrovascular disease (CVD), c pneumonia (PN), and d chronic lower respiratory disease (CLRD) mortality rates

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 (Fig. 2 and Fig. 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.

Fig. 2

Associations between air pollutant concentrations (□: CO, ○: SO2, △: NO2, ◇: O3, ⊠: PM10) and a ischemic heart disease (IHD), b cerebrovascular disease (CVD), c pneumonia (PN), and d chronic lower respiratory disease (CLRD) mortality rates in the capital (blue) or non-capital (red) areas

Fig. 3

Associations between air pollutant concentrations (□: CO, ○: SO2, △: NO2, ◇: O3, ⊠: PM10) and a ischemic heart disease (IHD), b cerebrovascular disease (CVD), c pneumonia (PN), and d chronic lower respiratory disease (CLRD) mortality rates in urban (magenta) or rural (green) 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 intervals

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 (Fig. 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 TNFα and IL-1β, 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.