기상청은 황사 감시 및 예보를 위하여 황사·연무통합예측모델 (ADAM3)을 매일 4회 (00, 06, 12, 18 UTC) 현업 운영하고 있다. ADMA3은 미국환경보호국에 서 개발된 대기화학 수송모델인 CMAQ (Community Multiscale Air Quality Modeling System)에 황사발원알고리즘을 통합하여 개발되었으며 황사 발원 모듈, 수송 및 확산모듈, 건·습성 침착모듈 등으로 구성되어 있다 (Park et al., 2010; Park and In, 2003). 모델 영역은 한반도, 일본과 황사 발원지를 포함하는 동아시아지역 (20~60°N, 70~150°E)으로 수평해상도는 25 km이고 연직층은 100 hPa까지 49개 층으로 구성되어 있으며 초기 시간에 따라 168시간 (00, 12UTC)과 72시간 (06, 18UTC)의 예측 시간을 갖는다 (NMC, 2021, 2020). 기상입력자료는 영국기상청의 통합모델인 UM (Unified Model)을 기반으로 기상청에서 생산되는 GDAPS (Global Data Assimilation and Prediction System) 자료를 이용하고, 지상 에어로졸 (PM₁₀) 농도의 실시간 예측오차를 반영하기 위해 앙상블 기반 3DAVR (Three dimensional VARiational data assimilation) 체계를 적용하고 있다 (NMC, 2020). ADAM3의 황사 발원지는 세계기상기구에 보고되는 종관기상실황 관측자료를 활용하여 황사가 보고된 영역으로 정의되었고, 황사 발원지의 토양 유형은 중국 토양도와 몽골 수문기상국의 토양정보를 이용하여 고비 (Gobi), 모래 (Sand), 황토 (Loess), 혼합 (Mixed)의 네 가지로 구분하고 있다 (Lee et al., 2012; Park et al., 2010; In and Park, 2002; Park, 2002). ADAM3에서 황사농도 (D)는 각 토양 유형별로 산정된 풍속 (WS₁₀)과 상대습도 (RH_1.5)의 각 임계 조건 (WS_10t, RH_1.5t)을 충족하면, 식 (1)의 황사 발원알고리즘을 통해 산정된다 (Part et al., 2010; Park and In, 2003).

여기에서, Rf와 PRf는 각각 식생 비율과 과거 강수에 의한 황사 발원 억제지수를 의미하고 u*4는 마찰속도를 의미한다. 황사 발원량은 마찰속도의 4제곱 (u*4)에 비례하여 산정되고 위성에서 관측된 NDVI (normalized difference vegetation index) 자료를 이용하여 시공간적으로 변화하는 식생 특성을 반영하며 (Cho et al., 2021; Park et al., 2010), 예측시점 이전 (3~6시간)에 황사 발원지에서 발생한 강수량을 고려하여 황사 발원량을 조절한다 (Lim et al., 2021). 즉, 식생이 증가할수록 강수량이 많을수록 황사 발원량이 감소하도록 설계되어 있다. ADAM3의 자세한 환경 설정은 표 1에서 확인 가능하다.

기상청의 황사 일기도와 위성자료를 활용하여 황사의 발원과 한반도로의 유입과정을 분석하였다. 황사 일기도는 GTS (Global Telecommunication System)를 통해 보고되는 육상 종관기상실황 전문자료 중 먼지현상을 일기도에 나타낸 것으로 1일 8회 3시간 간격으로 자료 활용이 가능하다. 위성자료는 GK2A (Geostationary Korea Multi-Purpose SATellite 2A) RGB 황사 영상을 이용하였다. 이는 GK2A에 탑재된 AMI (Advanced Meteorological Imager)에서 제공되는 11.2~10.5 μm 파장을 사용하여 황사를 탐지하는 위성영상으로 황사를 옅은 진분홍, 연분홍, 자홍, 보라색으로 표출하고 있으며 황사의 이동 분석 및 모델 평가 등에 활용할 수 있다 (Jang et al., 2022).

황사의 발원과 국내에 유입된 황사의 농도 및 영향 범위를 파악하기 위하여 황사감시기상탑이 설치되 어 있는 몽골 놈곤 (Nomgon) 지점과 중국에 위치한 황사 감시망 (운영장비: Tapered element oscillating microbalace 1400a 측정기)에서 관측한 PM₁₀ 농도자료를 이용하였고, 국내 기상청 지상 감시망 (운영장비: β-ray 부유분진측정기) 27개 지점에서 관측한 지상 PM₁₀ 농도 자료를 활용하였다. 또한, 백령도 및 서울에서 관측된 광학입자계수기 (Optical Particle Counter) 자료를 이용하여 한반도에 유입된 황사의 입자 크기별 부피농도를 분석하였다. 광학입자계수기는 직경 0.3~25 μm 입자의 산란광의 강도를 측정하여 입자의 크기별 개수를 측정하는 기기로 청정지역의 대기 중 부유하는 입자의 측정에 용이하다고 알려져 있다 (Chun et al., 1999). 마지막으로 황사 사례의 시·종 시각은 목측관측지점에서 황사 발생 시 황사 현상의 부호, 방향, 관측 시각 등을 기록한 황사 기사자료를 활용하였다. 본 연구에서 이용한 놈곤 지점과 중국 및 기상청의 황사 감시지점의 위치는 그림 2와 표 2를 참고하면 된다.

Fig. 2. 
ADAM3 domain (left) and geographical location of the employed stations. Detailed information of the dust monitoring tower in Nomgon of Mongolia (green triangle; A) and KMA-CMA joint monitoring stations (blue squares; a~g) and KMA (dark red circles; 1~27) are presented in Table 2.

황사 발원지에서의 기상 특성은 NCEP/NCAR R1 (National Centers for Environmental Prediction/National Center for Atmospheric Research Reanalysis 1; Kalnay et al., 1996) 재분석자료의 500 hPa 지위고도, 850 hPa 바람, 2 m 기온, 누적 강수량 자료를 이용하여 파악하였다. NCEP R1 자료는 약 2.5°×2.5°의 수평격자해상도와 6시간의 시간해상도를 가진다. 본 분석에는 기후기간과 아노말리 분석을 위해 1991년부터 2020년까지 자료와 2022년 12월의 일 및 월 평균 자료를 이용하였다.

기상청의 ADAM3은 황사 및 연무예측을 위해 매일 4회 (00, 06, 12, 18 UTC) 운영되고 황사 PM₁₀ 및 비황사 PM₁₀ 농도를 산출하고 있다. 2022년 12월 12일부터 14일까지 우리나라에서 관측된 고농도 황사 사례에 대하여 초기 시간별로 생산된 예측 황사 PM₁₀ 자료와 위성자료 (RGB 황사영상)를 비교하여 ADAM3에서의 황사 발생 및 이동 등의 공간분포를 정성적으로 평가하였다. 또한, ADAM3에서 산출된 황사 PM₁₀과 비황사 PM₁₀의 합인 총 PM₁₀의 평가를 위해서 지점 별 관측 PM₁₀ 농도와 입자 크기별 부피농도를 이용하여 비교하였다. 더불어, 초기 시간별로 갱신되는 총 PM₁₀의 예측성을 정량적으로 평가하기 위하여 평균편향오차 (Mean Bias Error, MBE), 평균제곱근오차 (Root Mean Square Error, RMSE)를 이용하였다. MBE (식 (2))와 RMSE (식 (3))는 예측값과 관측값의 차이를 나타내는 척도로 0에 가까울수록 예측값과 관측값의 수치가 서로 유사하다고 평가된다.

여기에서 N은 황사 발원지와 우리나라에서의 PM₁₀ 지상관측소 수이고, i는 예측 시간이며, f_i와 o_i는 각각 예측 시간 (i)에서의 예측값과 관측값을 의미한다.

기상청 목측관측지점을 기준으로 2022년 12월 12일 (최초 관측시점/지점: 1930 KST/백령도)부터 14일 (최종 관측시점/지점: 0115 KST/울산)에 관측된 황사의 발원과 이동을 분석하기 위해 그림 3에 황사 일기도와 GK2A RGB 황사영상을 나타내었다. 12월 11일 09 KST (Korea Standard Time) 500 hPa 일기도의 몽골 지역에 위치한 상승 기압골이 발달함에 따라 (그림 생략) 몽골 동부지역에 하층 저기압이 생성되었고 저기압 후면과 한랭한 고기압 전면에 기압경도력이 강화되었으며 몽골 중부와 중국 북부지역인 고비사막 및 내몽골 고원 일대에서 황사가 발원하였다. 황사 일기도에서 볼 수 있듯이 11일 12 KST와 12일 00 KST에는 등고도면이 조밀하여 강풍으로 인해 황사가 발원하기 좋은 조건이었다 (그림 3a, c). 또한, 황사 발원 시 850 hPa 고도에서 저기압 서쪽과 중심에 각각 온 도골과 온도능이 위치하였다 (그림 생략). Kim et al. (2022)의 연구에 따르면 강화된 온도경도를 가진 동서방향의 온도골이 남북으로 형성된 등고도선과 교차함에 따라 강한 경압불안정을 야기한다고 하였다. GK2A RGB 황사영상에서도 몽골 남부를 시작으로 황사 (분홍색)가 탐지되어 고비사막 및 내몽골 고원 전역에서 황사 (자홍색)가 발원한 것으로 확인되었다 (그림 3b, d).

Fig. 3. 
Surface weather charts for Asian dust (left), where the red weather symbols represent the Asian dust observed at surface synoptic observation stations, and GK2A dust RGB images (right), where pink color represents the Asian dust, at 1200 KST on 11 (a, b), at 0000 and 2100 KST on 12 (c~f), and at 1200 KST on 13 December 2022.

NCEP/NCAR R1 재분석자료를 이용하여 황사 발원지에서의 500 hPa 지위고도, 850 hPa 바람, 2 m 기온, 누적강수량의 평년 편차를 살펴보면 (그림 4), 황사 발원일 (12월 11일)에는 황사 발원 전 일주일 (12월 5일부터 11일)에 비해 몽골 남부 및 중국 북부지역의 황사 발원지에서 강한 저기압 편차가 형성되었고, 황사 발원지에서부터 한반도까지 평년보다 강한 북서기류가 형성되어 있는 것을 알 수 있다. 이는 황사 발원지에서 상승기류에 의해 대기상층으로 황사가 확산되고 북서풍을 타고 우리나라에 영향을 줄 수 있음을 의미한다 (그림 4a, d). 황사 발원지에서 기온은 황사 발원 전 일주일 (12월 5일부터 11일)에 비해 황사 발원일 (12월 11일)에 2~7°C의 강한 양의 편차를 형성하였며, 황사가 발원하기에 충분히 온난한 상태였던 것으로 보인다. 강수량은 평년과 비슷한 수준으로 나타났다 (그림 4b, e).

Fig. 4. 
NCEP/NCAR R1 500 hPa Geopotential height with 850 hPa wind (a, d), 2 m air temperature (b, e), and accumulate precipitation (c, f) anomalies for 5 to11 December 2022 (a~c) and 11 December 2022 (d~f). The climate period is from 1991 to 2020. Circles indicate the source region of the Asian dust.

그림 5는 황사 발원지인 몽골 남부에 위치한 놈 곤 지점에서 관측된 PM₁₀ 농도와 풍속을 나타낸 것이다. 놈곤에서는 11일 15 KST부터 12일 03 KST까지 9 m s^-1 이상의 강풍이 발생하였는데 이는 Park et al. (2010)이 제시한 12월의 황사 발원 풍속 임계값 (7.5 m s^-1) 이상의 풍속을 나타낸다. 강풍이 발생하기 시작한 후부터 지상 PM₁₀ 농도가 증가하기 시작하였고 11일 20 KST부터 1시간 평균 지상 PM₁₀ 농도가 300 μg m^-3 이상으로 관측되었으며, 최고 2,714 μg m^-3 (11일 23 KST)까지도 나타났다. 비슷한 시각에 중국 북부지역인 고비사막 및 내몽골 고원에 위치한 엘렌하오터와 우라터중치 지점에서도 12~15 m s^-1의 강풍과 함께 1시간 평균 지상 PM₁₀ 농도가 800 μg m^-3 이상의 매우 짙은 황사가 관측되었다 (그림 생략). 얼렌하오터와 우라터중치에서 1시간 평균 지상 PM₁₀ 농도는 각각 1,136 μg m^-3 (11일 22 KST)와 4,143 μg m^-3 (11일 23 KST)까지 발생하였다. 그 후 만주지역에 위치한 지점인 츠펑, 통랴오, 쓰핑 순으로 각각 1,083 μg m^-3 (12일 08 KST), 1,084 μg m^-3 (12일 13 KST), 890 μg m^-3 (12일 17 KST)의 PM₁₀ 최고농도가 나타났다. 11일부터 12일까지 내몽골 고원에서 발원한 황사는 대륙고기압의 확장과 함께 저기압 후면을 따라 발해만 방향으로 형성된 북서풍을 따라 남동진하였고 12일 15 KST에 화베이 지역에서 황사가 관측되었으며, 12일 21 KST경에 백령도를 시작으로 한반도 전역에 영향을 미쳤다 (그림 3e~h). 12일 21 KST와 13일 12 KST의 황사 일기도에서 만주지방과 동해에 위치해 있던 저기압과 고기압 사이의 기압경도력이 더 커지고 등고선이 더 조밀해진 것으로 보아 저기압 주위의 풍속이 증가하였고 이로 인해 저기압이 빠르게 북동진한 것으로 보인다. 12일 1930 KST 남한의 최북단 섬인 백령도에서부터 시작된 황사 (목측 기준)는 13일 00 KST경부터 서울과 서해안 일부 지역에 영향을 주었다. 그림 6에는 한반도 평균 (27개 지점 평균) PM₁₀ 농도와 300 μg m^-3 이상의 농도가 나타난 지점의 1시간 평균 PM₁₀ 농도를 나타내었다. 황사 영향을 받는 기간 동안 한반도는 평균적으로 최고 192 μg m^-3의 PM₁₀ 농도가 관측되었다 (그림 6a). 300 μg m^-3 이상의 PM₁₀ 농도는 12일 22 KST에 백령도에서 가장 먼저 관측되었고 이후에는 서울 (13일 04 KST), 속초·북춘천·군산·강화 (13일 11 KST), 전주 (13일 12 KST), 진도 (13일 15 KST) 순으로 관측되었다 (그림 6b~f). 백령도와 서울에서는 각각 472 μg m^-3 (13일 05 KST), 480 μg m^-3 (13일 11 KST)의 최고 PM₁₀ 농도를 보였다 (그림 6b~f).

Fig. 5. 
(a) PM₁₀ mass concentration (μg m^-3) and (b) 2 m wind speed (m s^-1) at the dust monitoring tower in Nomgon of Mongolia from 10 to 13 December 2022.

Fig. 6. 
Time series for PM₁₀ mass concentration at surface observation sites in South Korea from 11 to 14 December 2022.

그림 7은 황사 사례 중 최고 PM₁₀ 농도가 관측되었던 백령도와 서울에서 관측된 입자 크기별 부피농도 자료를 나타낸 것이다. 백령도의 경우, 지상 PM₁₀ 농도가 증가하기 시작한 12일 15 KST부터 21 KST까지는 입경이 0.6 μg 이하의 초미세먼지가 증가하였으며 (그림 7a), 이후 1시간 평균 PM₁₀ 농도가 472 μg m^-3 (13일 06 KST)으로 증가했고 13일 09 KST까지 200 μg m^-3 이상의 고농도 지상 PM₁₀이 관측되었다 (그림 6). 이때 증가한 PM₁₀ 농도는 입경이 약 4 μm (조대입자)에서 최대로 나타났다. 서울의 입자 크기별 부피농도 자료에서는 13일 00 KST에서 16 KST 사이에 조대입자 (직경 2.5 μm 이상) 영역의 부피농도가 높게 나타났으며 (그림 7b), 13일 11 KST에 1시간 평균 PM₁₀ 농도가 480 μg m^-3까지 관측되었다. 한반도 전역에 영향을 미친 황사는 마지막으로 14일 0115 KST에 울산에서 목측으로 관측되었다.

Fig. 7. 
Time series for particle size distributions observed by the optical particle counter at (a) Baengnyeongdo and (b) Seoul in 11~14 December 2022.

우리나라에 영향을 미친 황사를 정확하게 예측하기 위해서는 우선 모델이 발원지에서의 황사 (총 PM₁₀ 농도)를 얼마나 정확하게 모의했는지 평가해야 한다. 또한, 동일한 목표 시각에 대한 모델 황사농도 예측의 변화를 파악할 수 있도록 (Lee et al., 2012) 매 초기 시간별 갱신되는 예측정보의 정확도가 어떻게 변화하는지 예측경향성에 대한 평가가 필요하다. 그림 8은 발원지에서의 황사예측과 이동을 공간적으로 검증하기 위해 초기 시간별로 ADAM3에서 예측된 PM₁₀ 농도와 GK2A RGB 황사영상을 나타낸 것이다. 모델 초기 시간은 10일 09 KST부터 12시간 간격으로 나타내었고 (그림 8 세로방향) 위성자료와 비교된 예측 시각은 황사가 발원하기 시작하였던 12월 11일 18 KST부터 6시간 간격으로 나타내었다 (그림 8 가로방향).

Fig. 8. 
Spatial distributions of total PM₁₀ concentration in East Asia with four runs by initial date: 09 KST and 21 KST 10 December, and 09 KST and 21 KST 11 December 2022. The GK2A RGB dust images for the forecast date are shown on the top.

ADAM3은 10일 09 KST 모의에서 12일 18 KST (+21 h)부터 12일 06 KST (+33 h)까지 몽골 남부, 고비사막 및 내몽골 고원 부근 발원지에서 발원한 황사를 모의하였고 북서기류를 따라 동남진하는 황사를 모의하고 있다. 모델에서 모의된 황사 공간분포는 GK2A RGB 황사영상에서 탐지된 황사 (주황, 자색, 분홍 등)의 공간분포와 비슷하게 나타났다. 서로 다른 초기 시간으로부터 모의된 모든 예측장이 발원지에서 황사 발원과 이동 패턴을 비슷하게 모의하고 있어 모델이 본 황사 사례를 안정적으로 예측한 것을 볼 수 있다.

황사의 이동속도, 최고강도, 발원지의 먼지 발생 등의 정보들이 최신 자료로 갱신됨에 따라 (Lee et al., 2012) 다른 초기 시간으로부터 모의된 황사 강도와 패턴도 서로 다르게 예측된다. 이러한 차이는 각 지점별로 예측된 PM₁₀ 농도 시계열 분석을 통해 살펴볼 수 있다. 그림 9에는 고농도 황사가 관측된 지점 (얼렌하오터, 통랴오, 위서)과 발원지 풍하측에 위치한 지점 (후이민)의 초기 시간별로 모의된 PM₁₀ 농도의 시계열을 제시하였다. 그림에서 회색 음영은 관측에서 100 μg m^-3 이상의 PM₁₀ 농도가 관측된 시점을 나타낸 것이다. ADAM3의 10일 21 KST의 모의에서 11일과 12일에 각각 우라터중치와 통랴오에서 관측된 1,000 μg m^-3 이상의 매우 높은 PM₁₀ 농도를 상당히 과소모의하였지만 모델의 초기 시간이 갱신됨에 따라 PM₁₀ 농도를 높게 모의하는 것을 볼 수 있다 (그림 9a, b). 또한, 황사의 시작 및 종료시각을 유사하게 모의하였다 (그림 생략). 그림 9c의 위서 지점의 경우, 10일 21 KST의 모의에서 최고농도 황사를 관측보다 약 100 μg m^-3 과대모의 (12일 11 KST)하였으나, 시간이 지남에 따라 PM₁₀ 농도를 낮게 모의하여 11일 21 KST 모의에서는 관측값과 거의 유사한 강도 (약 144 μg m^-3)로 모의하였다. 선행 시간이 짧아질수록 관측에 가까운 PM₁₀ 농도를 모의하였지만 13일에 발생한 최고농도 시점을 관측보다 약 5시간 늦게 모의하였다. 마지막으로, 후이민 지점에서는 초기 시간별 모의된 PM₁₀ 농도의 차이를 볼 수 없었지만, 11일 21 KST 모의에서 예측된 최고농도가 204 μg m^-3으로 관측된 최고농도 (249 μg m^-3)와 비슷하게 모의하였다 (그림 9d).

Fig. 9. 
Time series of observed PM₁₀ concentration at (a) Erenhot, (b) Tongliao, (c) Yushe, and (d) Huimin in Asian dust source regions initialized from different simulation starting time. Gray shading indicates the range of the time when PM₁₀ concentrations above 100 μg m^-3 were observed.

한반도에서 예측을 검증하기 위해 한반도에 황사가 유입되기 전부터 황사 영향이 벗어난 시각까지 12시간 간격으로 예측된 황사 PM₁₀ 농도를 초기 시간별로 나타내었고, 관측과 비교를 위해 동일한 예측 시간의 GK2A RGB 황사영상을 함께 나타내었다 (그림 10). 4개의 서로 다른 초기 시간으로부터 모의된 황사 PM₁₀은 예측 시간이 서로 다름에도 불구하고 북서풍을 따라 요동반도와 (12일 12 KST) 서해로 (13일 00 KST) 유입된 황사가 우리나라의 서해와 남해 등에 영향을 미치고 (13일 12 KST) 14일 00 KST에 한반도의 영향에서 벗어난 것을 비슷하게 모의하고 있어 모델이 본 황사 사례를 안정적으로 예측한 것을 알 수 있다. 네 초기 시간 모두 12시간 후인 14일 00 KST 모의에서는 한반도의 예측농도가 약화되어 80 μg m^-3 이하로 나타났고 남해상의 황사는 동진하여 일본에 도달한 것으로 모의하였다. 위성에서 탐지된 황사는 중국 허베이성을 거쳐 (12일 12 KST) 12일 13 KST에 발해만을 통과하여 서해로 유입되었고, 모델에서는 황사가 허베이성을 거쳐 (12일 12 KST) 13일 00 KST에 발해만보다 남쪽으로 치우친 산둥반도 쪽으로 황사가 유입된 후 서해와 영향을 미쳐 약간의 차이를 보였지만 초기 시간별로 ADAM3에서 모의된 황사는 위성영상에서 탐지된 황사와 전반적으로 비슷한 분포를 나타내었다 (그림 10). 그러나, 초기 시간별 동일한 예측 시각에도 그 강도와 패턴은 다소 달라지는 것을 볼 수 있다. 한반도에 황사가 도달한 13일 00 KST에 대한 모의를 비교해 보면 초기 시간 11일 09 KST와 21 KST에는 백령도 부근에서 PM₁₀ 농도가 약 150 μg m^-3으로 모의되었지만 초기 시간 12일 21 KST에는 농도가 점점 증가하여 약 250 μg m^-3 이상이 모의되었다. 또한, 12일 14 KST에 예측된 황사의 경우, 예측 시간이 짧아질수록 (초기 시간 12일 21 KST) 한반도에서 빠져나간 황사영역대가 더욱 넓어진 것으로 모의되었다.

Fig. 10. 
Spatial distributions of total PM₁₀ concentration in South Korea with four runs by initial date: 09 KST and 21 KST 11 December, and 09 KST and 21 KST 12 December 2022. The GK2A RGB dust images for the forecast date are shown on the top.

ADAM3에서 예측된 황사농도에 대하여 지점별로 살펴보고자 관측에서 300 μg m^-3 이상의 고농도를 보였던 백령도, 서울 지점과 서해 및 남해안에 위치한 흑산도와 안면도의 관측 및 예측 PM₁₀ 농도 시계열을 그림 11에 제시하였다. 최고농도가 나타났던 백령도와 서울 지점의 경우 11일 09 KST 모의에서 황사의 시종 시각은 관측과 유사하게 모의하였지만 최고농도를 과소모의하였다. 예측 초기 시간이 갱신됨에 따라 모의농도가 약간 증가하였지만 가장 최근의 초기 시간인 12일 21 KST의 모의에서도 백령도와 서울의 예측된 최고농도는 각각 211 μg m^-3, 132 μg m^-3으로 나타나 관측에 비해 PM₁₀ 농도를 각각 450 μg m^-3, 350 μg m^-3가량 낮게 모의하여 여전히 과소모의하였다. 이러한 황사 과소모의는 고농도 황사가 발생했던 발원지역 (우라터중치, 엘렌하오터 등)에서의 과소모의가 원인일 수 있다. 상대적으로 낮은 PM₁₀ 농도가 관측되었던 흑산도와 안면도의 경우 초기 시간에 따라 예측된 PM₁₀ 농도의 차이가 크지 않았고, 12일 21 KST 모의에서 PM₁₀ 농도가 각각 168 μg m^-3, 166 μg m^-3으로 예측되어 관측 (흑산도: 261 μg m^-3, 안면도: 219 μg m^-3)과 약 100 μg m^-3 이하의 차이를 보였다. 한반도에서의 예측된 황사는 발원지에서의 황사발생뿐만 아니라 발원지에서 한반도까지의 기류 예측 등이 복합적으로 작용되기 때문에 발원지에 비하여 불확실성이 커지는 것으로 보인다.

Fig. 11. 
Time series of simulated and observed PM₁₀ concentration at (a) Baengnyeongdo, (b) Seoul, (c) Heuksando, and (d) Anmyeondo in South Korea, initialized from different simulation starting time.

초기 시간에 따라 변화하는 백령도와 서울에서 산출된 입자 크기별 부피농도를 그림 12에 나타내었다. 입자별 에어로졸 부피농도는 ADAM3에서 직접적으로 예측되지 않으므로 Whitby and McMurry (1997)의 방법에 따라 예측된 에어로졸 질량농도 및 수농도를 이용하여 산정하였다. 11일 09 KST에 모의된 백령도에서의 입자 크기별 부피농도는 12일 20 KST부터 13일 00 KST까지 입자 직경 0.8 μg 이하의 미세먼지가 증가하기 시작하였고 이후에는 3~10 μg의 입자 직경에서 최대 부피농도를 보여 (그림 12a), 초미세입자와 조대입자의 부피농도를 과소모의했지만 입자별 부피농도 분포는 관측과 유사하게 모의되었다 (그림 7a 참고). 서울에서도 11일 09 KST에 모의 결과에서 입자 직경별 부피농도를 관측보다 과소모의하였지만 조대입자 (5~10 μg) 부피농도가 증가하는 시각 (13일 00 KST)은 관측과 유사하게 모의하였다. 백령도와 서울 모두 초기 시간에 따른 입자별 부피농도 분포 결과는 비슷하게 나타났다.

Fig. 12. 
Time series for particle size distributions simulated at (a) Baengnyeongdo and (b) Seoul of South Korea, initialized from different simulation starting time.

황사가 발원하고 이동하여 한반도에 영향을 준 기간에 대하여 황사 발원지와 한반도에서의 MBE와 RMSE를 초기 시간에 따라 6시간 간격으로 표 3에 제시하였다. 초기 시간에 따라 72시간 (1~72 예측 시간)의 예측자료를 이용하였고 황사 발원지와 한반도 모든 지점에 대한 평균값을 나타내었다. 모든 초기 시간에 대한 평균 MBE는 황사 발원지와 한반도에서 각각 -16.50 μg m^-3과 -25.89 μg m^-3으로 나타났다. 선행 연구에 따르면 황사 발원지에서 ADAM3는 봄철에 PM₁₀ 농도를 관측에 비해 과대모의하는 것으로 알려져 있으나 (Ryoo et al., 2020), 본 황사 사례가 포함된 기간 동안에는 PM₁₀ 농도를 전반적으로 과소모의하였다. 모든 초기 시간에 대한 평균 RMSE는 황사 발원지와 우리나라에서 각각 137.76 μg m^-3, 35.95 μg m^-3 수준으로 평가되어 지상 PM₁₀ 농도가 높게 발생했던 황사 발원지에서의 오차가 더 크게 나타났다. 두 지역 모두 초기 시간에 따라 MBE의 차이가 거의 나타나지 않았고, RMSE는 최근의 초기 시간자료를 이용했을 때 작은 오차를 보이긴 했지만 그 차이는 작게 나타났다.