Development and Veriﬁcation of the Available Number of Water Intake Days in Ungauged Local Water Source Using the SWAT Model and Flow Recession Curves

: Climate change signiﬁcantly affects water supply availability due to changes in the magni-tude and seasonality of runoff and severe drought events. In the case of Korea, despite a high water supply ratio, more populations have continued to suffer from restricted regional water supplies. Though Korea enacted the Long-Term Comprehensive Water Resources Plan, a ﬁeld survey revealed that the regional government organizations limitedly utilized their drought-related data. These limitations present a need for a system that provides a more intuitive drought review, enabling a more prompt response. Thus, this study presents a rating curve for the available number of water intake days per ﬂow, and reviews and calibrates the Soil and Water Assessment Tool (SWAT) model mediators, and found that the coefﬁcient of determination, Nash–Sutcliffe efﬁciency (NSE), and percent bias (PBIAS) from 2007 to 2011 were at 0.92%, 0.84%, and 7.2%, respectively, which were “very good” levels. The ﬂow recession curve was proposed after calculating the daily long-term ﬂow and extracted the ﬂow recession trends during days without precipitation. In addition, the SWAT model’s ﬂow data enables the quantitative evaluations of the number of available water intake days without precipitation because of the high hit rate when comparing the available number of water intake days with the limited water supply period near the study watershed. Thus, this study can improve drought response and water resource management plans.


Introduction
Climate change threatens hydrologic systems, availability of water supply resources, and water supply capabilities [1,2]. In the case of Korea, the annual average precipitation is around 1277 mm, which is 1.6 times higher than the world average (807 mm), while its water intake, at 34%, is also among the highest in Organization for Economic Cooperation and Development (OECD) countries [3]. Moreover, most of the precipitation is concentrated in the monsoon season, and the steep incline of its rivers and the low water-holding capacity of the surface soil cause water to be discharged in high amounts, making it difficult to address water shortage and drought in Korea [4]. Furthermore, local precipitation, flow, and evapotranspiration rates in Korea have changed because of global warming and climate change, which have increased the possibility of extreme drought [5]. In 2018, Korea recorded a high water supply ratio of 99.2%. Nevertheless, the number of people subject to local water supply restrictions increased over the years, from 6000 in 2011, 44,000 in 2013, to 122,000 in 2015 [6].
This study selected the Hancheon Basin, which serves as the Yecheon Water Intake Station's water source in Yecheon-gun, Gyeongsangbuk-do. The Geographic Information System (GIS) data and metrological data were collected in the watershed to operate the SWAT model. The validation and calibration of the SWAT model parameters were done by dividing the watershed into several areas using the flow observation points of the flood control office as exit points and defining the sub-watersheds to include the water source studied in this paper. Subsequently, the resulting long-term flow data were used to generate an available number of water intake days curve at the water source, which was then validated by comparing it with actual drought cases. The purpose of this study is to draw a flow recession curve quantitatively at a local water source and use the curve to develop an available number of water intake days curve for drought control. The overall flow of the study is represented in Figure 1.
Water 2021, 13,1511 This study selected the Hancheon Basin, which serves as the Yecheon Water I Station's water source in Yecheon-gun, Gyeongsangbuk-do. The Geographic Inform System (GIS) data and metrological data were collected in the watershed to opera SWAT model. The validation and calibration of the SWAT model parameters were by dividing the watershed into several areas using the flow observation points of the control office as exit points and defining the sub-watersheds to include the water s studied in this paper. Subsequently, the resulting long-term flow data were used to erate an available number of water intake days curve at the water source, which was validated by comparing it with actual drought cases. The purpose of this study is to a flow recession curve quantitatively at a local water source and use the curve to de an available number of water intake days curve for drought control. The overall fl the study is represented in Figure 1.

Description of SWAT
The SWAT is a program developed by Dr. Jeff Arnold of the United States Department of Agriculture Agricultural Research Service (USDA ARS) for analyzing long-term hydrologic and water quality changes in large and complex watersheds with diverse soil properties and land uses [14]. The SWAT model is a semi-distributed long-term flow model that is more efficient as a dry season model compared with the Tank model and is more reliable than other semi-distributed models, such as the Precipitation-Runoff Modeling System (PRMS) model and the Hydrological Simulation Program-Fortran (HSPF) model [45,46]. Hence, this model was selected for this study as it has also been tested for its applicability to various areas in Korea and universal applicability to studies on surface flow, riverbed, lake, and underground water.
The logic of the SWAT model incorporates the characteristics of various models developed at the USDA ARS. Based on the Simulator for Water Resources in Rural Basins (SWRRB) model developed for nonpoint pollutant analysis, the SWAT model combines the logic of the Routing Outputs to Outlet (ROTO) model developed for tracking river channel flow and sediments [47,48]. Given how the first model's developed in the early 1990s, the model has been upgraded from SWAT 94.2 to 2009 through various bug fixes and optimizations, and this study used the latest version, SWAT 2012. In the SWAT model, water balance is the most important element of all interactions in the watershed. For accurate hydrology and water quality prediction, the model's hydrologic circulation should be identical to the watershed's natural phenomenon [49]. The simulation logic for hydrologic circulation in the SWAT model consists largely of the land phase and the water or routing phase. The SWAT model controls the loads on water, sediments, nutrients, and main channels of each sub-watershed using the land phase logic and uses the water or routing phase to define the movement of water and sediments from the river to the exit through the river network. The SWAT model's indicator phase defines the precipitation, evaporation and transpiration, surface flow, lateral flow, baseflow, and groundwater flow using the Water Balance Equation (1). In other words, when precipitation reaches the surface, it is either blocked by plants or falls on the surface to permeate into the soil layer or be discharged. In the short term, the discharged water gathers rapidly to form a small river, and the water that permeated into the soil layer either evaporates or permeates deeper underground to return as surface water. The SWAT model simulates the hydrologic cycle by hydrologic response unit (HRU).
where SW t is the final soil water content (mm), SW 0 is the initial soil water content on day i (mm), t is the time (days), R day is the amount of precipitation on day i (mm), Q sur f is the amount of surface runoff on day i (mm), E a is the amount of evapotranspiration on day i (mm), ω seep is the amount of water entering the vadose zone from the soil profile on day i (mm), and Q gw is the amount of return flow on day i (mm).

Calibration and Validation of SWAT
The SWAT model and other precipitation runoff models are used to estimate past events, to predict future events, or in cases where observational data are not sufficient. The model requires parameter calibration based on comparison with actual observations. In addition, statistics assessment is mainly used to validate the calibrated parameters and assess how well they represent the actual measurements. The statistical indicators used for the calibration and validation of the hydrologic models were analyzed differently depending on their users. Moriasi et al. [50] established three performance assessment criteria for statistical indicators based on the findings of various previous studies: coefficient of determination, Nash-Sutcliffe efficiency (NSE), and percent bias (PBIAS).
The coefficient of determination (R 2 ) compares the suitability between the true and predicted values of a hydrologic model. It is calculated as shown in Equation (2), and a value closer to 1.0 represents a more optimal value. Moriasi et al. [50] proposed model findings regarding flow (including streamflow, surface flow, and baseflow), sediment, and nitrogen (as a water quality indicator) as criteria for daily, monthly, and yearly statistical performance assessments. The assessment criteria related to flow were presented as very good (R 2 > 0.85), good (0.75 < R 2 ≤ 0.85), satisfactory (0.60 < R 2 ≤ 0.75), and unsatisfactory (R 2 ≤ 0.60).
where Q i is observed flow (m 3 /s), S i is simulated flow (m 3 /s), and Q and S are the average flow of each analysis period derived from the observed and simulated, respectively. Nash-Sutcliffe efficiency (NSE) compares the reflection of tendencies by the model's true and predicted values. It is calculated as shown in Equation (3), and a value closer to 1.0 represents a more optimal value. Moriasi et al. [50] presented the assessment criteria related to flow as very good (NSE > 0.80), good (0.70 < NSE ≤ 0.80), satisfactory (0.50 < NSE ≤ 0.70), and unsatisfactory (NSE ≤ 0.50).
where Q i is the observed flow (m 3 /s), S i is simulated flow (m 3 /s), and Q is the average observed flow of each analysis period. Percent bias (PBIAS) compares the sums of the true and predicted values of a hydrologic model. It is calculated as shown in Equation (4), and a value closer to 0 represents a more optimal value. Moriasi et al. [50] presented the assessment criteria related to flow as very good (PBIAS ≤ ±5), good (±5 < PBIAS ≤ ±10), satisfactory (±10 < PBIAS ≤ ±15) and unsatisfactory (PBIAS ≥ ±15).
where Q i is the observed flow (m 3 /s), and S i is simulated flow (m 3 /s ).

Description of Flow Recession Curve
The flow recession curve is used widely in hydrologic analysis. After precipitation, surface flow caused by the precipitation in the watershed decreases, along with a base flow consisting mainly of the reservoir water permeated or stored underground, as shown in Equation (5) [31].
where Q 0 is the initial flow (m 3 /s), Q is the flow after t hours (m 3 /s), and k is the ratio of Q after an hour at Q 0 , and t is time.
Kullman [39] suggested that during a prolonged period with no or little precipitation, the flow within the watershed declines, as shown in Figure 2. Based on his findings, this study used the daily flow hydrologic curve calibrated for the peak and base flow sections by comparing them with the actual measurements of the SWAT model and assumed that the non-precipitation date is the starting point of the initial flow. In addition, as suggested by Bako and Hunt [37], and Vogel and Kroll [38], recession segments were selected from the daily long-term flow hydrologic curve of the calibrated and validated SWAT model, and each recession segment was analyzed statistically to define the recession curve equation.

Description of Available Number of Water Intake Days
This paper defined the number of water intake days as the number of days of available flow that can be taken at the water source. This study used and enhanced the recession curve concept proposed by Kullman to calculate the number of water intake days and identify the periods during which the number of days between rainfalls exceeds a certain threshold [39]. The average recession curve of individual recession events was defined ( Figure 3) and normalized, and the master recession curve was proposed by extending the curve to the point where the end of the curve meets the drought flow (0.354 m 3 /s).

Description of Available Number of Water Intake Days
This paper defined the number of water intake days as the number of days of available flow that can be taken at the water source. This study used and enhanced the recession curve concept proposed by Kullman to calculate the number of water intake days and identify the periods during which the number of days between rainfalls exceeds a certain threshold [39]. The average recession curve of individual recession events was defined ( Figure 3) and normalized, and the master recession curve was proposed by extending the curve to the point where the end of the curve meets the drought flow (0.354 m 3 /s).

Description of Available Number of Water Intake Days
This paper defined the number of water intake days as the number of days of available flow that can be taken at the water source. This study used and enhanced the recession curve concept proposed by Kullman to calculate the number of water intake days and identify the periods during which the number of days between rainfalls exceeds a certain threshold [39]. The average recession curve of individual recession events was defined ( Figure 3) and normalized, and the master recession curve was proposed by extending the curve to the point where the end of the curve meets the drought flow (0.354 m 3 /s).  In addition, in accordance with the findings of Tallaksen [34] that recession curves progress at different speeds on account of seasonal weather changes, the seasonal recession curves of the study area were defined. Seasons were defined as spring, summer, fall, and winter, and, in accordance with the methods for seasonal meteorological division proposed by Lee [51], the twelve months were grouped into the respective seasons: spring (March to May), summer (June to August), autumn (September to November), and winter (December to February).
The water source's minimum water intake was defined as the flow until the lowest level that can be taken from the river (local water source) to the water intake station through the intake gate (that is, drought flow). This paper calculated the number of water-intake days in the following steps.
First, this study used the semi-distributed SWAT model to divide sub-watersheds, including the study's watershed of interest, and calculated daily long-term flow by calibrating and validating the parameters.
Second, the daily long-term flow was used to define a flow recession event with no precipitation and single-out recession segments.
Third, the selected recession events were listed to align the recession curve equations, which were used to calculate the master recession curve equations for the entire period and the respective seasons.
Lastly, considering the drought flows and the water demand for the study water source calculated through the flow analysis using the SWAT model, the available number of water intake days curves was developed.

Selection of the Study Watershed
This study's watershed of interest is the Hancheon watershed, which is located between 128 • 37 -128 • 55 E, 36 • 20 -36 • 28 N in the Naeseong River section of the Nakdong River System (Figure 4). Its total area is 192.4 km 2 , 29.7 km flow channel extension, an average width of 6.5 km, and a 0.2 watershed shape factor. The watershed at the Yecheon Water Intake Station is located near 432 Baekjeon-ri, Yecheon-gun, Gyeongsangbuk-do, Korea. The station takes riverbed water through concrete collection conduits at the Hancheon River, which serves as the main water source. The Yecheon Water Intake Station, having a capacity of 10,800 m 3 per day, supplies to the Yeocheon-eup and Yucheon-myeon areas and the Homyeong-myeon area through the Homyeong Reservoir. In the Yecheon-gun area, the water supply population size is 17,447 persons, and the daily water supply per person stands at 310 L per capita per day (lpcd). The water demand, calculated by multiplying the water supply population with the daily water supply per person, is 4409 m 3 per day.

Meteorological Data
The framework plan for the Hancheon River and other high-level plans were reviewed to collect meteorological data. The Korea Meteorological Administration (KMA) data regarding the study watershed were reviewed through GIS analysis, and the Yeongju and Mungyeong Weather Stations were selected, from which a 13-year data on daily temperature, precipitation, wind, and humidity were collected (1 January 2006 to 1 January 2018). The analysis found that the average highest temperature in the Yecheon Watershed was 30.1 • C, and the lowest temperature was −7.6 • C. The annual average precipitation and the highest precipitation (in July) were 1274.5 mm and 318.4 mm, respectively. The average humidity was 65%, and the average wind speed was 2.5 m/s (Table 1).

Meteorological Data
The framework plan for the Hancheon River and other high-level plans were reviewed to collect meteorological data. The Korea Meteorological Administration (KMA) data regarding the study watershed were reviewed through GIS analysis, and the Yeongju and Mungyeong Weather Stations were selected, from which a 13-year data on daily temperature, precipitation, wind, and humidity were collected (1 January 2006 to 1 January 2018). The analysis found that the average highest temperature in the Yecheon Watershed was 30.1 °C, and the lowest temperature was −7.6 °C. The annual average precipitation and the highest precipitation (in July) were 1274.5 mm and 318.4 mm, respectively. The average humidity was 65%, and the average wind speed was 2.5 m/s (Table 1).

Geographical Information System (GIS) Data
This study put together a soil map, a land use map, and topographical data, such as digital elevation model (DEM) data, to operate the SWAT model. The highest, lowest, and average DEM of the Yecheon Watershed derived through the watershed division were 1100 m, 79 m, and 381 m, respectively (Figure 5a). The soil map plays a crucial role in the surface flow, evaporation, and circulation of soil water and underground water in the SWAT model. This study used the 1:25,000 scale soil map provided by the Rural Development Administration (RDA) to build input data for the SWAT model in the GIS format (Figure 5b). In addition, the soil properties, encroachments, effective soil depths, and drainage grades derived from the RDA's Soil Environment Information System [52] were matched with the soil information, and the soil moisture and field capacity were calculated using the soil property estimation method proposed by Saxton et al. [53]. In the soil map developed for the Yecheon Watershed, the soil series distribution rates were 59.8% for SAMGAG (mountain and forest soil, granite, and granite-Gneiss), 14.2% for SEOGTO (upland soil, granite, and granite-gneiss), and 7.5% for SACHON (paddy soil, granite, and granite-gneiss). As for the land use map, this study used the 1:25,000 land cover map provided by the ME Environmental Geographical Information System [54], and the ME's land use classifications were reorganized in alignment with the SWAT model input codes. The reorganization results are shown in Figure 5c. In the land use map, mixed forests accounted for the highest percentage (68.6%), followed by crop fields (26.4%), residential areas (2.6%), other green lands (1.5%), marshes (0.4%), and basins (0.5%). The inclines of the watershed were calculated from the DEM. Inclines were categorized into three groups: 0-30; 30-60; and higher than 60. The percentage of the 30-60 inclines was the highest at 51%, followed by 0-30 (37%), and higher than 60 (12%). The incline map, soil map, and land use map were used to calculate the SWAT model's hydrologic response units (HRUs), which resulted in 2,274 HRUs (Figure 5d). Water 2021, 13, 1511 10 of 24

Flow Data
The Yecheon Water Level Gauging Station of the MOE Nakdong River Flood Control Office is located 2 km downstream from the Yecheon Water-Intake Station, at around 310-19 (Yecheongyo Bridge) Nambon-ri, Yecheon-eup, Yecheon-gun, Gyeongsangbuk-do, Korea. The Yecheon Water Level Gauging Station measures water level with T/M equipment. The station holds daily data from January 1963 to 2020. Despite the long period of measurement, the station lacks flow measurement data. Data before 2007 are not available because of the opening of weir gates. For this reason, this study selected the daily flow data calculated with the 2007 water level flow rating curve equation proposed in the Hydrologic Annual Report, except for the data on the days when the weir gates were open.
Even though a water level gauging station in the downstream section of the Hancheon River exists, this study conducted on-site measurement of flow rates to determine drought flow at the water-intake source, calculate the drought reference values (drought water level and drought flow), and collect the water level and flow data used for the validation of the SWAT model. Specifically, depths and flow rates were calculated based on the number of verticals. Flow rates were calculated using the midsection method, which calculates the flow rates and sectional areas of the segments, other than the first and the last segments, and sums up the respective segments' flow rates. Table 2 shows the flow measurement results between December 2016 and April 2017. Measurements were taken two to four times at certain intervals depending on the situation on-site.

Division of Subbasins
The DEM is a crucial element in determining the shape of a watershed and water movement direction. This study extracted the 7111 (intermediate contour: 2 m) and the 7114 (index contour: 10 m) from the 1:5000 digital evaluation map from the National Geographic Information Institute (NGII) and converted them to 5 m resolution GIS Raster files. The converted DEM data were then used to analyze the flow direction and accumulation of the river to create a river network. Then, 24 sub-watersheds were identified by setting two flow observation points: the Yecheon Water-Intake Station and the MOE Water Level Gauging Station (Yecheon Bridge) ( Figure 6).

Calculation of Long-Term Flow Using the SWAT Model
This study reviewed the 10-year water level and flow data at the MOE Yecheon Water Level Gauging Station to acquire reliable daily flow data and determined whether a water level flow rating curve exists. For this study, the 2007 water level data were determined to be available for removing artificial factors, such as the opening of weir gates, and the water level flow rating curve equation was used to calculate daily flow rates. To calibrate the flow rate parameter, as shown in Figure 7, the year 2006 was selected as the warm-up period, and the calibration period (2007)   Three parameters affecting the surface, interim, and baseflow, respectively, were selected from previous studies and the SWAT manual for parameter calibration (CN2, ESCO, CANMX, SOL_K, SLSOIL, LAT_TIME, GW_DELAY, GWQMN, ALPHA_BF). After analyzing each parameter's sensitivity, five of them were selected for the study, as shown in Table 3 (CN2, ESCO, GW_DELAY, GWQMN, ALPHA_BF). Then, the peak segment of the flow curve was calibrated closer to the actual measurements. The flow rates in the early conditions were overestimated above the actual measurements. The CN2 value, which represents the SCS flow curve indicator in AMC-2, was reduced by 25%, and the soil evaporation compensation factor (ESCO) was lowered from the initial 0.95 to 0.1 to address the issue. Using the calibrated parameters, the flow at the recession segment in the flow curve was adjusted closer to the actual measurement. To that end, the delay time for aquifer recharge was extended from 31 to 80 days. In addition, the base flow recession constant was raised from 0.048 day to 0.8 day, and the threshold water level in a shallow aquifer for base flow was raised from 1000 mm to 2000 mm by 50 mm increments.  Three parameters affecting the surface, interim, and baseflow, respectively, were selected from previous studies and the SWAT manual for parameter calibration (CN2, ESCO, CANMX, SOL_K, SLSOIL, LAT_TIME, GW_DELAY, GWQMN, ALPHA_BF). After analyzing each parameter's sensitivity, five of them were selected for the study, as shown in Table 3 (CN2, ESCO, GW_DELAY, GWQMN, ALPHA_BF). Then, the peak segment of the flow curve was calibrated closer to the actual measurements. The flow rates in the early conditions were overestimated above the actual measurements. The CN2 value, which represents the SCS flow curve indicator in AMC-2, was reduced by 25%, and the soil evaporation compensation factor (ESCO) was lowered from the initial 0.95 to 0.1 to address the issue. Using the calibrated parameters, the flow at the recession segment in the flow curve was adjusted closer to the actual measurement. To that end, the delay time for aquifer recharge was extended from 31 to 80 days. In addition, the base flow recession constant was raised from 0.048 day to 0.8 day, and the threshold water level in a shallow aquifer for base flow was raised from 1000 mm to 2000 mm by 50 mm increments.
The calibration resulted in that the daily flow hydrologic curves and the scatterplot graphs are shown in Figure 8a,b, which follow the overall actual measurements' tendencies.   The determination coefficient for comparing the true and predicted values of the SWAT model in 2007 (R 2 ) was 0.89, which is "very good" according to the statistical performance assessment criteria proposed by Moriasi et al. [50]. The value of NSE, which compares the reflection of tendencies between the true and predicted values of the SWAT model, was 0.88, which is "very good." In addition, PBIAS, which compares the sums of the true and predicted values of the SWAT model, was -3.1%, which is within 5% of the total and classified as "very good." The findings confirm that the parameter calibration of the SWAT model was carried out appropriately.
Then Figure 9e,f show the daily flow hydrologic curves and scatterplot graphs between 2007 and 2011. Throughout the calibration/validation period, the average determination coefficient (R 2 ) was 0.92 and NSE was 0.84, which are "very good" under the statistical performance assessment criteria. PBIAS was 7.24%, which is "good" under the statistical performance assessment criteria. The average annual precipitation was 1255 mm, and the total flow was 697 mm or 56% of the precipitation. The evapotranspiration was 495 mm or 39% of the precipitation. In the SWAT model, the total flow is the sum of the direct and indirect flow. The annual average of direct flow was 532 mm or 42% of the precipitation, and the annual average indirect flow was 165 mm or 13% of the precipitation.
A comparison of the flow measurements between December 2016 and May 2017 at the Yecheon Water Intake Station verified that the absolute values and tendencies are similar. In addition, an analysis of the water flow in the study watershed (Yecheon Water-Intake Station) using the long-term (10 years, 2007-2017) flow derived from the SWAT model, the abundant flow rate was 3.728 m 3 /s, the normal flow was 1.521 m 3 /s, the low flow rate was 0.778 m 3 /s, and the drought flow rate was 0.354 m 3 /s (Table 4). Figure 10 shows the flow comparison between the Yecheon Water Intake Station and the Yecheon Bridge.

Calculation Results of Water Supply Available Days Curve
Based on the daily long-term flow data from the calibrated and validated SWAT model, the recession events with 15 or more days of recession without precipitation were identified using the recession curve concept proposed by Kullman [39]. To that end, this study used the daily long-term flow calculations at the Yecheon Water Intake Station between 2007 and 2017, which resulted in 43 recession events. The 43 events were distributed across the four seasons as proposed by Tallaksen [34]. Five events were identified between March and May, which were used to define the master recession curve equation for spring and the determination coefficient ( = 0.99). The master recession curve equation for each season was calculated using the same method, as shown in Figures 11 and  12. Table 5 lists the determination coefficient and the available number of water-intake days curve equation for each season. The calculation found that the curves are steeper in spring and summer and slower in autumn and winter. The overall determination coefficient was 0.98. By season, the winter determination coefficient was the highest at 0.91, followed by summer (0.93), autumn (0.94), and spring (0.99).

Calculation Results of Water Supply Available Days Curve
Based on the daily long-term flow data from the calibrated and validated SWAT model, the recession events with 15 or more days of recession without precipitation were identified using the recession curve concept proposed by Kullman [39]. To that end, this study used the daily long-term flow calculations at the Yecheon Water Intake Station between 2007 and 2017, which resulted in 43 recession events. The 43 events were distributed across the four seasons as proposed by Tallaksen [34]. Five events were identified between March and May, which were used to define the master recession curve equation for spring and the determination coefficient (R 2 = 0.99). The master recession curve equation for each season was calculated using the same method, as shown in Figures 11 and 12. Table 5 lists the determination coefficient and the available number of water-intake days curve equation for each season. The calculation found that the curves are steeper in spring and summer and slower in autumn and winter. The overall determination coefficient was 0.98. By season, the winter determination coefficient was the highest at 0.91, followed by summer (0.93), autumn (0.94), and spring (0.99).

Comparison of Drought Cases and Applicability Evaluation
Based on the water supply restriction in the Nakdong River Watershed, which includes the Yecheon Watershed, this study determined the number of days of water supply restriction in May and June 2015 ( Figure 13, Table 6). Then, the findings were used to compare the time series of the available number of water intake derived from the available number of water-intake days curves. The findings are as shown in Figure 14. All water supply restriction cases began in the May 31-June 18 period, when the available number of water-intake days in the Yecheon Watershed dropped below the dry season level. The findings indicate a high level of accuracy, which suggests a high level of applicability.

Comparison of Drought Cases and Applicability Evaluation
Based on the water supply restriction in the Nakdong River Watershed, which includes the Yecheon Watershed, this study determined the number of days of water supply restriction in May and June 2015 ( Figure 13, Table 6). Then, the findings were used to compare the time series of the available number of water intake derived from the available number of water-intake days curves. The findings are as shown in Figure 14. All water supply restriction cases began in the May 31-June 18 period, when the available number of water-intake days in the Yecheon Watershed dropped below the dry season level. The findings indicate a high level of accuracy, which suggests a high level of applicability.    Giran -myeon  June 18-26  17  32  2  Jaesan-myeon  June 15-22  8  20  3  Jaesan-myeon  June 15-22  8  20  4 Myeongho-myeon June 14-26 3 10 5 Myeongho-myeon June 14-26 8 20 6 Myeongho-myeon June 11-22  11  18  7 Seokpo-myeon June 6-21 14 21 8 Jaesan-myeon May 30-June 22 18 40 Figure 14. Comparison of water supply available days and start date of water restriction events of Yecheon Basin (1 May to 1 July 2015).

Conclusions
The purpose of this study was to provide information for drought response and water resource planning aimed at effectively preventing the inconvenience and financial damage caused by water supply restrictions at local water sources on account of droughts and water shortage. Visits to local waterwork sites revealed that the personnel in the field are in need of the data on the number of days that water can be taken from the water source in case of no precipitation. Therefore, this study sought to develop a methodology for determining the available number of water-intake days quantitatively to ensure the local waterworks' resilience against droughts. The details and findings of this study are summarized below.
First, the SWAT model was used to calibrate and validate the daily flow hydrologic curve, which resulted in the determination coefficients for the total period ( ), NSE, and PBIAS at 0.70, 0.54, and 12.9, respectively, which fall under the "satisfactory" category based on the statistical performance assessment criteria proposed by Moriasi et al. [50]. The water balance of the study watershed was analyzed based on these findings, and the

Conclusions
The purpose of this study was to provide information for drought response and water resource planning aimed at effectively preventing the inconvenience and financial damage caused by water supply restrictions at local water sources on account of droughts and water shortage. Visits to local waterwork sites revealed that the personnel in the field are in need of the data on the number of days that water can be taken from the water source in case of no precipitation. Therefore, this study sought to develop a methodology for determining the available number of water-intake days quantitatively to ensure the local waterworks' resilience against droughts. The details and findings of this study are summarized below.
First, the SWAT model was used to calibrate and validate the daily flow hydrologic curve, which resulted in the determination coefficients for the total period (R 2 ), NSE, and PBIAS at 0.70, 0.54, and 12.9, respectively, which fall under the "satisfactory" category based on the statistical performance assessment criteria proposed by Moriasi et al. [50]. The water balance of the study watershed was analyzed based on these findings, and the average annual precipitation was 1255 mm, while the total flow was 697.4 mm or 55.6% of the precipitation. The evapotranspiration was 495.2 mm or 39.4% of the precipitation.
Second, based on the daily long-term flow at the MOE's local water source for the 2007-2017 period calibrated and validated using the SWAT model, the recession curve concept proposed by Kullman [39] was used to identify the periods of the recession events with 15 or more days of recession without precipitation. Based on the 43 precipitation events and the effect of seasonal weather proposed by Tallaksen [34], the available number of water intake days curves were calculated, accounting for the daily average water demands and drought flow rates throughout the total period and each season.
Third, this study applied the available number of water intake days curve equation developed for the local water source to actual cases of water supply restrictions in the study area and found that the starting dates of water supply restriction and the available numbers of water intake days matched the periods in which the water level dropped below the dry season level. In addition, this study proposed practical and straightforward concepts that can be used by local water source staff to respond to droughts.
In this paper, a practical and straightforward indicator for drought response in case of no rainfall used at the water intake station was developed, which was not be presented in previous studies. However, this study needs to be supplemented by further research on applying the methodology developed in this study to various areas and cases. Once optimized for different areas, the methodology is expected to be applicable to drought response and water resource planning to effectively prevent the inconvenience and financial damage caused by water supply restriction at local water sources on account of droughts and water shortage.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.