Home

공기괴 역궤적 분석은 대기 중 화학 관측 물질의 농도 상승과 관측 물질의 주요 배출원을 이해하는 데 유용하게 활용되고 있다. 역궤적 분석의 대표적인 도구로 단일 궤적 추정 모델이 꾸준히 사용되고 있는데 이는 단일 궤적 추정모델이 대기 경계층의 난류를 표현함에 있어 오차를 유발할 수 있는 단점에도 불구하고, 계산 속도가 빠르고 운영이 비교적 쉽다는 장점이 있기 때문이다. 그리고 HYSPLIT 모델(HYbrid Single-Particle Lagrangian Integrated Trajectory model, www.arl.noaa.gov/HYSPLIT.php)의 경우는 미 국립 해양 대기청(NOAA)에서 제공하는 최적화된 기상장을 무료로 사용할 수 있고, 정기적인 수정과 보완을 수행함에 따라 비교적 정확하게 역궤적을 모사할 수 있다. 단일 역궤적 모델의 경우, 다수의 역궤적을 통계적인 방법으로 결합하여 하나의 궤적이 초래할 수 있는 오차를 줄일 수 있기 때문에, 이러한 역궤적의 통계 분석방법은, 입자 확산 모델[e.g., FLEXPART 모델(Stohl et al., 2005; Kim et al., 2013)]과 함께, 화학물질의 배출 지역을 추정하기 위한 연구에 최근 활발히 응용되고 있다.

공기괴 역궤적의 통계 분석을 통한 관측 물질들의 잠재 배출원 지역을 추정하고자 할 때 이용되는 방법으로는 군집 분석(cluster analysis)이 있다(Moody and Samson, 1989; Harris and Kahl, 1990; Kim et al., 2008). 군집 분석은 관측 지점에 도달하는 역궤적의 이동경로와 속도를 기준으로 역궤적을 비슷한 몇 개의 군집으로 나누어, 각 군집에 포함된 역궤적의 시간과 상응하는 시간에 관측된 물질 농도의 상관성을 비교하여 특이한 농도를 보였던 시기의 공기괴 이동경로와 배출원 지역의 방향과 같은 대략적 정보를 얻는 방법이다. 이 방법은 역궤적 군집과 오염 농도(pollution event concentration)를 연관시킬 수는 있지만 배출원들의 지역적 위치를 추정하기는 어려운 단점이 있다.

다른 한가지 방법은 Seibert et al. (1994)와 Stohl(1996)이 개발한 역궤적의 공기괴 체류시간(residence time) (Ashbaugh, 1983)을 관측 농도에 적용하여 배출원 지역을 찾아내는 통계 분석방법이다. 오늘날 가장 활발하게 사용되고 있는 역궤적 통계 분석 방법이며(Charron et al., 2000; Aalto et al., 2002; Li et al., 2010), 지역별 배출량(예컨대, GtC per unit time for CO₂)을 계산하는 것은 아니지만, 통계적으로 누적된 격자별 농도가 그 격자의 배출강도(source strength)의 기준이 될 수 있다는 관점에서 배출원 지역을 결정하는 방법이다. 물론 상승농도가 관측된 시기에 도달한 역궤적의 경우, 통과한 모든 격자들에서의 각 체류시간이 상승 농도에 대한 그 격자의 기여도로 환산되기 때문에 배출원 부근에서의 갑작스러운 농도 변화는 다소 과소평가 될 수 있다는 단점도 있다. 그럼에도 불구하고 복잡한 전구 inversion 모델이나 입자 확산 모델에 비해 계산 시간이 짧을 뿐만 아니라, 이전 연구들에서 다양한 화학 성분들의 지역 배출원들을 신뢰성 있게 찾아냄으로써, 특히 현재 새롭게 연구되는 관심 물질들의 일차적인 배출원 지역을 찾는 데 많이 활용되고 있다(Reimann et al., 2004, 2008).

따라서 본 기술 논문에서는 이 역궤적 통계 분석법에 대해서 상세히 소개하고, 이 방법을 안면도 기후 변화감시센터(36.53^oN, 126.32^oE)에서 연속 관측한 이산화탄소(CO²²)의 농도에 적용하여 이산화탄소 오염농도에 기여한 잠재 배출원 지역을 규명하는 연구의 기초를 만들고자 한다. 다음 이어지는 2.1절, 2.2절과 2.3절에서는 HYSPLIT 모델을 이용한 역궤적 계산, 상승농도에 대응하는 시각에 안면도에 도달한 HYSPLIT 역궤적들의 체류시간 통계, 이산화탄소 상승 농도와 체류시간의 결합을 통한 잠재적 지역(격자별) 배출 강도 추정 방법에 대해서, 각각 서술하였다.

본 연구에서는 앞서 기술된 역궤적 모델 통계 분석법을 안면도 기후변화감시센터에서 비분산형적외선 분광기(NDIR)를 이용해 연속 관측한 대기 중 이산화탄소 농도의 1시간 평균자료 중 2010년 일년간의 자료에 적용하였다. 관측 자료로부터 이산화탄소의 지역 배경대기 값(regional background concentrations)을 결정하기 위해, 현재 Advanced Global Atmospheric Gases Experiment (AGAGE)와 같은 대기 질 연속 관측 및 모니터링을 위한 국제 협력 연구진들에서, 긴 수명의 온실기체들의 배경 값 결정에 광범위하게 응용되고 있는 통계적인 방법을 사용하였다(O’Doherty et al., 2001). 이 통계 방법은 안면도 지역의 배경대기 농도가 정규분포를 이룬다는 가정 하에, 이 정규분포의 평균에서 ± 2.5σ 이상 벗어나는 값들을 반복적으로 제거하여 배경농도를 구하는 것이다. 이산화탄소와 같이 일주기와 계절 주기의 자연 변동성에, 인위적 요인에 의한 오염 농도가 추가되는 경우, 위와 같은 통계적인 접근이 오류를 유발할 수 있다. 하지만 본 기술 분석에서 사용하는 농도는 상대적인 증가 농도(오염농도에서 배경농도를 뺀 값)이기 때문에 배경농도의 가능한 과대/과소 평가는 절대 배출강도에는 영향을 줄 수 있으나, 배출강도의 상대적 분포(즉, 배출원 지역의 gradient)에는 영향을 주지 않는다. 따라서, 배경농도 산출과 관련된 복합적인 논의는 본 기술 논문의 범위를 벗어나므로 자세히 언급 하지 않았다. 하지만 향후 정량적인 이산화탄소 배출량 산정을 위해서는 정확한 배경농도 산출이 중요하기 때문에, 현재 안면도와 같은 지역급 관측점들에서의 자연적 온실기체 배경농도 결정에 대한 다양한 연구가 진행중이며, 본 연구진도 이에 대한 추가적인 연구 결과 발표를 준비 중이다.

Figure 3은 안면도 기후변화감시센터에서 연속 관측한 2010년 이산화탄소 시계열 자료를 보여준다. AGAGE 통계법에 따라 결정된 배경농도와 이로부터 분리된 상승농도, 즉 오염농도가 제시되었다. 이 오염 농도와 상응하는 시기의 역궤적 체류시간 결과를 결합하여 얻은 배출강도 식 (4)는 Fig. 4의 동북아시아 지도 위에 ppm 단위로 도시하였다. 그림에서 안면도 기후변화 감시센터 위치는 검정색 원으로 표시되었다. 본 분석은 인위적인 오염농도에 적용하기 때문에, 해양을 통과하는 역궤적의 경로는 고려하지 않았다. 결과에서 보듯이, 높은 배출 강도를 보이며 주요 배출원으로 추정되는 지역은 중국 북경 인근의 산업화 지역과 북부 국경 지역이다. 특히 중국 북부지역은 겨울철 안면도에 도달하는 역궤적의 주요 경로(Fig. 2참조)로서, 이 지역의 난방을 위한 석탄 연료의 사용을 나타내고 있다. 이와 같은 배출원 지역 추정의 정확성은 높은 배출강도를 보이는 지역의 산업구조 및 인구 밀도 등을 모두 고려하여 검증될 필요가 있다. 현재의 분석 결과에 따르면, 안면도에서 관측되는 이산화탄소 오염농도의 주요 배출원은 중국 북경과 북부 국경 지역으로 판단된다. 한편 중국 주요 산업 지역인 양쯔 강 주변의 동남부 해안 지역의 배출강도가 낮은 것은, 안면도에 도달하는 공기괴들이 중국의 중북부 지역에 한정되어 있어 안면도의 오염농도에 대한 영향을 미치지는 않지만, 동남부 해안 지역이 이산화탄소의 배출원이 될 가능성은 배제할 수 없음을 보여준다. 따라서, 이와 같이 가능한 주요 배출원 지역들의 포괄적인 분석을 위해, 앞으로 지역적으로 고르게 분포된 가능한 많은 관측소들에 대한 연속관측지원과 양질의 자료의 확보가 중요하다.

Fig. 3. 
Hourly mean concentrations of CO₂ observed at AMY in 2010. The baseline concentrations, which were determined by the AGAGE statistical method, and the enhancements of polluted air masses are denoted in red and blue, respectively.

Fig. 4. 
Potential source regions of CO₂ pollutions observed at AMY in 2010.

본 기술 논문은, 대기 중 화학 성분들의 연속 관측 자료로부터 가능한 오염 배출원 지역을 추정하기 위해 역궤적 모델을 통계적으로 분석하는 방법을 소개하였다. 본 방법에 따라 결정된 공기괴 역궤적의 지역별, 계절별 체류시간 분포는 관측소의 종관 기상 패턴과 관측소에 영향을 주는 주요 공기괴의 이동경로를 이해하는 데 도움을 준다. 2010년도 안면도 기후변화감시센터에 도달한 역궤적들의 체류시간을 계절별로 계산한 결과, 안면도에 도달한 공기괴는 대부분의 계절에서 중국 동북부 지역의 영향을 받고 있으며, 여름철에는 깨끗한 북 태평양기단의 영향하에 있는 것으로 나타났다. 지역(격자)별 배출강도는 산출된 체류시간과 2010년 안면도에서 연속 관측한 이산화탄소의 농도 분포 중 오염농도를 결합하여 추정하였다. 그결과, 안면도의 이산화탄소 오염농도는 주로 중국 북경인근과 중국 북부 국경지역의 영향에 의한 것으로 나타났다. 본 연구는, 공기괴 역궤적 모델의 통계 분석이 다양한 배출원을 가지고 있는 이산화탄소에 대해서도 장거리 배출원 지역을 추정하는데 유용하게 활용될 수 있음을 보임에 따라, 이산화탄소뿐 아니라 인위적 생성원을 갖는 주요 화학물질들의 오염 배출원의 가능한 지역적 위치를 추정하는 데 응용될 수 있음을 제시한다.