X-Net-Based Radar Data Assimilation Study over the Seoul Metropolitan Area

: This study investigates the ability of the high-resolution Weather Research and Forecasting (WRF) model to simulate summer precipitation with assimilation of X-band radar network data (X-Net) over the Seoul metropolitan area. Numerical data assimilation (DA) experiments with X-Net (S- and X-band Doppler radar) radial velocity and reﬂectivity data for three events of convective systems along the Changma front are conducted. In addition to the conventional assimilation of radar data, which focuses on assimilating the radial velocity and reﬂectivity of precipitation echoes, this study assimilates null-echoes and analyzes the e ﬀ ect of null-echo data assimilation on short-term quantitative precipitation forecasting (QPF). A null-echo is deﬁned as a region with non-precipitation echoes within the radar observation range. The model removes excessive humidity and four types of hydrometeors (wet and dry snow, graupel, and rain) based on the radar reﬂectivity by using a three-dimensional variational (3D-Var) data assimilation technique within the WRFDA system. Some procedures for preprocessing radar reﬂectivity data and using null-echoes in this assimilation are discussed. Numerical experiments with conventional radar DA over-predicted the precipitation. However, experiments with additional null-echo information removed excessive water vapor and hydrometeors and suppressed erroneous model precipitation. The results of statistical model veriﬁcation showed improvements in the analysis and objective forecast scores, reducing the amount of over-predicted precipitation. An analysis of a contoured frequency by altitude diagram (CFAD) and time–height cross-sections showed that increased hydrometeors throughout the data assimilation period enhanced precipitation formation, and reﬂectivity under the melting layer was simulated similarly to the observations during the peak precipitation times. In addition, overestimated hydrometeors were reduced through null-echo data assimilation. This study was conducted by using S- and X-band radar with AWS data. Radar reﬂectivity was ﬁrstly partitioned to corresponding hydrometeor types based on the background model temperature. The classiﬁed reﬂectivity was assimilated into the model by converting this value into a mixing ratio of four hydrometeors with the reﬂectivity mixing ratio relationship. The numerical experiment for the initial and boundary conditions without data assimilation, the CTRL, was used as a measure of the degree of improvement of radar data assimilation. The traditional radar data assimilation method, referred to as the DA experiment, only applied precipitation radar echoes in the assimilation. The DA_NP experiment assimilated both weather radar echoes and null-echoes to suppress and correct erroneous model precipitation. data characteristics development of reﬂectivity the precipitation peak was not simulated in the CTRL, strong reﬂectivity below layer simulated in DA and DA_NP, to the observations at the peak predictability improved through null-echo data


Introduction
Heavy rainfall is frequent over the Korean Peninsula during the summer, causing property damage and many human casualties. A total of 740 million dollars in property damage, including loss of drainage and inundation, was caused by heavy rainfall from 2 to 11 July and on 14 July 2017 [1]. Weather forecasts are directly related to the lives of individuals and organizations; thus, the importance of predicting heavy rainfall is increasing. However, meso-β-scale heavy rainfall has a radius of 10 to 100 km and a duration of tens of minutes to hours and is, thus, very difficult to predict because of wide variations in precipitation intensity and precipitation area depending on the region [2,3]. Therefore, research is needed to improve the predictive accuracy by assimilating observations with high spatio-temporal resolution into a numerical weather prediction (NWP) model.    In all three events considered, every 30 min, cycling of data assimilation is performed for 3 h, totaling six assimilation cycles. In each assimilation cycle, observations from radar reflectivity, radial velocity, and AWS observations are used.

WRF 3D-Var Assimilation System
The 3D-Var data assimilation method in WRF system is capable of assimilating all types of conventional observations, as well as radar data. The 3D-Var system combines the observations with background information on the model state and uses a linearized forecast model to ensure that the observations are given a dynamically realistic and consistent analysis field [10]. In this study, we used an incremental WRF 3D-Var data assimilation system [24,25]. The 3D-Var incremental system minimizes a scalar term called the cost function (or objective function, J), defined as a function of the analysis increment relative to the first gauss (background), using a linearized observation operator. The minimizations are subjected to the constraint of observation uncertainty, the difference between the observations, and the analysis projected to the observation space using the observations operator (H).
The terms J o and J b are the cost functions obtained from observation and background term, whereas R is the observation error covariance matrix. The term d = y o − H(x b ) is the innovation vector, which measures the departure of the observation y o and the background x b , while H is the linearization of the nonlinear observation operator H. The term ν is the control variable, where ν = U (x − x b ), and U is the decomposition of the background error covariance B under the constrain B = UU T . The covariance matrix for the regional background error statistics was calculated using National Meteorological Center (NMC) method [26] by taking 24 h and 12 h forecasts during the summer period (June, July, and August) for the chosen domain with five control variables. These are the stream function, unbalanced temperature, unbalanced surface pressure, unbalanced velocity potential, and pseudo-relative humidity. The minimization of cost is performed using the conjugate gradient method.

Radar Reflectivity Observation Operator
The radar reflectivity was partitioned into the reflectivity of each hydrometeor type based on the model background temperature by using the hydrometeor classification method and then converted to the hydrometeor mixing ratio.
The observed reflectivity (Z o ) was converted from dBZ to mm 6 ·m −3 , which is the unit for input reflectance (Z e ) and is expressed as Z o = 10 log 10 Z e .
Remote Sens. 2020, 12, 893 5 of 24 Z e can be expressed as Z e = Z r + Z ds + Z ws + Z g , because it is a volume average that is observed by several hydrometeors, such as rain (r), dry snow (ds), wet snow (ws), and graupel (g) [27][28][29].
For the precipitation echo data assimilation (Z o > −15 dBZ), Wang et al. (2013) classified hydrometeors by using the model's temperature field (T (K)). Rain exists in a grid with a temperature of T ≥ 5 • C, and a grid temperature of −5 • C < T < 5 • C assumes that rain, wet snow, hail, and dry snow can coexist (Equations (4)- (7)). The α in Equations (5)-(6) represents a value of zero at −5 • C with α = 1 at 5 • C, and it varies linearly between zero and one with the model temperature (Equation (8)). (4) α= The reflectivity of the hydrometeors was converted into the mixing ratio (kg·kg −1 ) of each hydrometeor by using the equation of the reflectivity-mixing ratio relationship. The hydrometeor mixing ratio was then used as an indirect assimilation method to assimilate reflectivity into the model [30].
where ρ a is the density (kg·m −3 ) of air. To create an environment in which convective clouds are actively maintained, the water vapor mixing ratio was nudged as the saturated water vapor mixing ratio when the observed reflectivity was greater than 30 dBZ. The saturated water vapor mixing ratio was calculated by using the Clausius-Clapeyron equation (e (hPa)) for water, and the water vapor saturation mixing ratio (q s ) was calculated as follows: where L is 2.5 × 10 6 J·kg −1 by heat of evaporation, R v represents the gas constant of water vapor (461.51 J·kg −1 ·K −1 ), is the ratio of the gas constant of dry air to the gas constant of water vapor, and P (hPa) is the pressure of the model. For the null-echo data assimilation with the observed Z o < −15 dBZ, we substituted Z e into Equations (9)-(12) with −15 dBZ when the model background reflectivity (Z b ) exceeded −15 dBZ and the corresponding value for each hydrometeor mixing ratio (q x ), which is expressed as Remote Sens. 2020, 12, 893 6 of 24 We assumed that only liquid clouds existed when the temperature of the model background was larger than or equal to 5 • C. In this criterion, the relative humidity for water vapor in the model (RH b ) was reduced to 95% when close to saturation. Assuming mostly solid cloud particles when the temperature was less than or equal to −5 • C, only the super saturated water vapor mixing ratio was removed. If the temperature of the background was between −5 • C and 5 • C, the two preceding criteria were mixed and added linearly according to the temperature of the background. In other words, no precipitation except clouds remained in the atmosphere.
The saturation mixing ratio for ice was obtained by using Equation (13) for ice (L is 2.83 × 10 6 J·kg −1 , the latent heat of sublimation), and the water vapor pressure for ice was calculated with Equation (14). A schematic diagram for assimilating radar reflectivity with null-echo information is shown in Figure 2.

Doppler Radial Velocity Observation Operator
To assimilate the Doppler radar radial velocity (Vr), the wind variable of the model is converted to the Doppler radial velocity.
where u, v, and w (m•s −1 ) are the wind components of the model, x, y, and z are the locations of the radars, xi, yi, and zi are the locations of the radar observation data, ri (m) is the distance between the radar and the observation point, and VT (m•s −1 ) is the terminal velocity. Sun and Crook (1998) terminal velocity equation was used to relate rainwater mixing ratio to VT (Equation (22)) [31].
Where qr is the rainwater mixing ratio (kg•kg −1 ), and a is correctional term calculated using Equation

Doppler Radial Velocity Observation Operator
To assimilate the Doppler radar radial velocity (V r ), the wind variable of the model is converted to the Doppler radial velocity.
where u, v, and w (m·s −1 ) are the wind components of the model, x, y, and z are the locations of the radars, x i , y i , and z i are the locations of the radar observation data, r i (m) is the distance between the Remote Sens. 2020, 12, 893 7 of 24 radar and the observation point, and V T (m·s −1 ) is the terminal velocity. Sun and Crook (1998) terminal velocity equation was used to relate rainwater mixing ratio to V T (Equation (22)) [31].
where q r is the rainwater mixing ratio (kg·kg −1 ), and a is correctional term calculated using Equation (23).
where p 0 is ground pressure, and p is state pressure (Pa). More detailed information on how the operator of Equation (21) is connected to the control variables can be found in Xiao et al. (2005) [22].

Model Configuration and Experimental Design
This study was conducted by using the WRF version 3.9.1 model, which was developed by the National Center for Atmospheric Research (NCAR) [32]. As shown in Figure 3, this model consisted of three nesting domains with resolutions of 9 km, 3 km, and 1 km, and it was centered over the Seoul metropolitan area, where heavy rainfall occurred. Cumulus parameterization was applied to all areas because the Multiscale Kain-Fritsch Scheme suitable for high-resolution models was used [33]. The WRF Double Moment 6 class (WDM6) scheme was employed for cloud microphysics processes, and the planetary boundary layer (PBL) used the Yonsei university (YSU) scheme [34,35]. A rapid radiative transfer model (RRTM) was selected for long-wave processes, and the Dudhia short-wave scheme was selected for short-wave processes [36,37]. Surface layer processes used the revised MM5 Monin-Obukhov scheme, and land surface processes utilized the Unified Noah land surface model [38,39]. The initial and boundary conditions were based on the NCEP GFS 0.25 Degree Global Forecast Grids Historical Archive, which has a resolution of 0.25 • and was provided by the National Centers for Environmental Prediction/National Centers for Atmospheric Research (NCEP/NCAR). The model configurations are summarized in Table 2.
Remote Sens. 2020, 12, x FOR PEER REVIEW 8 of 26 scheme was selected for short-wave processes [36,37]. Surface layer processes used the revised MM5 Monin-Obukhov scheme, and land surface processes utilized the Unified Noah land surface model [38,39]. The initial and boundary conditions were based on the NCEP GFS 0.25 Degree Global Forecast Grids Historical Archive, which has a resolution of 0.25° and was provided by the National Centers for Environmental Prediction/National Centers for Atmospheric Research (NCEP/NCAR). The model configurations are summarized in Table 2.    Three experiments were performed to compare the effects of data assimilation and null-echo data assimilation ( Table 3). The first experiment, the CTRL experiment, generated a 12 h forecast by using only the initial and boundary conditions without performing any data assimilation. The CTRL was set as a measure of the degree of improvement for the data assimilation experiments. The second experiment, referred to as DA, involved data assimilation as performed in the majority of studies with experiments that assimilate AWS data with radar radial velocities and precipitation echoes. The third experiment, DA_NP, further assimilated null-echo information into the model that already included precipitation echoes from DA. These two experiments (DA and DA_NP) assimilated observations into domains 2 and 3 for 3 h with cycling intervals of 30 min totaling six times, and a 9 h forecast field was created. This avoids the need for thinning the input data and allows more radar and AWS data to be used in the model since the data are interpolated in the pre-processing stage. Furthermore, it can reduce the moisture inconsistency between domain 2 and domain 3 that can lead to forecast degradation when two-way nesting is applied. Table 3. Experimental design.

CTRL
No data assimilation DA AWS + radar radial velocity + radar reflectivity DA_NP AWS + radar radial velocity + radar reflectivity + null echoes

Description of the Cases
Three cases of heavy rainfall over the Seoul metropolitan area were selected to understand the predictability of NWP model forecasts with radar DA. These three cases were typical summertime convective systems along the Changma front in Korea. The forecast periods and characteristics of these cases are summarized in Table 4.

Accuracy Verification
For an objective comparison of the forecasts, we considered a couple of evaluation parameters as described below. Quantitative error verification was performed by using the cumulative precipitation from AWS data and the model simulations. The mean absolute error (MAE) and root-mean-square error (RMSE) can be defined as in Equations (24)-(26).
where P i is the predicted precipitation (mm), O i is the observed precipitation from AWS (mm), and N is the number of horizontal grids in the model.  Verification was performed to evaluate the occurrence of model precipitation compared to the AWS observations. A contingency matrix for skill score calculations is shown in Table 5. Accuracy determines whether the same forecast or non-forecast occurred, as expressed in Equation (27). Values closer to one indicate higher accuracy.
The critical success index (CSI) is an indicator that focuses on accurately predicting precipitation. The ratio of precipitation prediction (Equation (28)) is sensitive to A and becomes lower when B and C occur. A perfect prediction of precipitation has a value of one.
The equitable threat score (ETS) is similar to the CSI but removes the contribution from S f (Equation (29)), which is an accidental hit, with higher accuracy approaching one (Equation (30)).

Space-and Time-Averaged Comparison Fields (a) Contoured Frequency by Altitude Diagram (CFAD)
A CFAD is a contour plot that displays the frequency distribution of reflectivity in an area of detectable echo volume at each vertical height. Many studies showed that CFADs are a convenient tool to examine the characteristics of storms, especially in terms of their temporally averaged evolution [40][41][42][43]. Min et al. (2015) showed that CFADs can assess the characteristics of microphysics schemes in simulating summer monsoon and convective precipitation using radar observations [43].
A CFAD begins with a histogram calculation that uses a constant reflectivity bin width within a constant vertical volume. The CFAD ignores the horizontal echo structure and summarizes the frequency distribution information in a single two-dimensional plot. In this study, the radar and simulated reflectivity of the CFAD was −16 dBZ or greater, with a bin width of 0.5 dBZ and vertical intervals of 250 m. The histograms were normalized by the total number of points at each vertical level, and the results were converted into percentage values. The average reflectivity (dBZ1) and the average of the reflectivity factor in linear units (dBZ2, in mm −6 ·m −3 ) for each height were calculated as follows: where N is the total number of points at each height (level). The linear unit reflectivity factor was used to better identify the location of the bright band (BB). This phenomenon is caused by melting ice particles through the 0 • C isotherm layer. In addition, we calculated the cumulative reflectivity frequencies of the 25th, 50th, and 75th percentiles, which were included in each of the CFAD plots.
(b) Time-Height Cross-Sections (THCS) A time-height cross-section diagram is useful to identify the duration of precipitation and the various echo tops of the precipitating system [43]. This type of diagram is used to augment CFADs, which do not reveal the time evolution characteristics of reflectivity (i.e., precipitation) that are related to phase errors when simulating any storm. We calculated the THCS of the simulated reflectivity with a 1 h temporal resolution, whereas the radar observations had 20 min resolution (temporal). Additionally, all the observed datasets were vertically interpolated to 250 m from their native polar and model sigma vertical coordinates to the same vertical resolution. The range was limited to an approximately 100,000-m radius. The melting level in the time-height cross-section is represented by a dashed line that stretches horizontally with higher reflectivity below the neighboring area.

Increment of the Analysis Field
The predicted accuracy of a numerical forecast model depends on the fidelity of the analysis field. The horizontal increments of the DA and DA_NP analysis fields at the last cycle time were analyzed to understand changes in the analysis fields following the assimilation experiments. Figures 5-7 show the radar, CTRL, DA, and DA_NP reflectivity (dBZ) and increments of the water vapor mixing ratio (g·kg −1 ) for DA and DA_NP at 2.0 km above sea level. Figure 5 shows the results for 0900 UTC 2 July 2017. Reflectivity greater than 30 dBZ was observed in the West Sea and inland areas of the Korean Peninsula. The CTRL simulated reflectivity greater than 30 dBZ over the West Sea, but its area was smaller than that of the observations. DA and DA_NP, which included data assimilation, simulated greater reflectivity over the West Sea than the CTRL and were closer to the observations. The water vapor mixing ratio increased where the reflectivity was located in DA and DA_NP because of the effect of radar data assimilation. No water vapor reduction occurred in DA_NP because no over-simulated reflectivity was observed in the region where the reflectivity was below −15 dBZ. On 2100 UTC 22 July 2017, reflectivity greater than 30 dBZ was observed along the demilitarized zone (DMZ), and DA and DA_NP simulated reflectivity above 30 dBZ in the same area but the CTRL never simulated reflectivity ( Figure 6). The water vapor mixing ratio in DA and DA_NP increased in the area where the precipitation echo was located. Reflectivity below −15 dBZ was dominant in central Gyeonggi Province and Seoul. However, DA overestimated echoes above 30 dBZ. Because of the null-echo data assimilation, the water vapor mixing ratio of DA_NP was reduced to eliminate over-simulated reflectivity. At 0900 UTC 28 August 2018, over 30 dBZ of a well-defined precipitation echo along the frontal boundary (the Changma front) was observed in the West Sea (Figure 7). The CTRL did not simulate the arc-shaped echoes that developed and migrated over the West Sea. In contrast, DA and DA_NP simulated these arc-shaped echoes. However, echoes that did not exist in the radar observations were simulated in southern Gyeonggi because echoes were over-predicted in DA and DA_NP. In DA_NP, the water vapor mixing ratio was reduced after assimilating non-precipitation echo data.
along the frontal boundary (the Changma front) was observed in the West Sea (Figure 7). The CTRL did not simulate the arc-shaped echoes that developed and migrated over the West Sea. In contrast, DA and DA_NP simulated these arc-shaped echoes. However, echoes that did not exist in the radar observations were simulated in southern Gyeonggi because echoes were over-predicted in DA and DA_NP. In DA_NP, the water vapor mixing ratio was reduced after assimilating non-precipitation echo data.

Distribution of the Cumulative Precipitation
The distribution of the cumulative precipitation for each case is shown and analyzed for Cases 1, 2, and 3. Figure 8 shows the 9 h cumulative precipitation distribution of the AWSs, CTRL, DA, and DA_NP for Case 1. During this 9 h period, a narrow and elongated precipitation band of more than 100 mm was observed from the western to eastern metropolitan area. The CTRL simulated some

Distribution of the Cumulative Precipitation
The distribution of the cumulative precipitation for each case is shown and analyzed for Cases 1, 2, and 3. Figure 8 shows the 9 h cumulative precipitation distribution of the AWSs, CTRL, DA, and DA_NP for Case 1. During this 9 h period, a narrow and elongated precipitation band of more than 100 mm was observed from the western to eastern metropolitan area. The CTRL simulated some precipitation but underestimated the amount and could not accurately simulate heavy rainfall. DA simulated heavy rain over 110 mm from the west to east and simulated similar precipitation patterns to the AWSs. However, the precipitation areas were further north than those of the AWSs and did not simulate heavy rain over 110 mm to the east of Gyeonggi. DA_NP reduced over-precipitation over 110 mm in the northeastern portion of Gyeonggi, but the precipitation area was still farther north compared to the observations, and this model did not simulate heavy precipitation to the east of Gyeonggi. Figure 9 shows the 9 h cumulative precipitation distribution for Case 2. The convective system moved from the northwest to the southeast, dropping more than 100 mm of heavy rain. The CTRL experiment failed to simulate precipitation in the metropolitan area. DA simulated precipitation in the Seoul metropolitan area, and its observed precipitation distribution was more similar to that of the AWSs compared to the CTRL. However, the precipitation amount increased, overestimating the intensity by more than 160 mm to the south of Gyeonggi. In contrast, DA_NP simulated over 100 mm of precipitation in the Seoul area, similarly to the observations, and removed the over-predicted precipitation from DA. Figure 10 shows the cumulative precipitation distribution for Case 3. As the convective system moved eastward, more than 100 mm of precipitation occurred in the metropolitan area. The CTRL did not correctly simulate the convective system and simulated less than 10 mm of precipitation in the Seoul metropolitan area. DA simulated more than 180 mm for the major precipitation structure but overestimated the amount in the Seoul area by 50 mm compared to the observations. With DA_NP, the null-echo data assimilation simulated approximately 90 mm of precipitation, similarly to the observations in the Seoul metropolitan area. In the DA and DA_NP experiments, there are areas of 200 mm or greater precipitation to the West Coast that were overestimated in the model. The over-prediction occurred immediately after the cycling, which was caused by excessive amount of water vapor and hydrometeor contents due to radar data assimilation.
for Case 3. As the convective system moved eastward, more than 100 mm of precipitation occurred in the metropolitan area. The CTRL did not correctly simulate the convective system and simulated less than 10 mm of precipitation in the Seoul metropolitan area. DA simulated more than 180 mm for the major precipitation structure but overestimated the amount in the Seoul area by 50 mm compared to the observations. With DA_NP, the null-echo data assimilation simulated approximately 90 mm of precipitation, similarly to the observations in the Seoul metropolitan area. In the DA and DA_NP experiments, there are areas of 200 mm or greater precipitation to the West Coast that were overestimated in the model. The over-prediction occurred immediately after the cycling, which was caused by excessive amount of water vapor and hydrometeor contents due to radar data assimilation.

Model Verification
Accuracy verification of the rainfall prediction was performed by using quantitative error verification and categorical precipitation classification ( Table 5). The AWSs in South Korea have an effective radius of approximately 10 km 2 . Thus, the average value of nine model precipitation grid points was compared to the closest AWS location in domain 3.

Model Verification
Accuracy verification of the rainfall prediction was performed by using quantitative error verification and categorical precipitation classification ( Table 5). The AWSs in South Korea have an effective radius of approximately 10 km 2 . Thus, the average value of nine model precipitation grid points was compared to the closest AWS location in domain 3. Table 6 shows the model bias, MAE, and RMSE that were calculated by using the accumulated precipitation of the AWSs and forecast fields for each case. In all three cases (Cases 1-3), the CTRL showed negative bias after underestimating the precipitation location and intensity. However, the precipitation in the DA and DA_NP experiments was simulated through radar data assimilation, and the model bias, MAE, and RMSE were reduced compared to those of the CTRL. DA_NP showed the lowest error in all three cases (Cases 1-3), with excessively simulated water vapor and hydrometeors removed through null-echo data assimilation, simulating less precipitation than DA. The results of averaging the errors in all cases showed that the bias error was the greatest in the CTRL (−31.7 mm), and the precipitation in DA formed through data assimilation, but the bias of 6.9 mm indicates overestimated precipitation. DA_NP showed the lowest error rates, i.e., −1.2 mm for the bias, 20.9 mm for the MAE, and 31.2 mm for the RMSE, because of a reduction in the overestimated rainfall amount compared to DA.
The categories were also verified according to the occurrence of precipitation. The categorical validation calculated the accuracy, CSI, and ETS as verification statistics depending on whether the observational data from the AWSs and model experiments detected precipitation. The accuracy, CSI, and ETS were obtained for the three cases and averaged over time ( Figure 11). The effect of data assimilation continued for up to 6 h, and all three verification scores showed higher accuracy for DA and DA_NP compared to the CTRL. In DA_NP, the overestimated precipitation was reduced through null-echo data assimilation, showing the highest predicted accuracy. After 6 h, the verification scores in the experiment in which data assimilation was performed decreased because the overestimated rainfall amount increased the false alarm ratio (not shown).  Vertical profile verification of 3 h prediction was performed by using Osan radiosondes. Figure  12 shows the biases of the water vapor mixing ratio, temperature, and wind components of the models compared to the Osan radiosonde data. No radiosonde data were available for Case 2, so only Case 1 and Case 3 were validated. In both cases, an underestimated water vapor mixing ratio was Vertical profile verification of 3 h prediction was performed by using Osan radiosondes. Figure 12 shows the biases of the water vapor mixing ratio, temperature, and wind components of the models compared to the Osan radiosonde data. No radiosonde data were available for Case 2, so only Case 1 and Case 3 were validated. In both cases, an underestimated water vapor mixing ratio was obvious in the CTRL from 600 hPa to 800 hPa. DA overestimated the water vapor mixing ratio, while DA_NP reduced the over-predicted mixing ratio by assimilating the null-echo data assimilation. Radar data assimilation improved the temperature profiles and showed lower error in DA and DA_NP compared to the CTRL. The effect of radial velocity data assimilation also improved the wind forecasts compared to the radiosonde data.

CFADs
The CFADs of Cases 2 and 3 were analyzed only because Case 1 did not show significant effects from data assimilation. CFAD analysis was performed to analyze the microphysical development of precipitation. Both the radars and the models used reflectivity information within 100 km of the KWK location. Figure 13 shows a cumulative 9 h CFAD for the KWK radar site and the CTRL, DA, and DA_NP experiments. The CFAD of KWK shows that ice particles fell and grew from the upper level as the reflectivity increased with lower altitude. A strong reflectivity peak appeared at an altitude of 4 km, which is the location of the BB, and the growth of particles in the troposphere was adequately demonstrated. On the other hand, the CTRL hardly simulated precipitation; thus, no apparent microphysical structure was visible. In DA and DA_NP, ice particles fell and grew from the upper level, showing the highest frequency around 4.5 km below the melting layer at 5.5 km (dashed line). In the experiment in which data assimilation was performed, the reflectivity value was higher than the observed value for KWK, and the 50% line was shifted to the right. These results were caused by the double moment characteristics of WDM6, which employs a double moment only for liquid water droplets and, thus, tends to create a sharp increase in radar reflectivity below the melting level of 0

CFADs
The CFADs of Cases 2 and 3 were analyzed only because Case 1 did not show significant effects from data assimilation. CFAD analysis was performed to analyze the microphysical development of precipitation. Both the radars and the models used reflectivity information within 100 km of the KWK location. Figure 13 shows a cumulative 9 h CFAD for the KWK radar site and the CTRL, DA, and DA_NP experiments. The CFAD of KWK shows that ice particles fell and grew from the upper level as the reflectivity increased with lower altitude. A strong reflectivity peak appeared at an altitude of 4 km, which is the location of the BB, and the growth of particles in the troposphere was adequately demonstrated. On the other hand, the CTRL hardly simulated precipitation; thus, no apparent microphysical structure was visible. In DA and DA_NP, ice particles fell and grew from the upper level, showing the highest frequency around 4.5 km below the melting layer at 5.5 km (dashed line). In the experiment in which data assimilation was performed, the reflectivity value was higher than the observed value for KWK, and the 50% line was shifted to the right. These results were caused by the double moment characteristics of WDM6, which employs a double moment only for liquid water droplets and, thus, tends to create a sharp increase in radar reflectivity below the melting level of 0 • C [42]. Thus, this model tends to overestimate the reflectivity compared to radar observations. The average reflectivity below the melting layer in DA_NP decreased by 5 dBZ compared to DA, showing a similar value to that of the observations. This phenomenon resulted from a reduction in water vapor and hydrometeors from the null-echo data assimilation. Figure 14 is the accumulated 9 h CFAD for Case 3. The CTRL simulation again failed to simulate the convective system; thus, the structure of hydrometeors growing from the top level did not appear. The data assimilation in DA and DA_NP increased the reflectivity from the upper level, showing the growth of hydrometeors, with the greatest reflectivity below the melting layer, a similar pattern to the observations. However, Case 3 also overestimated the reflectivity because of the inaccuracy of the microphysics scheme in simulating convective storms and the assimilation method of radar data, as noted in Case 2. DA_NP reduced the reflectivity of the 50% cumulative line below the melting layer by eliminating the mixing ratio of excessively simulated water vapor and hydrometeors in the model.

Time-Height Cross-Sections (THCS)
The characteristics and development of reflectivity (precipitation) over time were analyzed through THCSs. Figure 15 shows the THCSs of the KWK radar site and the CTRL, DA, and DA_NP experiments. Precipitation was observed in the KWK radar throughout the forecast period, with a reflectivity peak of about 30 dBZ near 4 km from 2230 UTC 22 July to 0030 UTC 23 July. In the CTRL, reflectivity was not found until 0300 UTC 23 July because precipitation was not simulated by the model. In the DA and DA_NP experiments, the THCSs were similar to that of the KWK reflectivity, but strong reflectivity occurred early in the prediction during the cycling period. The hydrometeors and humidity were added to existing moisture field of the model, and a larger amount of precipitation was simulated at the beginning of the prediction, simulating a precipitation peak from 2200 UTC to 2300 UTC 22 July. DA_NP showed the most similar reflectivity compared to KWK by reducing the over-simulated reflectivity through null-echo data assimilation. Null-echo data assimilation removed the excessive precipitation during the forecast period by removing the water vapor in the atmosphere, which reduced the mean reflectivity during the 2-6 h forecast period. However, it was not effective in controlling the strong reflectivity that occurred immediately after the cycling period from 2200 UTC to 2300 UTC 22 July. Figure 16 is the THCS of Case 3. Reflectivity was observed throughout the simulated period in KWK, but the CTRL only simulated the reflectivity in the beginning and middle of the period, with the amount less than that of KWK. DA and DA_NP overestimated the reflectivity but simulated the precipitation peak, similarly to the observations. DA_NP in Case 3 also showed lower reflectivity than DA, reducing the over-simulated precipitation from 1400 UTC to 1800 UTC 23 July.

Time-Height Cross-Sections (THCS)
The characteristics and development of reflectivity (precipitation) over time were analyzed through THCSs. Figure 15 shows the THCSs of the KWK radar site and the CTRL, DA, and DA_NP experiments. Precipitation was observed in the KWK radar throughout the forecast period, with a reflectivity peak of about 30 dBZ near 4 km from 2230 UTC 22 July to 0030 UTC 23 July. In the CTRL, reflectivity was not found until 0300 UTC 23 July because precipitation was not simulated by the model. In the DA and DA_NP experiments, the THCSs were similar to that of the KWK reflectivity, but strong reflectivity occurred early in the prediction during the cycling period. The hydrometeors and humidity were added to existing moisture field of the model, and a larger amount of precipitation was simulated at the beginning of the prediction, simulating a precipitation peak from 2200 UTC to 2300 UTC 22 July. DA_NP showed the most similar reflectivity compared to KWK by reducing the over-simulated reflectivity through null-echo data assimilation. Null-echo data assimilation removed the excessive precipitation during the forecast period by removing the water vapor in the atmosphere, which reduced the mean reflectivity during the 2-6 h forecast period. However, it was not effective in controlling the strong reflectivity that occurred immediately after the cycling period from 2200 UTC to 2300 UTC 22 July. Figure 16 is the THCS of Case 3. Reflectivity was observed throughout the simulated period in KWK, but the CTRL only simulated the reflectivity in the beginning and middle of the period, with the amount less than that of KWK. DA and DA_NP overestimated the reflectivity but simulated the precipitation peak, similarly to the observations. DA_NP in Case 3 also showed lower reflectivity than DA, reducing the over-simulated precipitation from 1400 UTC to 1800 UTC 23 July.

Summary and Conclusions
This study investigated the ability of the high-resolution WRF model to simulate summer precipitation with dual polarized X-Net data over the Seoul metropolitan area. DA experiments with X-Net (S-and X-band Doppler radar) radial velocity and reflectivity data for convective summer monsoon events were evaluated. The conventional assimilation of radar data, which involved the assimilation of reflectivity and radial velocities from precipitation echoes, and the assimilation of null-echoes and its effects on short term QPF were studied. Three heavy rain events were simulated and analyzed in detail by comparing the cumulative precipitation, CFADs, THCSs, and verification statistics.

Summary and Conclusions
This study investigated the ability of the high-resolution WRF model to simulate summer precipitation with dual polarized X-Net data over the Seoul metropolitan area. DA experiments with X-Net (S-and X-band Doppler radar) radial velocity and reflectivity data for convective summer monsoon events were evaluated. The conventional assimilation of radar data, which involved the assimilation of reflectivity and radial velocities from precipitation echoes, and the assimilation of null-echoes and its effects on short term QPF were studied. Three heavy rain events were simulated and analyzed in detail by comparing the cumulative precipitation, CFADs, THCSs, and verification statistics.

Summary and Conclusions
This study investigated the ability of the high-resolution WRF model to simulate summer precipitation with dual polarized X-Net data over the Seoul metropolitan area. DA experiments with X-Net (S-and X-band Doppler radar) radial velocity and reflectivity data for convective summer monsoon events were evaluated. The conventional assimilation of radar data, which involved the assimilation of reflectivity and radial velocities from precipitation echoes, and the assimilation of null-echoes and its effects on short term QPF were studied. Three heavy rain events were simulated and analyzed in detail by comparing the cumulative precipitation, CFADs, THCSs, and verification statistics. This study was conducted by using S-and X-band radar with AWS data. Radar reflectivity was firstly partitioned to corresponding hydrometeor types based on the background model temperature. The classified reflectivity was assimilated into the model by converting this value into a mixing ratio of four hydrometeors with the reflectivity mixing ratio relationship. The numerical experiment for the initial and boundary conditions without data assimilation, the CTRL, was used as a measure of the degree of improvement of radar data assimilation. The traditional radar data assimilation method, referred to as the DA experiment, only applied precipitation radar echoes in the assimilation. The DA_NP experiment assimilated both weather radar echoes and null-echoes to suppress and correct erroneous model precipitation.
The results of the horizontal incremental analysis showed that the water vapor mixing ratio and hydrometeors increased in areas where radar reflectivity was observed for both DA and DA_NP. In DA_NP, null-echo information was added, and the over-simulated water vapor and hydrometeors were removed. Accuracy evaluation with the total cumulative precipitation resulted in reduced errors for DA and DA_NP compared to the CTRL. However, over-predicted precipitation is a known problem when radar data are assimilated. DA_NP showed the highest accuracy by reducing over-simulated precipitation through null-echo data assimilation. This effect was also shown in the verification statistics, with the highest score from DA_NP up to 6 h when the effect of data assimilation was sustained. CFAD analysis was performed to analyze the microphysical structure of precipitation. The data assimilation simulated the process of hydrometeors falling and growing from the upper atmosphere. The CFADs also showed a decrease in reflectivity because of the reduction in over-predicted hydrometeors in DA_NP. The characteristics and development of reflectivity (precipitation) over time were also analyzed with THCSs. Although the precipitation peak was not simulated in the CTRL, strong reflectivity below the melting layer was simulated in DA and DA_NP, similarly to the observations at the peak precipitation times. Thus, precipitation predictability was improved through null-echo data assimilation. This study demonstrated the applicability of null-echo data assimilation to simulate convective storms when radar data assimilation is utilized. The application of this radar data assimilation method can increase the predictive accuracy.
When performing radar data assimilation, problems with over-simulating reflectivity in the model and causing heavy precipitation can be found early in the forecast period. Further study is needed to evaluate the detailed performance of how microphysics and convective schemes play a role in assimilation and reduce the dynamical imbalance from discontinuity that is generated by radar data assimilation.
Author Contributions: J.-W.L. wrote the first draft, assisted with data curation, and performed formal analysis, investigation, software coding, and visualization of the research article. K.-H.M. conceptualized the paper, supervised and administered the project, performed formal analysis, provided resources and funding, and reviewed and edited the paper. Y.-H.L. and G.L. supported with funding, provided methods, and reviewed and edited the paper. All authors have read and agreed to the published version of the manuscript.