Home

본 연구에서는 국가태풍센터에서 북서태평양 태풍강도 예측 가이던스 모델로 사용하고 있는 11개 모델들에 대한 예측성능을 평가하였다. 사용된 가이던스 모델들은 7개의 전구모델, 4개의 지역모델로 구성되어 있다. 각 모델의 설명, 자료기간, 공간해상도, 자료 생산 간격 및 주기 등은 Table 1에 기술되어 있다.

11개 모델 중에 우리나라에서 운용 중인 모델은 UM, UM_R, TWRF, KWRF, KEPS이다. UM은 영국기상청의 통합수치모델(Unified Model, UM)을 도입하여 개발한 기상청의 전지구 자료동화예측시스템(Global Data Assimilation and Prediction System, GDAPS)으로 2010년부터 기상청의 현업예측모델로 사용되고 있다(Lee et al., 2011). UM은 전구모델과 지역모델로 구성되어 있다. UM 전구모델은 4차원 변분법을 사용하며 25 km의 공간해상도와 70개의 연직 층을 가지며, 하루 4번 예측정보를 생산한다. 자료기간은 2010년부터 2014년까지이다. 지역모델인 UM_R(UM Regional, UM_R)은 12 km의 공간해상도와 70개 연직 층을 가지는 기상청의지역예보시스템(Regional Data Assimilation and Prediction System, RDAPS)으로 하루 4번의 예측자료를 생산한다. 자료기간은 총 4년(2011~2014년)이다. KWRF는 WRF (Weather Research and Forecasting) 모델을 기반으로 한 기상청의 지역예보시스템 중의 하나이다. KWRF는 10 km의 고해상도 격자와 고급역학을 적용한 고정밀 모델로, 특히 초기 분석장 품질 향상을 위해 3차원 변분법을 이용한 자체분석을 수행할 뿐 아니라 태풍 초기장 개선을 위한 보거싱(Bogusing) 기법을 개발하여 적용하고 있다(Park et al., 2008). 하루 4번 예측자료를 생산하며, 자료기간은 총 4년(2011~2014년)이다. TWRF는 WRF를 기반으로 한 태풍전용 예측모델이다. 15 km의 공간해상도를 가지며 하루 2번 예측을 수행한다. 2010년부터 2014년까지 5년간의 자료가 존재하지만 최대풍속(Maximum Wind Speed, MWS) 자료는 2013년부터 포함되었다. KEPS(Korea Ensemble Prediction System)는 2011년부터 기상청에서 운영해 오고 있는 앙상블 모델로 4차원 변분법을 사용하여 전구모델로부터 초기장을 생성한다. 40 km의 공간해상도와 70개의 연직 층을 갖는다. 초기 섭동은 Ensemble Transform Kalman Filter (ETKF)를 사용하여 생산되고 총 24개의 앙상블 멤버로 구성되어 있다. 하루 2번 예측정보를 생산한다.

GFS는 미국 NCEP (National Centers for Environmental Prediction) 산하 EMC (Environmental Modeling Center)에서 개발되어 NWS (National Weather Service)에서 운영하고 있는 전지구 예측시스템(Global Forecast System)이다(NCEP, 2003). 수평해상도는 35 km이며, SAS (Simplified Arakawa-Schubert) 구름 모수화, 3차원 변분 GSI (Gridpoint Statistical Interpolation) 자료동화, 그리고 Vortex Relocation 방식 초기화 방법을 사용한다. 하루 2회 예측자료를 생산하며, 자료기간은 2008년부터 2014년까지 7년이다.

HWRF (Hurricane WRF)는 미국 NCEP에서 운용 중인 WRF 기반의 지역 허리케인 예측모델이다. GFS 전지구 예측장을 초기 및 경계 값으로 사용하며, 향상된 보거싱 방법과 물리과정 모수화 방법을 사용한다. 태풍 중심을 따라 이동하는 둥지 격자 시스템을 사용하여 태풍 중심 부근에서 3 km까지 고해상도로 예측이 가능하다(Vijay et al., 2014). NOGAPS (Navy Operational Global Atmospheric Prediction System)는 미 해군에서 개발 · 운영 중인 55 km 수평해상도와 30개의 연직 층으로 구성 된 전지구 모델이다. 3차원 자료동화 방법을 사용하며, 자료기간은 2008년부터 2014년까지이다.

ECMWF (European Centre for Medium-Range Weather Forecast)는 유럽 28개국을 회원국으로 구성되어 있는 독립적인 정부간(intergovernmental) 조직에서 운용 중인 전지구 예측모델이다. ECMWF는 최신의 역학코어(Dynamical core), 자료동화(4차원 변분), 물리과정 모수화를 사용하며, 가장 진로예측 성능이 우수한 전지구 모델 중의 하나로 알려져 있다. 25 km 수평격자와 91개의 연직 층을 가지며 하루 2번 예측정보를 생산한다. 자료기간은 2008년부터 2014년까지이다.

JGSM (Japan Global Spectral Model)은 일본 기상청의 전지구 예측모델로 하루 4번 예측정보를 생산한다. 수평해상도는 18 km이며, 4차원 변분 자료동화 방법을 사용한다(Nakagawa, 2009). TEPS (Typhoon Ensemble Prediction System)는 일본 기상청의 전지구 앙상블 모델로 수평해상도는 56 km이며, 11개의 앙상블 멤버로 구성되어 있다(Yamaguchi and Komori, 2009).

11개의 태풍 가이던스 모델들의 예측결과는 열대폭풍(Tropical Storm, TS)급 이상의 태풍에 대해 RSMCTokyo에서 제공하는 6시간 간격의 최대풍속, 중심기압, 그리고 태풍중심 위치자료와 비교하였다. Figure 1은 2012년 15호 태풍 볼라벤에 대해 모델예측 결과와 RSMC-Tokyo의 최대풍속을 비교한 예시이다.

Fig. 1. 
Comparisons of TC intensity (MWS) between model predictions and RSMC-Tokto best track for Typhoon Bolaven (1215) using (a) GFS, (b) JGSM, (c) HWRF, and (d) KWRF. The bold black line indicates the RSMC-Tokyo best track. Thin lines indicate the model predictions. Filled circles represent the issue time of the prediction.

최대풍속의 정의에 있어 RSMC-Tokyo에서는 세계기상기구(World Meteorological Organization, WMO) 권고 기준인 10분 평균 최대풍속을 사용하지만 JTWC는 1분 평균풍속을 사용한다. Simiu and Scanlan (1978)에 따르면 10분 평균풍속은 1분 평균풍속의 약 88%에 해당되며 풍속이 증가할수록 그 차이는 더 커지는 것으로 보고하고 있다. 따라서 1분 평균풍속을 사용하여 초기장을 생성하고 편차조정을 하는 모델들(GFS, HWRF 등)을 10분 평균풍속 기반의 RSMC-Tokyo 자료를 사용하여 성능평가를 할 경우 상대적으로 태풍강도를 과대 모의하는 경향이 커진다. 이러한 문제점을 해결하기 위하여 본 연구에서는 1분 평균풍속 기반의 모델들에 대해서는 RSMC-Tokyo와 JTWC의 관계식을 이용하여 최대풍속을 보정하여 성능을 평가하였다. Figure 2는 2008년부터 2013년까지 JTWC와 RSMC-Tokyo 자료의 최대풍속 산포도를 나타내며, 식 (1)은 이 자료를 바탕으로 구한 두 기관의 최대풍속 선형 회귀식을 나타낸다.

Fig. 2. 
A scatter plot of MWS between RSMC-Tokyo (10-min average) and JTWC (1-min average) best track data during 2008~2013. A thick line represents the linear regression between two data sets.

여기서 최대풍속이 21.65 m s^-1 이상인 경우에만 회귀식을 적용한 이유는 구해진 회귀식이 RSMC의 최대풍속이 21.65 m s^-1 이상일 때부터 JTWC가 더 큰 값을 가지기 때문이다.

대부분의 가이던스 모델들은 일 2회 예측자료를 생산하지만, 일부 모델들의 경우는 일 4회 자료를 생성한다. 본 연구에서는 동일한 조건으로 통계분석을 하기 위하여 일 2회(0000 UTC와 1200 UTC) 예측 결과만을 이용하였다. 분석은 2008년부터 2014년까지 기간에 대해 실시되었지만 모든 모델의 자료가 존재하는 2013년부터 2014년에 대해 대부분의 상세 통계분석이 이루어졌고 나머지 시기에는 연도별 오차변화 경향만이 조사되었다.

본 연구에서 태풍강도 가이던스 모델들의 예측정확도 분석은 주로 최대풍속(m s^-1) 자료를 이용하여 이루어졌다. 태풍 진로예측결과와 태풍 강도의 또 다른 지표인 태풍의 최저중심기압(hPa)은 최대풍속 자료와의 비교 · 분석을 위해서만 제한적으로 사용되었다. 통계분석은 모델과 최적경로자료 사이의 평균절대오차(MAE), 편차(Bias), 그리고 상관계수(Correlation Coefficient, R)를 이용하여 수행되었다(식 (2)와 (3) 참고).

여기서 F_i와 O_i는 각각 모델값과 관측값(최적경로 자료)을 각각 의미하며, n는 자료의 총 개수를 의미한다.

본 연구에서는 먼저 최근 2년간(2013~2014년) 상기 통계분석 항목을 예보시간별, 월별, 태풍 강도별, 급격한 강도 변화 시, 한반도 영향 태풍의 경우, 태풍 강도변화(ΔMWS) 기준, 이동경로별, 그리고 해역별로 나누어 평가하였다. 예보시간별 평가는 모든 모델들이 공통적으로 자료가 존재하는 72 시간까지만 이루어졌다. 사용된 11개 수치모델들은 전구모델(UM, JGSM, TEPS, NOGAPS, GFS, ECMWF, KEPS)과 해상도가 상대적으로 높은 지역모델(UM_R, TWRF, KWRF, HWRF)로 구분하여 분석하였다. 또한 현업 이용 시에 참고할 수 있도록 분석 항목마다 모델의 성능을 순위로 나타내었고, 모델 정확도에 따라 예측성능이 상대적으로 우수한 상위그룹(평균값보다 0.5 표준편차 이상), 평균보다 예측성능이 낮은 하위그룹(평균값보다 0.5 표준편차 이하), 그리고 이 두 중간에 해당에 해당되는 평균그룹으로 분류하였다. 모델 성능의 연변화 경향은 모델별로 2008년부터 2014년까지 자료가 존재하는 기간에 대해 48 h 예측 값을 이용하여 실시하였다. 이 분석에서는 강도뿐만 아니라 진로에 대한 성능 평가를 함께 실시하였다. 태풍 강도 등급별 평가는 국가태풍센터에서 사용하고 있는 태풍 강도 기준에 따라 약(17 m s^-1=MWS < 25m s^-1, weak), 중(25 m s^-1=MWS < 33m s^-1, medium), 강(33 m s^-1=MWS < 44m s^-1, strong), 매우 강(44 m s^-1 = MWS, very strong)으로 나누어 분석하였다.

Elsberry et al. (2014)는 “태풍의 급격한 강도변화”를 24시간 동안 최대풍속이 30 knot 변화하였을 때로 정의하여 사용하였다. 그러나 이 정의는 JTWC를 기준으로 한 1분 평균풍속을 사용한다. 그러므로 급격한 강도변화 크기에 따른 분석에서는 앞서 정의 한 RSMC-Tokyo와 JTWC 최대풍속 관계식(식 (1))을 적용하여 24시간 동안 강도의 증가 또는 약화가 약 12.81 m s^-1 이상인 경우를 급격한 강화(Rapid Intensification, RI)와 급격한 약화(Rapid Weakening, RW)로 정의하여 사용하였다.

한반도 영향태풍에 대한 분석은 국가태풍센터에서 발간한 2013년과 2014년 태풍분석보고서에 분류된 한반도 영향 태풍 6개, 리피(1304), 콩레이(1315), 다나스(1324), 너구리(1408), 할룽(1411), 봉퐁(1419)에 대해 실시하였다. 이 태풍들은 모두 태풍백서에서 정의하는 비상구역(28°N 북쪽, 132°E 서쪽)에 진입하여 한반도에 영향을 주었다(Fig. 3). 태풍 리피(1304)는 2013년 6월 18일에 0000 UTC에 발생하여 6월 21일 0000 UTC에 온대저기압으로 변질되었으며, 6월 20일 0000 UTC부터 비상구역에 진입하며 한반도에 영향을 주었다. 태풍 콩레이(1315)는 같은 해 8월 26일 0600 UTC에 발생하여 8월 30일 0600 UTC에 소멸하였고, 8월 30일 0300 UTC부터는 비상구역에 진입하였다. 태풍 다나스(1324)의 경우 10월 4일 0600 UTC에 발생하여 10월 9일 0000 UTC에 온대저기압으로 변질되었으며 10월 8일 0300 UTC부터 비상구역에 진입하며 한반도에 영향을 주었다. 2014년에 발생한 8호 태풍 너구리의 경우 2014년 7월 3일 0000 UTC에 발생하여 7월 11일 0000 UTC에 소멸하였으며 7월 8일 1800 UTC부터 비상구역에 진입하였다. 태풍 할룽(1411)은 7월 29일 0000 UTC에 발생하여 8월 8일 0600 UTC에 비상구역에 진입하였으며 8월 11일 0000 UTC에 온대저기압으로 변질 되었다. 마지막 태풍 봉퐁(1419)는 10월 3일 0000 UTC에 발생하여 10월 14일 0000 UTC에 온대저기압으로 변질되었으며, 10월 12일 0000 UTC 이후부터 비상구역에 진입하며 한반도에 영향을 주었다. 2014년 12호 태풍 나크리는 한반도 영향 태풍으로 분류되었지만 수명이 너무 짧아 본 연구에서는 분석에서 제외하였다.

Fig. 3. 
Track for typhoons that influenced the Korean peninsula during 2013~2014. Gray box represents a typhoon-emergency region defined by KMA.

태풍 진로별 모델의 성능평가는 Kim et al. (2011)에서 사용된 Fuzzy c-Means Clustering Method (FCM)을 이용하여 태풍진로를 6개의 군집으로 분류한 후 수행하였다. 군집분류에 사용된 자료는 1982부터 2014년까지 북서태평양에서 발생한 모든 태풍을 6시간마다 최대 5일(120 h) 간격으로 만든 총 15377개의 태풍진로이다. 5일 간격으로 태풍 진로를 구성한 것은 태풍의 진로가 긴 경우에 이동하는 중간에 환경이 달라져 모델의 예측성능에 차이가 발생하는 것을 구별하여 평가할 수 있게 하기 위함이다.

연도별로 가이던스 모델의 성능변화는 2008년부터 2014년까지 48시간 예측된 최대풍속(MWS), 중심기압(MCP), 태풍중심위치(Track)의 MAE를 이용하여 조사하였다(Fig. 5). 그림에서 세모(△)-실선 그래프는 전구모델의 성능변화를 나타내고 동그라미(○)-점선 그래프는 지역모델을 나타낸다. 최대풍속 평가에서 최근 가장 큰 향상을 보인 모델은 NOGAPS(2008년 13.1m s^-1에서 2014년 3.9 m s^-1로)와 GFS (2008년 12.6 m s^-1에서 2014년 3.7 m s^-1로)로 나타났다(Fig. 5a). 최근 들어 정확도가 도리어 낮아진 모델도 있었다. ECMWF는 2008년에 오차가 8.5 m s^-1이었으나 해마다 계속 오차가 증가하여 2014년에는 16.2 m s^-1까지 높아졌다. JGSM, UM, KEPS도 최근 들어 성능이 저하된 모델로 분류된다. 전반적으로 지역모델의 성능은 전구모델보다 평균적으로 높았고 시간이 지남에 따라 성능이 좋아지는 경향도 더 뚜렷이 나타났다.

Fig. 5. 
Comparisons of the changing model performance in terms of the MAE for (a) maximum wind speed, (b) minimum central pressure, and (c) track forecasts at the lead time of 48 h during 2008~2014. Solid lines with triangle represent global models. Dashed lines with circle represent regional models.

태풍 중심기압의 예측성능 변화경향은 최대풍속의 결과와 대체로 유사하게 나타났다(Fig. 5b). 그러나 한 가지 특징적인 것은 최근 최대풍속 예측성능이 낮아지고 있는 ECMWF가 중심기압 예측성능은 반대로 향상되고 있다는 사실이다. 실제로 ECMWF는 2008년에 모든 모델 중에 가장 큰 오차(21.8 hPa)를 보였지만 2014년에는 세 번째로 작은 오차(12.6 hPa)를 보였다. 이것은 일반적으로 알려진 태풍의 최대풍속과 중심기압의 선형적인 관계식이 적용되지 않는 결과로 추후 이 원인에 대한 추가적인 분석이 요구된다.

연도별로 가이던스 모델의 예측진로를 평가한 결과에서는 TWRF를 제외한 대부분의 모델들이 성능이 향상되거나 유지한 것으로 나타났다(Fig. 5c). 진로오차가 가장 많이 개선된 모델은 GFS로 나타났다. GFS는 2008년에 48시간 진로오차가 약 265 km에 달하였으나 해마다 정확도가 향상되어 2014년에는 오차가 96 km까지 낮아져 11개 모델 중에 가장 좋은 성능을 보였다. 최근(2012~2014년) 들어 진로오차가 평균적으로 가장 적은 모델은 ECMWF, GFS, HWRF으로 나타났다. ECMWF는 진로오차가 가장 낮은 모델로 알려져 있지만 2014년에는 GFS, HWRF 다음으로 좋은 성능을 보였다. 전구모델과 지역모델을 구분하여 비교하면, 태풍 진로는 강도예측 결과와는 달리 전구모델이 평균적으로 지역모델보다 예측 정확도가 높고 성능개선도 더 크게 이루어진 것으로 나타났다.

가이던스 모델의 계절별 예측성능은 표본수의 분포에 따라 총 4개의 기간으로 나누어 분석하였다(Table 3). 즉, 상대적으로 표본수가 적은 4~7월(봄~초여름)은 묶고, 표본수가 많은 여름철 8월과 9월은 따로 나누고, 10~11월(가을철)은 묶어 성능평가를 실시하였다. 계절별로 평균오차를 비교하면, 10~11월에 예측오차가 가장 크게 나타나는 것이 특징이다. 이것은 10~11월에 북서태평양에서 강한 태풍이 가장 많이 발생하는 것과 관련이 깊다. 실제로 계절별로 발생한 태풍의 최대풍속 평균값을 비교하면 10~11월이 47.4 m s^-1로 4개의 기간 중 가장 큰 값을 보였다. 따라서 강한 태풍에 대해 최대풍속을 크게 과소모의 하였던 전구모델들(UM과 KEPS 등)에 의해 이 시기의 예측오차가 더욱 커진 것으로 사료된다. 전구모델 중에서는 TEPS와 GFS 그리고 지역모델 중에서는 KWRF와 HWRF가 모든 계절에서 좋은 예측성능을 보였다. 반면 전구모델 중에 ECMWF, KEPS, UM 그리고 지역모델 중에서 TWRF는 대부분의 시기에서 하위그룹에 속하였다. 특정 계절에만 좋은 성능을 보이는 모델도 있었다. JGSM은 8월과 10~11월에, UM_R은 8월에, 그리고 TWRF은 10~11월에만 상위그룹의 예측성능을 보였다.

강도별 가이던스 모델들의 예측성능은 모든 예측시간(0~72시간)의 절대오차 평균값을 이용하여 비교하였다(Table 4). 태풍강도의 구분은 국가태풍센터에서 발행한 태풍백서에 표시된 분류방법을 따랐다. 분석결과, 전반적으로 태풍의 강도가 약할수록 가이던스 모델들의 오차가 작게 나타났으며 강도가 강해질수록 오차가 커지는 경향을 보였다. 이러한 현상은 가이던스 모델들에서, 특히 전구모델에서, 전반적으로 나타나는 과소모의 경향 때문에 나타난 것으로 분석된다. TWRF는 예외적으로 태풍 강도가 강해질수록 오차가 감소하는 경향을 보였다.

지역모델 중에 KWRF와 HWRF는 약한(weak 급) 태풍을 제외한 모든 강도의 태풍에서 예측성능이 우수하였다. 특히, 가장 높은 공간해상도(3 km)를 가진 HWRF는 매우 강한(very strong) 태풍에 대한 강도 예측성능이 탁월하였다(평균오차 3.8 m s^-1). 두 번째로 공간해상도가 높은 KWRF(10 km)의 경우에도 강한 태풍에 대한 예측성능이 우수하였다. 전구모델 중에는 TEPS와 GFS가 모든 강도 군에서 6.6 m s^-1 이하의 오차를 보이며 상위그룹에 속하였다. 이 두 모델은 공간해상도가 각각 60 km와 약 1°로 지역모델에 비해 해상도는 낮았음에도 강한 태풍을 잘 모의하여 그 의미가 더욱 크다. 한편, UM과 KEPS는 매우 강한 태풍에 대한 오차가 가장 큰 모델(모두 25 m s^-1 이상)로 나타났다.

태풍의 강도가 급격하게 변하는 경우에 대한 예측 능력은 전체적인 수치모델의 성능을 평가하는 중요한 기준으로 사용될 수 있다(Knaff et al., 2005). Table 5는 급격한 강도 변화 예측성능에 대한 평가결과를 나타낸다. 전구모델 중에는 GFS 그리고 지역모델 중에는 HWRF가 급격한 강도의 강화와 약화를 모두 잘 예측하였다. 반면, UM_R는 RW 경우 그리고 ECMWF와 KWRF는 RI의 경우에 상위그룹의 좋은 성능을 보였다. RI와 RW의 오차를 비교하면, 대부분의 가이던스 모델들(ECMWF, KWRF, HWRF 제외)이 RI의 경우에 오차가 더 커지는 것을 알 수 있다. 이것은 앞서 기술한 바와 같이 많은 수치모델들이 태풍강도를 과소 모의하는 경향이 크기 때문에 나타난 결과로 사료된다.

우리나라는 평균적으로 매년 3~4개의 태풍이 직 · 간접적으로 영향을 미친다(Lee et al., 1992). 본 연구에서는 2013년과 2014년에 우리나라에 영향을 주었던 6개의 태풍에 대해 72시간까지 모든 예측결과를 종합하여 모델별 성능분석을 실시하였다(Table 6). 6개의 태풍은 1304와 1315를 제외하면 모두 매우 강하게(very strong 급) 발달한 태풍이었다.

6개의 모든 태풍에 대한 평균적인 예측성능은 JGSM, TEPS, NOGAPS, GFS, HWRF가 상위그룹에 속하였다. 이 중 GFS와 HWRF는 모든 개별 태풍예측에서 상위그룹에 속하는 성능을 보여주었고, KWRF 또한 태풍 봉퐁(1419) 사례를 제외 한 모든 태풍에서 상위그룹에 속하였다. 반면 UM, ECMWF, KEPS의 경우 평균적인 예측성능이 하위그룹에 속할 뿐 아니라 대부분의 개별 분석에서도 하위그룹이었다. 일본 현업 모델인 JGSM과 TEPS의 경우 태풍 강도가 약했던 태풍(1304, 1315) 경우 더 우수한 예측성능을 보였다.

대부분의 가이던스 모델들은 초기시간부터 큰 오차를 보인다. 그러나 현업예보의 측면에서 보면 초기값은 이미 알고 있기 때문에 초기값의 성능평가는 의미가 없고 그보다 초기시간 이후의 강도변화를 잘 예측하는 것이 더 중요하다. Table 7은 초기값으로부터 태풍의 최대풍속변화(ΔMWS) 예측성능을 72시간까지 모든 예측결과를 종합하여 평가한 결과를 나타낸다.

분석결과, 강도변화 예측에 대한 상관도와 절대오차 분석에서 모두 상위그룹의 좋은 성능을 보인 모델은 GFS, TEPS, HWRF로 나타났다. HWRF는 강도 절대값 뿐 아니라 변화값 평가에서도 가장 좋은 예측성능을 보였다. 반면, 절대값의 상관도 분석에서 0.07의 낮은 값으로 하위그룹에 속했던 UM과 KEPS는 변화값의 평가에서도 0.1 이하의 낮은 상관도를 보여 최하위 성능을 보였다. ECMWF와 KWRF도 변화값에 대한 오차 평가에서는 가장 낮은 성능을 보였다. 한편, JGSM의 경우 초기값 예측 정확도가 높아 절대값 예측에서는 0.79의 상관도를 보이며 상위권 성능을 보였지만(Fig. 4 참고) 강도 변화값에 대한 평가에서는 상관도가 0.4로 낮아져 하위그룹에 속하였다. 편차 분석에서는 대부분의 모델들이 양의 값을 보여 강도변화를 실제보다 더 과대 모의하는 것으로 나타났다. 이러한 결과는 절대값 평가에서 많은 모델들이 과소 모의하는 경향과 상반된 것으로 대부분의 모델들이 과소 모의된 초기값을 보상하기 위해 강도변화를 과대하게 예측하는 경향이 있는 것으로 해석할 수 있다.

현재 기상청에서는 동적 데이터베이스를 기반으로 하는 태풍진로 예측시스템인 DYTRAP (Dynamic database Typhoon Track Prediction)을 사용하고 있다(Lee et al., 2011). 이 방법은 예측하고자 하는 태풍의 이동경로와 가장 유사한 과거 태풍을 선택하고 이 태풍들에 대해 과거 수치모델의 예측성능 결과를 바탕으로 앙상블 평균을 구하는 방식을 사용한다(Jun et al., 2015). 이렇게 유사경로 태풍에 대해 모델 예측성능을 평가하여 예측에 활용하는 이유는 태풍의 강도가 이동경로 상의 해양 및 대기의 환경에 의해 결정되고 모델별로 이러한 환경을 고려하는 성능에 차이가 나기 때문이다. 본 연구에서는 태풍의 이동경로별로 가이던스 모델의 예측성능을 평가하기 위하여 1982년부터 2014년까지 최적경로 자료를 이용하여 태풍의 이동경로를 6개의 군집으로 분류한 뒤 군집별로 오차분석을 실시하였다(Fig. 6).

Fig. 6. 
Six track patterns over the western North Pacific during the TC season for the period of 1982~2014 based on the Fuzzy c-mean clustering method. The thick lines represent the mean tracks of each cluster. Numbers represent the sampling number of each cluster.

분류된 군집별 태풍의 이동경로 특성을 요약하면 다음과 같다. Cluster1(이하, C1)은 평균진로가 남중국해를 지나 서진하여 중국 내륙에 상륙하는 저위도 태풍의 이동경로를 나타낸다. C2는 대만을 지나 동아시아에 상륙하는 태풍의 경로를 나타낸다. C3은 주로 일본 동쪽 해안을 따라 태평양을 향해 북동진하는 진로를 나타낸다. C4는 C3과 비슷한 경향을 보이지만 좀 더 북서쪽으로 치우쳐 일본 내륙을 직접적으로 관통하는 진로를 나타낸다. C5는 140°E 이동에서 발생하여 일본 동쪽의 먼 해양을 통과하여 북동쪽으로 이동하는 경로를 나타낸다. C6는 열대해역 동쪽에서 발생하여 북서진하는 저위도 이동경로를 나타낸다. 본 연구에서는 연구기간 동안에 자료가 적었던 군집 C5와 C6을 제외한 4개의 군집에 대해서 강도 예측성능을 평가하였다.

분석결과, 전구모델 중에는 TEPS와 GFS 그리고 지역모델 중에는 HWRF가 모든 군집에서 상위그룹에 속하는 예측성능을 보여주었다. UM과 KEPS는 모든 군집에서 하위그룹의 성능을 보였다. 특정 경로의 군집에서만 좋은 성능을 보이는 모델도 있었다. JGSM은 C2, NOGAPS는 C3, 그리고 UM_R은 C4에서 상위그룹의 좋은 성능을 보였다. 평균적으로 가장 오차가 큰 진로군집은 C3이었고 가장 오차가 작은 진로 형태는 C1이었다.

수치모델들의 해역별 예측성능을 비교하기 위해 48 h 최대풍속의 절대오차를 경위도 10° 간격으로 구간을 나누어 분석하였다(Fig. 7). 분석결과, 대부분의 모델들은 태풍이 가장 많이 발생하고 발달하는 지역(10~20°N, 130~140°E)에서 오차가 가장 컸다. 또한 대부분 저위도 지역에서 모델의 오차가 크고 위도가 높아질수록 오차가 작아지는 경향이 보인다. 이것은 고위도 해역에서 대부분의 태풍은 약화되는 경향이 있어 강도예측이 상대적으로 쉽고 저위도에서는 일반적으로 태풍이 강하게 발달하기 때문에 강도에 비례하여 오차도 커지기 때문에 발생하는 현상으로 사료된다. 해역별로 성능 차이가 많이 나는 모델은 UM, ECMWF, KEPS 등이며, 반대로 해역별 차이가 적은 모델은 NOGAPS, TWRF, HWRF 등으로 나타났다.

Fig. 7. 
Spatial distributions of MAE averaged over 10° × 10° grid for MWS predicted from 11 guidance models at the lead time of 48 h during 2008~2014.