Institute of Health & Environment Seoul National University
[ Article ]
The Korean Journal of Public Health - Vol. 53, No. 1, pp.17-26
ISSN: 1225-6315 (Print)
Print publication date Mar 2016
DOI: https://doi.org/10.17262/KJPH.2016.03.53.1.17

대기오염이 천식으로 인한 입원에 미치는 영향 (서울, 상하이를 대상으로)

전웅1 ; 김호1, *
1서울대학교 보건대학원 보건학과, 서울대학교 보건환경연구소
Associations of Air Pollution on Asthma Hospitalization in Seoul and Shanghai
Woong Jeon1 ; Ho Kim1, *
1Department of Biostatistics and Epidemiology, Graduated School of Public Health & Institute of Health and Environment, Seoul NationalUniversity

Correspondence to: *Ho Kim, Ph.D. ( hokim@snu.ac.kr., 02-880-2702) Department of Biostatistics and Epidemiology, School of Public Health, Seoul National University, 1 Gwanak-Ro, Gwanak-gu, Seoul 151-752, South Korea.

Abstract

Objectives

This study intends to analyze the associations of air pollution on asthma hospitalization and healthcare utilization in Seoul and Shanghai using national health insurance data (NHI-Korea health insurance, SHIB-Shanghai health insurance).

Methods

This data were analyzed with a time series design using generalized additive model.

Results

In the season-specific analysis, stronger associations were observed in the fall season for PM10 in Seoul, fall seasons for SO2, NO2 in Shanghai. In two-pollutant models, the associations between hospital admission for asthma and PM10 were strong and statistically significant in Seoul regardless of adjustment. Also, the associations between hospital admission for asthma and SO2, NO2 were strong after adjustment for PM10 in Shanghai. Lastly this study shows that there is a strong relationship on PM10 in Seoul, SO2 and NO2 in shanghai with hospital admission for asthma.

Conclusion

This study shows that outdoor air pollution was associated with increased risk of asthma hospitalization in Seoul and Shanghai. In China especially Shanghai, Gaseous air pollutants (SO2, NO2) seemed to have larger effects on asthma hospitalization compared to PM10. In contrast, Seoul, Korea, PM10 had much higher risk on asthma hospitalization than that of SO2 and NO2. These results suggested that associations of air pollution on asthma exacerbations should be considered in asthma related policy decisions and resource distributions of public health.

Keywords:

Asthma, Seoul, Shanghai, Pollution, Hospitalization, Generalized Additive model

Introduction

세계보건기구(WHO,2007)는 세계적으로 대기오염으로 인한 사망자수가 약3백만명이고, 전체사망자(55백만명)의 약 5%를 차지하는 것으로 추정하고 있다. 한편 대기오염과 그에 따른 대기 오염 물질의 증가에 따라 인구집단에게 미치는 영향은 크며 다수의 연구를 통해 알려져 왔다(Chen et al., 2004). 1952년 12월 5일~12월 9일까지의 런던스모그가 대기오염에 따른 사망의 영향을 보여준 대표적인 사건으로 사망자수가 평상시 보다 약 2.6배인 4000명 이상이 사망 했으며, 1953년 2월까지 8000여명의 사망자를 보여주었다. 그리고 대기오염이 건강에 미치는 영향은 연령과 성별, 주거지, 기상상황, 경제수준 등에 따라 다양하게 나타난다(Kan et al.,2008). 또한 특정 그룹에 있어서는 그 영향이 더욱 두드러지게 나타나며(WHO,2004), 현재 호흡기 질환을 가지고 있을 경우에는 그 영향에 더욱 취약하다(Goldberg et al.,2001).

특히, 호흡기 질환 중 전세계적으로 3억명에 가까운 사람이 천식으로 인한 영향을 받고 있는 것으로 나타났으며(Masoli et al.,2004), 사회 경제적 부담이 급증하고 있으며, 산업화되고 선진화된 지역에서 특히 유병률이 높은 것을 알 수 있다(To et al.,2012). 그리고 대기오염의 위험이 높은 중국에서는 천식증상을 겪는 전체인구가 3,000만명에 이르며(Chen et al.,2013), 한국에서는 성인의 질병부담 3위의 질환이다. 대기오염에 따른 천식의 영향을 정량적으로 분석한 기존연구에서는 건강영향(Health Outcome)으로 천식악화, 천식으로 인한 병원 입원 등을 통해 알아보고 있다.

다른 나라의 선행 연구를 살펴볼 경우 일본의 도쿄에서는 응급실 환자 중 1% 정도가 천식악화로 입원한 경우였으며, 낮은 기온이 이에 추가적인 영향을 보이는 것을 알 수 있었으며(Abe et al., 2009), 미국에서는 메디케이드(Medicaid) 자료를 이용해 살펴 보았을 때 PM10 (particulate matter less than 10 mm in aerodynamic diameter)과 SO2 (sulfur dioxide) 천식악화에 영향을 미쳤고(Jaffe et al.,2003), 유럽에서는 성인에서는 NO2 (nitrogen dioxide) 소아에서는 SO2 가 천식입원에 영향을 주었으며(Sunyer et al.,2000), 대만의 자료에서는 SO2, NO2, CO 성분이 천식 입원에 영향을 미친 것으로 나타났다(Yang et al.,2007).

본 연구에는 대기오염 물질에 따른 천식의 입원에 미치는 건강영향을 알아보는 것이 주 목적이며, 국가간의 비교를 통해 대기오염에 대한 나라별 급성 천식의 취약성을 살펴보는 것과 대기오염 물질의 지연효과(lag effect)를 통해 천식의 악화의 시계열적 고찰과, 계절에 따른 영향의 차이를 살펴보는 것으로서도 그 중요성이 있다.


Methods

1.Data

본 연구의 연구대상은 한국의 수도(Capital)이자, 가장 큰 대도시인 서울(Seoul), 중국의 상업의 중심지이며 대표적인 대도시로 상하이(Shanghai)를 분석지역으로 한다.

각 도시 별로 수집된 자료의 기간은 2001년 1월1일부터 2006년 12월31일까지(서울), 2005년 1월1일부터 2011년 12월31일까지(상하이)이다.

또한 각 지역별로 자료의 출처는 아래와 같다.

한국(서울)의 데이터는 한국 기상청의 기상자료(Korea Meteorological office), 한국보건환경연구소의 대기오염자료 (Research Institute of Public Health and Environment in Korea), 국민건강보험의 청구자료(National Health Insurance Claim, Seoul, 2001-2006년)의 자료를 활용할 것이며, 중국의 데이터는 중국기상자료(China Meteorological Data Sharing Service System), 상하이 환경모니터링 센터의 대기오염자료(SEMC-Shanghai Environmental Monitoring Center),상하이건강보험자료 (Shanghai Health Insurance, nine urban districts area in shanghai, 2005-2011년)를 활용하여 분석을 하였다.

2개 이상의 측정소가 있는 경우에는 산술평균을 내어 사용하였다.

Comparisons of the study areas

2. Statistical Analysis

일별 천식 입원환자수(Daily Asthma Hospitalization), 대기오염자료(Air pollution data), 기상자료(Meteorological data)는 날짜별로 연결할 것이며, 시계열 분석(Time-series analysis)이 필요하다. 그에 따라 본 연구에서는 GAM(Generalized Additive Model : 일반화 가법 모형) 모형을 사용한다.

GAM 모형은 비선형 함수(Non-linear function)를 포함한 회귀 모형(Regression model)을 의미하며 적합하고자 하는 각 변수 별 함수들의 합의 모양인 형식으로 이해를 해 볼 수 있다. Natural Cubic Splines 함수는 국소화된 여러 개의 3차 함수가 연결된 함수로, 시계열적 변화를 가진 대상을 설명하기 위한 함수로 널리 이용되고,(Peng, Dominici et al. 2006). 비선형적인 관계를 시계열 분석에서 반영하기 위하여 NS(Natural Cubic Splines) 를 사용한다.

다만 포아송 분포의 경우 평균과 분산이 같은 특성이 존재하며, 분석 시 분산이 평균보다 커지는 과산포(Overdispersion)의 문제가 있을 수 있다. 따라서 이를 보완하기 위해 실제 분석에서는 과산포을 가정한 Quasi-Possion 분포를 사용한다. (Dominici et al., 2004)

모델 설정은 (Bell, Ebisu et al. 2008)을 참고하였다.

연구모형은 아래와 같다.

logEYtc=α+Pollutanttc+nshumid,df+nstemp,df+factordow+nssn,df

Variables and Characters of the model

또한 공변량(Covariates) 으로서 잠재적인 비선형적(non-linear)관계,교호작용(interaction)의 보정을 하기 위해 기상조건 중 일 평균 온도(℃,daily mean temperature)는 자유도(degree of freedom)6으로, 일 평균 습도(%, daily mean humidity)는 자유도3으로 보정하였다. 또한 요일을 보정(Day of Week)하고, 측정되지 않은 장기적인 시계열 보정을 1년당 7의 자유도를 반영하여 보완 하였다. 한편 분석 시 대기오염 변수(PM10, SO2, NO2)에 대해서는 지연효과(Lag0~5:단일지연값,Lag01:이동평균값)를 사용하여 분석하였다. Lag(지연효과)란 관심 변수의 관측값을 일정기간 앞으로 당기는 것으로 이해할 수 있다. 지연효과를 사용하는 주된 이유는 대기오염 노출에 대한 건강 상태의 변화에 일정 시간이 필요하기 때문인 것으로 이해할 수 있다. 본 결과 기술 시에는 대기오염의 단일지연효과(lag 0~5)와 2일 이동평균값(lag01) 중 가장 사망에 큰 영향을 준 지연효과를 중심으로 살펴보기로 한다. 그리고 천식의 입원에 미치는 영향이 계절별로 어떤 영향을 미치는가와 대기오염 변수간의 영향을 보기 위해서 교호 작용을 고려하여 분석을 실행하며, 이는 계절에 따른 건강영향의 이질성과 대기오염 성분의 상호 관련성을 보여줄 것이다.

모든 분석은 R version 3.02 와 SAS 9.4를 이용하였다.

[Figure1]

Smoothing plots of asthma hospital admission in Seoul(2001-2005) and Shanghai(2006- 2011)


Results

1. Descriptive statistics

[Table3],[Table4]에서는 각각 서울(Seoul)과 상하이(Shanghai)의 기술통계량(Descriptive statistics)을 보여주고 있다. 서울과 상하이의 천식입원횟수는 총 기간(서울 5년 : 2001-2006년, 상하이 6년 : 2005-2011년) 동안 45,541회와 15,678회로 각각 나타났다. 상하이의 천식 입원횟수가 적은 이유로는 상하이의 의료보험자료는 상하이의 9개의 시내중심지역(Urbanized area, 약 7백만명)에 국한되어 다소 실제와 차이를 보이는 것으로 나타났다.

Summary statistics of daily asthmatic hospitalization admission (Seoul N=45541),air pollutant concentrations and weather conditions in Seoul(2001-2006)

Summary statistics of daily asthmatic hospitalization admission (Shanghai N=15,678), air pollutant concentrations and weather conditions in Shanghai(2005-2011)

서울(한국)은 일 평균 20.8회의 천식에 따른 입원이 나타났으며, 천식에 따른 입원횟수는 여름(17.4회)이 다른 계절에 비하여 적게 나타났고, 봄(24.회)이 가장 많았으며, 봄,겨울,가을,여름 순인것으로 나타났다. 또한 연구기간 동안 평균적인 PM10 농도는 65㎍/㎥ 이었으며, SO2 와 NO2는 각각 13㎍/㎥와 68㎍/㎥인 것으로 분석되었다. 일 평균온도는 12℃, 상대습도는 61%으로 나타났다. 그리고 성별과 연령을 구분 지어 살펴볼 경우 남성이 천식에 따른 입원을 하는 경우는 일 평균 9.8회, 여성은 11회를 차지했으며, 연령별로 볼 경우 큰 차이는 없지만 20-65세는 10.5회, 65세이상은 10.3회 인 것으로 나타났다.

상하이(중국)는 일 평균 5.5회의 천식 입원횟수를 보였고, 한국과 마찬가지로 여름(4.8회)이 가장 적은 천식입원횟수를 보였고, 겨울(6.8회)이 가장 많았다. 그 순서는 계절별로 살펴볼 경우 겨울,가을,봄,여름 순서로 나타났다. PM10농도는 87㎍/㎥, SO2 농도는 45㎍/㎥ 그리고 NO2농도는 59㎍/㎥으로 나타났고, 일 평균온도와 상대습도는 17℃, 69%로 각각 나타났다. 그리고 성별과 연령을 구분 지어 살펴볼 경우 남성이 천식에 따른 입원을 하는 경우는 일 평균 3.0회, 여성도 3.0회를 차지했으며, 연령별로 볼 경우 큰 차이는 없지만 18-65세는 3.0회, 65세이상은 2.0회 인 것으로 나타났다.

한편 서울과 상하이의 대기오염물질의 농도를 비교하기 위해 24시간 평균농도를 비교해 보면 [Figure2]에서 보는 것처럼 PM10 농도의 경우에는 서울(65.5㎍/㎥), 상하이(87.9㎍/㎥)로 상하이가 높고, SO2 의 농도는 서울(13.2㎍/㎥), 상하이(45㎍/㎥)로 상하이가 서울에 비해 3배보다 더 높음을 알 수 있다. 한편 NO2의 농도는 서울(68.4㎍/㎥), 상하이(59.8㎍/㎥)로 서울이 조금 더 높은 농도를 나타낸다. 기상조건[Figure3]에서는,평균온도는 17.5℃(상하이), 12.9℃(서울)로 상하이가 높으며, 상대습도도

69.5%(상하이), 61.8%(서울)로 높아, 두 기상조건 모두 상하이가 높게 나타난다.

[Table3],[Table4]에서는 각각 서울(Seoul)과 상하이(Shanghai)의 기술통계량(Descriptive statistics)을 보여주고 있다. 서울과 상하이의 천식입원횟수는 총 기간(서울 5년 : 2001-2006년, 상하이 6년 : 2005-2011년) 동안 45,541회와 15,678회로 각각 나타났다. 상하이의 천식 입원횟수가 적은 이유로는 상하이의 의료보험자료는 상하이의 9개의 시내중심지역(Urbanized area, 약 7백만명)에 국한되어 다소 실제와 차이를 보이는 것으로 나타났다.

[Figure2]

Comparisons of Air pollution concentrations (24-h average) in Seoul(2001-2005) and Shanghai(2006-2011)

[Figure3]

Comparisons of Meteorological Measures (24-h average) in Seoul(2001-2005) and Shanghai(2006-2011)

2. Correlation coefficients

서울(한국)와 상하이(중국) 두 지역 모두에서 세 가지 대기오염물질(PM10, SO2, NO2)과 일별 천식에 의한 입원 수와 높은 상관관계를 보였다[Table5]. 서울의 경우 천식입원과 PM10은 0.252, SO2 는 0.140, NO2는 0.183으로 나타났으며, 상하이의 경우는 각각 순서대로 0.0176, 0.0396, 0.1013으로 값이 나왔다.

한편 온도와 상대습도는 일별 천식에 의한 입원 수와 상관 관계가 없는 것으로 보이며 오히려 역의 상관관계가 관찰된다. 대기오염 물질과 일 평균 온도 및 상대습도도 특별한 상관성이 없다. 하지만 상대습도와 일 평균 온도는 일반적인 사실과 마찬가지로 서로간에 상관 관계가 높게 나타난다. 서울의 경우 0.406, 상하이의 경우 0.161의 값을 가진다.

대기오염물질의 피어슨 상관계수 (Pearson correlation coefficients)를 자세히 살펴보면 서울의 경우에는 PM10과 SO2 가 0.367, PM10과 NO2가 0.385, SO2 와 NO2가 0.753으로 상대적으로 SO2 와 NO2가 높은 상관관계를 보였다. 상하이에서는 PM10과 SO2 는 0.627, PM10과 NO2는 0.660, SO2 와 NO2는 0.703으로 서울에 비해 상하이는 대기오염물질이 서로 높은 상관성을 보이는 것으로 나타났다.

Pearson correlation coefficients between daily air pollutant concentrations and weather conditions in Seoul(2001-2006) and Shanghai(2005-2011)

3. The pollutants concentrations against asthma hospital admission increase

[Figure4]는 대기오염물질 중 PM10과 천식에 따른 병원입원의 상대위험도(Relative Risk)를 나타낸 것이며, [Figure5]는 SO2 와 천식에 따른 병원입원의 상대위험도를 나타낸 것 이고, [Figure6]는 NO2 와 천식에 따른 병원입원의 상대위험도를 나타낸 것이다. 세 그림 모두, 실선은 일별 천식으로 인한 입원의 상대위험도의 추정값을 나타내며, 점선은 95% 신뢰구간을 나타내준다.

우선, [Figure4]의 PM10 농도와 천식에 따른 입원의 관계를 볼 경우 PM10의 농도가 증가할 경우, 천식에 대한 상대위험도가 서울에서는 거의 선형적으로 증가하지만, 상하이에서는 이런 관계를 보이지 못하고, 오히려 고농도에서 떨어지는 형태를 보이고 있다.

[Figure4]

Smoothing plots of PM10 against asthma hospital admission. X-axis is the pollutants concentrations (㎍/㎥) (2day moving average lag, Lag01) (Left:Seoul, Right:Shanghai)

다음으로 [Figure5]에서는 대기오염물질 중 SO2 의 농도와 천식입원의 상대위험의 추정 값을 살펴 볼 경우, 서울에서는 일정한 관계를 보여주지 못하고 있지만, 상하이에서는 농도가 증가함에 따라 천식의 입원에 대한 상대위험도가 급격히 증가하는 형태를 보여준다.

[Figure5]

Smoothing plots of SO2 against asthma hospital admission.

마지막으로 [Figure6]는 NO2 농도와 천식입원에 대한 상대위험도의 형태를 나타낸다. SO2 의 경우와 유사한 그림을 나타내며, 서울의 경우 NO2의 농도와 천식에 따른 입원의 관계는 보이지 않지만, 상하이의 경우 선형적으로 증가하는 형태를 보인다.

[Figure6]

Smoothing plots of NO2 against asthma hospital admission.

4. The pollutants concentrations against asthma hospital admission increase by lag effect

[Table6], [Table7], [Table8]는 단일지연효과(Lag0-Lag5)와 이동평균법(lag01:2-day moving average lag)의 대기오염물질(PM10, SO2, NO2) 노출 농도의 10단위(10㎍/㎥) 증가에 따른, 천식 입원의 증가를 서울과 상하이에서 각각 확인해 본 것 이다.

통계적으로 유의한 관계는 PM10의 경우([Table6])에는 서울에서는 lag0(0.540, 95% CI: 0.333, 0.747), lag01(0.572, 95%CI: 0.312, 0.832), lag1(0.248, 95% CI : 0.016, 0.479)에서 나타났으며 이는 PM10 이 10㎍/㎥ 증가함에 따라, 0.57% 천식 입원이 증가함을 나타낸다(lag01기준). 상하이에서는 lag5(0.324, 95% CI : 0.044, 0.694) 에서만 유의한 값을 보였으며, PM10 이 10㎍/㎥ 증가함에 따라, 0.32% 천식 입원이 증가한다.

Percent change of risk (mean and 95% confidence interval) in daily asthma hospital admission per PM10 10 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

SO2 의 경우([Table7])는 PM10의 경우와 다르게 서울에서는 lag4(3.521, 95% CI : 0.661, 6.461)에서 유의하였으며, SO2 가 10㎍/㎥ 증가함에 따라, 3.52% 천식 입원이 증가한다. 상하이에서는 lag0(1.347, 95%CI: 0.422, 2.280), lag01(1.499, 95%CI: 0.447, 2.562), lag1(0.881, 95% CI : 0, 1.770) 에서 통계적으로 유의한 결과가 있었으며 SO2 이 10㎍/㎥ 증가함에 따라, 1.49% 천식 입원이 증가함을 의미한다.

Percent change of risk in daily asthma hospital admission per SO2 10 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

NO2의 경우([Table8])는 SO2 와 유사한 결과를 보인다. 서울은 lag4(0.603, 95% CI : 0.014, 1.224)로서, NO2 이 10㎍/㎥ 증가함에 따라, 0.6% 천식 입원이 증가함을 의미한다. 상하이는 lag0(1.772, 95% CI : 0.733, 2.822), lag01(2.462, 95%CI : 1.240, 3.697), lag1(1.871, 95%CI : 0.817, 2.937)에서 유의하며, NO2 이 10㎍/㎥ 증가함에 따라, 2.46% 천식 입원이 증가한다.

Percent change of risk in daily asthma hospital admission per NO2 10 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

[Figure7]

Relative risk (mean and 95% confidence interval) in daily asthma hospital admission per PM10 , SO2 ,NO210 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011)

3가지 대기오염물질의 지연효과를 고려하고, 주로 가장 큰 결과값을 보인 lag01(2-day moving average lag) 를 사용하여, 다른 분석에서도 사용하였다. 이는 천식의 입원에 있어서 급성(immediate, acute)효과가 장기 노출(long-term exposure)에 비해 그 영향이 크다는 것을 말해주고 있다.

5. The pollutants concentrations against asthma hospital admission increase by season effect

[Table9]은 각각의 대기오염물질에 따른 계절별 효과를 추정한것이다. 이동평균법(2-day moving average lag:lag01)을 적용하여 대기오염물질(PM10, SO2, NO2) 노출 농도의 10단위(10㎍/㎥) 증가에 따른, 천식 입원의 퍼센트 증가를 계절별로 서울과 상하이에서 각각 확인해 본 것 이다.

Percent change of risk (mean and 95% confidence interval) in daily asthma hospital admission per 10 ㎍/㎥ change of exposure in air pollutants by season in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

[Figure8]

Relative risk in daily asthma hospital admission per PM10, SO2, NO2 10 ㎍/㎥ change of exposure by season in Seoul(2001-2005) and Shanghai(2006-2011) (Left:Seoul, Right:Shanghai)

각 계절은 봄(3월-5월), 여름(6월-8월), 가을(9월-11월), 겨울(12월-2월)으로 서울과 상하이 두 지역 모두 구분하였다. 계절별 분석을 하였을 경우, 서울에서는 대기오염물질 중 PM10이 가을에 거의 유의한 값(0.702, 95% CI : -0.227, 1.639)을 나타냈고 PM10 이 10㎍/㎥ 증가함에 따라, 0.7% 천식 입원이 증가함을 의미한다. 상하이에서는 SO2 (2.279, 95% CI : 0.169, 4.434)와 NO2 (1.993, 95%CI : 0.139, 4.171)에 있어서 다른 계절에 비해 가을에 차이가 크고, 유의한 것으로 나타났으며 이 결과는 SO2 이 10㎍/㎥ 증가함에 따라, 2.27% 천식 입원이 증가하고, NO2 이 10㎍/㎥ 증가함에 따라, 1.99% 천식 입원이 증가함을 보여준다.

또한 두 지역 모두 SO2 와 NO2는 계절별로 천식으로 인한 입원에 대한 차이가 나타났지만, PM10은 계절에 따라 차이가 두드러지지 않았다.

6. The pollutants concentrations against asthma hospital admission increase by two pollutant mixed model

[Table10]은 두 대기오염물질(PM10+ SO2, PM10+ NO2, SO2 + NO2, SO2 + PM10, NO2+ PM10, NO2+ SO2:2-day moving average lag)을 동시에 고려하여 혼합모델을 만들어 천식에 대한 입원과의 영향을 살펴본 것이다. 우선 PM10 의 경우는 PM10 단독모델, PM10을 SO2 로 보정, PM10 을 NO2으로 보정해서 세가지로 살펴보았으며, SO2의 경우는 SO2 단독모델, SO2 을 NO2 로 보정, SO2 를 PM10 으로 보정하여 살펴보았다. 마지막으로 NO2 도 다른 것과 마찬가지로 NO2 단독으로 한 모델, NO2 를 PM10으로 보정한 모델, NO2을 SO2 로 보정한 모델을 살펴보았다.

Percent change of risk (mean and 95% confidence interval) in daily asthma hospital admission per 10 ㎍/㎥ change of exposure in air pollutants with single and two-pollutant models in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

서울의 경우에는 천식에 따른 입원과 대기오염 물질 중 PM10 의 경우에는 단독모델(0.572, 95% CI : 0.312, 0.832)에서도, PM10 이 10㎍/㎥ 증가함에 따라, 0.57%의 천식증가를 보여줬고, 이는 혼합모델인 SO2 (0.625, 95%CI : 0.363, 0.888)와 NO2 (0.597, 95% CI 0.335, 0.859)의 보정에도 큰 차이를 보이지 않았다. 한편 SO2 는 단독모델(-1.713, 95% CI : -5.391, 2.107)와 혼합모델(PM10으로 보정: -3.657, NO2으로 보정:-2.486)에서 유의한 값을 보이지 않았다. NO2 에서는 단일모델(-0.147, 95% CI : -1.015, 0.728)에서는 유의한 값을 나타내지 못했고, 혼합모델 중 SO2 로 보정한 경우 (0.257, 95% CI :-0.964, 1.494) 만 다소 유의미한 값을 나타냈고, PM10으로 보정한 경우는 (-0.500, 95% CI: -1.373, 0.380) 통계적으로 유의한 값을 보여주지 않았다.

상하이의 경우에는 SO2의 단일모델 (1.499, 95% CI : 0.447, 2.562) 은 통계적으로 유의한 값을 나타냈지만, PM10으로 보정하였을 경우(2.519, 95% CI : 1.122, 3.935) 더 큰 영향을 보여주였고, 이는 SO2 가 10㎍/㎥ 증가함에 따라, 2.51% 천식 입원이 증가함을 보여준다. 또한 이 결과는 PM10의 효과를 배제한 SO2 자체의 효과를 의미하며 보정 전에 비해 1.66배정도 높은 수치이다. NO2 단일모델의 경우(2.462, 95% CI : 1.240, 3.697)가 통계적으로 유의한 값을 나타냈으며, NO2가 10㎍/㎥ 증가함에 따라, 2.46% 천식 입원이 증가함을 의미한다. 한편 NO2를 PM10으로 보정한 혼합모델(4.395, 95% CI : 2.702, 6.116)의 경우 영향이 커졌고, 이 결과는 PM10의 효과를 배제한 NO2 자체의 효과를 의미하며 보정 전에 비해 1.76배정도 높은 수치로서, NO2가 10㎍/㎥ 증가함에 따라, 4.4% 천식 입원이 증가함을 의미한다. 하지만 SO2 로 보정한 경우 단일모델과 차이가 거의 나타나지 않았다(2.582, 95% CI : 0.780 ,4.415). 한편 PM10은 단독모델(0.013, 95% CI :0.425, 0.453)과 혼합모델(SO2 로 보정: -0.659, NO2으로 보정:-1.035)에서 유의한 값을 보이지 않았다.


Conclusion

본 연구는 대기 오염에 따른 천식 입원의 증가를 서울(Seoul)과 상하이(Shanghai)에서 살펴본 것이다. 상하이에서 대기오염물질 중 가스 성 물질(SO2, NO2)이 PM10보다 상대적으로 천식의 입원에 더 큰 위험 요인인 것으로 나타났다. 이와 반대로 서울에서는 PM10이 천식의 입원에 SO2와 NO2보다 더 큰 영향을 미치는 것으로 보여진다. 이러한 결과는 각 지역에서 나타나는 각 대기오염 물질의 농도로 해석될 수 있다. 앞서 설명 되어진 것처럼 SO2의 농도는 상하이가 서울에 비해 3배정도 더 많으며, 서울의 PM10농도는 다른 물질(SO2, NO2)에 큰 영향을 받지않고 상대적으로 고농도이다. 이와 더불어 본 연구자료의 NO2는 서울이 상하이 보다 약간 고농도이지만 2015년 NASA의 NO2 대기오염농도에서 서울과 상하이는 공동 세계4위로 높은 대기오염을 나타냈고, 다른 연구자료도 마찬가지로 서울이 높은 자동차보급률과 함께 NO2 가 높은 것은 합리적인 것으로 보였다(Kang TK et al., 2013).

한편 대기오염물질의 농도에 따른 경보시스템도 이에 영향을 미칠 수 있다. 대기오염물질의 급증에 따른 경보를 발령하는 중국의 상하이의 경우 실외 활동을 자제하게 되고, 가시적이며 민감한 대기오염물질인 고농도의 실외 PM10에는 노출이 실제로는 적고 오히려 지속적으로 노출된 SO2와 NO2가 큰 영향을 줄 수 있다. 이에 반해 한국에서는 대기오염예보가 있지만 중국의 경우에 비해서 외부활동을 제한하는 정도가 덜 하다고 보여지며, SO2와 NO2에 비해 상대적으로 높은 PM10 농도의 단기간 노출이 천식의 악화를 초래하고 폐기능을 떨어뜨릴 수 있다고 볼 수 있다(Kan, H. et al., 2014).

그리고 천식의 입원에 미치는 건강 영향은 지연 효과(lag effect)로 해석해 볼 수 있다. 서울에서는 PM10의 경우에는 lag0, lag1, lag01(이동평균법)가 유의하게 나타났으며, SO2와 NO2는 lag4가 통계적으로 유의한 값을 보여주었다. 상하이에서는 PM10의 경우는 lag5, SO2와 NO2는 lag0, lag1, lag01 에서 유의한 결과치를 나타냈다. 가장 두드러지게 관련성을 살펴볼 수 있었던 지연효과는 lag01로서 다른 모델에도 lag01를 적용하여 결과값을 살펴보았다. 다시 말하자면 대기오염물질(PM10, SO2, NO2)이 천식에 미치는 영향이 주로 짧은 지연효과를 나타냈고 이는, 급성적이고 즉각적인 대기오염물질의 농도증가에 따른 것임을 암시하고 있다.

또한 본 연구에서 천식의 입원의 증가에 미치는 중요한 요소로 계절도 들 수 있었다. 계절별 영향에서 천식은 다른 유사한 연구들에서 가을의 선선한 날씨와 일교차 및 환절기에 따른 천식 환자가 급증 하는 것과 마찬가지로 서울과 상하이 두 지역 모두에서 따뜻한 계절(여름)에 비하여 차가운 계절(가을, 겨울), 특히 가을의 대기오염과 천식의 입원이 강한 영향을 보여주었다. 이는 차가운 공기가 기도를 왕복하여, 점막이 건조해지면서 호흡기 보호층이 약화되어 기관지를 자극할 뿐만 아니라, 환절기가 되어 갑자기 떨어진 기온에 몸이 적응을 하지 못하고 면역력이 떨어져 생긴다고 말 할 수 있으며, 가을철 식물의 씨앗과 꽃가루 등에 의한 알러지 반응으로 해석 될 수도 있다(Kim SH et al, 2015).

본 연구의 제한 점은 다음과 같다. 우선 연구가 개인의 특성을 고려하기 어려운 생태적 디자인으로 설계되었기에, 천식의 급성악화가 대기오염이라는 원인과 결과를 직접적으로 유추하기에는 한계점을 지니고 있지만, 관련성을 추정할 수 있다. 또한 다수의 대기오염물질의 측정장소를 단순 평균하여 사용하였으므로, 대기오염의 개인별 노출의 정도가 다르다는 것을 고려하지 못한다. 그리고 다수의 연구결과에서 고려된 대기오염물질인 O3 및 PM2.5 는 중국 자료 수집의 한계로 인해서 본 연구에서 함께 분석되지 못했다.

Acknowledgments

이 논문은 서울대학교 국제협력본부 및 보건환경연구소의 지원을 통해 수행된 연구임

References

  • Abe, T., et al , The relationship of short-term air pollution and weather to ED visits for asthma in Japan, The American journal of emergency medicine, (2009), 27, p153-9. [https://doi.org/10.1016/j.ajem.2008.01.013]
  • Bell, ML., Ebisu, K., Peng, RD., et al , Seasonal and regional short-term effects of fine particles on hospital admissions in 202 US counties, 1999–2005, American journal of epidemiology, (2008), 168, p1301-10. [https://doi.org/10.1093/aje/kwn252]
  • Cai, J., Zhao, A., Zhao, J., et al , Acute effects of air pollution on asthma hospitalization in Shanghai, China, Environmental pollution, (2014), 191, p139-44. [https://doi.org/10.1016/j.envpol.2014.04.028]
  • Chen, B., Hong, C., Kan, H., Exposures and health outcomes from outdoor air pollutants in China, Toxicology, (2004), 198, p291-300. [https://doi.org/10.1016/j.tox.2004.02.005]
  • Chen, ZH., Wang, PL., Shen, HH., Asthma research in China: A five‐year review, Respirology, (2013), 18, p10-9. [https://doi.org/10.1111/resp.12196]
  • Dominici, F., Peng, RD., Bell, ML., et al , Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases, Jama, (2006), 295, p1127-34. [https://doi.org/10.1001/jama.295.10.1127]
  • Goldberg, MS., Burnett, RT., Bailar, J. 3rd, et al , Identification of persons with cardiorespiratory conditions who are at risk of dying from the acute effects of ambient air particles, Environmental health perspectives, (2001), 109, p487..
  • Jaffe, DH., Singer, ME., Rimm, AA., Air pollution and emergency department visits for asthma among Ohio Medicaid recipients, 1991–1996, Environmental Research, (2003), 91, p21-8. [https://doi.org/10.1016/S0013-9351(02)00004-X]
  • Kan, H., et al , Season, sex, age, and education as modifiers of the effects of outdoor air pollution on daily mortality in Shanghai, China, Environ Health Perspect, (2008), 116, p1183-8. [https://doi.org/10.1289/ehp.10851]
  • Kim, SH., et al , Effects of particulate matter in ambient air on the development and control of asthma, Allergy Asthma Respir Dis, (2015, Sep), 3(5), p313-319. [https://doi.org/10.4168/aard.2015.3.5.313]
  • Masoli, M., Fabian, D., Holt, S., Beasley, R., The global burden of asthma: executive summary of the GINA Dissemination Committee report, Allergy, (2004), 59, p469-78. [https://doi.org/10.1111/j.1398-9995.2004.00526.x]
  • McCreanor, J., Cullinan, P., Nieuwenhuijsen, MJ., et al , Respiratory effects of exposure to diesel traffic in persons with asthma, New England Journal of Medicine, (2007), 357, p2348-58. [https://doi.org/10.1056/NEJMoa071535]
  • Peng, RD., Dominici, F., Louis, TA., Model choice in time series studies of air pollution and mortality, Journal of the Royal Statistical Society: Series A (Statistics in Society), (2006), 169, p179-203. [https://doi.org/10.1111/j.1467-985X.2006.00410.x]
  • Sunyer, J., Spix, C., Quenel, P., et al , Urban air pollution and emergency admissions for asthma in four European cities: the APHEA Project, Thorax, (1997), 52, p760-5. [https://doi.org/10.1136/thx.52.9.760]
  • To, T., Stanojevic, S., Moores, G., et al , Global asthma prevalence in adults: findings from the cross-sectional world health survey, BMC public health, (2012), 12, p204. [https://doi.org/10.1186/1471-2458-12-204]
  • World Health Organization, (2004), Health aspects of air pollution: results from the WHO project" Systematic review of health aspects of air pollution in Europe.
  • Yang, C-Y., et al , Air pollution and hospital admissions for asthma in a subtropical city: Taipei, Taiwan, Journal of Toxicology and Environmental Health, Part A, (2007), 70, p111-7. [https://doi.org/10.1080/15287390600755059]
  • 강택구, 외 , (2013), 한·중대기오염 저감 관리 비교와 협력방안 보고서, KIEP 중국종합연구, p13-56.

[Figure1]

[Figure1]
Smoothing plots of asthma hospital admission in Seoul(2001-2005) and Shanghai(2006- 2011)

[Figure2]

[Figure2]
Comparisons of Air pollution concentrations (24-h average) in Seoul(2001-2005) and Shanghai(2006-2011)

[Figure3]

[Figure3]
Comparisons of Meteorological Measures (24-h average) in Seoul(2001-2005) and Shanghai(2006-2011)

[Figure4]

[Figure4]
Smoothing plots of PM10 against asthma hospital admission. X-axis is the pollutants concentrations (㎍/㎥) (2day moving average lag, Lag01) (Left:Seoul, Right:Shanghai)

[Figure5]

[Figure5]
Smoothing plots of SO2 against asthma hospital admission.

[Figure6]

[Figure6]
Smoothing plots of NO2 against asthma hospital admission.

[Figure7]

[Figure7]
Relative risk (mean and 95% confidence interval) in daily asthma hospital admission per PM10 , SO2 ,NO210 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011)

[Figure8]

[Figure8]
Relative risk in daily asthma hospital admission per PM10, SO2, NO2 10 ㎍/㎥ change of exposure by season in Seoul(2001-2005) and Shanghai(2006-2011) (Left:Seoul, Right:Shanghai)

[Table1]

Comparisons of the study areas

Seoul Shanghai
Study period January 1, 2001- December 31, 2006 (2191 days) January 1, 2005 - December 31,2011 (2922 days)
Total Asthma Hospitalization 45,541 15,678
Study area 605.18 km2 (Total area) 279 km2 (9 Urbanized area)
Total area 605.18 km2 6,340.5km²
Study population 10,297,004 Around 7,000,000
Total population 10,297,004 (2005) 13,602,600 (2005)
Latitude 37° 34′ 31° 14′
Longitude 126° 59′ 121° 29′

[Table2]

Variables and Characters of the model

Variable Character
Daily number of Asthma
hospitalization
Continuous
Asthma (J45): ICD-10
Seoul : ≥20years old
Shanghai : ≥18years old
PM10
NO2
SO2
Continuous unit: ㎍/㎥
2-day moving average : lag01
Day of week
factor(DOW)
Categorical
Temperature
ns(temp,df=6)
Continuous
Humidity
ns(humid,df=3)
Continuous
Timetrend
ns(sn,df=year*7)
Continuous
Seoul : 6years
Shanghai : 7years

[Table3]

Summary statistics of daily asthmatic hospitalization admission (Seoul N=45541),air pollutant concentrations and weather conditions in Seoul(2001-2006)

Seoul Mean SD Min P(25) P(50) P(75) Max
Daily Asthma hospitalization 20.8 11.1 2.0 13.0 18.0 26.0 84.0
Season Spring(Mar-May) 24.3 11.9 4.0 16.0 22.0 32.0 83.0
             Summer(Jun-Aug) 17.4 8.6 3.0 11.0 16.0 22.0 57.0
             Fall(Sep-Nov) 18.6 10.3 2.0 12.0 16.0 23.0 84.0
             Winter(Dec-Feb) 22.8 12.0 3.0 14.0 20.0 29.0 76.0
Sex   Male 9.8 5.7 0.0 6.0 9.0 13.0 41.0
           Female 11.0 6.5 0.0 6.0 10.0 14.0 58.0
Age  18-65 10.5 5.9 0.0 6.0 9.0 13.0 50.0
          ≥65 10.3 6.3 0.0 6.0 9.0 14.0 44.0
Air pollution Concentrations (24-h average)
   PM10 (㎍/㎥) 65.5 50.2 10.0 39.0 57.0 80.0 983.0
   SO2 (㎍/㎥) 13.2 5.6 5.1 8.9 11.8 16.2 43.2
   NO2(㎍/㎥) 68.4 23.0 20.9 50.9 66.0 84.3 170.1
Meteorological measures (24-h average)
   Temperature(℃) 12.9 10.1 -15.7 4.1 14.5 21.8 30.4
   Relative Humidity(%) 61.8 14.7 21.1 50.8 61.8 72.4 96.0

[Table4]

Summary statistics of daily asthmatic hospitalization admission (Shanghai N=15,678), air pollutant concentrations and weather conditions in Shanghai(2005-2011)

Shanghai Mean SD Min P(25) P(50) P(75) Max
Daily Asthma hospitalization 5.5 3.2 1.0 3.0 5.0 7.0 29.0
Season Spring(Mar-May) 5.0 3.0 1.0 3.0 5.0 7.0 18.0
            Summer(Jun-Aug) 4.8 2.6 1.0 3.0 5.0 6.0 15.0
            Fall(Sep-Nov) 5.8 3.0 1.0 3.0 6.0 8.0 18.0
            Winter(Dec-Feb) 6.5 3.9 1.0 4.0 6.0 9.0 29.0
Sex   Male 3.0 2.0 0.0 1.0 2.0 4.0 18.0
           Female 3.0 2.0 0.0 1.0 2.0 4.0 13.0
Age  20-65 3.0 3.0 0.0 2.0 3.0 5.0 20.0
          ≥65 2.0 2.0 0.0 1.0 2.0 3.0 12.0
Air pollution Concentrations (24-h average)
   PM10 (㎍/㎥) 87.9 55.9 7.0 50.0 74.4 110.0 643.0
   SO2 (㎍/㎥) 45.0 28.4 5.0 23.8 38.0 60.0 223.0
   NO2(㎍/㎥) 59.8 22.7 11.2 43.4 56.4 72.4 167.0
Meteorological measures (24-h average)
   Temperature(℃) 17.5 9.1 -3.4 9.6 18.7 25.2 35.7
   Relative Humidity(%) 69.5 12.3 23.0 62.0 70.0 79.0 95.0

[Table5]

Pearson correlation coefficients between daily air pollutant concentrations and weather conditions in Seoul(2001-2006) and Shanghai(2005-2011)

Seoul
SO2 NO2 Temperature Relative Humidity Asthma Hospitalization
PM10 0.367 0.385 -0.079 -0.073 0.252
SO2 0.753 -0.429 -0.229 0.140
NO2 -0.237 -0.187 0.183
Temperature 0.406 -0.194
Relative Humidity -0.126
Shanghai
SO2 NO2 Temperature Relative Humidity Asthma Hospitalization
PM10 0.627 0.660 -0.185 -0.285 0.0176
SO2 0.703 -0.304 -0.352 0.0396
NO2 -0.352 -0.180 0.1013
Temperature 0.161 -0.1799
Relative Humidity -0.0222

[Table6]

Percent change of risk (mean and 95% confidence interval) in daily asthma hospital admission per PM10 10 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

PM10 Seoul Shanghai
Mean 95% CI Mean 95% CI
lag0 0.540* 0.333 0.747 -0.040 -0.421 0.341
lag01 0.572* 0.312 0.832 0.013 -0.425 0.453
lag1 0.248* 0.016 0.479 0.056 -0.317 0.430
lag2 0.041 -0.193 0.275 -0.404 -0.779 -0.027
lag3 0.158 -0.068 0.385 -0.500 -0.875 -0.123
lag4 0.139 -0.080 0.359 -0.480 -0.853 -0.105
lag5 -0.046 -0.288 0.196 0.324* 0.044 0.694

[Table7]

Percent change of risk in daily asthma hospital admission per SO2 10 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

SO2 Seoul Shanghai
Mean 95% CI Mean 95% CI
lag0 -1.121 -4.218 2.077 1.347* 0.422 2.280
lag01 -1.713 -5.391 2.107 1.499* 0.447 2.562
lag1 -1.302 -4.398 1.894 0.881* 0.000 1.770
lag2 0.446 -2.472 3.451 0.031 -0.841 0.910
lag3 2.133 -0.739 5.089 -0.275 -1.141 0.599
lag4 3.521* 0.661 6.461 -0.284 -1.140 0.580
lag5 1.096 -1.676 3.945 0.790 -0.065 1.652

[Table8]

Percent change of risk in daily asthma hospital admission per NO2 10 ㎍/㎥ change of exposure in air pollutants in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

NO2 Seoul Shanghai
Mean 95% CI Mean 95% CI
lag0 0.182 -0.520 0.890 1.772* 0.733 2.822
lag01 -0.147 -1.015 0.728 2.462* 1.240 3.697
lag1 -0.410 -1.136 0.322 1.871* 0.817 2.937
lag2 0.062 -0.585 0.713 0.633 -0.404 1.681
lag3 0.427 -0.198 1.056 0.046 -0.978 1.080
lag4 0.603* 0.014 1.224 -0.038 -1.051 0.985
lag5 1.000 0.997 1.002 1.003 1.000 1.007

[Table9]

Percent change of risk (mean and 95% confidence interval) in daily asthma hospital admission per 10 ㎍/㎥ change of exposure in air pollutants by season in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

Seoul Shanhai
Mean 95% CI Mean 95% CI
PM10 spring 0.553 -1.341 2.485 -0.337 -2.342 1.709
summer 0.208 -2.176 2.650 0.068 -2.429 2.629
fall 0.702 -0.227 1.639 0.125 -0.749 1.006
winter 0.370 -1.804 2.592 0.260 -1.714 2.273
SO2 spring -2.089 -17.687 16.464 1.981 -3.083 7.309
summer -4.243 -24.029 20.695 0.742 -4.705 6.500
fall 6.926 -1.256 15.786 2.279* 0.169 4.434
winter -3.849 -19.002 14.138 1.210 -3.208 5.830
NO2 spring -0.769 -4.266 2.856 3.405 -2.044 9.157
summer 0.556 -3.553 4.841 -1.254 -6.568 4.364
fall 0.488 -1.113 2.114 1.993* 0.139 4.171
winter -0.589 -4.210 3.169 4.496 -0.511 9.755

[Table10]

Percent change of risk (mean and 95% confidence interval) in daily asthma hospital admission per 10 ㎍/㎥ change of exposure in air pollutants with single and two-pollutant models in Seoul(2001-2005) and Shanghai(2006-2011) (*:statistically significant)

Seoul Shanghai
Mean 95% CI Mean 95% CI
PM10 without adjustment 0.572* 0.312 0.832 0.013 -0.425 0.453
adjusted for SO2 0.625* 0.363 0.888 -0.659 -1.244 -0.070
adjusted for NO2 0.597* 0.335 0.859 -1.035 -1.644 -0.422
SO2 without adjustment -1.713 -5.391 2.107 1.499* 0.447 2.562
adjusted for PM10 -3.657 -7.322 0.153 2.519* 1.122 3.935
adjusted for NO2 -2.486 -7.573 2.881 -0.139 -1.666 1.413
NO2 without adjustment -0.147 -1.015 0.728 2.462* 1.240 3.697
adjusted for PM10 -0.500 -1.373 0.380 4.395* 2.702 6.116
adjusted for SO2 0.257 -0.964 1.494 2.582* 0.780 4.415