Comparative Application of Rain Gauge, Ground-and Space-Borne Radar Precipitation Products for Flood Simulations in a Dam Watershed in South Korea

: This study presents a comparative analysis of ﬂood simulations using rain gauge, ground-and space-borne radar precipitation products. The objectives are to assess the effectiveness of two radar-based data sources, namely the Radar-AWS Rainrates (RAR) and Integrated Multi-satellite Retrievals for GPM (IMERG), in a dam watershed with gauge observations, and explore the modeling feasibility of integrating the half-hourly IMERG satellite precipitation in regions with ungauged or limited observational area. Two types of HEC-HMS models were developed, considering areal-averaged and spatially distributed gridded data simulations utilizing eight selected storm events. The ﬁndings indicate that the RAR data, although slightly underestimate precipitation compared to the gauge measurements, accurately reproduce hydrographs without requiring parameter adjustments (Nash–Sutcliffe efﬁciency, E NS , 0.863; coefﬁcient of determination, R 2 , 0.873; and percent bias, PBIAS, 7.49%). On the other hand, ﬂood simulations using the IMERG data exhibit lower model efﬁciency and correlation, suggesting potential limitations in ungauged watersheds. Nevertheless, with available discharge data, the calibrated model using IMERG shows prospects for utilization (E NS 0.776, R 2 0.787, and PBIAS 7.15%). Overall, this research offers insights into ﬂood simulations using various precipitation products, emphasizing the signiﬁcance of reliable discharge data for accurate hydrological modeling and the need for further evaluation of the IMERG product.


Introduction
Precipitation is a key meteorological forcing variable in hydrologic modeling for flood simulation. Conventionally, areal precipitation data from rain gauges, estimated using an averaged method, have been adopted as an adequate input for simulating the rainfall-runoff process in watersheds due to its relatively simple application [1]. Recently, remotely sensed precipitation data obtained from ground-and space-borne radar systems have become available, offering high spatiotemporal resolution and wide coverage, particularly in datascarce or ungauged regions [2][3][4]. These advancements contribute to improved accuracy and provide insights into hydrologic modeling for flood simulation and forecasting. In particular, the utilization of satellite-based (i.e., space-borne radar) precipitation data has shown significant potential, especially in ungauged regions and developing countries where ground-based measurements are limited. Notably, the Global Precipitation Measurement (GPM) mission, developed as a specialized observation satellite, has played a crucial role in providing valuable scientific progress and data for areas lacking comprehensive in situ measurements [5][6][7].
The Global Precipitation Measurement (GPM) Core Observatory was launched in February 2014 to provide next-generation global rain and snow observations in near realtime. With many years' worth of calibrated and validated precipitation estimates, GPM's data are being used for scientific studies [5,6]. The Integrated Multi-satellitE Retrievals for GPM (IMERG) precipitation products were released in early 2015 [8]. This high-resolution methods that rely on gauge observations or precipitation estimates from weather radars, such as those conducted by Li et al. [10] and Gilewski and Nawalany [32]. Therefore, in order to expand the application of GPM precipitation data to ungauged watersheds, particularly in flood analysis, it is essential to conduct more comparative studies that assess the level of accuracy in the IMERG product-based flood simulation modeling by comparing it with those of well-established precipitation estimates (i.e., rain gauge observation and ground-borne radar precipitation).
In South Korea, the availability of a robust in situ data infrastructure for rainfall observation, which includes numerous rain gauges and a ground-based radar system for precipitation estimation, provides significant advantages for evaluating the precision of GPM IMERG space-borne radar precipitation products [39,40]. Moreover, South Korea's mountainous topography, covering approximately 65% of the land, introduces an additional dimension to the assessment; the spatial distribution of precipitation during rainfall storm events exhibits high variability, especially in dam watersheds that represent typical mountainous terrain. This intricate pattern of rainfall and the resulting changes in the runoff outflow show complicated characteristics, enabling a comprehensive evaluation that captures the influence of the unique characteristics (i.e., amount and spatial variability) from different types of precipitation products in flood simulation. However, there have been only a few precipitation validation studies conducted in South Korea using GPM IMERG data [39][40][41][42], and none of them have directly applied it to flood hydrological analysis. Kim et al. [39] and Nguyen et al. [40] demonstrated the proof of GPM ground validation studies using the radar network and gauge observations, respectively. Kim et al. [41] assessed precipitation products from GPM using gauge-based precipitation data of Far East Asia during the pre-monsoon and monsoon seasons. Lee et al. [42] implemented a validation of GPM IMERG by comparing it to multiple gauge-based analysis products (e.g., HRPPs and CMORPH) over East Asia. Therefore, it is necessary to conduct a study that explores the application of GPM IMERG precipitation data in hydrologic flood simulation, compares the modeling results with rain gauge observations and ground-based radar precipitation estimates in data-rich regions. With that, the suitability of space-borne satellite radar precipitation data in flood simulations also can be investigated.
Hence, this study aims to perform a comparative application of rain gauge, groundand space-borne radar precipitation products for flood simulations, specifically focusing on evaluating the effectiveness of the half-hourly GPM IMERG final run product using the HEC-HMS model in the Yongdam dam watershed-a representative dam area characterized by mountainous terrain and an abundance of monitoring stations (water level, discharge, etc.) and data in South Korea. By addressing the existing research gaps, this study endeavors to implement a rigorous scientific investigation and analyze the comparative modeling results in regions where observational data are relatively abundant. Furthermore, this study seeks to evaluate the effectiveness of incorporating satellite precipitation data (i.e., GPM IMERG; space-borne radar product) into flood hydrological modeling simulations in regions with limited observational data.

Study Area
The target watershed of this study is the Yongdam Dam, located in the upper reaches of the Geum River in South Korea, between 35 • 35 N to 36 • 00 N and 127 • 20 E to 127 • 45 E. It encompasses Chuncheongnam-do, Jeollabuk-do, and Gyeongsangnam-do, including the counties of Muju, Jinan, and Jangsu. The study area exhibits a temperate climate, and its precipitation patterns are significantly influenced by the East Asian monsoon system, with notable rainfall during the summer months. The watershed area covers 930 km 2 , which accounts for approximately 9.5% of the total Geum River basin. The Yongdam Dam watershed mainly consists of eight sub-watersheds, five inflowing streams, and the Yongdam Dam, which has a height of 70 m, a length of 498 m, and a total storage capacity of 815 million m 3 . Also, it has a total of 22 operational observation stations for the production Water 2023, 15, 2898 4 of 24 of hydro-meteorological data, including water level/flow, precipitation, soil moisture, and evapotranspiration ( Figure 1). watershed mainly consists of eight sub-watersheds, five inflowing streams, and the Yongdam Dam, which has a height of 70 m, a length of 498 m, and a total storage capacity of 815 million m 3 . Also, it has a total of 22 operational observation stations for the production of hydro-meteorological data, including water level/flow, precipitation, soil moisture, and evapotranspiration ( Figure 1).

Land Surface Data
The hydrologic model, which is necessary for flood simulations, is typically constructed based on land surface data (geo-spatial information) such as a Digital Elevation Model (DEM), land cover, and soil data. In this study, a 1-arcsecond DEM with a resolution of approximately 30 m was obtained from the National Geographic Information Institute's Land Information Platform [43]. For land cover and soil data, the 1:25,000 scale subdivision land cover map and the detailed soil map provided by the Ministry of Environment's Environmental Geographic Information Service [44] and the National Institute of Agricultural Sciences' Korean Soil Information System [45] were utilized, respectively. Figure 2 illustrates the secondary datasets of land cover classification, Hydrologic Soil Groups (HSGs), and Curve Numbers processed from the original land surface data for the hydrologic model development. It can be observed that the Yongdam Dam watershed is predominantly composed of mountainous areas and agricultural areas, accounting for 70.6% and 25.9%, respectively.

Land Surface Data
The hydrologic model, which is necessary for flood simulations, is typically constructed based on land surface data (geo-spatial information) such as a Digital Elevation Model (DEM), land cover, and soil data. In this study, a 1-arcsecond DEM with a resolution of approximately 30 m was obtained from the National Geographic Information Institute's Land Information Platform [43]. For land cover and soil data, the 1:25,000 scale subdivision land cover map and the detailed soil map provided by the Ministry of Environment's Environmental Geographic Information Service [44] and the National Institute of Agricultural Sciences' Korean Soil Information System [45] were utilized, respectively. Figure 2 illustrates the secondary datasets of land cover classification, Hydrologic Soil Groups (HSGs), and Curve Numbers processed from the original land surface data for the hydrologic model development. It can be observed that the Yongdam Dam watershed is predominantly composed of mountainous areas and agricultural areas, accounting for 70.6% and 25.9%, respectively. watershed mainly consists of eight sub-watersheds, five inflowing streams, and the Yongdam Dam, which has a height of 70 m, a length of 498 m, and a total storage capacity of 815 million m 3 . Also, it has a total of 22 operational observation stations for the production of hydro-meteorological data, including water level/flow, precipitation, soil moisture, and evapotranspiration ( Figure 1).

Land Surface Data
The hydrologic model, which is necessary for flood simulations, is typically constructed based on land surface data (geo-spatial information) such as a Digital Elevation Model (DEM), land cover, and soil data. In this study, a 1-arcsecond DEM with a resolution of approximately 30 m was obtained from the National Geographic Information Institute's Land Information Platform [43]. For land cover and soil data, the 1:25,000 scale subdivision land cover map and the detailed soil map provided by the Ministry of Environment's Environmental Geographic Information Service [44] and the National Institute of Agricultural Sciences' Korean Soil Information System [45] were utilized, respectively. Figure 2 illustrates the secondary datasets of land cover classification, Hydrologic Soil Groups (HSGs), and Curve Numbers processed from the original land surface data for the hydrologic model development. It can be observed that the Yongdam Dam watershed is predominantly composed of mountainous areas and agricultural areas, accounting for 70.6% and 25.9%, respectively.

Gauged Data
Hourly precipitation and discharge observations were obtained from seven rain gauges and six water level/flow stations, respectively, over the Youngdam Dam watershed, provided by the Ministry of Environment's Water Resources Management Information System [46]. The rainfall amounts measured with these gauges were utilized as input data in hydrological modeling to enable a comparison study with other precipitation products. Specifically, the Thiessen polygon method was employed to calculate the basin average precipitation. Discharge data were collected to facilitate the calibration and validation of the constructed hydrological model by comparing them with the simulated flow outputs. Table 1 provides the location information for each observation station.

Ground-Borne Radar Precipitation Product (RAR)
The Radar-AWS Rainrates (RAR) is a data product developed by the Korea Meteorological Administration (KMA) with the objective of utilizing precipitation forecasts. It provides quantitative estimates of real-time rainfall by applying the Z-R relationship, which is derived from radar reflectivity (Z) and rain gauge measurements (R). The RAR synthesis used in this study incorporates the algorithm of RAR version 2.0, which was enhanced and implemented from October 2015 onwards [47]. This version includes the Local Gauge Correction (LGC) technique, which is an enhancement of the Inverse Distance Weighting (IDW) method. The LGC approach applies weights that account for discrepancies between measurements from Automated Weather Stations (AWSs) and values derived from radar, allowing for the correction of ground radar-based gridded precipitation [48]. Detailed methodology regarding this technique can be found in Lee et al. [49], and Figure 3 illustrates the comprehensive process of the RAR algorithm, incorporating the Local Gauge Correction method [47,50].
Additionally, Table 2 presents specific details regarding the RAR data, including file names, resolution, data structure, and storage format. As indicated in the table, RAR data is primarily provided in a binary format, which requires conversion into a human-readable format (e.g., ASCII) for applications such as flood analysis and hydrological modeling. The RAR data are produced at 10 min intervals and represent a grid-based precipitation synthesis with a spatial resolution of 1 km. For reference, the "Lambert Conformal Conic projection" is adopted as a Projected Coordinate System (PCS) for RAR, which minimizes distortion when representing areas, and it finds extensive usage in hydrological modeling and watershed studies.  Additionally, Table 2 presents specific details regarding the RAR data, including file names, resolution, data structure, and storage format. As indicated in the table, RAR data is primarily provided in a binary format, which requires conversion into a human-readable format (e.g., ASCII) for applications such as flood analysis and hydrological modeling. The RAR data are produced at 10 min intervals and represent a grid-based precipitation synthesis with a spatial resolution of 1 km. For reference, the "Lambert Conformal Conic projection" is adopted as a Projected Coordinate System (PCS) for RAR, which minimizes distortion when representing areas, and it finds extensive usage in hydrological modeling and watershed studies.   The space-borne radar precipitation data utilized in this study is derived from the Integrated Multi-satellitE Retrievals for GPM (IMERG) product, which is produced globally by the NASA Goddard Space Flight Center. IMERG combines data from the core Global Precipitation Measurement (GPM) satellite and a constellation of around 10 supporting satellites equipped with Microwave Imagers [5]. This dataset provides gridded precipitation estimates with a resolution of approximately 10 km (a grid size of 0.1 • × 0.1 • ) within the latitude range of 60 • N to 60 • S in the WGS84 coordinate system. The temporal resolution ranges from as frequent as every 30 min to daily and monthly intervals, with options available for 3 h, 3 days, and 7 days aggregations to meet various research requirements.
IMERG offers data in different stages, namely "Early", "Late", and "Final" runs, catering to various application needs. Each run has a distinct latency period, with durations ranging from a few hours for rapid response applications like flash flood observations to several months for research purposes; approximately 4-5 h (suitable for flash flood observations), 14-15 h (agricultural applications), and around 3 months (research data) [9]. These runs facilitate the timely availability of precipitation information for diverse user requirements. The data production system is depicted in Figure 4, providing insights into the data and runs involved. Map information unsigned char -

Space-Borne Radar Precipitation Product (IMERG)
The space-borne radar precipitation data utilized in this study is derived from the Integrated Multi-satellitE Retrievals for GPM (IMERG) product, which is produced globally by the NASA Goddard Space Flight Center. IMERG combines data from the core Global Precipitation Measurement (GPM) satellite and a constellation of around 10 supporting satellites equipped with Microwave Imagers [5]. This dataset provides gridded precipitation estimates with a resolution of approximately 10 km (a grid size of 0.1° × 0.1°) within the latitude range of 60° N to 60° S in the WGS84 coordinate system. The temporal resolution ranges from as frequent as every 30 min to daily and monthly intervals, with options available for 3 h, 3 days, and 7 days aggregations to meet various research requirements.
IMERG offers data in different stages, namely "Early", "Late", and "Final" runs, catering to various application needs. Each run has a distinct latency period, with durations ranging from a few hours for rapid response applications like flash flood observations to several months for research purposes; approximately 4-5 h (suitable for flash flood observations), 14-15 h (agricultural applications), and around 3 months (research data) [9]. These runs facilitate the timely availability of precipitation information for diverse user requirements. The data production system is depicted in Figure 4, providing insights into the data and runs involved. The IMERG product continually evolves, with version 06 being the most recent iteration. It incorporates advancements and updates to improve the accuracy and reliability of the precipitation estimates. Notably, the integration of the Tropical Rainfall Measuring Mission (TRMM) satellite-based TMPA with IMERG has further enhanced the dataset, The IMERG product continually evolves, with version 06 being the most recent iteration. It incorporates advancements and updates to improve the accuracy and reliability of the precipitation estimates. Notably, the integration of the Tropical Rainfall Measuring Mission (TRMM) satellite-based TMPA with IMERG has further enhanced the dataset, offering a comprehensive satellite-based precipitation record from 2000 to present [8]. This integration enables the analysis of long-term precipitation patterns and trends for a wide range of research applications, including hydrological modeling. The specifics of the IMERG data can be found in Table 3. The data are primarily provided in ".HDF5" or ".nc4" file formats, and the grid images are recorded in the WGS84 geographic coordinate system, allowing for seamless integration with other geospatial data and enabling various geoprocessing tasks such as coordinate projection. These characteristics make the IMERG dataset suitable for hydrological modeling studies and watershed analyses.

Hydrologic Model
HEC-HMS (Hydrologic Engineering Center's Hydrologic Modeling System) is a widely used hydrological model developed by the US Army Corps of Engineers (US-ACE) for applications such as flood analysis and simulation [51]. It is primarily classified as a lumped hydrologic model as it predominantly utilizes relatively simple methods such as the SCS Curve Number, unit hydrograph, and Muskingum routing for estimating effective rainfall, watershed transformation, and river routing calculations, respectively. Although it mainly employs point-based data (i.e., rain gauge observation), it is possible to construct a model that incorporates gridded radar-based rainfall data by utilizing the tools of HEC-GeoHMS (Geospatial Hydrologic Modeling Extension) [52] and the ModClark (modified Clark) method [1,53,54]. Figure 5 depicts the schematic diagram of the ModClark watershed transformation method, which is a modification of Clark's unit hydrograph method [55] to enable the simulation of spatially distributed grid-based rainfall data. In this approach, individual unit hydrographs are computed for each grid cell, taking into account the transitions between cells. The determination of the unit hydrograph for each grid cell involves the consideration of the time of concentration, which relates to the travel time of water through the watershed, as described in Equation (1). These unit hydrographs undergo storage attenuation processes across the entire watershed, as indicated with Equation (2), resulting in the generation of individual cell hydrographs at the watershed outlet. The storage attenuation can be attributed to the overall retention characteristics of the watershed. By convolving these individual cell hydrographs with excess rainfall for each grid cell, the direct rainfall-runoff hydrograph for the entire watershed In this approach, individual unit hydrographs are computed for each grid cell, taking into account the transitions between cells. The determination of the unit hydrograph for each grid cell involves the consideration of the time of concentration, which relates to the travel time of water through the watershed, as described in Equation (1). These unit hydrographs undergo storage attenuation processes across the entire watershed, as indicated with Equation (2), resulting in the generation of individual cell hydrographs at the watershed outlet. The storage attenuation can be attributed to the overall retention characteristics of the watershed. By convolving these individual cell hydrographs with excess rainfall for each grid cell, the direct rainfall-runoff hydrograph for the entire watershed can be obtained. [

Flood Events Selection
For the flood simulations in this study, one or two relatively independent flood events per year from 2014 to 2018 that occurred in the Youngdam Dam watershed were selected. The selection of independent events is necessary due to the adoption of the event-based SCS Curve Number method in the HEC-HMS model for estimating effective rainfall and avoiding inaccurate results that can arise from the use of continuous flood events [4]. Table 4 presents the observed rainfall amounts from seven rain gauges and provides areal average precipitation of the watershed, peak discharge, and total runoff for a total of eight selected flood events. In the years 2015-2017, the selection of relatively smaller-scale floods was inevitable because of prolonged drought conditions. The largest flood event among the selected cases occurred in August 2018, with a total precipitation of approximately 300 mm and a discharge of 2200 m 3 /s. It is important to note that the rainfall and streamflow data used in this study were obtained at an hourly time resolution and were utilized for performing the HEC-HMS model simulation and calibration.

Radar Precipitation Products Processing
Generally, spatially distributed grid-based rainfall data are represented as images in a specific coordinate system and file format. The RAR data utilize the Lambert Conformal Conic projection as its projected coordinate system, while the IMERG data use the WGS84 geographic coordinate system. Therefore, for the purpose of HEC-HMS hydrological modeling, several processing steps are required, including (1) data format conversion, (2) geo-referencing for coordinate system transformation, and (3) DSS file generation. In this study, a total of 2 or 3 data processing steps were performed for each image file. NCL (NCAR Command Language) and Python-based script programs were developed and employed for this purpose. Table 5 outlines the overall data processing steps for spatially distributed grid-based rainfall data, specifically focusing on the procedures for RAR and IMERG datasets. The RAR dataset adopts the "Albers Equal-Area" map projection, while the IMERG dataset utilizes the "ITRF2000" coordinate system, both as the Standard Hydrologic Grid (SHG) for the HEC-HMS model. In the final stage of data processing, the "asc2dssGrid.exe" file provided by the USACE's HEC-GridUtil is utilized to extract ".dss" files, which are directly used as precipitation input data for the HEC-HMS model. Note: 1 The same Python script in ArcGIS (ArcPy) is utilized for both processing steps.

HEC-HMS Model Development
To perform flood simulations using ground-and space-borne radar-based spatially distributed gridded precipitation products, the HEC-HMS model needs to be constructed using ArcGIS applications such as Arc Hydro [56] and HEC-GeoHMS [52]. This allows the incorporation of the ModClark method, which is a grid-based rainfall-runoff analysis approach, into the HEC-HMS model. In this study, in addition to the ModClark model, a model based on Clark's watershed routing method, which has traditionally been applied as one of the lumped hydrologic modeling approaches, was also developed simultaneously. This enables the comparison of the flood simulations using the areal-averaged precipitation through spatial interpolation with Thiessen polygons with point data. For the construction of each model, the delineation of stream network and subbasins from the DEM was performed in Arc Hydro by setting a threshold of 3% of the total watershed area of the Youngdam Dam watershed. The specific calculations for runoff depth, watershed transformation, baseflow, and river reach routing in each hydrologic modeling process were configured using the methods described in Table 6 within the HEC-GeoHMS options.

Flood Simulation and Evaluation
The flood simulations of the constructed models were performed using the hourly gauged and radar-based rainfall data for the selected eight flood events. For each model, the model parameters were adjusted to derive simulated hydrographs that closely matched with the discharges at the six flow measurement stations within the Youngdam Dam watershed. In this study, unlike the typical calibration process of using a single parameter set for rainfall-runoff prediction models, each parameter was individually adjusted for model verification. By implementing this, the approach aimed to assess the effectiveness and accuracy of the spatially distributed gridded rainfall data and the ModClark method in analyzing the flood characteristics of the mountainous dam watershed, compared to the conventional method that uses gauge points for deriving the areal average rainfall. Equations (3)-(5) introduced in the study are statistical indicators used to compare and evaluate the simulated hydrographs against the observed discharge values. These indicators include Nash-Sutcliffe efficiency (E NS ), coefficient of determination (R 2 ), and percent bias (PBIAS (%)) to assess the agreement between the observed and simulated discharge values.
where O i and S i represent the observed and simulated discharge, and O and S represent the mean values of the observed and simulated discharge, respectively.

Results and Discussion
3.1. Precipitation Estimates 3.1.1. Amounts Table 7 compares the amounts of areal average precipitation total derived from seven rain gauges in the Yongdam Dam watershed with ground-and space-borne radar data. The areal-averaged values were calculated using Thiessen polygons for the gauged data, while the radar-based gridded data utilized the arithmetic mean over the grids. The table also presents the maximum and minimum values of the radar grid for each flood event. As evident from the comparison of total amounts in the table, precipitation based on ground-borne radar tends to slightly underestimate compared to the gauged values. On the other hand, space-borne radar precipitation includes some instances of overestimation in relation to the gauge observations. This trend is also apparent in Figure 6, which shows all hourly data for the flood events used in this study. For RAR, the total sum is 689.6 mm, and for IMERG, it is 816.2 mm, resulting in underestimations of 28.2% and 15.1%, respectively, in comparison with the gauged total of 960.9 mm. The correlation coefficients (R 2 ) further provide insight into the relationship between the datasets. RAR illustrates a high correlation with the observations, with values of 0.86, while IMERG exhibits a significant difference in correlation, with a value of 0.46. Therefore, the relevant model parameter values were adjusted during the flood simulations to reduce the discrepancies between the IMERG outputs and the actual observations, and to incorporate them into the analysis.
As evident from the comparison of total amounts in the table, precipitation based on ground-borne radar tends to slightly underestimate compared to the gauged values. On the other hand, space-borne radar precipitation includes some instances of overestimation in relation to the gauge observations. This trend is also apparent in Figure 6, which shows all hourly data for the flood events used in this study. For RAR, the total sum is 689.6 mm, and for IMERG, it is 816.2 mm, resulting in underestimations of 28.2% and 15.1%, respectively, in comparison with the gauged total of 960.9 mm. The correlation coefficients (R 2 ) further provide insight into the relationship between the datasets. RAR illustrates a high correlation with the observations, with values of 0.86, while IMERG exhibits a significant difference in correlation, with a value of 0.46. Therefore, the relevant model parameter values were adjusted during the flood simulations to reduce the discrepancies between the IMERG outputs and the actual observations, and to incorporate them into the analysis. Figure 6. Scatter plots comparing the areal-averaged hourly RAR and IMERG products with gauged data, respectively. Figure 6. Scatter plots comparing the areal-averaged hourly RAR and IMERG products with gauged data, respectively.

Spatial Variability
The spatial distributions of the total cumulative gridded precipitation from groundand space-borne radar for eight flood events, extracted using the aforementioned data processing methods, are depicted in Figure 7. Figure 7a displays the 1 km × 1 km grid-based RAR precipitation total for each individual flood event. The varying precipitation values across different grids indicate that even within the same watershed, there can be regional differences in the rainfall amounts for the same event. While several factors may contribute to this, the similarity in distribution among six out of the eight events suggests that it is primarily influenced by the topographic characteristics of mountainous areas. Therefore, if the areal average precipitation data from rain gauges are used for flood simulations, the lack of deliberation for this spatial distribution may lead to inaccurate analysis results or difficulties in model calibration. On the other hand, Figure 7b represents the spatial distribution of IMERG precipitation data at an approximate resolution of 10 km (0.1 • × 0.1 • grid). Although there are differences in the total precipitation across grids, the spatial resolution of this IMERG dataset is less dense compared to the RAR data, which may pose some limitations for related analyses. Nevertheless, it is deemed to be highly useful for flood analysis in ungauged or global river basins where rainfall data are unavailable. The spatial distributions of the total cumulative gridded precipitation from groundand space-borne radar for eight flood events, extracted using the aforementioned data processing methods, are depicted in Figure 7.  Figure 7a displays the 1 km × 1 km grid-based RAR precipitation total for each individual flood event. The varying precipitation values across different grids indicate that even within the same watershed, there can be regional differences in the rainfall amounts for the same event. While several factors may contribute to this, the similarity in distribution among six out of the eight events suggests that it is primarily influenced by the topographic characteristics of mountainous areas. Therefore, if the areal average precipitation data from rain gauges are used for flood simulations, the lack of deliberation for this spatial distribution may lead to inaccurate analysis results or difficulties in model calibration. On the other hand, Figure 7b represents the spatial distribution of IMERG precipitation data at an approximate resolution of 10 km (0.1° × 0.1° grid). Although there are differences in the total precipitation across grids, the spatial resolution of this IMERG dataset is

Model Development
Two types of HEC-HMS models were constructed with different hydrological modeling approaches as presented in Figure 8. The first model, which applies Clark's method, can simulate areal-averaged precipitation based on rain gauge observations, and it consisted of a total of 25 subbasins. The ModClark model, designed to utilize grid-based RAR and IMERG precipitation, includes 1337 and 76 analysis units of 1 and 10 km grid cells, respectively. In the ModClark model, the calculated flow lengths and the gridded Curve Number for each grid were incorporated as essential information for model development.
In addition, discharge data from the gauging points of six upstream rivers and watershed outlet were also inputted to assess and verify the flood analysis results against the model simulated outputs.

Flood Simulations
In order to investigate the results of flood simulations in relation to model calibration, the initial values of key parameters for each HEC-HMS model are provided in Table 8.
Based on these values, flood simulations were performed using hourly areal-averaged rain gauges and gridded precipitation data (two radar products). Figure 9 represents whether or not the initial parameter values were adjusted for each model after the required calibrations using the eight selected storm events.

Model Development
Two types of HEC-HMS models were constructed with different hydrological modeling approaches as presented in Figure 8. The first model, which applies Clark's method, can simulate areal-averaged precipitation based on rain gauge observations, and it consisted of a total of 25 subbasins. The ModClark model, designed to utilize grid-based RAR and IMERG precipitation, includes 1337 and 76 analysis units of 1 and 10 km grid cells, respectively. In the ModClark model, the calculated flow lengths and the gridded Curve Number for each grid were incorporated as essential information for model development.
In addition, discharge data from the gauging points of six upstream rivers and watershed outlet were also inputted to assess and verify the flood analysis results against the model simulated outputs.

Flood Simulations
In order to investigate the results of flood simulations in relation to model calibration, the initial values of key parameters for each HEC-HMS model are provided in Table 8. Based on these values, flood simulations were performed using hourly areal-averaged rain gauges and gridded precipitation data (two radar products). Figure 9 represents whether or not the initial parameter values were adjusted for each model after the required calibrations using the eight selected storm events. Upon examining the results, it is evident that the models utilizing gridded precipitation, particularly RAR data, required relatively less adjustment of parameters compared to simulations using areal-averaged precipitation. This can be attributed to the fact that in the processes of rainfall-runoff transformation (time of concentration, Tc, and storage coefficient, R, parameters used) and river routing (Muskingum K and X parameters used), which are associated with flow travel times, simulations were achieved without the need for adjustments of parameter values, yielding results similar to the observed hydrographs. In contrast, simulations using rain gauge observations and IMERG gridded precipitation, where the spatial distribution of precipitation across the entire watershed was not sufficiently considered, necessitated adjustments of parameter values related to travel times in almost all cases.
This indicates the increased efficiency of flood discharge analysis when utilizing the detailed grid resolution of spatially distributed precipitation, which allows for the Upon examining the results, it is evident that the models utilizing gridded precipitation, particularly RAR data, required relatively less adjustment of parameters compared to simulations using areal-averaged precipitation. This can be attributed to the fact that in the processes of rainfall-runoff transformation (time of concentration, T c, and storage coefficient, R, parameters used) and river routing (Muskingum K and X parameters used), which are associated with flow travel times, simulations were achieved without the need for adjustments of parameter values, yielding results similar to the observed hydrographs. In contrast, simulations using rain gauge observations and IMERG gridded precipitation, where the spatial distribution of precipitation across the entire watershed was not sufficiently considered, necessitated adjustments of parameter values related to travel times in almost all cases.
This indicates the increased efficiency of flood discharge analysis when utilizing the detailed grid resolution of spatially distributed precipitation, which allows for the incorporation of localized rainfall patterns and amounts into model simulations. Particularly, in watersheds dominated by mountainous terrains, such as dam watersheds, where there are significant variations in rainfall across different regions, the introduction of data that considers spatial distribution, such as radar precipitation, is crucial. However, it should be noted that gridded precipitation data at the grid-cell level still tends to underestimate or overestimate the actual observed values. Therefore, simulations using such data require adjustments of parameter values associated with precipitation quantities, and the IMERG product exhibited relatively larger differences in this aspect. Figure 10 shows the simulated hydrographs, which provide an overall flood analysis simulation result for this study. With a few exceptions, the majority of the simulated hydrographs closely resemble the observed hydrographs. In terms of the statistical comparison of the simulation results (Table 9), the model simulations using gauged data exhibit an average Nash-Sutcliffe efficiency (E NS ) of 0.895, a coefficient of determination (R 2 ) of 0.906, and a percent bias (PBIAS) of 7.42%. For the ground-borne radar precipitation used model, the E NS is 0.863, R 2 is 0.873, and PBIAS is 7.49%. Although there are some differences, as mentioned earlier, the simulations of RAR precipitation data perform well without the need for parameter adjustment, reaffirming the utility of detailed grid-based precipitation data for flood discharge simulations. Additionally, the space-borne radar data simulations also demonstrate relatively satisfactory results with an E NS of 0.776, an R 2 of 0.787, and a PBIAS of 7.15%, despite the lower correlation with the observed values. Note: * Average: ENS and R 2 arithmetic mean; PBIAS arithmetic mean of absolute value. Figure 10. Hydrographs of flood simulation for the developed models (# stands for event number); X-axis represents simulation time (hours) and Y-axis represents discharge (m 3 /sec). Table 9. Statistical results of flood simulation for the developed models. In addition, Figure 11 and Table 10 compare the hydrographs and evaluation statistics of the simulated results for each water level/flow gauge location. In the comparison with individual flow discharge observations, the models using RAR ground-borne radarbased gridded precipitation exhibit good agreement with the observed values, with an average E NS ranging from 0.53 to 0.92, without the need for adjustment of parameter values related to travel times. However, for the models using areal mean of the gauged data and IMERG space-borne gridded radar precipitation, the average E NS ranges from 0.66 to 0.87 and 0.43 to 0.84, respectively. In the case of IMERG, the precipitation grid is relatively larger compared to RAR, which indicates that the spatial information of the applied precipitation in each watershed is not sufficiently reflected, leading to unsatisfactory results. For reference, Figure 11 illustrates only the hydrographs for storm event #1 adopted in this study, while the hydrograph results for the remaining seven events (Figures A1-A7) have been included in Appendix A. and IMERG space-borne gridded radar precipitation, the average ENS ranges from 0.66 to 0.87 and 0.43 to 0.84, respectively. In the case of IMERG, the precipitation grid is relatively larger compared to RAR, which indicates that the spatial information of the applied precipitation in each watershed is not sufficiently reflected, leading to unsatisfactory results. For reference, Figure 11 illustrates only the hydrographs for storm event #1 adopted in this study, while the hydrograph results for the remaining seven events (Figures A1-A7) have been included in Appendix A. Figure 11. Hydrographs of flood simulation for the developed models at each flow gauge; X-axis represents simulation time (hours) and Y-axis represents discharge (m 3 /sec). The dashed line represents the aggregated values of individually measured discharges, accounting for 93% of the observed discharge at the outlet point (storm event #1).   Additionally, from Table 9, it can be inferred that relatively small catchments may experience a greater incidence of missing or inaccurate discharge data compared to the measured data at the outlet point, which poses challenges in calibrating the models. However, it is crucial to acknowledge that the same occurrence can also arise due to errors in the precipitation data used in this study.

Ground-Borne Radar Data Simulation
Consequently, the flood analysis utilizing various precipitation datasets with diverse production sources in this study revealed simulation results that are not significantly different from the major findings of previously discussed studies in the introduction; notably less accurate for the satellite-based data (i.e., space-borne radar precipitation). However, considering the primary objective of this research, which is to evaluate the usability of satellite-based precipitation products for flood applications, the anticipation of suboptimal simulation results was intrinsic. Nonetheless, for watersheds with available observed discharge data, this study performed an investigation to improve the usability of satellite precipitation data. This was achieved through hydrological model calibration, with the aim of determining the extent to which it can be effectively employed for accurate flood analysis.

Conclusions
The present study conducted a comparative analysis of flood simulations by utilizing rain gauge, ground-and space-borne radar precipitation products. The primary objective was to assess the efficacy of incorporating two radar-based gridded data sources, namely RAR and GPM IMERG, in a dam watershed where gauge observations are relatively abundant. Additionally, this study aims to seek the feasibility of utilizing the half-hourly IMERG product in regions characterized by ungauged or limited availability of observed data. To achieve these objectives, data processing programs and two types of HEC-HMS models were developed for areal-averaged and spatially distributed gridded data based on flood simulations, focusing on the selected eight rainfall storm events. The key findings are as follows: (1) Overall, the Korea Meteorological Administration's ground-borne weather radar data (RAR) exhibited a slight underestimation of precipitation values compared to gauge point observations, with a difference of 28.2% across all eight flood events (R 2 0.86).
Consequently, adjustments to related model parameter values were necessary for flood analysis using the RAR data. However, due to the spatial distribution advantages of radar rainfall, aspects related to rainfall-runoff transformation and river channel routing, such as travel time, could be simulated to generate hydrographs that closely resembled the observed discharge without the need to adjust the related parameter values; (2) For space-borne (i.e., satellite-based) IMERG precipitation data, the total observed amount was underestimated by 15.1% across the eight flood events. However, the correlation coefficient (R 2 ) was 0.46, indicating significant differences from the gauged data. Even with adjustments to parameter values, flood simulations using the IMERG product demonstrated relatively lower correlation and model efficiency compared to the observations. This indicates a significant limitation in using the half-hourly IMERG data for flood modeling in ungauged watersheds. Despite this constraint, it is still inferred that in locations where some discharge data are available, the utilization of the model through verification and calibration is feasible; (3) When comparing the flood simulation results using the conventional method based on rain gauge observations with those using weather radar and satellite-based precipitation data, the models utilizing both radar data sources exhibited an average Nash-Sutcliffe efficiency (E NS ) of 0.863 and 0.776, an R 2 of 0.873 and 0.787, and a percent bias (PBIAS) of 7.49% and 7.15%, respectively. The model using the arealaveraged values showed an E NS of 0.895, an R 2 of 0.906, and a PBIAS of 7.42%. Although there were some differences, simulations using the RAR data demonstrated relatively satisfactory performance without adjusting parameter values, confirming their utility and efficiency; (4) Despite varying ratios ranging from 26.0% to 82.2% depending on antecedent rainfall conditions, the analysis of the eight selected flood events revealed the characteristic of watershed flood discharge with an average runoff ratio of 52.5%.
Finally, the aforementioned research results could be compared due to the availability of reliable discharge data obtained through extensive measurements in the upstream reaches of the watershed. Most of all, the presence of observational data is crucial in assessing the effectiveness of any hydrological model. Consequently, the utilization of hourly discharge data from multiple gauge observation points in this simulation of flood discharge can be considered as a noteworthy case study in the modeling of watershed flood management. In addition, although the GPM IMERG half-hourly product has not yet been established as a reliable dataset for hydrological modeling through this research, it is essential to apply it in diverse regions to evaluate its performance with additional case studies. However, as stated earlier, if discharge measurement data are accessible for the rivers within the specific watershed of interest, it is feasible to calibrate the model using this data, allowing for the integration of the IMERG dataset into flood discharge forecasting with sufficient reliability.
The hydrological modeling using the radar-based gridded rainfall data developed in this study is planned to be implemented in an ongoing research project for estimating inflow to ungauged dam reservoirs. Through this application, further research aims to facilitate the integration of hydrological models into ungauged and data-limited dam reservoir waterbody monitoring studies, utilizing other satellite data (optical or SAR; Synthetic Aperture Radar imageries). This mutual integration will be used as a framework to complement the current research efforts and establish linkages with hydrological modeling for inflow estimation in ungauged watersheds. The integration will also involve analyzing and estimating water levels, surface areas, and storage capacities obtained from satellite imagery data, thereby enhancing the accuracy and completeness of the analyses and estimations.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Acknowledgments: Hydrometeorological data support from the Korea Water Resources Corporation (K-water) for this work is gratefully acknowledged.

Conflicts of Interest:
The author declares no conflict of interest.