안전한 항공기 운항을 위한 현업 전지구예보모델 기반 깊은 대류 예측 지수: Part 2. 계절별 최적화 및 사례 분석

2.1 1년간 자료 기반 최적화의 한계

앞서 언급한 것처럼 1년 기간에 대한 최적화가AU와 TSS를 동시에 증가시키면서 예측 성능을 향상시켰다. 하지만 대류 현상은 계절별로 예측 성능을 파악할 필요가 있다. 이를 위해, 예측 기간에 따른 ACI_NoOpt와 ACI_YrOpt의 ROC curve를 장기간이 아닌 계절별로 나누어 살펴보았다(Fig. 1). 여기서 겨울철은 12월, 1월, 그리고 2월(DJF), 봄철은 3~5월(MAM), 여름철은 6~8월(JJA), 그리고 가을철은 9~11월(SON)로 정의되었다. 겨울을 제외한 대부분의 계절에서는 기존의 ACI_NoOpt보다 ACI_YrOpt의 곡선이 상대적으로 더 좌측 상단으로 확장되어 성능이 향상된 모습을 보여준다. 하지만 여기서 중요한 점은 DJF 기간에서 ACI_YrOpt의 예측 성능이 오히려 현저하게 저하되었다는 것이다. ACI_NoOpt의 6시간, 24시간, 42시간 예측 기간에 대한 AUC는 0.816, 0.811, 그리고 0.739인 반면, ACI_YrOpt의 AUC는 각각 0.769, 0.733, 그리고 0.712로 매우 크게 감소하였다.

Fig. 1.

ROC curves of ACINoOpt (da shed l ine) and ACIYrOpt (solid line) of 6-hour (black), 24-hour (red), and 42-hour (blue) forecast in (a) DJF, (b) MAM, (c) JJA, and (d) SON. The circle, triangle, and square marker indicate the ACI of 0.25, 0.50, and 0.75, respectively.

DJF 기간에서만 1년 자료 최적화의 결과가 다른 양상을 보인 이유는 해당 시기의 관측된 대류 발생 빈도와 큰 관련이 있다. 장기간에 대한 최적의 가중치 조합은 계절과 상관 없이 전체 기간에 대한 대류를 예측하는 데 유리하도록 설정된다. 그러므로 설정된 깊은 대류 발생 빈도가 작은 특정 계절에서는 최적화의 효과가 작을 수 밖에 없다. 검증에 사용된 ETH 관측장에서 15 dBZ 이상의 대류가 관측된 전체 격자수를 계절별로 비교해 보았다(Table 1). DJF, MAM, JJA, 그리고 SON 기간에 대류가 발생한 격자수는 181,797개, 324,186개, 423,795개, 306,471개로, 다른 계절과 비교했을 때 DJF 시기의 관측 수는 다른 계절 대비 대략 절반 정도였으며 전체 격자 수의 고작 14.71%였다. 물론 계절에 따라 사용된 수치 예보모델의 예측 성능도 달라지기 때문에, 단순히 대류 발생 빈도에 따라 선형적으로 최적화의 효과가 증가하는 것은 아니다. 하지만 특정 계절의 매우 낮은 대류 발생 빈도는 이와 관련된 입력 변수의 강도 및 빈도 또한 낮다는 것을 의미하기 때문에 1년 기간을 대상으로 가중치를 최적화하거나 M F를 구축하는 과정에서 이러한 계절적 특성이 고려되는 데 한계가 있다.

Table 1.

Total number of grids where convection with reflectivity over 15 dBZ is observed in radar.

DJF를 제외한 모든 계절에서 ACI_YrOpt의 전체 예측기간에 대한 평균 AUC는 기존의 ACI_NoOpt와 비교하여 모두 증가하였다(Table 2). MAM과 JJA 기간에는 기존 대비 0.906%, 그리고 0.877% 향상되었으며, 특히 SON 의 경우 1.639%가 향상되었다. TSS 값도 해당 계절에서 각각 1.829%, 0.957%, 그리고 4.490%의 향상을 보였다. 하지만 DJF 기간의 경우, 평균 AUC와 TSS의 변화율은 각각 -4.359%와 -31.604%로 모두 매우 크게 감소한 것을 알 수 있다. 이로 인해 전체 계절에 대한 평균 AUC와 TSS들을 비교해 보면 오히려 1년 기간 최적화 이후 감소하는 결과를 보였다. 즉, 대부분의 계절에서는 예측 성능이 향상 되었지만, DJF 기간에 큰 폭의 성능 저하로 인해 전체 계절 평균 예측 성능은 오히려 감소하였다. 이것이 계절에 따라 발생 빈도 및 강도의 특성이 크게 달라지는 대류 현상을 예측하는데 있어서 예측 지수를 1년 단위 기간에 대해 최적화 했을 때 나타나는 한계이며, 계절 변동성을 고려하여 ACI를 개발하고 계절에 따라 검증해야 하는 이유이다.

Table 2.

Forecast period averaged AUC and TSS of ACINoOpt and ACIYrOpt with rate of change in each season.

2.2 계절별 최적화

한반도의 경우 여름철에는 북태평양 고기압이 세력을 확장하면서 그 가장자리를 따라 많은 양의 수증기가 빈번하게 유입된다. 또한, 고기압의 확장으로 인한 지면 가열에 더불어 발해만에서 발달하여 동진하는 상층 기압골 시스템은 종종 찬공기를 남하 시켜 대기를 열적으로 불안정한 상태로 만들기도 한다. 뿐만 아니라, 여름철에 대륙 고기압과 북태평양 고기압과 대치하여 한반도를 중심으로 강한 온위 경도를 이루어 전선을 발달 시킬 수 있다. 이러한 수증기, 불안정도, 그리고 강제 상승과 같은 기상학적 요인들은 여름철을 중심으로 강한 대류를 발생시키는 원인이 된다(Doswell, 1987; Schultz and Schumacher, 1999). 이러한 대류의 계절적 특성은 ACI의 입력 변수들의 계절별 확률 누적 분포에서도 잘 나타난다(Fig. 2). 해당 그림은 각 M F이 정의되는 위험 범위를 나타낸다. CAPE와 APCP의 위험 범위는 JJA 기간에 특히 넓고 높게 나타났다. 이는 대류 특성에 따라 두 변수 모두 여름철에는 강도가 크고 빈번하게 발생하기 때문이다. 반면에 DJF 기간에는 다른 기간에 비해 전반적으로 매우 낮은 값을 보였다. OLR의 경우, 계절별로 위험 범위는 비슷한 폭으로 나타났지만, 강한 대류가 빈번한 JJA 때의 범위는 가장 낮은 구간에 위치해 있다. 이처럼 대류 뿐만 아니라 이와 관련된 입력 변수들도 계절에 따른 강도 및 빈도의 차이는 크게 나타난다. 본 연구에서는 이러한 입력 변수들의 계절적 변동성을 고려하여 ACI 산출에 사용되는 M F를 각 계절별로 구축하였다(Table 3). MF를 정의하는 방법은 파트 1에 자세히 설명되어 있다. 이전의 ACI_NoOpt와 ACI_YrOpt에 사용된 선형 MF들은 모두 6차 방정식으로 이루어져 있다. 하지만 DJF 기간에 대한 APCP의 M F은 6차 방정식으로 정의하였을 때, 위범 범위 내의 높은 구간에서 CDF와 잘 맞지 않는 것을 확인하였다(Fig. 3). 본 연구에서는 더 정교한 7차 회기방정식과의 비교를 통해 겨울철 MF_APCP는 정성적으로 잘 맞으면서 상대적으로 평균 제곱근 오차(Root Mean Square Error, RMSE)도 낮았던 7차 방정식으로 정의하였다.

Fig. 2.

Variability on seasonal MF (CAPE, APCP; 1σ ~ 99%, OLR; 1% ~ -1σ).

Table 3.

Coefficients of MF coefficients in each season.

Fig. 3.

Cumulative distribution of Accumulated Precipitation (APCP) between +1 standard deviation and 99th percentile value in DJF. The blue and black line indicate the 7th and 6th-order polynomial function with formula below. The Root Mean Square Error (RMSE) of each polynomial function is also noted.

계절별로 구축된 각 변수의 M F는 범위가 다르기 때문에, 모든 기간에 동일한 가중치를 적용시킬 수 없다. 이를 해결하기 위해, 계절별로 모든 가중치 조합에 대한 AUC 값을 비교하여 그 시기에 맞는 최적의 가중치 조합을 선정하였다(Table 4). 자세히 말하면, 파트 1에서 전체 1년 기간에 대한 최적화 과정과 동일한 방식을 사용하여, 특정 기간 마다 모든 가중치 조합에 대한 AUC 값들을 비교하여 최대 값을 갖는 조합을 최적의 가중치 조합으로 정의하였다. 흥미롭게도 계절마다 최대 AUC를 갖는 상수 조합은 다르게 나타났다. DJF 시기에는 CAPE와 OLR의 가중치가 작고 APCP은 매우 높은 가중치를 갖는다. 반면, MAM 때는 CAPE의 비중이 더 높았으며, JJA와 SON 시기에는 OLR에 더 높은 가중치가 할당되었다. 이러한 결과는 AUC라는 통계치의 특성 및 사용된 모델이 모사한 강수 및 대류운 모수화의 계절 변동성과 관련이 깊다. 앞서 설명한 것처럼, 본 연구에서는 가중치를 최적화하는 과정에서 최종적인 AUC 값만 비교하여 그 조합을 선정하였다. 여기서 AUC는 고려하는 계절마다 큰 차이를 보일 뿐만 아니라 최적화 과정에서 다양한 요인에 따라 증가될 수 있다. 실제로, 파트 1의 1년 기간에서 보았을 때는 APCP를 거의 고려하지 않고 CAPE와 OLR을 주로 고려한 조합이 장기간 전체에 대한 PODY는 감소시켰지만 POFD를 더 감소시켜 줌으로써 AUC를 증가시킬 수 있었다. 하지만 뒤에서 다룰 계절별로 보이는 AUC가 증가되는 방식은 다양하다. 예를 들면, POFD가 증가하더라도 PODY가 더 큰 폭으로 증가할 때 AUC가 증가할 수 있다. 이렇듯 본 연구에서 사용하는 모델에서는 특정 시기마다 AUC가 증가되는 방식이 달라질 뿐만 아니라, AUC를 증가시키는 데 기여하는 변수들과 기여하는 정도 또한 달라진다. 그러므로 특정 계절에 대한 최적화 이후 가중치가 높게 할당된 변수는 해당 계절에서 보여진 AUC가 증가되는 방식에 더 많이 기여한 변수라고 이해할 수 있다. 이렇게 구축된 계절별 MF와 최적의 가중치 조합을 통해 각 계절에 맞는 ACI를 계산할 수 있으며, 이를 계절별 최적화된 ACI(Seasonal Optimized ACI, ACI_SnOpt)라고 정의하였다(Eq. 1).

Table 4.

Optimized weighting coefficients in each season.

3.1 ACI_SnOpt와 ACI_YrOpt의 예측 성능 비교

본 연구에서는 관측자료의 15 dBZ ETH와 비교를 통해 ACI_YrOpt와 ACI_SnOpt의 예측 성능을 계절별로 검증하였다. 즉, ACI_SnOpt을 이루는 계절별 ACI는 최적화 과정에서 정의된 계절별 기간에 대해 검증이 수행되었으며, ACI_YrOpt는 1년 전체 기간에 대해 최적화가 이루어졌지만, 마찬가지로 계절별로 예측 성능이 검증되었다. 각 계절마다 6시간, 24시간, 그리고 42시간 예측에 대한 ACI_YrOpt와 ACI_SnOpt의 ROC curve를 비교해보았다(Fig. 4). 먼저, 1년 기간 최적화로 인해 예측 성능이 가장 감소했던 DJF기간을 보면, 계절에 맞는 최적화를 진행했을 때 각 예측 기간의 AUC가 0.809, 0.803, 그리고 0.728로 기존과 대비하여 큰 폭으로 증가한 것을 알 수 있다. 또한, MAM의 경우, 전반적으로 ACI_SnOpt의 곡선이 기존보다 좌측 상단으로 약간 넓어지는 결과를 보였다. 6시간과 24시간 예측에서는 AUC가 각각 0.716과 0.675로 기존보다 약간 증가하였지만, 장기간인 42시간 예측에서는 0.636으로 기존의 0.637과 비교하여 약간 감소 또는 비슷한 결과를 보여준다. 즉, 단기(short-term)에 대한 예측 성능이 계절 최적화를 통해 상대적으로 더 크게 향상되었다. JJA 기간에는 단기(6, 24시간) 예측은 비슷하지만 42시간 예측에서 성능이 소폭 감소하였다. 6시간 예측 기간에 대한 AUC는 증가하였지만, 24시간부터 오히려 AUC가 감소하는 모습을 보여준다. 마지막으로 SON 기간의 ROC curve 변화는 M AM 기간과 비슷한 결과를 보여준다. 대부분의 예측 기간에서 곡선은 좌측 상단으로 넓어지며 AUC가 소폭 증가하였으며, 중기(medium range)보다 단기(short-term) 예측에 있어서 계절별 최적화를 통한 성능 향상이 뚜렷하였다.

Fig. 4.

ROC curves of ACIYrOpt (dashed line) and ACISnOpt (solid line) of 6-hour (black), 24-hour (red), and 42-hour (blue) forecast in (a) DJF, (b) MAM, (c) JJA, and (d) SON. The circle, triangle, and square marker indicate the ACI of 0.25, 0.50, and 0.75, respectively.

3.2 계절별 예측 기간에 따른 ACI_SnOpt와 ACI_YrOpt의 AUS와 TSS 비교

예측 기간에 따른 ACI_YrOpt와 ACI_SnOpt의 AUC와 TSS 값을 보면, DJF 때는 모든 예측 기간에서 두 통계치 모두 증가한 것을 알 수 있다(Fig. 5a). 여기서 TSS는 모든 예측 기간에 대해 동일한 증가폭을 보여주고 있으며, AUC는 상대적으로 짧은 예측 기간에서 더 큰 폭으로 증가하였다. 전체 예측 기간에 대한 평균 AUC는 기존의 0.746에서 0.772로 3.485% 증가하였고, TSS는 최적화 이후 0.321로 기존 대비 121.379%나 증가하였다(Table 5). 또한, MAM 기간의 최적화를 통해 TSS는 모든 예측 시간에서 증가하였다(Fig. 5b). 하지만 증가폭은 중기예측보다는 단기예측에서 크게 나타났다. 뿐만 아니라, 앞서 ROC curve를 통해 분석한 것처럼 AUC 또한 단기 예측에서는 증가하지만 36시간 이상 예측 기간에서는 소폭 감소하였다. 해당 기간 ACI_SnOpt의 평균 AUC와 TSS는 각각 0.670, 그리고 0.180으로 기존 대비 0.299%와 7.784%의 향상을 보였다. 계절별 최적화 이후, JJA 기간에는 대부분의 예측 시간대에서 예측 성능이 다소 감소하였다(Fig. 5c). 특히, 짧은 예측 시간대에서는 큰 차이가 없지만, 예측 시간이 길어지면서 AUC와 TSS 모두 감소하는 모습을 보여준다. 기존 ACI_YrOpt의 평균 AUC와 TSS는 0.690과 0.211인 반면 ACI_SnOpt는 0.686과 0.196으로 각각 0.580%, 그리고 7.109% 예측 성능이 저하되는 결과를 보였다. SON의 경우, 대부분의 예측 시간에서 성능이 향상되었다(Fig. 5d). 특히, MAM 기간과 마찬가지로 단기예측 결과에서 AUC와 TSS 모두 큰 폭으로 증가하였으며, 중기예측에 대한 성능은 상대적으로 증가폭이 작았다. 전체 예측 시간에 대해 AUC와 TSS는 0.403%, 그리고 15.102% 향상되었다. 정리하면, 입력 변수의 계절적 변동성을 고려하여 M F와 가중치를 최적한 이후, JJA를 제외한 모든 계절에서 전반적으로 예측 성능이 향상되었다. 특히, 1년 기간 자료기반 최적화 과정에서 큰 성능 저하를 보였던 DJF 시기에서는 성능이 큰 폭으로 향상되었다. 또한, MAM과 SON 기간의 경우 긴 시간 예측보다 짧은 시간의 예측 결과에서 성능이 상대적으로 크게 향상되었다. 반면에, JJA 시기에는 오히려 계절에 따른 최적화 이후 예측 성능이 전반적으로 감소하였으며, 이러한 성능 저하는 예측 시간이 증가하면서 커지는 경향이 나타났다.

Table 5.

Seasonal averaged AUC and TSS of ACINoOpt, ACIYrOpt, and ACISnOpt.

Fig. 5.

AUC (black) and TSS (grey) values of ACIYrOpt (dashed line) and ACISnOpt (solid line) for each forecast time in (a) DJF, (b) MAM, (c) JJA, and (d) SON.

3.3 계절별 ACI_YrOpt와 ACI_SnOpt 사이의 PODY와 POFD 변화 분석

최적화 이후 각 계절에서 나타난 AUC의 변화는 ACI_SnOpt와 ACI_YrOpt의 PODY와 POFD의 차이로 설명될 수 있다. 본 연구에서는 ACI를 0.01 단위로 나누어 각각에 대응되는 PODY와 POFD로 ROC 곡선을 분석하였기 때문에 특정 예측 기간 마다 총 101개의 점들이 존재한다. 또한 해당 점들은 ACI_SnOpt와 ACI_YrOpt의 ROC 곡선에서 각 ACI 값의 PODY와 POFD의 변화와 동일하기 때문에, Fig. 6에서 각 점들의 y-축 값은 ACI_SnOpt의 PODY 값과 ACI_YrOpt의 PODY 값의 차이와 동일하고 x-축은 두 ACI의 POFD 값의 차이와 동일하다. 그러므로 점들이 분포하는 방향은 곧 계절별 최적화 이후 ROC 곡선이 움직이는 방향으로 이해할 수 있다. 해당 그림에서 점선 상단에 위치한 점들이 많다는 의미는 ROC 곡선이 주로 좌측 상단으로 이동하여 곡선이 이루는 면적, 즉 AUC가 증가할 수 있다는 것을 의미한다. 본 연구에서는 편의를 위해 점선 상단을 (+)로 표현하였다. 그러므로 곡선의 하단 영역을 의미하는 (-)영역의 점보다 (+)영역의 점이 더 많을수록, 그리고 그 거리가 멀수록 AUC 증가에 기여할 수 있다. 해당 그림에 대한 해석 방법은 파트 1에도 자세히 설명되어 있다.

Fig. 6.

Scatter plot of change of PODY and POFD between ACIYrOpt and ACISnOpt (ACISnOpt - ACIYrOpt) for 06-hour (black), 24-hour (red), 42-hour (blue) forecast time in (a) DJF, (b) MAM, (c) JJA, and (d) SON. And circle, triangle, square marker indicates the ACI of 0.25, 0.50, and 0.75 for each forecast time.

먼저, 계절별 최적화를 통해 예측 성능이 가장 뚜렷하게 향상되었던 DJF 기간에는 대부분의 점이 1사분면에 위치한 것을 알 수 있다. 이는 PODY와 POFD가 모두 최적화를 통해 증가하였다는 의미이다. 그럼에도 불구하고 AUC가 크게 증가할 수 있었던 이유는 PODY의 증가폭이 더 컸기 때문이다. 실제로 6시간 예측 결과에서 ROC curve를 좌측 상단으로 옮기는 (+)점은 75개, 그리고 우측 하단으로 옮기는 (-)구역의 점은 26개로, 그 차이가 49개였다. 이러한 결과는 24시간과 42시간 예측 시간에서도 유사했다. 하지만 예측 기간이 짧을수록 점들의 위치는 원점으로부터 멀리 위치하여, 단기간 예측 결과에서 PODY의 증가폭이 더 큰 것을 보여준다. 이것이 42시간 예측 결과에서 (+)점과 (-)점의 개수 차이가 짧은 예측 기간 때보다 컸음에도 불구하고 오히려 단기간 예측에서 AUC가 더 크게 증가했던 이유이다(Fig. 4a과 Fig. 5a). 겨울을 중심으로 최적화한 이후, 전반적으로 POFD가 증가했기 때문에, 기존의 ACI_YrOpt보다 ACI_SnOpt가 더 넓은 영역을 예보한다는 것을 알 수 있다. 이러한 통계적 특징은 뒤에 이어질 사례 분석을 통해 다시 확인할 수 있다.

MAM기간에도 PODY와 POFD 모두 증가하였다. 하지만, AUC를 증가시킬 수 있는 (+)지역에 더 많은 점들이 위치해 있다(Fig. 6b). 즉, 적중률과 오탐률이 모두 증가하였지만 전반적으로 적중률이 더 크게 증가하여 예측 성능이 향상된 경우이다. 예측 기간이 짧을수록 원점을 중심으로 더 멀리 위치해 있기 때문에, 단 기간 예측에서 PODY가 더 크게 증가한 것을 알 수 있다. 뿐만 아니라, 6시간 예측 결과에서 (+)점과(-)점의 개수 차이는 42시간 예측 결과보다 더 크다. 이러한 요인들로 인해 Fig. 4b와 Fig. 5b에서 본 것처럼 장기간 보다 단기간 예측에서 예측 성능이 더 향상될 수 있었다. DJF시기와 마찬가지로 최적화 이후 기존보다 POFD가 증가하였기 때문에 ACI_SnOpt가 ACI_YrOpt보다 더 넓은 영역을 예측하는 것을 알 수 있다.

계절별 최적화 과정에서 유일하게 AUC가 감소하였던 JJA 때는 다른 계절과 달리 PODY와 POFD가 모두 감소하였다(Fig. 6c). 특히, 3사분면의 상단의 (+)점보다 하단의 (-)점이 더 많이 존재했으며, 이는 전반적으로 오탐률이 감소하는 정도보다 적중률의 감소폭이 더 컸다는 의미이다. 예측 기간별 차이를 살펴보면, PODY의 감소폭은 전반적으로 비슷하였지만, POFD는 42시간 예측에서 더 적게 감소하였다. 실제로 (+)와 (-)점의 개수 차이는 6시간 예측 기간 때는 13개였지만, 42시간 예측에서는 37개로 급격하게 증가하였다. 이로 인해 예측 기간이 길어질수록 AUC는 더 크게 감소하였다(Fig. 4b와 Fig. 5b).

SON 시기의 최적화 결과는 DJF와 M AM 때의 결과와 유사했다. 마찬가지로 대부분의 ACI점들은 1사분면에 위치해 있었으며, PODY가 POFD보다 더 큰 폭으로 증가하면서 AUC가 증가하였다. 또한, 짧은 예측 기간에서 적중률이 더 많이 증가하였다. 1사분면 상단의 (+)점과 하단의 (-)점의 개수 차이도 단기간 예측에서 제일 높았다. 이것이 가을 시기의 최적화 과정에서 예측 기간이 짧을수록 AUC가 더 많이 증가할 수 있었던 이유이다(Fig. 4d와 Fig. 5d).

앞서 설명한 것처럼 계절별 최적화를 통해 대부분의 계절에서 예측 성능이 향상된 것을 확인하였다. 하지만 여름철의 특성을 고려한 이후, 기존과 대비하여 ACI_SnOpt가 과소 예측하면서 오히려 성능이 소폭 감소하는 결과가 발생했다. 이는 위험 범위를 기반으로 구축된 M F과 관련이 깊다. 여름철의 경우 대류의 발생 강도와 빈도가 매우 높기 때문에, 다른 계절에 비해 CAPE와 APCP의 위험 범위는 매우 높고 OLR은 매우 낮은 값의 구간에서 설정되었다(Fig. 2) 특히, JJA 기간의 위험 범위는 파트 1에서 설명된 1년 기간에 대한 위험 범위와 차이가 크기 때문에, 같은 입력 변수의 값이라도 1년 기간의 M F보다 여름철을 중심으로 설정된 M F에서 더 낮은 잠재도 값이 산출될 수 밖에 없다. 이러한 요인으로 인해 ACI_SnOpt는 기존의 ACI_YrOpt와 비교하여 과소 예보하면서 AUC가 소폭 감소하였다. 그럼에도 불구하고 다른 계절에서는 모두 예측 성능이 향상되면서 전체 계절 평균 AUC는 0.719로 가장 높았다(Table 5). 평균 TSS도 유사한 결과를 보였다. ACI_NoOpt의 계절 전체 평균 TSS는 0.207이지만, ACI_YrOpt는 0.195로 오히려 감소하였다. 이는 마찬가지로 다른 계절 평균값은 증가하였지만, 겨울철의 급격한 감소가 원인이었다. 이후, 계절 변동성을 고려하면서 DJF 시기의 TSS는 급격히 증가하였고 M AM과 SON 시기에도 TSS는 각각 0.180와 0.282로 ACI_NoOpt와 ACI_YrOpt에 비교하여 예측 성능이 꾸준히 향상되었다. 여름철에는 TSS 또한 소폭 감소함에도 불구하고 다른 계절의 큰 향상으로 인해 전체 계절 평균 TSS는 0.245로 가장 높은 결과를 보여주었다. 즉, 정리하면 계절 평균 AUC와 TSS를 비교하였을 때, 본 연구를 통해 개발된 다양한 형태의 ACI 중에서 계절 변동성을 고려한 ACI_SnOpt가 가장 높은 예측 성능을 보였다. 특히, ACI_SnOpt의 AUC와 TSS는 기존의 ACI_YrOpt와 대비하여 각각 0.983%, 그리고 25.641% 향상되었다.

4.1 사례1 - 2021.12.16.0000 UTC

첫 번째 사례는 2021년 12월 16일 0000 UTC에 제주도 및 남부 지역을 중심으로 발생한 집중호우 사례이다. 당시 종관 기압장을 보면, 한반도를 중심으로 서해에는 저기압, 그리고 남동쪽에는 고기압이 위치해 있다(Fig. 7a). 이러한 기압 배치에서 남해 부근에는 남서풍 기류가 발달하면서, 동중국해에 위치한 지상 저기압이 빠르게 한반도를 향해 이동하였다. 이로 인해 제주도를 포함한 남부 내륙에는 넓은 지역에 강수가 발생한 것을 시간당 강수량에 대한 레이더 Plan Position Indicator (PPI) 영상에서도 확인할 수 있다(Fig. 7b). 또한, 한반도로 이동하는 지상 저기압과 고기압 사이로 강한 기압 경도력이 발생하여 저기압 전면에는 강한 하층 제트가 발달하였다. 이 강한 남풍으로 인해 많은 양의 수증기가 한반도 남해안으로 유입되었다(Fig. 7c). 이 하층제트 뿐만 아니라 서해의 저기압성 흐름으로 인한 북풍, 그리고 내륙의 지형적 요인으로 인해 한반도 남해안에는 국지적으로 강한 하층 수렴대가 형성되었다(Fig. 7d).

Fig. 7.

(a) Analysis synoptic map at the surface, (b) Hourly precipitation of radar Plan Position Indicator (PPI) image, KIM-GDAPS Auxiliary analysis chart of (c) 850 hPa Moisture flux, and (d) 850 hPa Convergence at 0000 UTC 16 December 2021 from the Korea Meteorological Administration.

해당 사례에서 15 dBZ ETH를 보면, 강수가 발생한 남해 및 남부 일부 지역에 대류가 존재하는 것을 알 수 있다(Fig. 8a). 내륙의 ETH는 FL200까지 나타났으며, 강수량이 높았던 남해 일부 지역에서는 FL250까지 대류가 발생하였다. 본 연구에서 최종적으로 예측하고자 하는 FL250 이상의 15 dBZ ETH 영역을 그림에서 검은색 선으로 표현하였고, 이를 깊은 대류 영역(deep convective area)으로 정의하였다. 반대로 예측장에서 설명하는 깊은 대류 영역은 그림에서 보라색으로 표현하였다. ACI_YrOpt의 24시간 예측 결과에서는 깊은 대류의 공간 분포가 관측과 매우 다르게 나타났다. 남해에는 FL250 이상의 강한 대류를 예측하였지만, 실제로 관측되지 않았던 한반도 내륙의 중부 지역에도 FL200 정도의 추가적인 대류를 예측하였다(Fig. 8b). 이로 인해 전반적인 대류의 분포에서 관측과 많은 차이가 나타났다. 하지만 계절 변동성이 고려된 ACI_SnOpt의 24시간 예측 결과에서는 관측장과 유사하게 남해 및 남해안에 중심을 둔 강한 대류 시스템만 예측하면서, 중부 지역의 큰 오차를 상당히 개선시켰다(Fig. 8c). 하지만 남해 지역을 중심으로 보면, ACI_SnOpt가 오히려 ACI_YrOpt보다 대류를 더욱 과대 예측하는 것을 알 수 있다. ACI_SnOpt의 예측 결과에서는 FL300 이상 고도까지 발달한 대류가 관측보다 더 넓은 지역에 발달하였다. 이렇게 관측과 더불어 기존의 ACI_YrOpt보다 과대 모의하는 특징은 앞서 Fig. 6a 통해 이해한 통계적 특성과 동일하다. 앞선 통계적 분석 결과에 따르면, DJF 기간에 맞는 M F와 가중치 조합을 적용시키면서 전반적인 예측 성능은 크게 향상되었다. 이 과정에서 PODY 뿐만 아니라 POFD도 증가하였다. 본 사례에서도 마찬가지로 ACI_SnOpt 예측장에서는 보라색으로 표현된 예측된 깊은 대류 영역이 검은색의 관측장에서의 깊은 대류 영역을 모두 설명하고 있다. 반면에 ACI_YrOpt에서는 거의 설명이 되지 못했다. 이를 통해 PODY의 증가를 알 수 있다. 하지만, ACI_SnOpt예측장의 보라색 영역이 실제 깊은 대류 영역보다 넓은 영역을 예측하며 POFD 또한 증가한 것을 보여준다.

Fig. 8.

(a) 15 dBZ ETH, 24-hour forecast of (b) ACIYrOpt and (c) ACISnOpt at 0000 UTC 16 December 2021. Grey and black line indicate the area of 15 dBZ ETH over FL0 and the deep convective area (15 dBZ ETH ≥ FL250) in (a), that is also expressed in (b) and (c). Purple line indicates the deep convective area (≥ FL250) in forecast field of (b) and (c).

앞서 언급한 것처럼 본 연구에서 ROC 곡선을 사용한 최적화 및 검증은 15 dBZ ETH가 FL0을 넘는 영역에 대해서 수행되었다. 해당 영역은 Fig. 8에서 회색으로 표현되었다. 본 사례에서 나타난 중요한 특징 중 하나는 최적화 활용된 FL0을 넘는 영역뿐만 아니라 그 바깥 지역에서도 공간 분포의 개선이 뚜렷하게 나타났다는 것이다. 이러한 대류 공간 분포에 대한 개선 과정을 자세히 살펴보기 위해, 입력 변수로부터 산출된 최종 잠재도를 비교해보았다. 여기서 최종 잠재도는 M F로부터 산출된 값에 가중 상수가 곱해진 값을 의미하며, 그 합이 최종적인 ACI 결과와 동일하다. ACI_YrOpt와 ACI_SnOpt에서 산출된 CAPE의 최종 잠재도를 비교해보면, 전반적인 공간 분포는 유사하지만 ACI_YrOpt에서 잠재도가 더 좁은 영역에서 높은 값을 보인다(Figs. 9a, d). 이러한 차이는 가중치 조합과 멤버십 함수가 모두 다르기 때문이다. 각각의 MFCAPE 범위만 고려해본다면, ACI_SnOpt에서 설정된 범위가 상대적으로 더 낮은 값에서 정의되었기 때문에 계산 과정에서 산출되는 CAPE의 최종 잠재도는 더 높아야 한다. 그럼에도 불구하고 ACI_YrOpt가 높은 이유는 장기간에 대한 최적의 가중 상수가 0.45로, 계절에 맞는 가중치보다 더 높게 설정되었기 때문이다. APCP의 경우, ACI_YrOpt 계산 과정에서 가중치가 0.05로 매우 낮게 설정되면서 최종 잠재도가 매우 낮게 산출되었다(Fig. 9b). 하지만 계절 변동성을 고려한 ACI_SnOpt에서는 상대적으로 작은 값에서 MF가 구축되었을 뿐만 아니라, 가중치가 0.7로 조정되었기 때문에, ACI_SnOpt의 예측 결과에서 0.65 이상의 잠재도가 남해의 넓은 지역에 예측되었다(Fig. 9e). ACI_SnOpt의 예측 결과에서 대류의 공간 분포가 크게 개선된 이유는 OLR의 최종 잠재도의 변화 때문이다(Figs. 9c, d). 먼저, ACI_YrOpt 계산 과정에서 산출된 OLR 예측장을 보면 기본적으로 깊은 대류의 발생 위치를 잘 설명하지 못하는 것을 알 수 있다. 하지만 장기간에 대한 최적의 가중치는 0.5로 설정되었기 때문에 이러한 잘못된 예측 결과가 많이 반영되었다. 이로 인해 최종적인 ACI_YrOpt 예측 결과에서 관측에는 발생하지 않았던 한반도 중부 지역의 대류를 잘못 예측하였다. 하지만 계절별 최적화 이후, MF 범위는 유사하지만 가중치는 0.1로 조정되면서 OLR의 최종 잠재도는 ACI_SnOpt 계산 과정에서 거의 나타나지 않았다. 이를 통해 중부 지역에 잘못 예측된 대류 분포를 상당히 보정할 수 있었다.

Fig. 9.

(a) αMCAPE, (b) βMAPCP, (c) γMOLR of ACIYrOpt, and (d) αMCAPE, (e) βMAPCP, (f) γMOLR of ACISnOpt. Bracket refers to the range of MFs.

4.2 사례2 - 2021.05.28.0300 UTC

해당 사례는 중국의 산동 지역에 지면 가열로 인해 발생한 저기압이 발달과 동시에 동진하며 한반도를 관통하여 집중호우를 발생시킨 사례이다. 집중호우 발생 시간 기준 3시간 이전의 지상 분석 일기도에서는 상층의 뚜렷한 기압골과 강한 제트 기류가 존재한다(Fig. 10a). 하지만 지상 저기압의 중심은 기압골의 우측, 혹은 제트 기류 입구의 우측에 위치하지 않기 때문에, 상층 시스템으로 인한 비-지균 상승의 직접적인 영향은 있다고 보기 어렵다. 실제로 보조 분석도의 상층 발산에서는 상층 기압골과 제트 기류로 인한 넓은 발산 구역과 서해안에 지상 저기압으로 인한 국지적인 발산 구역은 뚜렷하게 구분된다(Fig. 10c). 서해안에 위치한 저기압은 한반도로 동진하며, 주변에는 강한 하층 제트가 발달하였다. 이로 인해 저기압 전면에는 선형 형태의 강한 강수대가 발달하였으며, 시간이 지나면서 더욱 발달하여 위성 영상에서는 뚜렷한 콤마 형태의 구름 구조로 강화되었다(Fig. 10b). 또한, 해당 사례가 발생한 시기에는 서해안에는 국지적으로 높은 불안정도가 존재한 것이 특징이었으며, 이는 쇼월터 안정도 지수(Showalter Stability Index)에 대한 보조 분석도에서 확인할 수 있다(Fig. 10d).

Fig. 10.

(a) Analysis synoptic map at the surface, (b) GK2A RGB image, KIM-GDAPS Auxiliary analysis chart of (c) 300 hPa Divergence and wind, and (d) 850-500 hPa Showalter Stability Index at 0000 UTC 28 May 2021 from the Korea Meteorological Administration.

해당 사례는 선형 강수대가 한반도를 지나는 과정에서 항공기가 대류를 피하지 못하고 관통하였던 사례이다. 이로 인해, 항공기는 CIT를 조우하였다. 항공기 관측 자료인 Aircraft Meteorological Data Relay(AMDAR)에 기록된 순간 최대 난류 강도는 0.21m^2/3s^-1였으며 난류강도는 중간 강도로써 상공 약 6.28 km에서 관측되었다. 발생 당시인 2021년 5월 28일 0300 UTC 때의 15 dBZ ETH를 보면 FL250 이상으로 깊은 대류의 중심에서 파란색 점으로 표시된 CIT가 발생한 것을 알 수 있다(Fig. 11a). 앞선 사례와 마찬가지로 관측장에서 FL250 이상의 깊은 대류 영역은 검은색 선으로 표현되었으며, 예측된 깊은 대류 영역은 예측장에서 보라색으로 표현되었다. ACI_YrOpt와 ACI_SnOpt의 21시간 예측 결과에서는 관측에서 나타난 대류의 전반적인 공간 분포가 잘 설명되었다(Figs. 11b와 c). 하지만 CIT 발생 지점을 중심으로 보면, ACI_YrOpt은 관측과 비교하여 대류 강도를 약하게 예측한 반면에 ACI_SnOpt 예측 결과에서는 ETH를 FL250 이상으로 대류의 강도를 관측과 유사하게 예측되었다. 뿐만 아니라, ACI_YrOpt 예측장에서 보라색으로 표시된 깊은 대류 영역은 검은색으로 표시된 관측 자료의 깊은 대류 영역을 모두 설명하지 못하였지만, ACI_SnOpt를 통해 예측된 깊은 대류 영역은 실제 영역을 모두 설명하였다. 이는 곧 PODY의 증가를 의미하며, 앞서 Fig. 6b를 통해 분석한 M AM 기간의 최적화 과정에서 보인 통계적 특성과 동일하다. DJF와 마찬가지로 해당 계절에는 POFD와 PODY가 모두 증가하였지만 PODY의 증가폭이 더욱 크면서 예측 성능이 향상되었다. 본 사례에서도 기존의 ACI_YrOpt보다 과대 예측하면서 관측된 대류영역을 모두 예측하긴 했지만, 관측되지 않은 영역 또한 더 많이 예측하면서 POFD도 증가하는 것을 보여준다.

Fig. 11.

(a) 15 dBZ ETH, 21-hour forecast of (b) ACIYrOpt and (c) ACISnOpt at 0300 UTC 28 May 2021. Grey and black line indicate the area of 15 dBZ ETH over FL0 and the deep convective area (15 dBZ ETH ≥ FL250) in (a), that is also expressed in (b) and (c). Purple line indicates the deep convective area (≥ FL250) in forecast field of (b) and (c). Here, the intensity of Convectively-Induced Turbulence is indicated by the dots with colors in (a). And the location of CIT encounter is indicated by the blue circle in (b) and (c).

앞서 해당 사례의 쇼월터 안정도 지수에 대한 보조 분석도에서 확인한 것처럼 대류 발달 당시 불안정도가 높은 환경이었다(Fig. 10d). 마찬가지로 ACI_YrOpt와 ACI_SnOpt에서 계산된 CAPE의 잠재도 모두 대류가 존재하였던 한반도 내륙에서 매우 강하게 나타났으며, 다른 변수와 비교하여 최종적인 ACI에 가장 많은 기여를 한 것을 알 수 있다(Figs. 12a, d). 계절별 최적화 전후 CAPE의 가중 상수는 크게 조정되지 않았다. 하지만 M F의 범위가 더 낮게 설정되면서 ACI_SnOpt에서 CAPE의 최종 잠재도는 기존의 ACI_YrOpt 보다 높게 나타났다. 이러한 변화는 해당 사례에서 PODY를 증가시키는 데 직접적인 기여를 하였다. APCP의 예측장에서도 선형 형태의 대류가 잘 나타났다. 하지만 ACI_YrOpt 계산 과정에서는 APCP의 최종 잠재도가 매우 낮게 고려되었다(Fig. 12b). 즉, 해당 변수의 좋은 예측 결과가 거의 반영되지 못했다. 이러한 이유는 장기간에 대한 최적화 과정을 거치면서 강수량에 대한 가중 상수는 0.05로 매우 낮게 설정되었기 때문이다. 하지만 계절적 특성이 고려되면서, 가중치가 0.15로 증가하였다. 이를 통해 ACI_SnOpt 계산 과정에서는 APCP의 잠재도가 더 높게 산출되면서, 좋은 예측 결과가 최종적인 ACI 결과에 더 많이 반영될 수 있었다(Fig. 12e). 이것이 ACI_SnOpt가 ACI_YrOpt보다 더 강한 대류를 예측하여 실제 관측과 더 유사해 질 수 있었던 원인이었다. 반면에 ACI_SnOpt와 ACI_YrOpt에서 OLR은 전반적으로 깊은 대류에 대한 공간 분포를 잘 설명하지 못하였다(Figs. 12c, f). 영동 지역을 포함한 동해안 일 부와 긴 형태의 상층운이 존재했던 남해 일부 지역에서 대류를 예측하였지만, CIT를 발생시킨 한반도 중부 지역의 깊은 대류를 전혀 예측하지 못하였다. 그러므로 깊은 대류를 예측하는 데 있어서 OLR은 적게 고려될수록 유리하다. 두 가지의 ACI를 비교해보면, ACI_SnOpt에서 더 작게 산출된 것을 알 수 있다. 계절별 최적화 전후 M F 범위의 차이는 크지 않지만, 기존 ACI_YrOpt의 OLR 가중 상수는 0.5인 반면 ACI_SnOpt에서는 0.35로 더 낮게 설정되었다. 이를 통해 불필요한 예측 결과를 더 적게 고려할 수 있었다.

Fig.12.

(a) αMCAPE, (b) βMAPCP, (c) γMOLR of ACIYrOpt, and (d) αMCAPE, (e) βMAPCP, (f) γMOLR of ACISnOpt. Bracket refers to the range of MFs.