Current Issue

본 연구에 사용된 전지구 접합대순환 모형은 PNU CGCM이다. PNU CGCM의 성분모형은 NCAR (National Center for Atmospheric Research) Community Climate Model version 3 (CCM3; Kiehl et al., 1996) AGCM, GFDL Modular Ocean Model version 3 (MOM3; Pacanowski and Griffies, 1998) OGCM, Los Alamos National Laboratory (LANL) Elastic-Viscous-Plastic Model (EVP) Sea-Ice Model (Hunke and Dukowicz, 1997) 그리고 Land Surface Model (LSM; Bonan, 1998)로 구성된다. PNU CGCM의 성분 모형 및 해상도에 대한 상세한 정보는 Table 1에 나타냈다.

PNU CGCM의 대기성분 초기장은 NCEP/DOE(National Centers for Environmental Prediction/Department of Energy Reanalysis2(이하 NCEP-R2)) 자료를 이용하여 생산하였고 지면 초기조건생산을 위하여 Atmospheric Model Intercomparison Project (AMIP)-type 적분을 이용하였다. AMIP-type 적분은 관측된 해양정보를 경계조건으로 처방하고 모형 내에서 대기 및 지면상태가 동화되도록 하는 것이다. 대기 모형에 관측된 대기 상태를 처음 한 번만 초기조건으로 입력하고 관측된 해수면온도를 경계조건으로 하여 연속적으로 적분을 수행하였다. 해양성분 초기장은 Global Ocean Data Assimilation System (GODAS)의 자료를 사용했으며, Variational Analysis using a Filter (VAF)로 40층 전체에 대하여 자료동화를 적용하였다. PNU CGCM에 사용되는 대기, 해양 및 지면 초기장 생산 방법은 Ahn and Lee (2015)와 Ahn et al. (2018a, 2018b)에서 자세히 설명하고 있으며, Fig. 2에 나타냈다.

지역기후모형은 WRF version 3.5 (Weather Research and Forecasting model; Skamarock et al., 2008)를 사용하였고, WRF는 5 km의 수평해상도를 가지며 한반도를 중심으로 동서방향 150개, 남북방향 245개의 격자로 구성된다(Fig. 1). WRF의 초기 및 경계조건은 전지구 접합대순환모형인 PNU CGCM을 이용하여 생산한 전지구 장기예측 결과를 매달 초기화하여 처방하였으므로 해수면 온도 및 심층토양온도(Deep Soil Temperature) 등은 한달 간격으로 업데이트된다. 또한 매달 준 실시간으로 한반도 지역에 대한 고해상도의 예측자료를 생산하기 위하여 둥지 격자(Nesting) 도메인을 적용하지 않고 T42 GCM 해상도에서 5 km로 규모 축소하도록 실험을 수행하였다(Ahn et al., 2018a). 사용된 물리 방안은 지면 과정과 지표층에 각각 Noah Land-Surface Model (Chen and Dudhia, 2001), Monin-Obukhov scheme (Paulson, 1970), 행성 경계층에 Yonsei University (YSU) scheme (Hong et al., 2006), 적운 모수화에 Kain-Fritsch scheme (Kain, 2004), 미세물리 과정에 WRF Single-Moment 3-class (WSM3) scheme(Hong et al., 2004), 그리고 단파, 장파복사 모수화에 각각 Dudhia scheme (Dudhia, 1989), Rapid Radiative Transfer Model (RRTM) scheme (Mlawer et al., 1997)을 사용하였다(Table 2).

Fig. 1. 
The locations of ASOS in South Korea and topography heights (in meters) of the WRF domain.

Fig. 2. 
Schematic diagram of PNU CGCM.

본 연구에서는 PNU CGCM으로 생산된 전구 Hindcast 자료를 WRF의 초기 및 경계조건으로 하여 남한 지역에 대한 역학적 규모축소를 수행하는 방식으로 남한 지역 상세 예측자료를 생산하였다(Ahn et al., 2018a). PNU CGCM을 이용한 Hindcast는 매월 초기화하여 8개월의 lead (lead0~lead7)로 적분하는 방식으로 생산하였고, 1989년부터 2018년까지 30년에 대하여 시간별(hourly) 자료(Hindcast)를 생산하였다. 1월 초기장으로부터 예측을 시작하는 적분을 01RUN, 2월 초기장으로부터 시작하는 적분을 02RUN 등으로 명명하기로 한다. WRF 경계조건은 PNU CGCM으로부터 생산되어 60분 간격으로 처방되고, 후처리를 통하여 일별(daily) 자료를 생산하였다. WRF를 이용하여 역학적 규모축소를 적용하는 lead2~lead7 중 여름철(6월~8월, JJA)에 해당하는 lead(예를 들어 01RUN의 경우 lead 5, lead 6, lead 7)에 대해서만 분석하였다(Fig. 3).

Fig. 3. 
Lead-times in 8-month lead Hindcast experiment. Left column indicates initialized month and top line indicates predicted month. Lead2~7 (dark grey shaded) are dynamically downscaled by WRF.

본 연구에서는 폭염의 예측성을 살펴보기 위하여 폭염일수(Number of Heatwave days, HWD), 폭염강도(Heatwave Intensity, HWI) 및 폭염주의보 발령일수(Heatwave Warning, HWW)의 폭염지수에 대해 분석하였다. 기상청에서 정의한 폭염의 기준은 관측된 일최고기온이 33^oC 이상인 날이다. 따라서 폭염일수는 JJA 동안(총 92일) 일 최고기온이 33^oC 이상인 날의 수, 폭염강도는 일 최고기온이 33^oC를 초과할 경우 초과한 만큼의 기온을 JJA 동안 누적한 것으로 정의한다. 예를 들어 일 최고기온이 36^oC인 경우 그날의 폭염강도는 3^oC이다. 폭염주의보 발령일수는 JJA 동안일 최고기온이 33^oC 이상인 날이 2일 이상 지속되는 날의 수로 정의한다.

일 최고기온 기후값(1989~2018년, 30년)의 남한 평균을 살펴보면 6월, 7월, 8월, JJA에서 72개 지점의 기상청 종관기상관측자료(Automated Synoptic Observing System, ASOS; 이하 관측)(Figs. 4a-d)의 경우 26.38^oC, 28.84^oC, 29.67^oC, 28.29^oC, WRF (Figs. 4m-p)의 경우 20.34^oC, 23.67^oC, 25.79^oC, 23.27^oC로 6.04^oC, 5.17^oC, 3.88^oC, 5.02^oC의 차이를 보이며 WRF가 과소 모의하는 것으로 나타난다. 이처럼 모형의 기후(climatology)는 관측의 기후와 다르기 때문에 WRF 모델 기후에 대한 폭염의 기준을 정의하였다. 관측의 일 최고기온 남한 평균(총 2760개)을 1989년부터 2018년까지 92일(JJA)동안 모두 나열했을 때, 기상청에서 폭염의 기준으로 정의한 33^oC는 상위 6.2%(2590번째)에 해당한다. WRF에서 모의한 지상으로부터 2m 고도에서의 일 최고기온을 관측과 마찬가지 방법으로 나열했을 때 상위 6.2%(2590번째)에 해당하는 기온은 RUN별(01RUN~04RUN)로 각각 27.46^oC, 27.76^oC, 27.82^oC, 27.90^oC이다. 따라서 WRF의 예측결과에서 폭염 기준은 각각 RUN별로 지상으로부터 2m 고도에서의 일 최고기온이 27.46^oC, 27.76^oC, 27.82^oC, 27.90^oC 이상인 날로 정의하였다. 기상청 폭염 기준과 WRF의 RUN별 폭염 기준을 평균한 27.74^oC의 차이는 약 5.26^oC로 이는 WRF와 관측의 일 최고기온 남한 평균 차이의 범위(3.88^oC~6.04^oC)에 포함된다. 본 연구에서는 기상청에서 정의한 폭염 기준인 33^o C에 맞추기 위하여 남한 평균 일 최고기온의 질량 밀도 함수(Probability Distribution Function) 분포의 상위 93.7^th 퍼센타일을 폭염 기준으로 정하였으며, 많은 선행연구에서도 폭염, 한파 등 극한 기온을 다룰 때 이와 같이 상위 90^th, 98^th, 99^th 퍼센타일 등을 이용하여 기준을 선정한다(Mehl and Tebaldi, 2004; Hajat et al., 2006; Kunkel et al., 2010; Zhang, 2011; Vautard et al., 2013; Yin and Sun, 2018).

Fig. 4. 
Spatial distribution of daily maximum temperature (^oC) derived from ASOS (a-d), NCEP-R2 (e-h), PNU CGCM (i-j) and WRF (m-p) in the WRF domain for period of Hindcast (1989~2018).

PNU CGCM-WRF Chain의 HWD, HWI, HWW에 대한 예측성을 검증하기 위하여 시간상관계수(Temporal Correlation Coefficient, TCC), 평균제곱근오차(Root Mean Square Error, RMSE) 및 Heidke 기술 점수(Heidke Skill Score, HSS), 적중률(Hit Rate, HR) 등의 검증 지수들을 활용하였다. TCC는 두 시계열간의 상관성을 나타내며 ±1에 가까울수록 상관성이 높고 0에 가까울수록 상관성이 없음을 의미한다. RMSE는 두 시계열간의 오차를 제곱하여 평균한 값의 제곱근으로서 상호간 편차를 의미한다. 관측과 예측의 RMSE를 구했을 경우, 0에 가까울수록 예측의 오차가 작으며, 값이 클수록 예측의 오차가 커짐을 의미한다. HSS는 당연히 올바르게 예측할 경우를 제거한 정확도를 의미한다. 예를 들면 사막에서 비가 오지 않을 확률처럼 당연한 경우는 예측성에서 배제한다. HSS는 1에 가까울수록 예측성이 좋고 0에 가까울수록 예측성이 없음을, 0 미만일 경우 우연히 맞출 확률이 예측성보다 높음을 의미한다. HR은 예측의 정확성을 백분율로 나타낸 것이다. 그 정도가 클수록 좋으며, 정확도가 완벽할 경우 1을 나타낸다. HSS와 HR은 매년의 폭염 지수들(HWD, HWI, HWW)의 편차(Anomaly)를 시계열 분포로 나타냈을 때, 시계열의 ± 0.43σ(Standard deviation)을 기준으로 하여, 분포를 Above normal, Normal, Below normal로 나누어 관측과 비교평가하였다(Wilks, 1995). 또한, HR은 3분위로 나누어 무작위의 수로 예측했을 때 실제와 일치할 확률이 약 0.33이므로 0.33 이상의 점수를 나타내는 경우 예측성이 있다고 판단하였다. 검증 지수를 구하는 수식은 Table 3에 나타냈다. M, O는 모델값, 관측값을 의미하며 M, O는 모델값 평균, 관측값 평균을, n은 값의 수, s'_M, s'_O는 모델과 관측의 표준편차를 의미한다. 검증을 위하여 사용한 자료는 72개 지점의 기상청 종관기상관측자료(Automated Synoptic Observing System)와 NCEP-R2에서 일단위로 처리하여 제공하는 지상으로부터 2m 고도에서의 일 최고기온 자료를 사용하였다(Fig. 1). 기간은 Hindcast 자료와 마찬가지로 1989년부터 2018년까지 총 30년이다.

Verification method	Define
TCC (Temporal Correlation Coefficient)	CORR=∑MiOi-nM-O-ns'Ms'O (Perfect: 1)
RMSE (Root Mean Square Error)	RMSE=∑i=1nMi-Oi2n (Perfect: 0)
HR (Hit Rate) and HSS (Heidke Skill Score)	Hit Rate=A+F+KP (Perfect: 1) Heidke Skill Score=A+F+K-C1P-C1 (Perfect: 1) C1 = (M×D) + (N×H) + (O×L)÷(P×P)

PNU CGCM-WRF Chain의 폭염 예측성 검증에 앞서 일 최고기온에 대한 PNU CGCM과 WRF의 모의 성능을 ASOS 및 NCEP-R2와 비교하여 검증하였다. Figure 4는 전체 기간(1989~2018) 동안 ASOS, NCEP-R2, PNU CGCM, 그리고 WRF의 여름철 일 최고기온의 기후값을 나타낸다. 6월, 7월, 8월 그리고 JJA에 대하여 월별 분석을 수행하였고, lead2의 적분 결과를 월별로 나타내었다. 먼저 ASOS의 남한 평균 일 최고 기온(Figs. 4a-d) 은 각각 26.38^oC, 28.84^oC, 29.67^oC로 6월에서 8월로 진행될수록 높아진다. 지점별 일 최고기온을 살펴보면, 동해안을 따라 태백산맥이 위치한 곳은 주변보다 기온이 낮게 나타나는 특성이 나타난다. PNU CGCM에서 모의한 여름철 일 최고기온 공간분포를 살펴보면(Figs. 4i-l) 6, 7월에 내륙지역에서 해양지역보다 높은 값을 보인다. 8월에는 대부분의 한반도 영역에서 24^oC 이상의 분포를 나타내고, 남부지방에서는 28^oC 이상의 분포를 나타내며 위도가 낮아질수록 온도가 증가하는 분포를 보인다. 6, 7, 8월 남한 지역(33.5-39^oN, 125-130^oE) 영역평균 일 최고기온은 각각 19.82^oC, 23.99^oC, 25.71^oC이며 ASOS와 마찬가지로 6월에서 8월로 진행될수록 높아진다. 이와 같은 특징은 NCEP-R2의 공간분포(Figs. 4e-h)에서도 매우 유사하게 나타난다. NCEP-R2에서 남한 지역의 여름철 평균 일 최고기온의 분포 역시 내륙지역에서 해양지역보다 높은 값을 보였으며, 월별 영역 평균 일 최고기온은 각각 20.72^oC, 24.03^oC, 25.50^oC로 나타났다. PNU CGCM은 전지구 규모의 모형이기 때문에 한반도 영역의 복잡한 지형효과를 모의하는 데에 한계가 있다. 이러한 한계를 극복하고자 WRF를 이용하여 역학적 규모 축소한 결과를 살펴보면(Figs. 4m-p) PNU CGCM과 NCEP-R2의 일 최고기온 공간분포보다 지형적 특성이 상세하게 나타난다. WRF의 여름철 한반도 일 최고기온 분포는 NCEP-R2 및 PNU CGCM의 결과와 마찬가지로 해양보다 내륙에서 더 높게 나타나고 위도가 낮아질수록 증가한다. 한반도 영역 평균 월별 일 최고기온은 각각 20.34^oC, 23.67^oC 25.79^oC로 NCEP-R2 및 PNU CGCM과 매우 유사하게 나타난다. 반면, 한반도 동쪽에 위치한 태백산맥과 개마고원에 의해 나타나는 한반도의 지형적 특성(동고서저)이 PNU CGCM과 NCEP-R2의 공간분포에서는 잘 나타나지 않았으나 WRF 결과에서는 잘 나타난다. 앞서 언급한 바와 같이 WRF는 ASOS와 6.04^o C, 5.17^oC, 3.88^oC의 차이를 보이며 일 최고기온을 과소 모의하는 특징이 나타나는데, 이는 WRF의 초기 및 경계조건으로 사용한 PNU CGCM의 일 최고기온 분포의 영향 때문인 것으로 사료된다.

Figure 5는 일 최고기온에 대한 ASOS와 WRF의 TCC, RMSE, HR 및 HSS를 나타낸 그림이다. 먼저 TCC를 살펴보면 -0.14~0.34의 범위로 나타나며 월별로 살펴보면 6월은 02RUN, 03RUN, 04RUN, 01RUN, 7월은 04RUN, 03RUN, 01RUN, 02RUN, 8월의 경우 04RUN, 01RUN, 02RUN, 03RUN의 순서로 높게 나타난다. 04RUN의 8월에서는 95% 이상의 신뢰수준으로 유의한 값이 나타난 반면, 01RUN의 6월과 03RUN의 8월에서는 음의 값이 나타난다. 01RUN과 04RUN의 경우에는 오히려 lead가 멀어질수록 TCC가 높게 나타난다. RMSE의 경우 0.70~1.32의 범위에서 나타나며, 6월은 04RUN, 02RUN, 03RUN, 01RUN, 7월은 04RUN, 03RUN, 02RUN, 01RUN, 8월의 경우 04RUN, 01RUN, 02RUN, 03RUN 순서로 낮게 나타난다. HR의 경우 전체적으로 0.27~0.43의 범위에서 나타나며 6월은 02RUN, 03RUN/04RUN, 01RUN, 7월은 04RUN, 01RUN, 02RUN/03RUN, 8월은 02RUN, 01RUN, 04RUN, 03RUN 순서로 높게 나타난다. HSS의 경우 -0.11~0.14 범위에서 나타나며 6월은 02RUN, 03RUN/04RUN, 01RUN, 7월은 04RUN, 01RUN, 02RUN/03RUN, 8월은 01RUN, 02RUN, 04RUN, 03RUN 순서로 높게 나타난다. TCC, RMSE, HR 그리고 HSS를 살펴본 결과, 전반적으로 여름철 일 최고기온 분포에 대한 예측성은 02RUN과 04RUN에서 비교적 높게 나타났다.

Fig. 5. 
The (a) Temporal Correlation Coefficient (TCC) and (b) Root Mean Squared Error (RMSE), (c) Hit Rate (HR), and (d) Heidke Skill Score (HSS) of WRF prediction result. The upper x-axis indicates initialized month, lower x-axis indicates predicted month. Filled (open) circle of (a) TCC indicates the values that are statistically significant at the 95% (90%) confidence level.

Figure 6은 여름철(June-August, JJA) 총 92일 동안의 HWD, HWI 및 HWW에 대한 WRF와 관측의 TCC 및 RMSE를 평균하여 나타낸 그림으로 RUN별(01RUN~04RUN) 72개 지점에 대한 폭염의 예측성을 살펴보고자 하였다. TCC는 1, RMSE는 0에 가까울수록 관측과 비슷함을 의미한다. 관측 72개 지점에 해당하는 RUN별 예측값은 모두 원으로 나타내었으며, 72개 지점을 평균한 남한지역 평균은 별 모양으로 나타내었다. 먼저 HWD, HWI 및 HWW의 평균 TCC를 RUN별로 살펴보면, 01RUN~04RUN은 각각 -0.25~0.28, -0.22~0.31, -0.15~0.42, -0.23~0.64 범위에서 나타나며, 남한 평균은 각각 0.14, 0.02, 0.15, 0.31이다. TCC가 90% (95%) 신뢰수준에서 통계적으로 유의하게 나타나는 지점의 수는 04RUN에서 14개(8개)로 가장 많고 03RUN은 3개(3개)이며, 02RUN과 01RUN에서는 각각 1개(0개)와 0개(0개)로 거의 나타나지 않는다. 95% 신뢰수준에서 통계적으로 유의하게 나타나는 해당 지점은 03RUN의 경우 동두천, 춘천, 철원, 04RUN의 경우 대관령, 대전, 전주, 마산, 보령, 부안, 정읍, 해남으로 RUN별로 다른 지점에서 나타난다. HWD, HWI 및 HWW의 평균 RMSE의 경우 01RUN~04RUN은 각각 2.77~18.99, 3.10~18.86, 3.15~19.61, 3.46~19.93 범위에서 나타나며, 남한 평균은 각각 7.95, 8.19, 8.73, 7.73이다. RMSE가 5 미만으로 낮게 나타난 지점은 01RUN~03RUN의 경우 속초, 흑산도, 장수, 봉화 총 4개, 04RUN의 경우 속초, 흑산도, 장수, 봉화, 대관령 총 5개이다. TCC가 높게 나타난 지점은 RUN별로 다르게 나타난 반면, RMSE는 모든 RUN에서 속초, 흑산도, 장수, 봉화 4개 지점에서 일관성 있게 상대적으로 낮게 나타났다.

Fig. 6. 
Average time Correlation Coefficient (TCC) and average Root Mean Square Error (RMSE) of Number of Heatwave days (HWD), Heatwave Intensity (HWI) and Heatwave Warning (HWW). Filled circle indicate the values of ASOS (72 stations) and star indicates the mean value of South Korea.

RUN별로 살펴보면 TCC는 04RUN, 03RUN, 01RUN, 02RUN 순으로 높게 나타났으며, RMSE는 04RUN, 01RUN, 02RUN, 03RUN 순으로 낮게 나타난 것으로 보아 HWD, HWI, HWW에 대한 예측성은 04RUN에서 상대적으로 높음을 알 수 있다. 04RUN에서 TCC가 가장 높고 RMSE가 가장 낮게 나타나는 것은 일 최고기온의 TCC 및 RMSE가 04RUN에서 가장 높게 나타난 결과(Fig. 5)와도 일치한다.

Figure 7은 HWD, HWI, HWW에 대한 WRF의 HSS를 나타낸 그림이다. HSS는 1에 가까울수록 예측이 정확함을 의미하며 0은 no skill을 나타낸다. HWD에 대한 RUN별 HSS (Figs. 7a-d)를 살펴보면, 남한 평균이 01RUN~04RUN에서 각각 0.02, 0.06, 0.02, 0.04로 02RUN과 04RUN에서 비교적 높게 나타났다. 01RUN(Fig. 7a)의 경우 HSS가 -0.41~0.26의 범위로 나타나며, 서울, 마산, 임실, 거제 등 지점에서는 0.21 이상의 비교적 높은 HSS가 나타나는 반면, 태백, 서귀포 등 일부 지점에서는 약 -0.25 이하의 강한 음의 값이 나타난다. 02RUN (Fig. 7b)에서는 HSS가 -0.41~0.26의 범위로 나타나며 0 이상의 HSS가 나타나는 지점이 다른 RUN들에 비해 상대적으로 많다. 울산, 의성, 태백 등 지점에서는 0.25 이상의 비교적 높은 HSS가 나타나고, 부산 및 태백 등 일부 지점에서는 -0.27 이하의 낮은 HSS가 나타나기도 한다. 03RUN (Fig. 7c)의 경우 HSS가 -0.41~0.31의 범위로 나타나며 춘천, 제주, 서울 에서는 0.30 이상의 비교적 높은 HSS가 나타나는 반면, 태백, 포항 등 지점에서는 -0.31 이하의 음의 값이 분포하기도 한다. 또한 HSS가 0 이하로 나타나는 지점이 03RUN에서 가장 많다. 04RUN(Fig. 7d)에서는 HSS가 -0.41~0.35의 범위로 나타나며 원주, 전주, 보은, 진주 등 일부 지점에서는 0.27 이상의 비교적 높은 HSS가 나타난다. 남한 전체 평균값은 02RUN이 가장 높게 나타나지만, 0.20 이상의 HSS가 나타나는 지점 수는 02RUN과 04RUN에서 동일하게 나타난다. HWI에 대한 RUN별 HSS (Figs. 7e-h)는 남한 평균이 01RUN~04RUN 에서 각각 0.02, 0.09, 0.03, 0.04로 02RUN에서 상대적으로 높게 나타났다. 01RUN (Fig. 7e)은 HSS가 -0.48~0.30의 범위로 나타나며, 고흥, 의성, 고산 등 지점에서는 0.24 이상의 비교적 높은 HSS가 나타나고, 태백과 장수에서는 -0.37 미만의 강한 음의 HSS가 나타난다. 02RUN (Fig. 7f)의 경우 HSS가 -0.48~0.68의 범위로 나타나며 HSS가 0 이상으로 나타나는 지점과 0.2 이상 높게 나타나는 지점 둘 다 다른 RUN들에 비해 많다. 태백지역에서는 0.67의 높은 HSS가 나타나기도 한다. 광주, 완도, 태백 등 지점에서는 0.35 이상의 비교적 높은 HSS가 나타나는 반면, 부산, 성산 등 일부 지점에서는 -0.27 이하인 음의 HSS가 나타나기도 한다. 03RUN (Fig. 7g)의 경우 HSS가 -0.48~0.44의 범위로 나타나며 0보다 큰 HSS가 나타나는 지점이 가장 작다. 제주, 임실, 부여 등 지점에서는 0.20 이상의 비교적 높은 HSS가 나타나고, 태백과 성산 지점에서는 -0.22 이하의 음의 HSS가 나타나기도 한다. 04RUN (Fig. 7h)의 경우 HSS가 -0.48~0.43의 범위로 나타나며 광주, 보령, 해남 등 지점에서는 0.3 이상의 비교적 높은 HSS가 나타나고, 태백과 장수에서는 -0.41 이하의 음의 HSS가 나타나기도 한다. HWW에 대한 RUN별 HSS (Figs. 7i-l)에서는 남한 평균이 01RUN~04RUN에서 각각 0.01, 0.01, 0.02, 0.01로 RUN별로 큰 차이가 없으며, HWI, HWW의 남한 평균에 비해 낮게 나타난다. 01RUN(Fig. 7i)의 경우 HSS가 -0.54~0.26의 범위에서 나타나며, 마산, 인천, 장흥 등 지점에서는 0.25 이상의 비교적 높은 HSS가 나타난다. 반면 여수, 태백 등 일부 지점에서는 -0.41 이하의 음의 HSS가 나타나기도 한다. 02RUN (Fig. 7j)의 경우 HSS가 -0.54~0.46의 범위에서 나타나며, HSS가 0보다 크게 나타나는 지점이 가장 작다. 마산, 제주 등 지점에서는 0.30 이상의 비교적 높은 HSS가 나타나는 반면, 여수 및 태백 지점에서는 -0.41 이하의 음의 HSS가 나타난다. 03RUN(Fig. 7k)의 경우 HSS가 -0.54~0.30의 범위에서 나타나며, 제주 및 충주 등 지점에서는 0.3 이상의 비교적높은 HSS가 나타나고 여수 및 태백 지점에서는 -0.41 이하의 음의 HSS가 나타난다. 04RUN (Fig. 7l)의 경우 HSS가 -0.54~0.42의 범위에서 나타나며, 충주 및 진주 등 지점에서는 0.3 이상의 비교적 높은 HSS가 나타나고, 여수, 장수 및 태백 지점에서는 -0.41 이하의 음의 HSS가 나타난다. HWI의 경우 임실, 고산, 부여 총 3개 지점, HWW의 경우 보은, 의성 총 2개 지점에서 HSS가 3개 이상의 RUN에서 일관되게 0.2 이상 높게 나타난다.

Fig. 7. 
Spatial distribution of Heidke Skill Score (HSS) for Number of Heatwave days (top), Heatwave Index (middle) and Heatwave Warning (bottom). The values at the top left of each figure are mean HSS.

Figure 8은 HWD, HWI, HWW에 대한 WRF의 HR을 나타낸 그림이다. HR은 HSS와 마찬가지로 1에 가까울수록 예측이 정확함을 의미하며, 0.33 이상이어야 예측성이 있다고 할 수 있다. HWD에 대한 RUN별 HR (Figs. 8a-d)을 살펴보면, 남한 평균이 01RUN~04RUN에서 각각 0.35, 0.37, 0.35, 0.36로 HSS와 마찬가지로 02RUN과 04RUN에서 상대적으로 높게 나타났다. 01RUN (Fig. 8a)의 경우 HR이 0.13~0.50의 범위에서 나타나며, 특히 거제, 마산, 서울, 임실 등지점에서는 0.47 이상의 비교적 높은 HR이 분포하였다. 02RUN (Fig. 8b)의 경우 HR이 0.13~0.80의 범위에서 나타나며 전반적으로 0.33 이상의 HR이 나타나는 지점이 01RUN~04RUN 중 가장 많으며 태백, 보은, 울산 등 지점에서는 0.50 이상의 높은 HR이 나타난다. 03RUN (Fig. 8c)의 경우 HR이 0.13~0.53의 범위에서 나타나며, 특히 춘천, 서울, 제주, 봉화 지점에서는 0.50 이상의 높은 HR이 나타난다. 04RUN (Fig. 8d)의 경우 HR이 0.13~0.57의 범위에서 나타나며, 0.4 이상의 HR이 나타나는 지점이 다른 RUN에 비해 가장 많다. 특히, 전주, 진주 등의 지점에서는 0.5 이상의 높은 HR이 나타난다. HWI에 대한 RUN별 HR(Figs. 8e-h)에서는 남한 평균이 01RUN~04RUN에서 각각 0.35, 0.39, 0.36, 0.37로 HSS와 마찬가지로 02RUN에서 상대적으로 높게 나타났다. 01RUN (Fig. 8e)의 경우 HR이 0.10~0.53의 범위에서 나타나고 0.33 이상의 HR이 나타나는 지점이 02RUN 다음으로 많으며, 고산, 고흥, 의성에서는 0.50 이상의 높은 HR이 나타난다. 02RUN (Fig. 8f)의 경우 HR이 0.10~0.80의 범위에서 나타나며, 0.4 이상의 HR이 나타나는 지점이 01RUN~04RUN 중에서 가장 많다. 특히, 광주, 장흥, 태백, 완도, 제주 등 지점에서는 0.57 이상의 높은 HR이 나타나기도 한다. 03RUN (Fig. 8g)의 경우 HR이 0.10~0.63의 범위에서 나타나며, 제주, 임실, 춘천, 보은 등 지점에서는 0.47 이상의 높은 HR이 나타난다. 04RUN (Fig. 8h)의 경우 HR이 0.33~0.63의 범위에서 나타나며, 광주, 보령, 해남 등 지점에서는 0.53 이상의 높은 HR이 나타난다. HWW에 대한 RUN별 HR (Figs. 8i-l)에서는 남한 평균이 01RUN~04RUN에서 각각 0.34, 0.34, 0.35, 0.34으로 HSS와 마찬가지로 RUN별로 큰 차이가 없다. 01RUN (Fig. 8i)의 경우 HR이 0.07~0.50 범위에서 나타나며 서울, 의성, 문경 등 지점에서는 0.47 이상의 높은 HR이 나타난다. 02RUN(Fig. 8j)의 경우 HR이 0.07~0.64 범위에서 나타나며 HR이 0.33 이상인 지점들이 다른 RUN들에 비해 비교적 작게 나타난다. 고흥, 마산, 제주 등에서는 0.53 이상의 높은 HR이 나타난다. 03RUN (Fig. 8k)의 경우 HR이 0.07~0.53 범위에서 나타나며, 서울, 충주, 제주 보은 등 지점에서는 0.53 이상의 높은 HR이 나타난다. 04RUN (Fig. 8l)의 경우 HR이 0.03~0.60 범위에서 나타나며 진주, 보은, 군산, 충주지점에서는 0.5 이상의 높은 HR이 나타난다. HWD의 경우 마산, 고산, 보은, 임실, 장수 총 5개 지점, HWI의 경우 서울, 수원, 대전, 마산, 고산, 양평, 부여, 임실, 문경, 합천 총 10개 지점, HWW의 경우 마산, 진주, 보은, 임실, 문경, 의성 총 6개 지점에서 HR이 3개 이상의 RUN에서 일관되게 0.43 이상 높게 나타난다.

Fig. 8. 
Same as Fig. 5 bur for Hit Rate (HR).

HSS와 HR을 통하여 HWD, HWI, HWW에 대한 지역별 예측성을 살펴본 결과 남한 평균이 HSS는 0.01~0.09 범위로 모두 0 이상, HR은 0.34~0.39의 범위로 모두 0.33이상의 예측성이 나타난다. RUN별로 살펴보면, HSS와 HR 모두02RUN, 04RUN, 01RUN/03RUN, HWI는 02RUN, 04RUN, 03RUN, 01RUN, HWW의 경우 03RUN, 01RUN/02RUN/04RUN 순으로 높게 나타났다. 따라서 HWD, HWI, HWW 모두 대체적으로 02RUN 예측성이 가장 높고 01RUN의 예측성이 가장 낮음을 알 수 있다. 이는 일 최고기온에 대한 HSS는 02RUN/04RUN, 01RUN, 03RUN, HR은 02RUN, 04RUN, 01RUN, 03RUN 순서로 높게 나타난 것과 비교적 일치한다(Fig. 5). HWI는 고산, 부여, 임실, HWW는 보은, 의성 지점의 경우 3개 이상의 RUN에서 일관되게 예측성이 높게 나타나고 HWD는 태백, 부산, 통영, 수원, 남해, HWI는 태백, 부산, 통영, 남해, 대관령, 속초 등, HWW는 태백, 부산, 통영, 대관령, 속초, 남해 등의 경우 3개 이상의 RUN에서 일관되게 예측성이 낮게 나타난다. 이처럼 부산, 남해와 같이 모든 RUN에서 예측성이 낮게 나타나기도 하지만 이외 대부분 지역은 지역별로 뚜렷한 특징이 나타나지 않으므로 예측성이 지형의 영향을 받기보다는 lead time에 따라 매번 다르게 나타나는 것을 알 수 있다. 즉예측성이 높고 낮은 지역이 특정 지역에만 국한되어 있지 않고 고르게 나타나는 것으로 보인다.

Figure 9는 1994년과 2018년과 같이 기록적인 폭염이 발생한 해에 대한 WRF의 폭염 예측성을 살펴보기 위하여 관측과 WRF의 RUN별 Hindcast 자료를 이용하여 평년(기후값), 1994년 그리고 2018년의 HWD, HWI, HWW를 나타낸 그림이다. 분석 지점은 남한 지역(전체 지점의 평균값, S.Korea)과 2018년 관측 시작 이후 111년만의 극값(일 최고기온이 39.6^oC)이 나타난 서울(Seoul) 그리고 분지 지형으로 우리나라에서 가장 더운 대구(Daegu)로 하였다. 먼저 지점별 관측 값을 비교해보면, 남한 평균과 서울에서는 HWD, HWI, HWW의 강도가 2018년, 1994년, 평년 순서로 크게 나타나므로 폭염이 2018년, 1994년 순서로 평년보다 많이 발생하였음을 알 수 있고, 대구에서는 1994년, 2018년, 평년 순서로 크게 나타나므로 폭염이 1994년, 2018년 순서로 평년보다 많이 발생하였음을 알 수 있다. HWD, HWI, HWW의 평년(이후 HWD_CLIM, HWI_CLIM, HWW_CLIM)은 남한 평균과 서울이 비슷하지만, 1994년과 2018년에 발생한 HWD, HWI, HWW(이후 HWD_1994, HWI_1994, HWW_1994, HWD_2018, HWI_2018, HWW_2018)는 서울이 남한 평균보다 크다. 대구에서 발생하는 HWD, HWI, HWW는 평년, 1994년, 2018년 모두 남한 평균보다 크게 나타난다. WRF의 RUN별 예측성능을 살펴보면, 남한 (Fig. 9a)의 경우 HWD_CLIM과 HWW_CLIM은 모든 RUN에서 관측보다 과대 모의하는 것으로 나타나며 그정도는 HWD_CLIM의 경우 04RUN, 03RUN, 01RUN, 02RUN, HWW_CLIM은 04RUN, 01RUN, 03RUN, 02RUN 순서로 크게 나타난다. HWD_CLIM의 경우 WRF 폭염의 기준을 관측과 마찬가지로 남한 평균 일 최고기온의 상위 6.2%에 해당하는 기준으로 동일하게 정의하였음에도 관측과 다르게 나타난다. 이는 HWD_CLIM이 모든 지점에 동일한 폭염 기준을 적용하여 구한 지점별 HWD를 평균한 값이기 때문이다. 먼저 모든 지점의 일 최고기온을 평균한 뒤, 폭염 기준을 적용하여 HWD를 구하면 관측과 모델의 HWD_CLIM은 5.7일로 같은 값이 나타난다. HWI_CLIM의 경우 04RUN은 마찬가지로 과대 모의하지만, 01RUN, 02RUN, 03RUN은 관측보다 과소 모의하는 것으로 나타난다. 03RUN과 04RUN의 경우 HWD_1994, HWI_1994, HWW_1994, HWD_2018, HWI_2018, HWW_2018모두 관측과 마찬가지로 평년보다 크지만, 평년과 1994년의 차이가 11.8~37.2, 평년과 2018년의 차이가 15.5~50.6 범위로 나타나는 관측에 비해 03RUN과 04RUN에서는 평년과 1994년의 차이가 3.0~5.5, 평년과 2018년의 차이가 3.4~9.5 범위로 그 차이가 매우 작다. 반면, 01RUN과 02RUN에서는 1994년과 2018년의 값이 평년과 비슷하거나 오히려 더 작게 나타난다. 다음으로 서울(Fig. 9b)에서 HWD_CLIM, HWI_CLIM, HWW_CLIM을 살펴보면, 모든 RUN에서 WRF가 관측보다 과대 모의하는 것으로 나타나며 그 정도는 HWD_CLIM과 HWW_CLIM의 경우 04RUN, 01RUN, 03RUN, 02RUN, HWI_CLIM의 경우 04RUN, 03RUN, 01RUN, 02RUN 순서로 크게 나타난다. 03RUN과 04RUN의 HWD_1994, HWI_1994, HWW_1994, HWD_2018, HWI_2018, HWW_2018과 02RUN의 HWD_2018, WHI_2018, HWW_2018은 관측과 마찬가지로 평년보다 크게 나타나지만, 평년과 1994년의 차이가 16.0~51.8, 평년과 2018년의 차이가 24.0~82.9 범위로 나타나는 관측에 비해 03RUN과 04RUN에서는 평년과 1994년의 차이가 1.9~23.3, 평년과 2018년의 차이가 1.9~31.1 범위로 그 차이가 매우 작다. 반면, 01RUN의 HWD_1994, HWI_1994, HWW_1994, HWD_2018, HWI_2018, HWW_2018과 02RUN의 HWD_1994 HWI_1994. HWW_1994는 평년보다 더 작게 나타난다. 마지막으로 대구(Fig. 9c)에서 HWD_CLIM과 HWW_CLIM은 모든 RUN에서 관측보다 과대 모의하는 것으로 나타나며 그 정도는 04RUN 03RUN, 01RUN, 02RUN 순서로 크게 나타난다. 반면, HWI_CLIM은 모든 RUN에서 과소 모의하는 것으로 나타났다. 01RUN의 HWD_2018, HWI_2018, HWW_2018, 03RUN의 HWD_1994, HWI_1994, HWW_1994, HWW_2018, 04RUN의 HWD_1994, HWI_1994, HWW_1994에서는 평년보다 크게 나타나지만, 평년과 1994년의 차이가 24.5~119.3, 평년과 2018년의 차이가 8.5~76.9 범위로 나타나는 관측에 비해 WRF에서는 1994년의 차이가 1.1~6.8, 평년과 2018년의 차이가 1.1~12.7 범위로 그 차이가 매우 작다. 그 외에는 오히려 평년보다 작게 나타난다. 또 2018년보다 1994년의 폭염이 더 강했던 관측과 달리 WRF는 1994년보다 2018년을 더 크게 모의한다. 평년, 1994년 그리고 2018년의 HWD, HWI, HWW가 대구, 서울, 남한 순서로 전반적으로 크게 나타나는 관측과 다르게 WRF의 경우 서울, 대구, 남한 순서로 크게 나타난다.

Fig. 9. 
Number of Heatwave days (HWD), Heatwave Intensity (HWI) and Heatwave Warning (HWW) during JJA in 1994, 2018 and climatology over South Korea, Seoul and Daegu.

HWD_CLIM, HWW_CLIM은 01RUN~04RUN에서 전체적으로 과대 모의하는 경향이 있었으며 HWI_CLIM은 과소 모의하는 경향이 나타났다. 과대 모의하는 정도는 04RUN에서 가장 크게 나타나는데, 이는 지점별 표준편차가 RUN별로 각각 HWD의 경우 11.13, 9.70, 12.44, 17.48, HWI의 경우 17.97, 16.09, 18.40, 25.59, HWW의 경우 7.65. 6.77. 9.16. 14.82로 HWD, HWI, HWW 모두 04RUN, 03RUN, 02RUN, 01RUN에서 크게 나타나는 것을 보아 지점별 변동성이 04RUN에서가장 크기 때문에 과대모의 되는 정도도 가장 큰 것으로 보인다. 평년, 1994년 그리고 2018년의 HWD, HWI, HWW가 대구, 서울, 남한 순서로 전반적으로 크게 나타나는 관측과 다르게 WRF의 경우 서울, 대구, 남한 순서로 크게 나타나는 것을 보아 WRF가 지역별 특성을 관측에 비해 잘 모의하지 못하는 것으로 판단된다. 1994년과 2018년 같이 폭염이 평년대비 많이 발생한 특정 해에 대해서는 남한평균과 서울의 경우 03RUN과 04RUN, 대구의 경우 03RUN이 가장 잘 모의하는 것으로 나타나지만, 평년과 강도 차이는 관측보다 작게 나타난다.