Estimation of the Total Amount of Enhanced Rainfall for a Cloud Seeding Experiment: Case Studies of Preventing Forest Fire, Drought, and Dust

. In this study, a method for verifying the efect of cloud seeding in the case of a mixture of natural and artifcial rainfall bands was proposed, and its applicability to each experimental case was evaluated. Water resources that could be secured through cloud seeding were also quantifed for the experiments on forest fre prevention, drought mitigation, and dust reduction in 2020. Data on numerical simulations, radar-derived rainfall, rain gauge-derived rainfall, and weather conditions were applied. Areas with seeding and nonseeding efects were classifed according to the numerical simulation results and wind system, and enhanced rainfall was determined by comparing the changes in rainfall between the two areas. Te amount of water resources was determined by considering the area of the seeding efect and rainfall density. As a result, 1.74 mm (4.75 million tons) of rainfall increased from the experiment on forest fre prevention, 0.84mm (1.30million tons) on drought mitigation, and 2.78mm (24.44 million tons) on dust reduction. Tus, an average rainfall of 1.0 mm could be achieved through the experiment. Tese results helped verify the pure seeding efect and achieve the experimental purpose.


Introduction
Water-related disasters occur frequently worldwide because of climate change, and their damage is rapidly increasing. In some areas, fash foods occur because of heavy rainfall and typhoons, resulting in food damage, power outages, and landslides [1][2][3][4][5]. In other areas, rainfall does not occur for a long period, resulting in drought problems and large forest fres due to arid climates [6][7][8][9]. Terefore, rainfall has been comprehensively investigated to secure available water resources and improve the atmospheric environment [10,11].
Various methods for rainfall management worldwide have been proposed. For example, the amount of water resources can be increased by storing rainfall on the ground with dams, reservoirs, and rain gutter or by growing cloud droplets. Rainfall storage has been addressed in water resources engineering [12][13][14], whereas the latter has been conducted through cloud seeding experiments in meteorology [15][16][17]. Cloud seeding enhances rainfall by seeding hygroscopic materials in clouds to identify the processes of cloud physics.
In Korea, since the National Institute of Meteorological Sciences (NIMS) introduced the aircraft in 2017, studies on weather modifcation have been conducted to mitigate drought, prevent forest fre, dissipate fog, and reduce dust. It plans to perform a cloud seeding experiment based on numerical simulation (NS) and applies various types of equipment that can observe rainfall and cloud particles to verify the seeding efect on a target area. Te NIMS also aims to develop a seeding strategy based on cloud characteristics by securing data from various experiments every year. Tus, research by the NIMS has signifcantly contributed to the feld of cloud physics [30].
Despite these studies, accurately analyzing the changes in the amount of rainfall through cloud seeding experiments is difcult [15,[31][32][33][34]. In most studies, seeding efects are determined by statistically analyzing long-term experimental cases. In the United States, Climax I and II experiments conducted to verify the seeding efect on a mountainous terrain in winter from 1960 to 1970 revealed that the rainfall was increased by 18% [35][36][37]. We reanalyzed this result and made corrections for a 10% increase in rainfall. Recently, Rasmussen et al. [38] statistically analyzed an experiment for 6 years (2008-2013) and confrmed that 1.5% of the enhanced rainfall (ER) can occur in the analysis regions after seeding. In Israel, various cloud seeding experiments were performed from 1961 to 1975 to statistically analyze long-term data; they confrmed that the rainfall was increased by approximately 15% in the northern part of the experimental area [39,40]. In the United Arab Emirates, Hosari et al. [41] analyzed heavy rainfall events in the years when seeding was not performed  and in the years when seeding was performed (2003-2019); they analyzed and confrmed the annual average of 23% ER in the seeding area. Tese studies have difculties securing long-term data; therefore, for rapid technological improvement, an analysis technique for each event is needed despite uncertainties.
A method that can separate the pure rainfall amount from a natural cloud system should be developed to verify the seeding efect in each case. Te simplest method is determining the amount of rainfall by comparing the characteristics of the seeding and nonseeding efect areas. Trough ER calculation, the usefulness of cloud seeding in securing water resources can be verifed. Furthermore, it can be helpful in planning economic and efcient experiments for various purposes, such as forest fre prevention and ecological protection in the future.
In this study, a method for calculating the efect of cloud seeding was developed by classifying the seeding and nonseeding efect areas. Changes in the rainfall system according to seeding materials were confrmed, and the amount of water resources that could be secured through the experiments for each case was calculated.
In Section 2, the characteristics of the data used in this study are summarized. In Section 3, the algorithm proposed for calculating the amount of ER that could be secured through a weather modifcation experiment is introduced. In Section 4, after the algorithm was applied to the three events in 2020, the seeding efect for each case was quantitatively determined. In Section 5, results are summarized, expected efects are discussed, and directions for future studies are presented.

Rain Gauge-Derived Rainfall and Radar-Derived Rainfall.
Rainfall afected by the seeding efect can be observed using ground observation and remote sensing equipment. Rain gauges are used to easily observe rainfall from the ground and record the rainfall intensity at 1-minute intervals. Tis value represents the average intensity of the rainfall observed during this period. Rainfall data calculated by accumulating this value for a certain period were analyzed in this study. Te cumulative time was determined to be within 3-6 h. Rain gauge observation data can easily capture the rainfall characteristics of the target area.
Remote sensing observation data can be used to analyze the characteristics and track the rainfall band afecting the surrounding area. A radar is a device that indicates the intensity of rainfall by analyzing the electrical signal power refected from the raindrops. In this study, radar-derived rainfall was applied to efectively analyze the seeding efect, which can appear spatially wide. Rainfall was calculated using radar refectivity, the temporal resolution of the data was 5 min, and the spatial resolution was 1 km × 1 km. Hybrid surface rainfall (HSR) radar refectivity [42], which is observed at an altitude closest to the ground, was applied to minimize the deviation between the rain gauge-and radarderived rainfall. Radar-derived rainfall was calculated using the following Z-R relationship, which is used to estimate the spatial distribution of rainfall by the Korea Meteorological Administration (KMA): where Z is the radar refectivity (mm 6 /m 3 ) and R is the rainfall (mm/h). Refectivity is a value indicating the number of water droplets per unit volume. Since the range of observations varies widely, data were stored in units of decibel relative to Z such as 10log (Z). Te KMA was used to provide a synthetic rainfall image in real time via equation (1).

Numerical Simulation for Weather Modifcation
Experiment. Te NS results of the efect of the seeding material are required to determine the seeding efect because directly tracking the seeding material difused into the atmosphere is practically difcult. In addition, numerical models are helpful in determining the optimal experimental area where the seeding efect is expected because the diffusion rate and direction can be initially determined according to the fight area and weather conditions through a preliminary simulation. Considering this advantage, the NIMS determines the target area based on the NS before the experiment [30]. Table 1 summarizes the numerical models in NIMS [43]. Te numerical model simulates the cloud growth and extinction process using the Morrison microphysics scheme. Te model also utilizes the Yonsei University (YSU) planetary boundary layer scheme, the NOAH land surface scheme, and the MM5 similarity surface layer scheme. In this study, cloud seeding experiments were analyzed using the numerical model with settings shown in Table 1.
Numerical models show that the growth and extinction of rainfall as seeding materials are spatially difused. Considering this phenomenon, the time when the ER is maximized after seeding can be determined. Te efect of natural rainfall in an area with a simulated ER may be checked by using other observational data in the same period. In this study, data from the NS that showed the changes in cumulative rainfall (CR) were used to verify the seeding efect. Changes in the CR indicated an increase or decrease in rainfall in areas expected to be afected by the seeding material. Tus, regions where changes are expected in cloud seeding experiments can be compared. Terefore, the spatial scale of a seeding material in a rainfall system may be determined, consequently allowing the analysis of the seeding efect in an area where a change in rainfall is expected even if it is not a target area.

Weather and Topography Conditions.
When the seeding material is difused into the target area, it is signifcantly afected by meteorological factors. Te weather conditions in the experimental area should be considered to determine the seeding efect quantitatively. Because the seeding material is spread directly from the aircraft, it is signifcantly afected by the wind system in the experimental area. Te wind direction and velocity infuencing difusion difer depending on the seeding altitude. However, verifying the seeding efect in the target area is difcult if the experiment is not conducted in the planned area because of fight limitations. Even if the aircraft spreads the seeding material at the seeding altitude as planned, rainfall may not occur in the target area, depending on the wind characteristics of the infow rainfall system. Te meteorological data observed from the aircraft after the cloud seeding experiment were analyzed. Weather charts and satellite images were also analyzed to fgure out the cloud characteristics. Tese images show the efect of the low pressure or cold front on clouds with a mixed seeding efect and the cloud distribution near the experimental regions. For the analysis of satellite images, infrared images corresponding to the 8.7 μm wavelength of the GK2A satellite were applied.
One of the meteorological factors to consider along with the wind system at the seeding altitude is temperature. Temperature afects seeding material determination used in the cloud seeding experiment. When the temperature of the seeding altitude is above 0°, CaCl 2 is applied; when the temperature is below 0°, AgI spreads to induce the enhancement of ice formation. Te updraft rate at the seeding altitude and liquid water content (LWC) also afect experimental results. Te stronger the vertical velocity of the cloud in which the seeding material is difused or the more the LWC, the more favorable the rainfall enhancement. Te condition of LWC in the seeded cloud was observed by the aircraft. In addition, the image of P-velocity, which represents the index of vertical velocity, was generated by a very short-range data assimilation and prediction system (VDAPS) for the analysis of temperature and updraft in this study. Te intensity of updraft according to the region, the contour line indicating the region with a similar temperature, the relatively cold regions (C) and the warm regions (W), and the wind speed are included in this image.
Te impact of topographical factors on the rapid growth of clouds and rainfall bands should be considered. However, interpreting orographic efects on cloud growth is not straightforward [44]. As clouds are more likely grow on a mountainous terrain than in urban areas, orographic effects should be excluded to verify the pure ER by the seeding material. In the case of the mountainous terrain around the target area, rain gauges with a diference of ≥300 m than the average altitude of the observation stations were excluded from the seeding efect analysis in this study.

Calculation of the Total Amount of the Enhanced Rainfall
Changes in clouds before and after seeding can be analyzed to verify the seeding efect, and regions afected by the seeding material can be compared with regions not assigned for the same cloud system. In the former case, a cloud can be stagnant or the characteristics do not change signifcantly over time; therefore, meeting the conditions when conducting the experiment on the Korean Peninsula is difcult. In this study, a method for estimating the amount of rainfall was proposed by comparing the area afected by the seeding material with the neighboring area where clouds with similar characteristics are distributed. Tis system can classify the seeding and nonseeding efect areas based on the wind system and calculate the total amount of water resources that can be secured through the cloud seeding experiment ( Figure 1). Te algorithm consists of four steps ( Figure 1): determination of analysis time (AT), determination of seeding efect box (SB) and nonseeding efect box (NB), calculation Advances in Meteorology of ER, and determination of the total amount of enhanced rainfall (TAER). In AT determination, the major period when the rainfall changed in the target area according to the difusion of the seeding material is determined. In this study, AT was identifed on the basis of the period of rainfall change near the target area in the NS. In addition, considering that the CR may be excessive depending on the rainfall characteristics fowing into the target area, AT was set to prevent exceeding a maximum of 6 h. Once AT has been determined, the areas afected and not afected by the seeding material should be classifed. First, the area of the major rainfall change during the AT is marked with a red-dotted-line box referred to as a model-simulated box (MB). Te seeding-afected area is determined on the basis of the area with the largest radar CR change in the MB. Considering the range in which the efective seeding efect can be verifed in the rainfall band, the seeding efect area is determined as a red box (50 km × 50 km) and referred to as the seeding efect box (SB). Te nonseeding-afected area is determined as a blue box (50 km × 50 km, such as the SB) according to the wind direction of the ER and natural rainfall band shown in the NS and referred to as the nonseeding efect box (NB). If the two directions are similar, the NB is determined on the basis of the center of the seeding line of the aircraft located on the windward side of the SB. If the rain gauge-derived analysis of the NB is insufcient, the location of the NB can be slightly corrected to include as much of the seeding line as possible. If the direction of the ER in the NS is diferent from that of the observed rainfall band, the NB is determined on the basis of the area with the greatest rainfall among the same natural rainfall systems without overlapping the MB.
Te ER can be calculated by comparing SB and NB. If two or more rain gauges are present where the CR exceeds

Decision of AT
Determine the AT based on significant enhanced rainfall in NS.

Determine SB and NB
Determine MB based on area of rainfall change in NS.
Determine SB centered on the point where the accumulated radar rainfall is the largest in MB during AT.
Simulated ER and the observed rainfall bands are same direction in AT.
Determine NB as not overlaping the MB, and the largest rainfall area in similar rainfall system.  0.1 mm during the AT in both areas, the cumulative average rainfall (CAR) of each area is calculated using the observation data. At this time, the altitude of the rain gauges in the NB does not difer by more than 300 m from the average altitude of the rain gauge in the SB. If the CAR is small or missing and if less than two rain gauges are available under the corresponding conditions, the radar-derived CAR estimated at each location of the rain gauge is analyzed. However, considering the estimated rainfall accuracy, the amount of the radar rainfall should be higher than 60% compared to that of the rain gauge rainfall. ER can be calculated as the diference between the CAR of the two areas. Tis technique is based on assumption that the seeding efect is mixed with the rainfall system fowing into the SB. Te calculated value is the average amount of rainfall that can be generated in the area afected by the cloud seeding experiment. If the estimated area where rainfall may be generated is applied to the calculated ER, the amount of water resources that can be secured through the cloud seeding experiment can be determined. In this study, this parameter was defned as the TAER. Te estimated rainfall area was calculated using NS, which is referred to as the efective rainfall area (ERA). Te ERA was determined by adding the areas where CR was expected to be 0.1 mm or more during the AT, and only the rainfall area of the land part, excluding the sea, was considered. Te TAER can be calculated by multiplying the amount of artifcial rainfall, efective rainfall area, and rainfall density for unit conversion, as shown in following equation:

Calculation of ER
where the unit of TAER is tons and R d is the rainfall density assumed to be 1 g/cm 3 . Te calculated TAER represents the amount of water resources that can be secured through a cloud seeding experiment.

Events and Weather Conditions.
In this study, three cases were selected from the cloud seeding experiments performed at NIMS in 2020, and the seeding efect was analyzed. Experimental cases on forest fre prevention, drought mitigation, and dust reduction were selected. Cloud seeding experiments were then performed in diferent areas according to each purpose. A research aircraft (Beachcraft King Air 350 HW) was used to conduct the experiment. For forest fre prevention, the experiment was conducted in Gangwon-do, and the Cloud Physics Observation Site (CPOS), a verifcation site for cloud seeding in NIMS, was set as the target area (Event 1). For drought mitigation, the experiment was conducted to increase the storage of Boryeong Dam in Chungnam Province (Event 2). For dust reduction, the experiment was carried out in the capital area, and the KMA observation site located in Seoul (capital) was set as the target area (Event 3). Te experimental results of the events and weather conditions observed by the aircraft during the experiment are summarized in Table 2. Te analyzed weather map, satellite image, and cloud image of the experimental area are shown in Figure 2.
In Table 2, stratiform clouds were observed in all three experiments conducted at an altitude of approximately 2 km for Events 1 and 2 and less than 1 km for Event 3. CaCl 2 was spread as the seeding material because of the temperature higher than 0°C at the seeding altitude. Te wind speed in Event 2 was the highest, and the wind direction was afected by the southwest wind in all three cases. Te vertical velocity was the same in Events 1 and 2, and Event 3 was the largest among the three cases. Te LWC at the seeding altitude was 0.2-0.4 g/m 3 in all three cases, confrming that the experiment was conducted on clouds containing sufcient water droplets. Te weather chart in Figure 2 shows that the cloud seeding experiment was performed on clouds afected by a low-pressure front in all the three experiments. Terefore, the cold front passing across the Korean Peninsula after seeding afects the rainfall development. Te satellite image shows that Event 1 had thicker clouds than the two other cases. Te image taken while seeding from the aircraft confrmed that the experiment was performed in the clouds. Figure 3 shows the changes in the updraft around the experimental area after seeding. In Table 2, Event 3, which shows a large updraft, can be afected by seeding material difusion and cloud growth. Figure 3 shows the P-velocity, which represents the index of vertical velocity at the 700 hPa altitude; it was generated by the VDAPS. Considering that the reaction time of the seeding material was diferent for each event, the fgures were shown at intervals of 1 h or 2 h. In Event 1, no signifcant change was observed in the updraft with time compared with that in the other events, but a strong temperature valley occurred along the east coast. In Event 2, the updraft was weak in the West Sea, but a strong updraft occurred between 14:00 and 15:00 in the center of the inland area. In Event 3, a strong updraft occurred in the northern area at 14:00 as the cold front passed. Tis updraft infuenced the deeper vertical transport of seeding material which afects the enhancement of precipitation.

Determination of Areas with Seeding
Efect and Nonseeding Efect. Te NS, rainfall data, and weather conditions were analyzed according to the algorithm shown in Figure 1 to verify the seeding efect in each experimental case. First, AT was determined on the basis of the change in the time of rainfall in the target area. In Event 1, the amount of rainfall after the experiment was not considerably higher than that of the two other cases; therefore, the analysis time was determined to be 6 h from 12:00 to 18:00. In Event 2, the movement speed of the rainfall fowing into the experimental area was relatively fast; thus, the analysis time was 3 h from 13:00 to 16:00 after the experiment. Event 3 was also decided to be 5 h from 13:00 to 18:00 in consideration of the movement velocity of the cloud. Figure 4 shows the time when SB and NB are determined using the results of the ER in the NS and rainfall data from the rain gauge and radar during the AT. Radar-derived rainfall was calculated by applying equation (1) to the HSR refectivity observed at the rain gauge location.

Advances in Meteorology
In the analysis result of Event 1 (Figure 4(a)), the cloud seeding experiment was performed in the western part of the CPOS, and the ER appeared widely over the East Sea because of the infuence of the southwest wind. Te SB was determined to include the eastern part of Gangwon-do and a part of the East Sea, where the major rainfall changes occurred. Te SB was determined to include the target area CPOS and the vicinity of Gangneung (GNG), where the radar-derived rainfall was large. Gangneung is a representative coastal city adjacent to the East Sea in central Gangwon-do. Te spatial distribution of the radar CR showed that the natural rainfall system moved southeast. Tis indicates that the wind direction in the coastline was northwest diferent from that in the seeding efect area due to the mountain ranges. Terefore, because the simulated ER and natural rainfall bands had diferent directions, the NB was determined in the northern region of the SB without the seeding efect among the regions with similar rainfall clouds to the SB. NB included Sokcho (SOC), a city in the northern part of the East Sea, and the northernmost region in South Korea.
In Event 2, the experiment was conducted over the west sea, southwest of Boryeong Dam (BRD), which is the target area, and the seeding efect was expected diagonally because of the infuence of the cold front accompanied by the southeast wind (Figure 4(b)). Although the area of the ER was smaller than that of Event 1, the SB was wider. SB included the Boryeong Dam (BRD) and Gunsan (GSN) areas. Gunsan is an urban area adjacent to the West Sea, near the border between the central and southern regions of South Korea. Because the movement direction of the natural rainfall band appeared to the northeast in the same way as in the direction of the simulated ER, unlike Event 1, NB was determined to include the windward seeding line. Considering that the cloud seeding experiment was performed in the West Sea, the NB was determined to include the analyzed rain gauge and seeding line as much as possible.
Te wind system of Event 3 appeared to be in the southwest similar to Event 2; because it was aimed at the KMA station in Seoul, the experiment was conducted in the West Sea (Figure 4(c)). Compared with the two other events, the widest regions were the ones where rainfall changes would occur because of the rapid movement of the cold front. Te SB, including the West Sea, the capital area, and Gangwon-do, was determined to be the largest. Te SB also included Ansan (ASN), a satellite city located in the northern part of Gyeonggi-do and Seoul, where the radar-derived rainfall was the strongest. NB was determined in the windward seeding line area as in Event 2, considering that the wind direction of the rainfall system and the simulated ER are the same. Table 3 summarizes the basic information of the analyzed rain gauge included in each box determined by the event. In addition to the location of the rain gauge, the altitude from the sea level and the CR of the rain gauge and radar during the AT are shown in Table 3. In Event 1, data from rain gauges distributed near the shoreline were analyzed in NB and SB to exclude the orographic efect on the rainfall cloud. Terefore, seven stations, including Gangneung (GNG), were selected for SB and four stations, including Sokcho (SOC), were chosen for NB. In Event 2, SB had 15 stations, including Boryeong Dam (BRD) and Gunsan (GSN), and NB had 6 stations. At Daecheonhang (DCH) in the SB and Galmaeyeo (GMY) in the NB, no rainfall was observed during the AT. In Event 3, in Figure 4(c), only stations with an elevation of less than 100 m above sea level were analyzed considering the excessive concentration of rain gauges around the capital area. Accordingly, 25 stations, including the KMA, were selected in the SB, and 5 stations were chosen for the NB. At the GDA station in the NB, no rainfall was observed in the AT.
In Table 3, the CR during AT was greater on the rain gauge than on the radar. At the SHA station in Event 3, the rain gauge-observed rainfall was at least three times greater than the radar-derived rainfall because the rainfall observation methods of rain gauge and radar are diferent; furthermore, rain gauge rainfall observations can be signifcantly infuenced by the surrounding weather condition. Although radar-derived rainfall may have a lower measurement accuracy than rain gauge-derived rainfall, radar rainfall is advantageous for analyzing the average characteristics of rainfall systems distributed in the atmosphere. Terefore, in this study, the changes in the radar average refectivity before and after seeding in the two areas were compared ( Figure 5).
In all three events, the increasing and decreasing trends of the average refectivity of SB were similar to those of NB ( Figure 5). Te locations of the two regions were diferent, and the time of increase in refectivity did not coincide because clouds afected the SB after they passed through the  In Event 2, the radar refectivity in both areas decreased. In the SB, refectivity was continuously observed; in the NB, it was observed until 15:30 only. At the time of analysis, refectivity decreased in the SB, but it temporarily increased by 5.02 dBZ from 8.81 dBZ to 13.83 dBZ at 14: 25-14:55. Tis relative increase is not well apparent from Figure 5(b) in the written time range. Presently, it is not scientifcally clear how to distinguish between natural rainfall and cloud 332 seeding-induced rainfall [21]. However, in Figure 5(b), the average refectivity of the efective time (green-shaded box) is 13.8 dBZ in the SB and 13.1 dBZ in the NB though it is an indirect verifcation.
In Event 3, refectivity was observed only before 17:00 in both areas. Unlike the two other events, the refectivity between the two areas before the AT largely difered but showed similar increasing and decreasing trends during the AT. Tis fnding showed that the characteristics of rainfall fowing into the two regions changed consistently (homogeneously) over time. In the AT, as in Event 1, a two-fold increase in radar refectivity was observed because of the efect of the seeding material. In the NB, refectivity changed signifcantly over time. Terefore, rainfall in both the areas should be analyzed to verify the seeding efect. Figure 5 illustrate that the seeding efect can be verifed using rain gauge-measured data. 14

Determination of the TAER.
Once the SB and NB are determined, the amount of ER can be calculated according to Figure 1. Te CAR of the two areas can be determined and the ER can be calculated using the rain gauge-and radar-derived rainfall data shown in Table 3. Table 4 shows the statistical data of the CR of the rain gauge-and radar-derived data in each box by event and the amount of ER calculated by the diference in average rainfall between the two areas.
In Table 4, rain gauge-and radar-derived rainfall is statistically analyzed using only the data observed for the CR of ≥0.1 mm (Figure 1). As shown in Table 4, the CAR of the radar was smaller than that of the rain gauge. Te maximum value and variance of radar rainfall were also smaller than those of rain gauge rainfall except in Event 1. In the NB of Event 1, the variance of the radar-derived rainfall was the smallest among all cases. In rain gauge-derived rainfall, the variance was the largest in the SB of Event 3 among all cases. Tese results showed that radar-derived rainfall has a smaller spatial deviation than rain gauge-derived rainfall [45]; however, the amount of rainfall is small [46,47] because the radar observes raindrops in the atmosphere, but as droplets fall, observation loss may occur in the rain gauge depending on weather conditions. Te amount of ER calculated using the diference in CAR between the two boxes was similar to that of the radar and rain gauge. Tis fnding indicated that the observed values might difer at individual points; however, the cloud characteristics distributed in the two areas were similar ( Figure 5). Of course, the ER calculated by this analysis method may vary depending on the location of the area. Te SB is determined based on the numerical simulation, but the NB is determined among regions where similar systems are distributed, so it cannot be an objective. Terefore, it is necessary to analyze the sensitivity of the CAR determined according to the location of the NB. In this study, as shown 10 Advances in Meteorology in Figure 6, the results of the sensitivity analysis are shown for each event. Here, the CAR that changes when the NB in Figure 4 is moved in all directions is compared. Te movement directions were considered as southwest, south, southeast, west, east, northwest, north, and northeast based on the NB location in Figure 4 (Basic). In addition, the movement range was set to 25 km, which corresponds to half the size of NB, so that the analysis rain gauge could be changed without overlapping with the SB as much as possible.
In Figure 6(a), showing the result of Event 1, the CAR of the rain gauge in the NB determined by the algorithm was medium compared to that calculated in other locations. Rainfall increased as the unafected area moved      12 Advances in Meteorology south and decreased as it moved north. As shown in Figure 4, this is because precipitation was stronger in the southern area of the NB. In addition, since rain gauges are rarely distributed in the northern area, the CR was afected by the small observed value, and when the NB moved northeast, the rain gauges were not included. In the case of the radar rainfall, it was found that the CAR did not change signifcantly because the precipitation system variation was not large around the NB.
In Events 2 and 3, the CAR of the rain gauge and the radar before the location of the NB was moved which was almost larger than when it was moved to another location (Figures 6(b)-6(c)). In the case of Event 2, when the location of the NB was moved to the west (southwest, west, and northwest), the CAR was not observed due to rarely distributed rain gauge on the sea and no precipitation detection. On the other hand, when the NB was moved inland where the rain gauges are densely distributed, the CAR was   Advances in Meteorology calculated to be larger. In Event 3, the rain rate was strong both in the capital area and in the sea, but the CAR was calculated to be large when moving the NB to the inland, where there are many rain gauges. Te radar rainfall also showed similar results to the rain gauge rainfall in both cases, but the change in the CAR according to the location of the NB was larger in Event 3 than in Event 2.
Comparing the analysis results, Event 1 showed an ideal result that was close to the median value among the amount for which the CAR could be calculated, and the other two cases had a relatively large CR. However, the amount of the ER calculated in this study is a value calculated by comparing the area afected by seeding with the area on the windward side. If another area is determined for the NB, an ideal result can be obtained, but comparing the seeding-afected area with the nonwindward area is not reasonable for the ER calculation. Terefore, the results of this study were relatively evaluated as the objective CAR in the NB necessary for seeding efect analysis. Based on this sensitivity analysis result, the ER was calculated using ground measurement data assumed to be a true value according to the algorithm shown in Figure 1, in this study. Table 5 shows the calculated ER for each case and TAER calculated by considering the ERA.
In Event 3, the ER was 2.78 mm, which was the largest among the three events and 38.2% compared with that in the SB (Table 5). In Event 2, the ER was the smallest at 0.84 mm (19.4% compared with that in the SB); in Event 1, the ER was 1.74 mm (39.6% compared with that in the SB). Tese results showed that the pure seeding efect can be separated when natural and artifcial rainfall is mixed. Te ERA was also estimated in terms of the amount of ER; in particular, in Event 3, rainfall occurred in a larger area than in other cases ( Figure 4). As a result, a large amount of water resources could be secured: 4.75 million tons for forest fre prevention, 1.30 million tons for drought mitigation, and 24.44 million tons for dust reduction. Tese fndings were signifcant because the seeding efect was quantitatively determined for each event. In addition, the expected amount of water resources can help achieve the experimental purposes. In Event 1, the risk of forest fres in dry areas could be lowered by increasing the humidity around the target area. In Event 2, preparations for drought could be made by increasing the amount of water stored in the Boryeong Dam basin. Te Boryeong Dam supplied an average of 200,000 tons of water per day as of 2020, and the amount that could be secured through the experiment was equivalent to approximately a week of water supply. Te water resources secured in Event 3 could be expected to reduce the concentration of dust in the capital area. Although fundamentally solving drought, forest fre, and dust is difcult, these results are signifcant because we can secure usable water resources to prepare for disasters through cloud seeding experiments. Te analysis results of 20 cloud seeding experiments performed in 2020 in Korea and the experiments corresponding to the three events analyzed in this study are presented in Table 6 as a reference in this study.
Eleven experiments were performed in 2020 (Table 6): six for forest fre prevention, two for drought mitigation, and three for dust reduction. Forest fre prevention experiments were conducted relatively frequently except during the rainy season (June to August) in summer. Te drought mitigation experiment was conducted in May, and the dust reduction experiment was performed in January and November. Tis fnding shows that the risk of forest fres and droughts was frequently observed in 2020 because of the dry climate, and dust problems in autumn were signifcant. If the events were classifed according to the seeding material, four experiments for AgI and seven for CaCl 2 were performed. In the drought mitigation experiment, AgI was not used because the experiments were conducted in spring when temperature was high.
Te amount of ER was at least 0.16 mm (January 29, 2020) to 3.49 mm (January 19, 2020) among the 11 experiments. Te ERA varied according to the NS of each case, and the TAER was determined according to the amount of the ER and the size of the ERA. For the minimum rainfall area (April 27, 2020), the total rainfall was as low as 950 tons. For the TAER, Event 3 showed the largest volume of all cases. Te ER for 11 cases was a total of 12.19 mm, confrming that an average of at least 1.0 mm of rainfall could be enhanced in one experiment. Tis analysis result was signifcant because it could quantitatively verify the efect of cloud seeding in the NS by using the ground measurement data. Fundamentally, it helps achieve experimental purposes by securing domestic water resources in Korea.

Conclusions
In this study, a method for verifying the efect of cloud seeding by analyzing the seeding efect and nonseeding efect areas was developed, and its applicability was evaluated. Numerical models, radar and rain gauge-derived rainfall data, and weather condition data were applied. After the seeding efect and nonseeding efect areas based on the numerical data and wind system of rainfall in the target area were determined, the amount of rainfall that increased in the seeding efect area was calculated. Te amount of water resources that could be secured through the experiment was also determined in this study. Te major results were summarized as follows: (1) Te new analysis method consists of four steps (determination of AT, SB, NB, ER, and TAER). Te AT was determined on the basis of the change in the period of rainfall in the target area as shown in the NS. SB was determined to be centered on the area where the radar-derived rainfall was large in the area with the major change in rainfall of the NS. NB was determined according to the consistency of the simulated EP and observed rainfall bands. When the two directions coincided, the seeding line was included as much as possible in the NB. Te ER was calculated by comparing the characteristics of the two areas, and the fnal TAER was determined by considering the ERA and rainfall density. (2) Tree experimental cases were selected to prevent forest fre, mitigate drought, and reduce dust, and SB and NB were determined on the basis of the wind system to verify the seeding efect. In Events 2 and 3, because the wind direction of rainfall was identical to that of the ER, SB was determined in the area with a strong radar-derived rainfall, and NB was determined at the windward seeding line. Conversely, in Event 1, because the directions of the simulated rainfall and natural system were diferent, the NB was determined in the area where a similar cloud was distributed. Te rain gauges distributed near the coastline were analyzed to exclude orographic efects. During AT, the CR was observed to be greater on the rain gauge than on the radar, and the characteristics of clouds distributed in the two areas were confrmed to be similar. (3) Te amount of ER was calculated for each case by comparing the rainfall in SB and NB, and the amount of water resources that could be secured through the experiment was determined. In Event 1, 1.74 mm of rainfall (4.75 million tons) was generated to prevent forest fres. Furthermore, 0.84 mm (1.30 million tons) for drought mitigation in Event 2 and 2.78 mm (24.44 million tons) for dust reduction in Event 3 were secured. Terefore, an average of 1.0 mm of ER could be secured through the experiments in Korea.
Tese results showed that the amount of ER was calculated for each case. However, data analyzed with a rain gauge varies depending on the SB and NB locations, an error might occur in ER calculation. Since this study determined the area with the high rain rate as the afected area under the assumption that seeding material is totally converted to artifcial rainfall, there is clearly room for error to be refected in the analysis results if enhanced rainfall does not occur. Terefore, in this study, the change in rainfall according to the location of the NB was compared through sensitivity analysis, and the calculated rainfall was evaluated as objective. Together, numerical model uncertainty could also afect the analysis results; therefore, SB was determined on the basis of the change in radar-derived rainfall in this study. A numerical model verifed in previous studies on cloud seeding in Korea [43,48] was used. In addition, the statistical analysis of long-term experimental cases, such as the existing method, can reduce the error in the calculation of the seeding efect, but a case-by-case analysis method is needed for rapid technological development. More scientifc improvement is possible for the next plan by immediately performing an analysis after the experiment.
Te results of this study are signifcant because they can verify the pure seeding efect by isolating natural rainfall from the clouds afected by the seeding material. Of course, since the method of this study cannot completely remove natural variation of precipitation, additional case analysis is required. However, the result of this study is signifcant because it can determine the seeding efect on a case-by-case compared to previous studies. In the short term, the forecasting performance of the numerical model can be verifed; in the long term, water resources should be supplemented through cloud seeding, and the damage caused by forest fres, droughts, and dust must be reduced.
In the future, this study will be used in research felds such as cloud microphysics processes, optimal experimental strategical planning, and improvement of the forecasting performance of numerical models. Furthermore, the results of this study are the efect of an experiment using one aircraft; therefore, increasing the seeding materials or applying several aircraft should be considered to obtain more water resources. Additional ground observation networks and various experimental cases should also be secured to determine the quantitative seeding efect.

Data Availability
Te data that support the fndings of this study are available from the corresponding author from NIMS, KMA, Korea.

Conflicts of Interest
Te authors declare that they have no conficts of interest.