MODIS는 미국 항공우주국(NASA)에서 계획된 지구관측 위성계획 (Earth Observing System)의 일환으로 1999년에 발사된 TERRA위성과 2002년에 발사된 AQUA위성에 탑재된 다파장 복사계 (multi spectral radiometer)이다. 본 연구에서 사용된 MODIS 자료는 위성이 관측하는 직하점에서 약 1 km의 공간해상도를 가지는 원시자료인 level 1 calibrated radiance 자료(코드명 MOD02, MYD02)이다. MODIS 자료는 현재 NASA MODIS 자료 검색 시스템인 LAADS Web (http://ladsweb.nascom.nasa.gov/data/)으로부터 획득하였다. MODIS 관측자료의 범위는 2010년 4월 15일부터 5월 12일 동안 아이슬랜드의 Eyjafjallajökull 화산이 위치한 지역을 포함하는 북위 60.63±12도, 서경 19.62±24도에 해당하는 영역이다(그림 1 참조).

Fig. 1. 
Map of the study area (48.63N~72.62N, 43.62W~4.38E). Eyjafjallajökull volcano is indicated by star.

본 연구에서 궤적분석을 위하여 사용된 모델은 미국 해양 대기청(NOAA)에서 제공하는 Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) 모델(버전 4.7) (Draxler and Rolph, 2003; Rolph, 2003)이다. HYSPLIT 모델의 전진궤적 (Forward trajectory) 수행을 위하여 필요한 기상모델 입력자료는 미국 National Centers for Environmental Prediction (NCEP)에서 제공되는 기상 예측 모델결과인 The Global Data Assimilation System (GDAS) 자료를 사용하였다. GDAS 자료는 현재 NOAA의 기상예측모델 자료 저장용 FTP서버 (ftp://arlftp.arlhq.noaa.gov/pub/archives/gdas1)에서 획득하였다. GDAS는 매일 3시간 간격으로 하루 4회(00, 06, 12, 18 UTC) 자료를 생산하고 있으며, 공간해상도는 경위도 1도의 등간격으로 경도 360개, 위도 181개의 격자단위로 구성되어 있다.

Combination of Satellite and Trajectory Analysis(CSTA)는 위성자료와 궤적분석자료를 이용하여 화산재 구름의 현재 및 미래의 위치정보를 분석하기 위하여 개발된 알고리즘이다. CSTA는 화산재 감시 정보 산출을 위한 위성자료 처리부분, 화산재 이동경로 예측을 위한 궤적분석 부분, 결과물의 가시화 부분으로 구성되어 있으며, 그림 2는 전체적인 자료처리 과정을 나타낸다.

Fig. 2. 
Flow chart for data processing in CSTA.

CSTA에서 첫 단계는 MODIS 위성자료로부터 위성관측시간과 복사량 자료를 이용하여 가시광 채널 자료는 반사도(reflectance)로, 적외선 채널 자료는 밝기온도(brightness temperature)의 물리량으로 전환하는 과정을 거친다. 그리고 이 물리량은 Lee et al. (2014)이 위성자료로부터 화산재의 산출을 위하여 개발된 하이브리드 알고리즘을 이용하여 화산재 마스크(ash mask)와 레이어 고도 (H_ash), 그리고 입자 광학두께(Aerosol Optical Thickness; AOT)와 유효반경(R_eff)값을 생산한다. 적외채널의 사용으로 인한 장점으로는, 우선 야간에도 관측이 가능한 점과, 적외채널에 민감하지 못한 미세입자성분의 영향을 최소화할 수 있는 점이 있다. 하이브리드 알고리즘은 화산재 탐지를 위하여 주로 3.7 μm, 10.8 μm, 12.0 μm을 가지는 적외선 채널의 밝기온도차를 이용하며, 화산재로 탐지된 픽셀은 화산재 마스크로 저장된다. 이를 위하여 위성 영상의 각 픽셀별 밝기온도 자료로부터 다음의 식 (1) (Prata, 1989)과 식(2) (Ellord, 2003)을 만족하는 테스트를 수행한다.

여기서 BT_3.7, BT₁₁, BT₁₂는 각각 3.7 μm, 10.8 μm, 12.0μm에서의 밝기온도값이다. 기본적으로, 화산재 탐지는 식(1)과 식(2)의 조건을 모두 만족시키는 화소 단위의 밝기온도차 기법을 사용한다. 그리고 여기서 탐지된 화산재 화소는 화산재 입자의 특성정보를 산출하기 위하여 다음의 과정을 거친다. 우선, 복사전달 모델을 이용하여 미리 계산된 조견표(Look-Up Table; LUT)를 이용하여 위성이 관측한 복사량과 조견표와 의 평균 제곱근 편차 (Root Mean Square Difference; RMSD)(식(3))를 구하여 RMSD값이 가장 작은 경우를 선택하여 화산재의 고도와 AOT를 역으로 산출할 수 있게 된다.

여기서, N은 식(3)에서 사용된 채널i의 총 수로서 본 연구에서는 3.7 μm, 10.8 μm, 12.0 μm의 세 개의 채널 자료(N=3)가 사용되었다. 위 첨자 calc와 obs는 각각 계산된 값과 관측된 값을 의미한다. 본 연구에서 사용된 하이브리드 알고리즘은 대형 화산 폭발사례에 대하여 적용되었으며, 다른 위성산출물(Envisat Absorbing Aerosol Index (AAI)와 Cloud Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) Vertical Feature Mask (VFM) 자료)과의 비교검증을 통하여 화산재 탐지 정확도가 검증되었다(Lee et al., 2014).

두 번째 단계에서는 이전의 위성자료 처리를 통하여 산출된 화산재 정보를 토대로 궤적분석을 수행한다. 이때 궤적분석 모델의 수행은 이전 단계에서의 MODIS 위성자료 분석결과로부터 화산재 구름의 위치정보 (위도, 경도, 고도)와 시간정보(위성관측시간)를 전달받아서 24시간 후의 등온위 (isentropic) 전진궤적 (forward trajectory)분석을 수행하였다. 모델링 결과에서 궤적의 최종 도착지는 현재의 화산재 구름이 이동하여 도착하는 지점으로 가정하였다.

마지막으로 가시화 단계에서는 위성자료로 분석된 화산재 정보와 궤적분석 결과를 동시에 화면에 표출 한다. 여기서 사용된 방법은 위치정보와 속성정보를 동시에 가지는 결과물로부터 구글어스(Google Earth^TM)에서 사용되는 Keyhole Markup Language (KML)파일로 코딩된다. 생성된 KML파일은 구글어스를 통하여 쉽게 화면상에 가시화할 수 있는 특성을 가지므로, 사용자가 화산이 폭발한 지역뿐만 아니라 화산재가 위치한 시공간적인 정보를 보다 쉽게 이해할 수 있다.

Fig. 3. 
MODIS color composite images (left) and volcanic ash pixels determined by CSTA (right) for (a, b) 15 April 2010 UTC11:35, (c, d) 6 May 2010 UTC11:55, and (e, f) 11 May 2010 UTC13:50, respectively.

화산재 예보를 위한 국제기구인 화산재 예보 센터(Volcanic Ash Advisory Centre; VAAC)의 Eyjafjallajökull 화산재 보고서에 의하면, 2010년 4월과 5월 동안의 기간 중에서 비교적 큰 규모의 폭발이 여러 번 발생한 것으로 기록되었다. 본 연구에서는 전체 연구기간 중 2010년 4월 14일~18일, 2010년 5월 5일~7일, 2010년 5월 10일~11일에 발생한 폭발사례를 3개의 이벤트 기간으로 설정하였다.

그림 3은 아이슬랜드의 Eyjafjallajökull 화산 사례에 대하여 CSTA를 적용하여 탐지된 화산재 분석결과 사례로서 분화구로부터 배출되고 있는 화산재 구름의 형태가 잘 나타나는 경우인 2010년 4월 15일 UTC11:35, 2010년 5월 6일 UTC11:55, 2010년 5월 11일 UTC13:50에 관측된 위성자료를 선택적으로 나타내었다. MODIS 위성자료로부터 산출된 화산재 마스크는 화산재에 해당하는 영역만이 강조된 결과로서 칼라는 적외채널 온도차를 나타낸다. 그리고 MODIS칼라합성영상(color composite image)은 가시광선 영역의 채널을 이용하여 인간의 눈으로 볼 수 있는 자연색에 가까운 색으로 재 합성하여 생성된 영상정보를 제공한다(Lee, 2013b). MODIS 칼라합성영상과 화산재 탐지 결과의 비교를 통하여 화산재 구름과 일반 기상학적 구름의 차이를 정성적으로 파악할 수 있으며, 위성 자료의 분석결과는 화산재 구름을 비교적 정확히 탐지되었음을 알 수 있다.

2010년 4월 15일~19일의 위성영상에 나타난 화산재는 편서풍을 타고 주로 동쪽으로 이동하고 있지만, 이후 2010년 4월 17일에 발생한 저기압의 후면을 따라 남쪽으로 방향을 바꾸어 이동하였다. 2010년 5월 5일~6일 기간 동안의 화산재 구름은 남진하고 있는 것을 확인할 수 있으며, 5월 10일~11일의 화산재구름 남쪽으로 약 2,000 km에 걸쳐 이동하고 있는 모습이 관측되었다. 이러한 화산재 구름의 장거리 이동 특성은 화산에서 배출된 화산재의 배출조건과 기상 조건에 따라 주변지역에 미치는 영향이 다르게 나타난다.

Fig. 4. 
MODIS retrieved AOT11 (left), layer height (middle), and effective radius (right)_for (a, b, c) 16 April 2010 UTC11:35, (d, e, f) 6 May 2010 UTC11:55, and (g, h, i) 11 May 2010 UTC13:50, respectively.

그림 4는 그림 3과 동일한 기간에 아이슬랜드의 Eyjafjallajökull 화산재 사례에 대하여 CSTA를 적용하여 산출된 화산재 AOT₁₁, H_ash, R_eff 결과이다. 2010년 4월 16일은 화산재의 이동방향이 동쪽을 향하고 있으며, MODIS 위성영상에서 확인된 평균값은 약 AOT₁₁=0.51±0.30, H_ash=1.96±1.19 km, R_eff=4.20±1.80 μm, 화산재 구름의 넓이가 약 17,497 km로 산출되었다. 2010년 5월 7일의 화산재 구름의 경우는 AOT₁₁값이 약 0.39±0.34로서 상대적으로 적으나, 화산재의 고도는 약 4km~5 km 부근(전체 평균 3.70±1.19)에서 정동-남향 방향으로 굽은 형태로 약 29,445km²의 비교적 넓은 영역에 분포되어 있음을 확인할 수 있다. 그리고 2010년 5월 10일의 경우는 화산재구름의 특성값은 각각 AOT₁₁=0.27±0.17, H_ash=3.04±1.16 km, R_eff =3.12±1.37 μm으로 산출되었다. 화산재의 AOT₁₁값은 가장 작게 나타났으며, 화산재 구름의 영역은 약 8,794 km²이었다. 분화구로부터 지속적으로 약하게 배출된 화산재가 정남 방향으로 이동하였음을 알 수 있다.

그림 5는 전체 연구기간 동안의 Eyjafjallajökull 화산재의 AOT₁₁, 고도, 유효반경의 시계열 분포를 나타낸다. 이 기간 동안의 평균 분석값은 AOT₁₁=0.32±0.11, H_ash=2.84±1.06 km, R_eff=3.32±0.60 μm, 화산재가 영향을 나타낸 영역은 평균 20,447±19,135km²으로 상당히 넓은 지역에 걸쳐 화산재 구름이 분포했던 것으로 나타났으며 화산재 구름이 분포했던 영역에 대한 일변 변화도 크게 나타났다.

Fig. 5. 
Time series of the MODIS retrieved AOT₁₁, H_ash, effective radius for Eyjafjallajökull volcanic ash observed from 15 April to 12 May 2010.

본 연구에서는 화산폭발 이후 대기중으로 배출된 화산재 구름의 예측을 위하여 위성자료로부터 실제 관측된 화산재 구름의 속성 정보를 입력자료로 하는 궤적분석모델링을 수행하였다. 그림 6은 화산재 구름의 예측을 위한 단계별 과정을 나타낸다. 우선, 이전의 MODIS 화산재 탐지 과정에서 획득한 화산재 구름의 위치정보로부터 9*9개의 화소로 구성된 크기를 가지는 격자 (grid)를 생성한다. 이때, 9*9개의 화소중 화산재로 탐지된 픽셀이 적어도 60개 이상(약 75%)인 경우만 사용한다. 그리고 각 격자의 중심점을 출발점으로 하는 전진궤적분석을 수행하여 각 궤적선이 도달하는 최종 지점으로부터 24시간 이후에 화산재 구름이 이동하게 될 지점을 예측한다. 즉, 전진궤적 분석결과는 화산재가 위치한 지점으로부터 향후 24시간 동안의 공기의 이동 궤적을 나타내므로, 현재의 MODIS 화산재 탐지 결과로부터 화산재의 분포양상을 확인 후, 다음날의 화산재의 분포양상을 예측할 수 있는 중요한 수단이 된다. 그리고 예측된 화산재 구름의 위치정보는 다음날 위성으로부터 실제 관측된 화산재 구름의 정보와 비교를 통하여 검증이 가능하다.

Fig. 6. 
Illustration of 24 hr Forecast for volcanic ash plume from satellite determined ash pixels and forward trajectory modelling.

그림 7은 2010년 4월 15일 MODIS 화산재 탐지결과 (푸른점)와 예측된 화산재의 위치정보 (붉은점)를 나타낸다. 주어진 현재시간의 위성관측자료로부터 24시간 이후의 전진궤적 분석결과를 통하여 얻은 화산재의 예측정보를 통하여 화산재 구름의 영향 범위를 예측할 수 있다. 이 경우, Eyjafjallajökull 화산재 구름은 동남쪽으로 이동하여 24시간 이후에는 영국동부와 네덜란드에 도달하는 것으로 나타났다. 따라서 화산재 영향 범위를 위성자료와 비교하면, 대상지역에 대한 화산재 경보와 같은 부분에서 응용될 수 있을 것이다.

화산재 예측결과의 정확도 검증을 위하여, 예측된 화산재의 위치와 위성으로 관측된 화산재의 위치가 동일한 지점을 나타내는 격자의 수(N_mat)와 전체 예측된 격자의 수 (N_for)의 비율(Forecast to Measurement Ratio; FMR)을 식(4)와 같이 계산하였다.

본 연구에서 Eyjafjallajökull 화산재 예측을 위하여 CSAT 기법을 사용하여 MODIS 영상자료를 분석한 결과, 전체 총 25일의 사례일 중에 대하여 FMR을 계산한 결과, 약 77.8%의 정확도를 가지는 것으로 확인하였다. 이러한 결과를 통하여 화산재 감시 및 예측에 있어 위성자료와 궤적분석자료의 통합 분석기법이 매우 효율적으로 사용될 수 있는 것으로 판단된다.

Fig. 7. 
(left) MODIS detected (blue area) and 24-hr forecasted (red aera) volcanic ash plume and (right) next-day MODIS ash image for (top) 16 April 2010, (middle) 6 May 2010, and (bottom) 11 May 2010, respectively.