전국 시군구 PM_2.5 연평균 농도 추정 방법 제안: 전국규모 국가대기오염측정망 측정 자료가 없는 2001~2014년을 중심으로

1. 서 론

미세먼지가 건강에 미치는 악영향은 오랜 동안 다양한 지역에서 수행된 역학 연구를 통해서 보고되었다. 그 중에서도 공기역학직경이 2.5 μm 이하의 미세먼지(PM_2.5)의 단·장기 노출은 사망 및 질병 이환의 위험을 증가시켰다(US EPA, 2009; WHO, 2003; Burnett et al., 2000). 특히, 북미와 유럽을 중심으로 한 대 규모 코호트 연구에서는 PM_2.5의 장기 노출이 건강에 미치는 영향을 보고했으며, 최근 다른 지역 국가 및 전 세계를 연구지역으로 PM_2.5 장기 노출의 사망 영향을 추정하는 연구가 진행되었다 (Cohen et al., 2017; Di et al., 2017; Beelen et al., 2014; Hoek et al., 2013; Lim et al., 2012).

코호트를 기반으로 PM_2.5 장기 노출의 건강 영향을 추정하는 역학 연구를 수행하기 위해서는, 코호트 추적 기간의 개인별 노출 값을 산출하는 노출 평가가 필수적이다. PM_2.5의 노출 평가는 지상의 대기오염 측정소에서 측정한 농도값을 이용해왔다. 그런데, PM_2.5에 대한 대규모 지역의 지속적인 측정은 상대적으로 최근 시기에 시작되었다. 북미와 서유럽에서도 전국규모의 지속적인 측정은 각각 1990년대 후반과 2000년대 중반 이후 시작하였다(US EPA, 2009; WHO, 2004). 그러나 PM_2.5 장기 노출이 건강에 미치는 영향을 평가하기 위한 자료원이 되는 많은 코호트들은 PM_2.5의 측정시기보다 앞서 구축되었다. 따라서, 측정 농도값이 없는 과거 시기가 코호트 연구 기간에 포함되어 있어 노출 평가에 어려움이 있다(Kim et al., 2017c).

우리나라에서도 전국 단위로 PM_2.5 장기 노출과 건강 간의 연관성을 평가하는 연구가 수행되지 못한 가장 큰 이유는, 과거 PM_2.5 측정 자료가 존재하지 않기 때문이다. 국내에서는 1995년부터 국가 및 지자체에서 운영하는 일반대기오염 측정망의 대기오염 측정소를 통해서 공기역학직경 10 μm 이하의 미세먼지(PM₁₀) 농도 측정 자료를 수집하였고, 2001년부터의 자료를 공개하고 있다(MOE, 2011). 이렇게 2001년부터 공개된 PM₁₀을 자료를 이용해서 2002년부터 시작된 전국 단위의 코호트 자료와 함께 PM₁₀ 장기 노출의 사망 영향을 추정한 국내 연구가 최근에 수행되었다(Kim et al., 2017b). 반면 PM_2.5는 2011년 대기환경기준이 신설되고 2015년 1월 1일부터 적용되어 2015년부터 대기환경기준물질에 추가됨에 따라 전국규모측정값이 공개되기 시작하였다(MOE, 2017). 때문에, 2015년 이전에 구축된 코호트에 적용하기에 제한점을 가진다. 또한, PM₁₀과 비교하여 전국규모 측정소 개수가 상대적으로 매우 적다.

국가에서 운영하는 광범위한 대기오염 측정소의 최근 도입으로 인한 과거 PM_2.5 농도 측정값의 부재는 우리나라뿐만 아니라 유럽 등 다른 국가의 미세먼지 건강 영향 평가에서도 공통되는 제한점이다. 이러한 제한점을 극복하기 위해 장기 노출의 건강 영향을 추정하는 많은 역학 연구에서 과거 PM_2.5 농도 예측방법을 개발하였다. 이들 연구에서는 PM_2.5와 연관성이 높은 다른 대기오염 물질들을 이용하거나, 시정거리, 배출 자료, 확산 모형과 광화학반응 모형, 통계모형, 인공위성 측정값, 물질 간의 비율 등을 이용하여 시공간적으로 다양한 범위의 PM_2.5 농도를 예측하는 모형과 방법을 개발하였다(Kim et al., 2017c; Li et al., 2017; Brauer et al., 2016; Beelen et al., 2014; Hogrefe et al., 2009; Paciorek et al., 2009; Ozkaynak et al., 1985).

이러한 PM_2.5 예측 방법 중, 많은 연구에서는 PM_2.5 농도와 입자가 큰 미세먼지 농도의 비율을 이용했다. 즉, PM_2.5 측정값이 없는 과거시기에 측정된 입자가 큰 미세먼지 농도값에 PM_2.5 측정값이 있는 시기 동안 계산된 비율을 곱해서, 과거 시기의 PM_2.5을 예측했다. 비율을 이용하는 방법은 적용법이 단순하지만, 쉽게 사용이 가능하고 상대적으로 높은 예측력을 보였다. 미국의 한 연구에서는 전국규모 측정망이 설치된 1999년 이전의 PM_2.5 농도를 예측하기 위해, 1999년 이후의 PM₁₀과 총부유입자(total suspended particle: TSP)의 비율과 PM_2.5와 PM₁₀의 비율을 이용하였다.이 연구에서는 1999~2000년의 PM_2.5/PM₁₀과 1989~1997년의 PM_2.5/TSP 비율을 산출하고 두 비율을 곱하여 PM_2.5/TSP를 계산하였다. 이렇게 산출된 두 비율과 1999년 이전에 측정된 TSP나 PM₁₀ 농도값을 직접 곱하여 과거 PM_2.5 농도를 예측하였다(Lall et al., 2004). 또 다른 미국의 연구에서는 PM_2.5이 측정된 시기 동안 계산한 PM_2.5와 PM₁₀의 비율을 이용한 시공 간 모형을 통해서 PM_2.5 측정값이 없는 시기의 비율을 예측한 후, PM₁₀ 농도값에 곱하여 PM_2.5 농도를 추정하였다(Yanosky et al., 2009). 전 세계 단위의 노출예측 연구에서는 PM_2.5 측정값이 없는 지역에서 PM₁₀ 측정값과 PM_2.5와 PM₁₀의 비율을 이용하여 PM_2.5 농도값을 예측하였다(Brauer et al., 2016). 뿐만 아니라 홍콩의 환경보호부(Environmental Protection Department)는 PM_2.5 농도를 예측하기 위한 가이드라인에 PM_2.5와 PM₁₀의 비율을 사용하는 방법을 직접 제시하고 있다(HKSAR-EPA, 2016).

우리나라에서도 장기간 노출된 미세먼지의 건강영향에 대한 관심이 증가하면서, 코호트 연구들이 증가하고 있다. 예를 들어, 국민건강보험공단이 전국민 대표성을 가지도록 설계해서 2014년에 배포한 건강보험 표본코호트는 2002~2014년의 추적자료를 포함하고 있다. 이러한 자료들을 이용하여 전국 단위 PM_2.5 장기 노출의 건강 영향을 평가하기 위해서는, 전국단위 측정값이 없는 과거 시기인 2015년 이전 전국의 PM_2.5 농도를 예측할 수 있는 방법의 개발이 필수적이다. 그러나, 우리나라에서는 이러한 예측 방법 개발을 위한 연구가 수행되지 않았다. 본 연구에서는 PM_2.5 장기 노출의 건강 영향을 평가하기 위해, 단순하지만 예측력이 높다고 보고된 PM₁₀과의 비율을 이용하는 방법을 적용하여 과거 PM_2.5 농도를 예측하고자 한다. 즉, 2001년부터 2014년까지 PM_2.5 측정값이 없는 전국 시군구의 과거 PM_2.5 연평균 농도를 예측하고자 한다. 예측의 공간적인 단위는, 건강보험, 사망 및 출산 자료 등 우리나라 정부기관에서 제공하는 대부분 건강자료의 주소지 공개 단위인 시군구로 설정하였다.

2. 연구 방법

2. 2 PM_2.5 농도 예측 방법

본 연구에서는 PM₁₀ 농도에 PM_2.5와 PM₁₀의 비율을 곱하는 기존 연구(HKSAR-EPA, 2016; Lall et al., 2004)를 바탕으로, 측정값이 없는 2015년 이전 2001~2014년의 전국 시군구에서 PM_2.5 연평균 농도를 예측하는 방법을 개발하였다 (식 1). PM₁₀ 농도(PM_10,st)는 선행 연구(Kim and Song, 2017; Song and Kim 2016)에서 개발된 방법으로 예측된 전국 시군구별 2001년부터 2014년까지의 PM₁₀ 연평균 농도값이었다. PM_2.5와 PM₁₀의 연평균 비율(Ratio_R,T)은 PM_2.5와 PM₁₀의 측정값을 이용하여 산출하였다. 또한, 시간과 공간에 따라 변하는 비율을 보정하기 위해, 보정요소(Adjustment Factor_R,2015)를 포함하였다.

$\[ {\mathrm{P}\mathrm{M}}_{2.5,\mathrm{s}\mathrm{t}}={\mathrm{P}\mathrm{M}}_{10,\mathrm{s}\mathrm{t}}\times {Ratio}_{R,T}\times {Adjustment factor}_{R,2015}\mathrm{ } \]$

(1)

t: year(2001 to 2014)
s: sigungu area (1 to 252)
R: region (1 to 16 metropolitan cities and provinces)
T: year when PM_2.5 is available (t or fixed year in 2015 or for 2001~2014 depending on the prediction approach shown in Table 1)

예측 방법에 이용된 전국 시군구 단위의 2002~2014년 연평균 PM₁₀ 농도는 선행 연구의 PM₁₀ 농도 예측모형을 통해 산출된 시군구 대푯값을 차용하였다. 선행 연구에서는 2001~2014년 PM₁₀ 측정 자료를 바탕으로 일반크리깅(Universal Kriging) 예측 모형을 이용해서 전국규모 점 예측 모형을 개발하였다(Kim and Song, 2017). 일반크리깅 모형에는 322개 지리변수로부터 부분최소제곱(Partial Least Squares)을 통해 추정한 3개 예측변수와 공간적인 상관관계를 포함하였다. 이 모형을 이용해서 전국 6만여 개 집계구내 최대거주지역 중심점에 PM₁₀ 연평균 농도를 예측한 후, 이들 예측값의 시군구별 평균값을 계산하였다(Song and Kim 2016). 2010년 기준으로 40%에 해당하는 시군구에 측정소가 존재하지 않고, 다른 대부분의 시군구도 1~2개에 불과한 측정소를 포함하는 점을 고려할 때, 집계구 중심점은 인구대표성 있는 시군구 대푯값을 제공할 것으로 기대되었다.

PM_2.5와 PM₁₀의 비율은 PM_2.5와 PM₁₀의 측정 자료가 동시에 있는 측정소의 PM_2.5와 PM₁₀ 연평균 농도를 이용하여 측정소별로 계산한 후 16개 시도(7개 특별·광역시와 9개 도) 단위로 연도별 산술 평균하여 산출하고자 하였다. 그러나 이용 가능한 PM_2.5 자료의 시공간적인 제한으로, 이 비율은 시간 또는 공간적으로 일정한 비율, 서울에 한정된 연도별 비율, 또는 1년(2015년)에 한정된 시도별 비율로 산출하였다. PM_2.5와 PM₁₀의 비율을 시군구가 아닌 16개 시도단위로 산출한 이유는, 2015년 기준 50% 이상의 시군구의 PM_2.5 측정 자료가 없어서 시군구 단위로 비율을 산출할 수 없었기 때문이었다.

본 연구에서는 시공간적으로 제한된 PM_2.5 자료를 이용해서 계산한 비율을 보정하고자, 5가지 예측 방법을 탐색하였다. 각 방법은 PM_2.5와 PM₁₀의 연평균 비율의 시공간적인 분포에 대한 가정 하에 계산된 PM_2.5와 PM₁₀의 비율(Ratio)과 이용된 보정요소(Adjustment factor)에 따라 표 1과 같이 정리할 수 있다.

Table 1.

Five approaches using ratios of PM2.5 to PM10 for predicting PM2.5 annual average concentrations for 2001~2014 in South Korea.

각 5개의 예측 방법에서 사용한 PM_2.5와 PM₁₀의 비율은 다음과 같다. 방법 1(A1: spatial average)과 2(A2: temporal average)는 2015년 전국 평균 비율(Ratio_{National,2015})과 2001~2014년 서울 평균 비율 (Ratio_{Seoul,2001-14})을 이용하였다. 방법 3 (A3: regional average)은 2015년 각 지역 비율(Ratio_R,2015)을 이용한 반면, 방법 4 (A4: annual average)는 서울의 각 연도별 비율(Ratio_Seoul,t)을 이용하였다. 방법 5 (A5: regional and annual average)는 방법 3 또는 방법 4에 보정요소를 추가함으로써 PM_2.5와 PM₁₀ 비율 분포의 시공간적 변이를 보정한 방법이다. 즉, 방법 5-1은 방법 3과 동일한 비율(Ratio_R,2015)을 사용하고 보정요소로 서울의 2015년 비율 대비 각 연도별 비율을 이용하였다. 방법 5-2는 방법 4와 동일한 비율(Ratio_Seoul,t)을 이용하되, 보정요소로 2015년 서울의 비율 대비 각 지역 비율을 적용하였다. 본 연구에서 이용한 측정자료의 시공간적 한계로 인하여 방법 5의 식은 두 가지 서로 다른 가정 적용 과정을 통해서 수식상 동일하게 유도할 수 있다. 즉, 방법 5-1은 방법 3의 비율과 서울의 연도별 비율을 시간 보정요소로 이용하여 방법 5의 결과식을 유도할 수 있다. 또한 방법 4의 비율과 2015년 각 지역 비율을 공간 보정요소로 사용한 방법 5-2 역시 방법 5의 결과식을 유도할 수 있다. A1-A4에서는 보정요소를 고려하지 않았다.

본 연구에서 제안한 5가지 PM_2.5 농도 예측 방법은 다음과 같은 가정 하에 고안되었다. A1은 과거 어느 해의 어느 지역이든 2015년 전국 비율의 총 평균한 값을 이용하여 예측하는 방법으로, 연평균 PM_2.5와 PM₁₀ 비율의 공간적 평균이 연도에 따라 크게 달라지지 않는다는 가정을 적용하였다. A2는 과거 어느 해의 어느 지역이든 2001~2014년 서울 비율의 총 평균한 값을 이용하여 예측하는 방법이고, 비율의 시간적 평균이 공간 즉, 지역에 따라 다르지 않다고 가정하였다. A3은 과거 어떤 해의 한 지역의 PM_2.5 농도를 2015년 동일 지역의 비율을 이용하여 예측하는 방법으로 지역의 비율이 연도에 따라 변하지 않는다는 가정이다. A4는 과거 한 해 어느 지역의 농도는 동일시기의 서울에서 계산된 비율을 이용하여 예측하는 방법이고 동일 연도의 비율은 공간 즉, 지역에 상관없이 일정하다는 가정을 적용하였다. 마지막으로, A5는 A3 또는 A4를 기준으로 시공간 보정요소를 이용하여 예측하는 방법으로 결과적으로 비율이 시공간적으로 달라진다고 가정하였다.

본 연구에서 개발된 PM_2.5 농도 예측 방법의 예측력은 2011년부터 2014년까지 부산광역시의 24개 도시대기 측정망의 PM_2.5 측정 자료를 이용하여 평가하였다. 5가지 예측 방법을 통해 2011년부터 2014년까지 24개 측정소별로 PM_2.5 연평균 농도 예측값(predicted PM_2.5) 농도를 산출하고, 이 값을 동일 측정소의 측정값으로부터 계산된 연평균(observed PM_2.5) 농도와 비교했다. 두 값의 비교를 위한 검정통계량으로 평균제곱근오차(root mean square error: RMSE)와 평균제곱오차(mean square error: MSE)를 이용한 결정계수(MSE-based r-squared value: R²)를 산출하였다(Kim and Song, 2017; Kim et al., 2016b; Keller et al., 2015).

$\[ \mathrm{R}\mathrm{M}\mathrm{S}\mathrm{E}=\sqrt{\frac{\sum {(predicted {PM}_{2.5}\mathrm{ }-\mathrm{ }observed\mathrm{ }{PM}_{2.5})}^2}{n}} \]$

(2)

$\[ {\mathrm{R}}^2=\frac{1-MSE \left(predicted\mathrm{ }{PM}_{2.5}\right)\mathrm{ }}{Var\left(observed {PM}_{2.5}\right)\mathrm{ }}=1-\frac{\sum {(predicted {PM}_{2.5}\mathrm{ }-\mathrm{ }observed\mathrm{ }{PM}_{2.5})}^2}{\sum {(observed\mathrm{ }{PM}_{2.5}\mathrm{ }-\mathrm{ }orbserved {PM}_{2.5})}^2} \]$

(3)

평균제곱근오차(RMSE)는 PM_2.5 연평균 농도와 동일한 단위를 갖고 있으므로 PM_2.5 연평균 농도의 평균이나 표준편차와 오차의 크기 비교가 가능하다(식 2). 평균제곱오차를 기반으로 한 결정계수(R²)는 회귀선을 기준으로 상관계수의 제곱으로 표현되는 회귀 결정계수와는 다르게 기울기가 1인 1 : 1 대응선을 기준으로 하기 때문에 측정값과 예측값 사이의 일치성을 비교하기에 더 적합하다(식 3). 또한 공간적 예측력을 검증하기 위해서 측정소별 2001~2014년 4년 평균을 이용하여 예측값과 측정값 사이의 공간 평균제곱근오차(spatial RMSE)와 공간 결정계수(spatial R²)를 계산하였다.

3. 연구 결과

3. 1 미세먼지 (PM_2.5와 PM₁₀) 측정값 및 비율

본 연구에서 PM_2.5와 PM₁₀의 연평균 농도를 계산하는데 이용한 도시대기 측정망의 측정소는 2015년을 기준으로 PM_2.5는 일반대기오염측정망의 107개 측정소 중에서 95개, PM₁₀은 일반대기오염측정망의 268개 측정소 중에서 215개 측정소에서 측정값을 제공하였다(그림 1). 2015년 도시대기 측정소 226개 중, 84개 측정소에서 PM_2.5와 PM₁₀의 측정값을 동시에 제공하여 PM_2.5와 PM₁₀의 비율을 산출하였다.

Fig. 1.

Map 84 sites for PM2.5/PM10 out of the 226 urban background sites of the air quality regulatory monitoring network used in this study for 2015 in South Korea.

전국 도시대기 측정망 측정소의 측정값으로부터 산출한 PM₁₀ 연평균 농도의 평균은 2001년에 61.39 μg/m³에서 2015년에 47.86 μg/m³으로 감소하는 경향을 나타냈다(표 2, 그림 2). 이러한 경향은 대부분의 지역에서 차이가 없었으며 서울에서 30% 이상의 가장 큰 감소 추세를 보였다. PM_2.5는 서울에서 2001년의 연평균 농도값의 평균이 47.82 μg/m³이었고 2015년에는 23.26 μg/m³으로 PM₁₀의 감소 추세와 비슷한 경향을 보였다. 그러나, 2005년에는 거의 변화가 없거나 약간 감소하는 경향을 보이는 PM₁₀과는 다르게, PM_2.5는 2005년 전후 년도와 비교하여 크게 높은 연평균 농도값을 보이고 있었다(그림 2). 2015년부터 PM_2.5의 측정이 전국 단위로 확대되었으나 충청북도, 충청남도, 전라북도, 경상북도에는 여전히 측정소가 없었고, 대전과 제주도에는 1개의 측정소만 있었다. 2015년 PM_2.5 연평균 농도의 평균은 대전에서 28.99 μg/m³로 가장 높게 나타났고 서울에서 23.26 μg/m³으로 가장 낮게 나타났다.

Table 2.

Descriptive statistics of the annual average concentrations of PM2.5 and PM10 for 2001~2015 in 16 metropolitan cities and provinces of South Korea.

Fig. 2.

Box plots of the annual average concentrations of (a) PM10 in South Korea, (b) PM10 in Seoul, and (c) PM2.5 in Seoul for 2001~2015 (N: the numbers of sites; Boxplot: upper and lower whiskers=1.5 * interquartile range (IQR) from 75% quantile and -1.5 * IQR from 25% quantile; upper and lower limits of boxes=75% and 25% quantile; solid line in the middle of box=median/ 50% quantile).

PM_2.5와 PM₁₀의 서울지역 연평균 비율은 2001~2014년 전 기간 평균이 0.53이었다(그림 3). 다른 해에 비해 이상적으로 PM_2.5 농도가 높았던 2005년을 제외하면, 2001년에 0.58로 가장 높았고 2008년에 0.47로 가장 낮았으나, 연도별 평균의 차이는 크지 않았다. 반면, 2015년 전국 단위의 PM_2.5와 PM₁₀의 비율의 시도별 평균은 측정소가 없었던 4개 도를 제외한 12개의 시도 중 대전에서 0.67로 가장 높았고 경기도에서 0.47로 가장 낮았다. 시도별 평균은 서울에서의 연도별 평균에 비해 상대적으로 큰 차이를 보였다.

Fig. 3.

Ratios of PM2.5 to PM10 at air quality regulatory monitoring sites (a) across 14 years for 2001~2014 in Seoul and (b) across 16 metropolitan cities and provinces for 2015 in South Korea (N: the numbers of sites; Boxplot: upper and lower whiskers=1.5 * interquartile range (IQR) from 75% quantile and -1.5 * IQR from 25% quantile; upper and lower limits of boxes=75% and 25% quantile; solid line in the middle of box=median/50% quantile).

3. 2 PM_2.5 농도 예측

본 연구에서 개발된 PM_2.5 농도 예측 방법을 부산광역시의 PM_2.5 측정 자료를 이용하여 검증한 결과, 평균제곱근오차는 8.91~21.73 μg/m³, 결정계수는 0.49~0.79였다(표 3). 공간 평균제곱근오차(spatial RMSE)는 5.82~19.66 μg/m³, 공간 결정계수(spatial R²)는 0.45~0.84였다.

Table 3.

Validation statistics of the predicted PM2.5 annual average concentrations from five prediction approaches using ratios of PM2.5 to PM10 across 24 sites in Busan for 2011~2014.

5개의 PM_2.5 농도 예측 방법에서는 전반적으로 과소추정하는 것으로 나타났다 (그림 4). 그러나, 비율의 분포가 시공간적으로 일정하지 않다고 가정하여 공간 보정 변수를 추가한 A5가 1 : 1 대응선에 가장 가깝게 다가가며, 가장 낮은 평균제곱근오차(8.91 μg/m³)와 가장 높은 결정계수 (0.79)를 보였다. 이러한 높은 예측력은 연도마다의 변이를 제거한, 공간 평균 제곱근오차와 공간결정계수에서도 일관적이었다(5.82 μg/m³; 0.84). 반면, 비율의 시간적 평균값이 공간에 따라 일정하다고 가정한 방법(A2)은 가장 높은 평균제곱근오차 및 공간 평균제곱근오차(21.73 μg/m³; 19.66 μg/m³)와 가장 낮은 결정계수 및 공간결정계수(0.49; 0.45)를 보였다. 그러나, 시간적인 분포가 공간에 따라 일정하다고 가정한 A4는 공간적인 분포나 공간평균이 시간에 따라 일정하다고 가정한 A1과 A3 보다 예측력이 높았다.

Fig. 4.

Scatter plots of observations and predictions of the annual average PM2.5 concentrations from five prediction approaches using ratios of PM2.5 to PM10 across 24 sites in Busan for 2011~2014.

A5를 이용해서 예측한 2001년부터 2014년 전국 시군구 단위의 PM_2.5 연평균 농도 예측값의 평균은 2001년에 29.66 μg/m³에서 2014년에 24.99 μg/m³으로 감소하는 경향을 나타냈다(그림 5). 서울과 경기도 및 광역시에서 PM_2.5 연평균 농도 예측값이 높게 나타났고, 2005년을 제외하고 대부분의 지역에서 매년 감소하는 경향을 예측값 지도에서 확인할 수 있었다.

Fig. 5.

Maps of the annual average concentrations of PM2.5 predicted from the prediction approach using spatiotemporally-adjusted ratios of PM2.5 to PM10 in (a) 2001, (b) 2007, and (c) 2014 in South Korea.

4. 고 찰

본 연구는 전국규모 PM_2.5 자료가 이용가능해진 2015년 이전 기간인 2001년부터 2014년까지의 전국시군구 단위 PM_2.5 연평균 농도를 예측하기 위해, 2001년부터 2014년 서울과 2015년 전국의 PM_2.5와 PM₁₀의 대기오염측정망 측정 자료를 바탕으로 계산한 PM_2.5와 PM₁₀의 비율을 이용하였다. 이때 비율의 시공간적인 변화를 고려하고자, 시간과 공간에 따라 일관된 비율부터 보정을 통해 변화를 반영한 비율까지 5가지 다른 예측 방법을 비교하였다. 예측력의 비교 결과, PM_2.5와 PM₁₀의 비율의 시공간적인 분포의 차이를 가정한 방법에서 예측력이 가장 높은 것을 확인할 수 있었다.

5. 결 론

본 연구의 결과는 우리나라의 전국 단위 PM_2.5 장기 노출의 건강 영향을 평가를 위해 PM_2.5와 PM₁₀의 비율을 이용하는 농도 예측 방법의 높은 활용 가능성을 제안한다. 우리나라에서 PM_2.5 장기 노출의 건강영향에 대한 연구는, 최근 표본코호트를 포함해서 전국적으로 사용 가능한 코호트가 확대됨에 따라 연구기회가 확대되어 왔다. 그러나, 전국규모 PM_2.5의 2015년 이전 자료 부재로 이러한 연구 기회는 2015년 이전 PM_2.5 자료 이용이 가능한 서울지역 연구에 제한되어 왔다(Kim et al., 2017a; Kim et al., 2016a). 본 연구에서는 비교적 간단하면서도 설명력이 높은 농도 예측 방법을 제시함으로써, 향후 PM_2.5와 관련된 오염분포, 건강영향 및 질병부담 등에 대한 후속 연구에서 손쉽게 적용될 수 있을 것으로 기대된다.

Acknowledgments

본 연구는 2013년 교육부의 재원으로 한국연구재단의 지원을 받아 수행된 기초연구사업 (과제번호: 2013R1A6A3A04059017, 2018R1A2B600460)과 국립암센터의 기관고유사업(NCC-1810220-01)의 일환으로 수행되었습니다.

References

• Beelen, R., Raaschou-Nielsen, O., Stafoggia, M., Andersen, Z.J., Weinmaryr, G., Hoffmann, B., Wolf, K., Samoli, E., Fischer, P., Nieuwenhuijsen, M., Vineis, P., Xun, W.W., Katsouyanni, K., Dimakopoulou, K., Oudin, A., Forsberg, B., Modig, L., Havulinna, A.S., Lanki, T., Turunen, A., Oftedal, B., Nystad, W., Nafstad, P., De Faire, U., Pedersen, N.L., Östenson, C.G., Fratiglioni, L., Penell, J., Korek, M., Pershagen, G., Eriksen, K.T., Overvad, K., Ellermann, T., Eeftens, M., Peeters, P.H., Meliefste, K., Wang, M., Bueno-de-Mesquita, B., Sugiri, D., Krämer, U., Heinrich, J., de Hoogh, K., Schindler, C., Schikowski, T., Adam, M., Phuleria, H., Vilier, A., Clavel-Chapelon, F., Declercq, C., Grioni, S., Krogh, V., Tsai, M.Y., Ricceri, F., Sacerdote, C., Galassi, C., Migliore, E., Ranzi, A., Cesaroni, G., Badaloni, C., Forastiere, F., Tamayo, I., Amiano, P., Dorronsoro, M., Matsoulis, M., Trichopoulou, A., Brunekreef, B., Hoek, G., (2014), Effects of long-term exposure to air pollution on natural-cause mortality: an analysis of 22 European cohorts within the multicentre ESCAPE project, Lancet, 383(9919), p785-795.
• Brauer, M., Freedman, G., Frostad, J., van Donkelaar, A., Martin, R.V., Dentener, F., van Dingenen, R., Estep, K., Amini, H., Apte, J.S., Balakrishnan, K., Barregard, L., Broday, D., Feigin, V., Ghosh, S., Hpoke, P.K., Knibbs, L.D., Kokubo, Y., Liu, Y., Ma, S., Morawska, L., Sangrador, J.L.T., Shaddick, G., Anderson, H.R., Vos, T., Forouzanfar, M.G., Vurnett, R.T., Cohen, A., (2016), Ambient air pollution exposure estimation for the Global Burden of Disease 2013, Environmental Science & Technology, 50, p79-88.
• Burnett, R.T., Brook, J., Dann, T., Delocla, C., Philips, O., Cakmak, S., Vincent, R., Goldberg, M.S., Krewski, D., (2000), Association between particulate- and gas-phase components of urban air pollution and daily mortality in eight Canadian cities, Inhalation Toxicology, 12(Suppl. 4), p15-39.
• Cohen, A., Brauer, M., Burnett, R., Anderson, H.R., Frostad, J., Estep, K., Balakrishnan, K., Brunekreef, B., Dandona, L., Dandona, R., Feigin, V., Freedman, G., Bubbell, B., Jobling, A., Kan, H., Knibbs, L., Liu, Y., Martin, R., Morawska, L., Pope, III, C.A., Shin, H., Straif, K., Shaddick, G., Thomas, M., van Dingenen, R., van Donkelaar, A., Vos, T., Murray, C.J.L., Forouzanfar, M.H., (2017), Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: an analysis of data from the Global Burden of Diseases Study 2015, Lancet, 389, p1907-18. [https://doi.org/10.1016/s0140-6736(17)30505-6]
• Di, Q., Wang, Y., Zanobetti, A., Wang, Y., Koutrakis, P., Choirat, C., Dominici, F., Schwartz, J.D., (2017), Air Pollution and Mortality in the Medicare Population, The New England Journal of Medicine, 376, p2513-2522.
• Eeftens, M., Tsai, M.Y., Ampe, C., Anwander, B., Beelen, R., Bellander, T., Cesaroni, G., Cirach, M., Cyrys, J., de Hoogh, K., De Nazelle, A., de Vocht, F., Declercq, C., Dédelé, A., Eriksen, K., Galassi, C., Gražulevičienė, R., Grivas, G., Heinrich, J., Hoomann, B., Iakovides, M., Ineichen, A., Katsouyanni, K., Korek, M., Krämer, U., Kuhlbusch, T., Lanki, T., Madsen, C., Meliefste, K., Mölter, A., Mosler, G., Nieuwenhuijsen, M., Lodenwening, M., Pennanen, A., Probst-Hensch, N., Quass, U., Raaschou-Nielsen, O., Ranzi, A., Stephanou, E., Sugiri, D., Udvardy, O., Vaskövi, É., Weinmayr, G., Brunekreef, B., Hoek, G., (2012), Spatial variation of PM_2.5, PM₁₀, PM_2.5 absorbance and PMcoarse concentrations between and within 20 European study areas and the relationship with NO₂ – Results of the ESCAPE project, Atmospheric Environment, 62, p303-317. [https://doi.org/10.1016/j.atmosenv.2012.08.038]
• Gehring, U., Heinrich, J., Krämer, U., Grote, V., Hochadel, M., Sugiri, D., Kraft, M., Rauchfuss, K., Eberwein, H.G., Wichmann, H.-E., (2006), Long-Term Exposure to Ambient Air Pollution and Cardiopulmonary Mortality in Women, Epidemiology, 17, p545-551. [https://doi.org/10.1097/01.ede.0000224541.38258.87]
• HKSAR-EPA, (2016), Guidelines on the Estimation of PM_2.5 for Air Quality Assessment in Hong Kong, http://www.epd.gov.hk/epd/english/environmentinhk/air/guide_ref/guide_aqa_model_g5.html (accessed on Sep. 5, 2018).
• Hoek, G., Krishnan, R.M., Beelen, R., Peters, A., Ostro, B., Brunekreef, B., Kaufman, J.D., (2013), Long-term air pollution exposure and cardio- respiratory mortality: a review, Environmental Health, 12(1), p43. [https://doi.org/10.1186/1476-069x-12-43]
• Hogrefe, C., Lynn, B., Goldberg, R., Rosenzweig, C., Zalewsky, E., Hao, W., Doraiswamy, P., Civerolo, K., Ku, J.Y., Sistla, G., Kinney, P.L., (2009), A combined model–observation approach to estimate historic gridded fields of PM_2.5 mass and species concentrations, Atmospheric Environment, 43, p2561-2570. [https://doi.org/10.1016/j.atmosenv.2009.02.031]
• Hystad, P., Demers, P.A., Johnson, K.C., Brook, J., van Donkelaar, A., Lamsal, L., Martin, R., Brauer, M., (2012), Spatiotemporal air pollution exposure assessment for a Canadian population-based lung cancer case-control study, Environmental Health, 11, p22. [https://doi.org/10.1186/1476-069x-11-22]
• Keller, J.P., Olives, C., Kim, S.Y., Sheppard, L., Sampson, P.D., Szpiro, A.A., Oron, A.P., Lindström, J., Vedal, S., Kaufman, J.D., (2015), A unified spatiotemporal modeling approach for prediction of multiple air pollutants in the Multi--Ethnic Study of Atherosclerosis and Air Pollution, Environmental Health Perspectives, 123(4), p301-309.
• Kim, E., Park, H., Hong, Y.C., Ha, M., Kim, Y., Kim, B.N., Kim, Y., Roh, Y.M., Lee, B.E., Ryu, J.M., Kim, B.M., Ha, E.H., (2014), Prenatal exposure to PM₁₀ and NO₂ and children’s neurodevelopment from birth to 24 months of age: Mothers and Children’s Environmental Health (MOCEH) study, Science of the Total Environment, 481, p439-445. [https://doi.org/10.1016/j.scitotenv.2014.01.107]
• Kim, H., Kim, J., Kim, S., Kang, S.H., Kim, H.J., Kim, H., Heo, J., Yi, S.M., Kim, K., Youn, T.J., Chae, I.H., (2017a), Cardiovascular Effects of Long-Term Exposure to Air Pollution: A Population-Based Study With 900 845 PersonYears of Follow-up, Journal of the American Heart Association, 6, pe007170. [https://doi.org/10.1161/jaha.117.007170]
• Kim, K.N., Lim, Y.H., Bae, H.J., Kim, M., Jung, K., Hong, Y.C., (2016a), Long-Term Fine Particulate Matter Exposure and Major Depressive Disorder in a Community-Based Urban Cohort, Environmental Health Perspectives, 124(10), p1547-1553. [https://doi.org/10.1289/ehp192]
• Kim, O., Kim, S.Y., Kim, H., (2017b), Association between Long-Term Exposure to Particulate Matter Air Pollution and Mortality in a South Korean National Cohort: Comparison across Different Exposure Assessment Approaches, International Journal of Environmental Research and Public Health, 14, p1103.
• Kim, S.Y., Sheppard, L., Bergen, S., Szpiro, A.A., Sampson, P.D., Kaufman, J.D., Vedal, S., (2016b), Prediction of fine particulate matter chemical components with a spatio-temporal model for the Multi-Ethnic Study of Atherosclerosis cohort, Journal of Exposure Science and Environmental Epidemiology, 26(5), p520-528. [https://doi.org/10.1038/jes.2016.29]
• Kim, S.Y., Olives, C., Sheppard, L., Sampson, P.D., Larson, T.V., Keller, J.P., Kaufman, J.D., (2017c), Historical prediction modeling approach for estimating long-term concentrations of PM_2.5 in cohort studies before the 1999 implementation of widespread monitoring, Environmental Health Perspectives, 125(1), p38-46.
• Kim, S.Y., Song, I., (2017), National-scale exposure prediction for long-term concentrations of particulate matter and nitrogen dioxide in South Korea, Environmental Pollution, 226, p21-29. [https://doi.org/10.1016/j.envpol.2017.03.056]
• Lall, R., Kendall, M., Ito, K., Thurston, G.D., (2004), Estimation of historical annual PM_2.5 exposures for health effects assessment, Atmospheric Environment, 38, p5217-5226. [https://doi.org/10.1016/j.atmosenv.2004.01.053]
• Li, L., Wu, A.H., Cheng, I., Chen, J., Wu, J., (2017), Spatiotemporal estimation of historical PM_2.5 concentrations using PM₁₀, meteorological variables, and spatial effect, Atmospheric Environment, 166, p182-191.
• Lim, S.S., Vos, T., Flaxman, A.D., Danaei, G., Shibuya, K., Adair-Rohani, H., Amann, M., Anderson, H.R., Andrews, K.G., Aryee, M., Atkinson, C., Bacchus, L.J., Bahalim, A.N., Balakrishnan, K., Balmes, J., Barker-Collo, S., Baxter, A., Bell, M.L., Blore, J.D., Blyth, F., Bonner, C., Borges, G., Bourne, R., Boussinesq, M., Brauer, M., Brooks, P., Bruce, N.G., Brunekreef, B., Bryan-Hancock, C., Bucello, C., Buchbinder, R., Bull, F., Burnett, R.T., Byers, T.E., Calabria, B., Carapetis, J., Carnahan, E., Chafe, Z., Charlson, F., Chen, H., Chen, J.S., Cheng, A.T., Child, J.C., Cohen, A., Colson, K.E., Cowie, B.C., Darby, S., Darling, S., Davis, A., Degenhardt, L., Dentener, F., Des Jarlais, D.C., Devries, K., Dherani, M., Ding, E.L., Dorsey, E.R., Driscoll, T., Edmond, K., Ali, S.E., Engell, R.E., Erwin, P.J., Fahimi, S., Falder, G., Farzadfar, F., Ferrari, A., Finucane, M.M., Flaxman, S., Fowkes, F.G., Freedman, G., Freeman, M.K., Gakidou, E., Ghosh, S., Giovannucci, E., Gmel, G., Graham, K., Grainger, R., Grant, B., Gunnell, D., Gutierrez, H.R., Hall, W., Hoek, H.W., Hogan, A., Hosgood, III, H.D., Hoy, D., Hu, H., Hubbell, B.J., Hutchings, S.J., Ibeanusi, S.E., Jacklyn, G.L., Jasrasaria, R., Jonas, J.B., Kan, H., Kanis, J.A., Kassebaum, N., Kawakami, N., Khang, Y.H., Khatibzadeh, S., Khoo, J.P., Kok, C., Laden, F., Lalloo, R., Lan, Q., Lathlean, T., Leasher, J.L., Leigh, J., Li, Y., Lin, J.K., Lipshultz, S.E., London, S., Lozano, R., Lu, Y., Mak, J., Malekzadeh, R., Mallinger, L., Marcenes, W., March, L., Marks, R., Martin, R., McGale, P., McGrath, J., Mehta, S., Mensah, G.A., Merriman, T.R., Micha, R., Michaud, C., Mishra, V., Mohd Hanafiah, K., Mokdad, A.A., Morawska, L., Mozaffarian, D., Murphy, T., Naghavi, M., Neal, B., Nelson, P.K., Nolla, J.M., Norman, R., Olives, C., Omer, S.B., Orchard, J., Osborne, R., Ostro, B., Page, A., Pandey, K.D., Parry, C.D., Passmore, E., Patra, J., Pearce, N., Pelizzari, P.M., Petzold, M., Phillips, M.R., Pope, D., Pope, III, C.A., Powles, J., Rao, M., Razavi, H., Rehfuess, E.A., Rehm, J.T., Ritz, B., Rivara, F.P., Roberts, T., Robinson, C., Rodriguez-Portales, J.A., Romieu, I., Room, R., Rosenfeld, L.C., Roy, A., Rushton, L., Salomon, J.A., Sampson, U., Sanchez-Riera, L., Sanman, E., Sapkota, A., Seedat, S., Shi, P., Shield, K., Shivakoti, R., Singh, G.M., Sleet, D.A., Smith, E., Smith, K.R., Stapelberg, N.J., Steenland, K., Stöckl, H., Stovner, L.J., Straif, K., Straney, L., Thurston, G.D., Tran, J.H., Van Dingenen, R., van Donkelaar, A., Veerman, J.L., Vijayakumar, L., Weintraub, R., Weissman, M.M., White, R.A., Whiteford, H., Wiersma, S.T., Wilkinson, J.D., Williams, H.C., Williams, W., Wilson, N., Woolf, A.D., Yip, P., Zielinski, J.M., Lopez, A.D., Murray, C.J., Ezzati, M., AlMazroa, M.A., Memish, Z.A., (2012), A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010, Lancet, 380(9859), p2224-2260.
• Ministry of environment (MOE), (2011), Guidelines for installation and management of national air quality monitoring networks, (in Korean).
• Ministry of environment (MOE), (2017), ANNUAL REPORT OF AIR QUALITY IN KOREA 2016, (in Korean).
• Ozkaynak, H., Schatz, A.D., Thurston, G.D., Isaacs, R.G., Husar, R.B., (1985), Relationships between aerosol extinction coefficients derived from airport visual range observations and alternative measures of airborne particle mass, Journal of the Air Pollution Control Association, 35, p1176-1185. [https://doi.org/10.1080/00022470.1985.10466020]
• Paciorek, C.J., Yanosky, J.D., Puett, R.C., Laden, F., Suh, H.H., (2009), Practical large-scale spatio-temporal modeling of particulate matter concentrations, The Annals of Applied Statistics, 3(1), p370-397. [https://doi.org/10.1214/08-aoas204]
• Querol, X., Alastuey, A., Ruiz, C.R., Artiñano, B., Hansson, A.C., Harrison, R.M., Burningh, E., ten Brink, H.M., Lutz, M., Bruckmann, P., Straehl, P., Schneider, J., (2004), Speciation and origin of PM₁₀ and PM_2.5 in selected European cities, Atmospheric Environment, 38, p6547-6555. [https://doi.org/10.1016/j.atmosenv.2004.08.037]
• Son, J.Y., Bell, M.L., Lee, J.T., (2011), Survival analysis to estimate the association between long-term exposure to different sizes of airborne particulate matter and risk of infant mortality using a birth cohort in Seoul, Korea, Environmental Health Perspectives, 119(5), p725-730.
• Song, C., He, J., Wu, L., Jin, T., Chen, X., Li, R., Ren, P., Zhang, L., Mao, H., (2017), Health burden attributable to ambient PM_2.5 in China, Environmental Pollution, 223, p575-586. [https://doi.org/10.1016/j.envpol.2017.01.060]
• Song, I., Kim, S.Y., (2016), Estimation of Representative Area-Level Concentrations of Particulate Matter (PM₁₀) in Seoul, Korea, Journal of the Korean Association of Geographic Information Studies, 19(4), p118-129, (in Korean with English abstract). [https://doi.org/10.11108/kagis.2016.19.4.118]
• U.S. Environmental Protection Agency (US EPA), (2009), Integrated Science Assessment (ISA) for Particulate Matter (Final Report, Dec 2009), U.S. Environmental Protection Agency, Washington, DC, EPA/600/R-08/139F, 2009, https://cfpub.epa.gov/ncea/risk/recordisplay.cfm?deid=216546 (accessed on Sep. 5, 2018).
• Wang, M., Beelen, R., Eeftens, M., Meliefeste, K., Hoek, G., Brunekreef, B., (2012), Systematic Evaluation of Land Use Regression Models for NO₂, Environmental Science & Technology, 46, p4481-4489. [https://doi.org/10.1021/es204183v]
• Weichenthal, S., Villeneuve, P.J., Burnett, R.T., van Donkelaar, A., Martin, R.V., Jones, R.R., DellaValle, C.T., Sandler, D.P., Ward, M.H., Hoppin, J.A., (2014), Long-term exposure to fine particulate matter: association with nonaccidental and cardiovascular mortality in the agricultural health study cohort, Environmental Health Perspectives, 122(6), p609-615. [https://doi.org/10.1289/ehp.1307277]
• World Health Organization (WHO), (2003), Health aspects of air pollution with particulate matter, ozone, and nitrogen dioxide, http://www.euro.who.int/__data/assets/pdf_file/0005/112199/E79097.pdf (accessed on Sep. 5, 2018).
• World Health Organization (WHO), (2004), Review of Methods for Monitoring of PM_2.5 and PM₁₀, http://www.euro.who.int/__data/assets/pdf_file/0020/143156/E85769.pdf?ua=1 (accessed on Sep. 5, 2018).
• Yanosky, J.D., Paciorek, C.J., Suh, H.H., (2009), Predicting Chronic Fine and Coarse Particulate Exposures Using Spatiotemporal Models for the Northeastern and Midwestern United States, Environmental Health Perspectives, 117, p522-529. [https://doi.org/10.1289/ehp.11692]
• Yi, S.J., Kim, H., Kim, S.Y., (2016), Exploration and Application of Regulatory PM₁₀ Measurement Data for Developing Long-term Prediction Models in South Korea, Journal of Korean Society for Atmospheric Environment, 32(1), p114-126, (in Korean with English abstract).