Next Article in Journal
Phytoplankton Time-Series in a LTER Site of the Adriatic Sea: Methodological Approach to Decipher Community Structure and Indicative Taxa
Previous Article in Journal
Tuning the Catalytic Water Oxidation Activity through Structural Modifications of High-Nuclearity Mn-oxo Clusters [Mn18M] (M = Sr2+, Mn2+)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

A Study to Suggest Monthly Baseflow Estimation Approach for the Long-Term Hydrologic Impact Analysis Models: A Case Study in South Korea

Department of Rural Construction Engineering, Kongju National University, 54 Daehak-ro, Yesan 32439, Korea
*
Author to whom correspondence should be addressed.
Water 2021, 13(15), 2043; https://doi.org/10.3390/w13152043
Submission received: 31 May 2021 / Revised: 22 June 2021 / Accepted: 26 July 2021 / Published: 27 July 2021

Abstract

:
Changes in both land use and rainfall patterns can lead to changes in the hydrologic behavior of the watershed. The long-term hydrologic impact analysis (L-THIA) model has been used to predict such changes and analyze the changes in mitigation scenarios. The model is simple as only a small amount of input data are required, but it can predict only the direct runoff and cannot determine the streamflow. This study, therefore, aimed to propose a method for predicting the monthly baseflow while maintaining the simplicity of the model. The monthly baseflows for 20 watersheds in South Korea were estimated under different land use conditions. Calibration of the monthly baseflow prediction method produced values for R2 and the Nash–Sutcliffe efficiency (NSE) within the ranges of 0.600–0.817 and 0.504–0.677, respectively; during validation, these values were in the ranges of 0.618–0.786 and 0.567–0.727, respectively. This indicates that the proposed method can reliably predict the monthly baseflow while maintaining the simplicity of the L-THIA model. The proposed model is expected to be applicable to all the various forms of the model.

1. Introduction

Changes in land use or rainfall pattern not only change the behavior of direct runoff and baseflow in the watershed, but also affect the occurrence of nonpoint source (NPS) pollution. Urbanization and industrialization have increased the ratio of impervious surfaces in watersheds. Therefore, it is necessary to analyze such changes in land use and investigate scenarios that may reduce the impacts of these changes. Hydrologic models are generally used for such analyses; Bieger et al. [1] assessed the impact of land use changes using the soil and water assessment tool (SWAT) [2,3], and three land use scenarios for forest, cropland, and orchard area changes were established. The result indicated that forest, cropland, and orchard area changes of −34.48%, +265.32%, and 204.51% led to surface flow increases of 46.1%. Guse et al. [4] analyzed the impact of spatially distributed five crop rotations using the SWAT model; nitrate loads were reduced with dynamic changes in agricultural crop rotations. Martin et al. [5] used the regional hydro-ecological simulation system model [6]; the study result indicated that high flows (highest 10 percentiles) increased by 37–88% and that low flows (lowest 10 percentiles) increased by 23–37% by land use changes. Additionally, Srivastava et al. [7] and Aghsaei et al. [8] reported that vegetation can provide a significant effect on hydrological components with the variations in the physical characteristics of the land surface, soil, and vegetation, which are the roughness, albedo, architectural resistance, infiltration capacity, leaf area index, root depth, and stomatal conductance.
The long-term hydrologic impact analysis (L-THIA) model has been in use for this purpose since 1994 [9]. The L-THIA model was first developed in the form of a spreadsheet in 1994 [9], followed by redevelopment so that it could be integrated with geographic information systems (GIS) [10,11]. However, the model is limited in terms of reflecting the various types of land use. To counter such limitations, Lim et al. [12] developed the L-THIA/NPS WWW model so that L-THIA could be used in the form of a spreadsheet in L-THIA/NPS GIS, which is based on the ArcView software. Bhaduri et al. [13] used the GIS-based L-THIA/NPS GIS model to predict changes in the direct runoff and pollution from nonpoint sources resulting from land use changes; they reported that the average annual direct runoff increased by 80%, while with pollution from lead, copper, and zinc increased by more than 50% as the urban and impervious area increased by 18% during the period 1973–1991. Wilson and Weng [14] analyzed changes in the direct runoff that were caused by land use changes using the L-THIA NPS model, which is an ArcHydro GIS extension. They reported a two-fold increase in the direct runoff in places where the residential land use increased by 37.3%, even though the precipitation increased by <30%. Liu et al. [15] considered best management practices (BMPs) such as the implementation of wet ponds, green roofs, and bioretention or the construction cost, annual maintenance cost, and interest rate of 12 techniques that belong to low impact development (LIDs) using the L-THIA-LID model [16]. They determined cost-effective strategies for 15 scenarios to reduce direct runoff, total nitrogen (TN), total phosphorus (TP), total suspended solids (TSS), Pb, biochemical oxygen demand (BOD), and chemical oxygen demand (COD). Eaton [17] analyzed a direct runoff reduction method using the L-THIA LID model to analyze green infrastructure screening and reported that direct runoff can be reduced by 12% by using bioretention and raingardens in the watershed. Li et al. [18] analyzed changes in the surface runoff caused by land use and rainfall changes using the ArcL-THIA 10.1 model that is based on ArcGIS 10.0 [19]. The authors calibrated the ArcL-THIA 10.1 model by using the baseflow filter program (BFLOW) [20] model to separate direct runoff; they reported that enhanced precipitation contributed more significantly to the observed changes compared to land use during the period 2005–2015.
Different methods for analyzing scenarios resulting from changes in land use or rainfall conditions, which also consider baseflow, have been proposed for the L-THIA model. Ahiablame et al. [21] analyzed the effects of rain barrel/cistern and porous pavement using the LTHIA-LID model in their study, which reflected the influence of baseflow. The annual baseflow regression equation, which includes the area of the watershed, annual precipitation, and baseflow index, was utilized for baseflow analysis in the L-THIA model. However, as this method can only predict the annual baseflow, it has limitations when used to determine monthly characteristics. Ryu et al. [22] improved the prediction process of the existing L-THIA in detail by using the asymptotic curve number (ACN) instead of CN for predicting the direct runoff. This made it possible to predict streamflow by adding modules for baseflow predictions and channel routing to overcome the most significant limitations of the existing L-THIA. However, this process involves an increase in the complexity of the model as three additional model parameters related to direct runoff prediction, four model parameters related to baseflow prediction, and three model parameters related to channel routing prediction were required, and the need for an optimization algorithm was suggested for calibrating the model.
The L-THIA model is based on spreadsheets [9,16] and GIS [10,11,12,13,14,19]. It has been in continuous use since its conception for analyzing direct runoff resulting from changes in land use or rainfall patterns as well as the effects of BMP and LID techniques. However, a comparison of the predictions made using this model with actual streamflow is required to improve the utility of this model. In addition, it is necessary to reflect the influence of baseflow on the analysis process based on this model. However, the model is generally used for the above purposes because of its simplicity, as it requires only input data to define CN and its computation process is not complicated. It is, therefore, necessary to maintain the current benefits of this model while predicting baseflow. This study aims to propose a method for predicting the baseflow while maintaining the simplicity of the current L-THIA model.

2. Materials and Methods

2.1. Decription of the Study Area

As the purpose of this study is the proposal of a method that can predict baseflow while maintaining the simplicity of the current L-THIA model, streamflow data describing the runoff in a watershed, the land use map, and monthly precipitation data were required as input into the L-THIA model.
The land use map (scale at 1:5000) was provided by the Environmental Geographic Information Service [23]; the monthly precipitation data were provided by the Korea Meteorological Administration [24]; the daily stream flow data were acquired from the Water Resources Management Information System (WAMIS) [25]. These data were used in this study. The watersheds selected for the study include the flow data measurement points operated by WAMIS as watershed outlets, and they are distributed such that the conditions in various regions of Korea are considered without any spatial overlap (Figure 1). A total of 20 watersheds were selected with areas ranging from 5694.8 to 155,805.9 ha (Table 1). The land use types were classified into urban, agriculture, forest, pasture, wetland, bare land, and water. Forest occupied the largest area in all the watersheds studied, followed by agriculture. The exception to this trend was Wsd-02 in which agriculture (2741.0 ha) represented the largest proportion (41.5%), followed by urban (2042.1 ha), which accounted for 30.9% of the watershed area (Table 1).
The analysis period was set to the five-year period from January 1st 2016 to December 31st 2020. The minimum monthly precipitation in each watershed ranged from 0.5 mm (Wsd-01) to 8.1 mm (Wsd-07), whereas the maximum monthly precipitation ranged from 454.5 mm (Wsd-17) to 822.0 mm (Wsd-14). The maximum precipitation was, therefore, 73 times (Wsd-07) to 1271 times (Wsd-01) higher than the minimum in each watershed, indicating a significant difference in monthly precipitation (Table 2).
The minimum daily stream flow ranged from 0.01 m3/s (Wsd-03, Wsd-04, Wsd-06, Wsd-07, Wsd-14, and Wsd-17) to 0.60 m3/s (Wsd-19), whereas the maximum daily stream flow ranged from 118.65 m3/s (Wsd-01) to 1160.00 m3/s (Wsd-20). The average daily stream flow ranged from 1.31 m3/s (Wsd-01) to 21.72 m3/s (Wsd-20) (Table 2). The difference between the minimum and maximum daily stream flow was significant in each watershed. Therefore, it was judged necessary to reflect monthly precipitation conditions in predicting the monthly baseflow.

2.2. Baseflow Separation

As the purpose of this study is to propose a method for predicting monthly baseflow, measurement of the actual baseflow is required for comparison. The measured streamflow is the sum of the surface runoff and the baseflow, which means that it is difficult to evaluate the monthly baseflow predictions method solely by comparison with the corresponding streamflow. Therefore, the monthly baseflow prediction method was examined by separating the baseflow from the streamflow in each watershed using the method suggested by Eckhardt (Equation (1)). In this baseflow separation method [26], the baseflow on day t (bt) is determined using the maximum value of the long-term ratio of baseflow to streamflow (BFImax), the filter parameter (a), the filtered baseflow at the t − 1 time step (bt−1), and the streamflow at t time step (Qt):
b t = 1 B F I m a x a b t 1 + 1 a B F I m a x Q t 1 a B F I m a x
After separating the baseflow (m3/s) from the streamflow (m3/s) in each watershed using Equation (1), the baseflow can be expressed (in m3). It reflects the area of each watershed as shown in Table 3. Overall, the minimum, maximum, and average values of the streamflow and baseflow showed a tendency to rise as the area of the watershed increased. The mean flow percentage (%), which is the ratio of baseflow to streamflow, ranged from 32.7% (Wsd-15) to 59.2% (Wsd-19), indicating significant differences among the mean flow percentages in each watershed despite the seemingly insignificant differences observed using the actual numbers. Wsd-15, with a relatively large watershed area, exhibited the minimum percentage of 32.7%; however, the mean flow percentage showed a tendency to slightly increase as the watershed area increased. This indicates that it is necessary to consider parameters related to the area under investigation when predicting the monthly baseflow.

2.3. Monthly Baseflow Estimation Approach

A remarkable method that can be used for baseflow prediction is the spreadsheet tool for the estimation of pollutant load (STEPL) [27] model. This model was proposed by the U. S. Environmental Protection Agency (U. S. EPA) to establish a total maximum daily load plan for the United States, and it can analyze the effects of more than 50 BMPs while also simulating the average annual runoff, sediment, nitrogen, phosphorus, and the BOD. In this model, prediction of the baseline is based on the annual precipitation, land use, and hydrologic soil group (HSG). For example, it predicts that the baseflow will correspond to 36, 24, 12, and 6% of the annual precipitation in an urban area when the HSG is A, B, C, and D, respectively. In this model, land uses other than urban include cropland, pasture, and forest, but the baseflow is estimated using the same method for all types of area. In other words, the model predicts that the baseflow will correspond to 45, 30, 15, and 7.5% of the annual precipitation for the HSGs A, B, C, and D, respectively. Since this method is a statistical approach with a simple prediction process, it was judged that the application of this method could also maintain model simplicity in the baseflow prediction process of the L-THIA model.
A basic formula is required to predict the monthly baseflow. Land use conditions need to be reflected because the baseflow is likely to vary alongside land use. It was also judged that monthly rainfall requires consideration because both the baseflow and streamflow are affected by precipitation. The basic formula for predicting the monthly baseflow is shown in Equation (2). For predicting the baseflow in month i, this equation considers the conditions under which precipitation occurs in month i, the coefficient for urban land (CURBN) and the area covered by urban land (AURBN), the coefficient for agriculture (CAGRL) and the area of agricultural land (AAGRL), the coefficient for forest (CFRST) and the area of forested land (AFRST), the coefficient for pasture (CPAST) and the area covered by pasture (APAST), the coefficient for wetland (CWTLD) and the area covered by wetland (AWTLD), the coefficient for bare land (CBARE) and the area covered by bare land (ABARE), and the coefficient for water (CWATR) and the area covered by water (AWATR).
B a s e f l o w i = P r e c i p i t a t i o n i × ( C U R B N × A U R B N + C A G R L × A A G R L + C F R S T × A F R S T + C P A S T × A P A S T + C W T L D × A W T L D + C B A R E × A B A R E + C W A T R × A W A T R )
Equation (2) can reflect the conditions surrounding both monthly precipitation and land use for predicting the baseflow, as the coefficients allow the degree of influence from each land use type to vary. However, the coefficients require definition, for which a genetic algorithm (GA) [28] was used in this case. GA, which is similar to the evolutionary process of nature, is used to obtain the optimized solution to a given problem. It is, therefore, useful in solving highly complex problems in the fields of business or engineering [29,30,31]. In this study, the selected watersheds were divided into two groups for determination of the seven coefficients. The first group included 11 watersheds. It consisted of the watersheds that were given even numbers: Wsd-01, Wsd-02, Wsd-04, Wsd-06, Wsd-08, Wsd-10, Wsd-12, Wsd-14, Wsd-16, Wsd-18, and Wsd-20, and it included the watersheds with the minimum and maximum areas. The second watershed group consisted of the remaining nine watersheds: Wsd-03, Wsd-05, Wsd-07, Wsd-09, Wsd-11, Wsd-13, Wsd-15, Wsd-17, and Wsd-19. After the coefficients of Equation (2) were defined through GA for the watersheds in the first group, Equation (2) was applied to the watersheds in the second group, along with a definition for all the coefficients, to examine the monthly baseflow prediction method. The processes described in Section 2.1–2.3 can be expressed as the following Figure 2.

3. Results and Discussion

3.1. Determination of Regression Model Coefficients

Optimal values for the seven coefficients in Equation (2) for predicting the monthly baseflow were determined by GA based on the monthly baseflow estimated for each watershed and the monthly flow that was distinguished using the Eckhardt filter equation (Equation (1)). The coefficients for which the smallest difference between the monthly baseflow separated using Equation (1) and the monthly baseflow estimated using Equation (2) was observed in each watershed were determined to be optimal. In this study, it was deemed necessary to measure the separated monthly baseflows and to develop the criteria for examining the validity of the estimated monthly baseflow. Duda et al. [32] reported that the estimated result is applicable when the R2 is higher than 0.65 and the difference is 45% or less. Skaggs et al. [33] mentioned that the result is applicable when the Nash–Sutcliffe efficiency (NSE) is higher than 0.50, and Wang et al. [34] suggested that it is applicable when the NSE is higher than 0.50, the R2 is higher than 0.60, and the PBIAS is ±15%. Moriasi et al. [35] deemed the result to be applicable when the R2 is higher than 0.60, NSE is higher than 0.50, and the PBIAS is less than 15%. In other words, there are various criteria for determining the applicability of the estimated result. The different criteria were summarized in this study, and the applicability of the monthly baseflow was estimated using a scatter plot with an NSE higher than 0.50 and an R2 higher than 0.60.
In terms of the optimized coefficients for the 11 watersheds in the first group, CURBN ranged from 0.00289 (Wsd-04) to 0.06314 (Wsd-04), CAGRL from 0.10076 (Wsd-01) to 0.49061 (Wsd-06), CFRST from 0.01256 (Wsd-02) to 0.28413 (Wsd-10), CPAST from 0.0.07913 (Wsd-06) to 0.31460 (Wsd-20), CWTLD from 0.05222 (Wsd-02) to 0.78599 (Wsd-16), CBARE from 0.03488 (Wsd-04) to 0.23603 (Wsd-16), and CWATR from 0.04919 (Wsd-04) to 0.38384 (Wsd-16) (Table 4). The optimized coefficients in each watershed were obtained when the difference between the separated and estimated baseflows describing baseflow was minimal. Since the purpose of this study is to propose a method for predicting the monthly baseflow in multiple watersheds rather than judging the accuracy of the monthly baseflow prediction for a specific watershed, the final values of the coefficients of Equation (2) were determined based on the average of each coefficient. The coefficients were, therefore, determined to be 0.04 for CURBN, 0.40 for CAGRL, 0.20 for CFRST, 0.18 for CPAST, 0.48 for CWTLD, 0.15 for CBARE, and 0.22 for CWATR.
In general, the contribution of land use to the baseflow can be considered to be related to the impervious surface ratio. The contribution of urban areas to the baseflow will be low because the impervious surface ratio is high. A CURBN of 0.04 was finally determined for urban areas based on the optimized results; this reflects the conditions of impermeability, as this coefficient is relatively lower than those obtained for other land uses. In contrast, the contribution of wetlands and reservoirs to the baseflow is high because of the constant infiltration of water. The value of 0.48 for CWTLD also appears to reflect this condition, as it is relatively high as compared to the coefficients for other types of land use.
It is noteworthy that the coefficient for agricultural land, CAGRL, had the second highest value after CWTLD. Agriculture in Korea is dominated by rice paddies, which are maintained in pond conditions during the rice cultivation period from May to October, resulting in a similar contribution to the baseflow as wetland. The contribution of agriculture to the baseflow should, therefore, be similar to that of wetland; this condition is sufficiently reflected in the coefficient for agriculture.
The monthly baseflow in the first watershed group was calculated again using Equation (2) by applying the finally determined coefficients, and the suitability of the monthly baseflow was determined by the values of R2, NSE, and the scatter plot. Both R2 and NSE showed applicable ranges with R2 ranging from 0.600 (Wsd-06) to 0.817 (Wsd-16) and NSE from 0.504 (Wsd-01) to 0.677 (Wsd-18 and Wsd-20) (Table 5). Based on the scatter plot, the estimated monthly baseflow tended to be slightly lower than the separated monthly baseflow for Wsd-06, Wsd-08, and Wsd-10; Wsd-20 showed a scattered tendency that was comparable to the other watersheds. Overall, however, there were no significant differences in the tendencies or values obtained via prediction and the separated monthly baseflow (Figure 3).
In the plot of flow duration curves, the estimated monthly baseflow did not capture the separated monthly baseflow in the dry-conditions (flow duration intervals from 60% to 90%) and the low-flow (flow duration intervals from 90% to 100%) regimes often; however, it does in the other flow regimes, which are the high-flow (flow duration intervals from 0–10%), the moist-conditions (flow duration intervals from 10% to 40%), and the mid-range flow (flow duration intervals from 40% to 60%) regimes (Figure 4).
Therefore, the application of optimized coefficients to each watershed will render the estimated monthly baseflow similar to the separated monthly baseflow. Thus, the application of the finally determined coefficients is expected to result in satisfactory predictions.

3.2. Validation of Regression Model Coefficients

When the model parameters that were determined for use for calibration were adjusted for the associated watersheds during calibration, the results derived for these watersheds could be considered satisfactory. However, it is necessary to apply these model parameters to watersheds other than those used in the calibration process to examine whether the model parameters were well calibrated and determine whether the estimated results are applicable. Therefore, in this study, the values determined for the coefficients that correspond to the model parameters in the calibration process were applied to the second group to determine whether the estimated monthly baseflow is also applicable in this group.
The estimated monthly baseflow was determined to be applicable because the R2 ranged from 0.618 (Wsd-09) to 0.786 (Wsd-17) and the NSE from 0.567 (Wsd-07) to 0.727 (Wsd-05) (Table 6). Based on the scatter plot, the estimated monthly baseflow tended to be slightly lower than the separated monthly baseflow for Wsd-03, Wsd-11, and Wsd-19. Overall, however, no significant differences were observed in either tendency or value for the estimated and separated monthly baseflow (Figure 5).
Similar to the results of calibration, the estimated monthly baseflow did not capture the separated monthly baseflow in the dry-conditions and the low-flow regimes often in the flow duration curve plots; however, it does in the other flow regimes, which are the high-flow, the moist-conditions, and the mid-range flow regimes (Figure 6). Based on flow duration curves in both calibration and validations processes, the estimated monthly baseflow fit to the separated monthly baseflow reasonably in the high-flow and the moist-conditions; however, it did not especially in the low-flow regime. This means that the monthly baseflow approach will be reasonable in the applications with the issues regarding the high flow and the moist conditions such as flooding or nonpoint source pollution analysis. However, caution needs to be exercised when the approach is used for any applications regarding low flow such as water supply simulations in drought.

4. Conclusions

This study aimed to propose a method for predicting the baseflow in order to improve the shortcomings of the current long-term hydrologic impact analysis (L-THIA) model, which can only predict surface runoff. To achieve this aim, the measured streamflow was separated into surface runoff and baseflow to determine whether the L-THIA model can predict the hydrologic behavior in a given watershed through comparison with the measured flow. The model is limited in that it cannot determine the influence of the baseflow because only the surface runoff can be estimated. It has been used until now because it requires only land use maps, soil maps, and precipitation data and its computation process is simple from a user point of view. Therefore, it is necessary to improve the current L-THIA model so that it can predict the baseflow while maintaining its simplicity. As such, a method to improve the limitations of the L-THIA model was proposed by applying the baseflow prediction method based on land use and annual precipitation, which is used in the spreadsheet tool for the estimation of pollutant load (STEPL) model.
To this end, 20 independent watersheds in South Korea were selected and divided into two groups. The first group of 11 watersheds was used to propose the monthly baseflow prediction method, and the second group of nine watersheds was used to examine the proposed method. Since the proposed method uses monthly precipitation data along with the areas and coefficients associated with seven different land use types, it does not exceed the range of data used in the current L-THIA model to predict surface runoff. The simplicity of the model can, therefore, be maintained. In addition, the coefficients represent the contribution of precipitation to the baseflow for each land use type. In the case of urban areas, which are generally considered to have a low contribution to the baseflow due to their high impervious surface ratio, the coefficient was determined to be lower than that of other land uses. Wetland is considered to have a high contribution to the baseflow because of the continuous infiltration that is associated with water storage. Therefore, it has a higher coefficient than the other land use types. In addition, ponded paddy fields, which are characteristic of the agriculture in Korea, are assumed to contribute to the baseflow in a similar manner as wetland during the cultivation period from May to October. Therefore, the coefficient for agricultural land was determined to be similar to the that for wetland. In other words, it appears that the proposed monthly baseflow prediction method can sufficiently reflect the conditions of the different land uses while reflecting seasonal conditions as it uses monthly precipitation.
A proposed step in use of the approach is to separate direct runoff and baseflow from measured streamflow first, to calibrate direct runoff estimated with curve numbers, and to calibrate baseflow estimated by the suggested method, in turn. Additionally, since the proposed method does not reflect the hydrologic behavior beneath the surface and can be seen as a statistical or probabilistic approach based on the characteristics of the watersheds in Korea, caution needs to be exercised when it is applied to watersheds without ponded paddy fields. Attention also needs to be paid to the application of this monthly baseflow prediction regression equation when the watershed area used in the derivation process is exceeded. In addition, the suggested regression was developed on a monthly basis; therefore, it will be limited to apply for daily or weekly baseflow estimations.

Author Contributions

Conceptualization, Y.-S.P.; methodology, H.-S.C.; investigation, M.-S.C.; data curation, M.-S.C.; writing—original draft preparation, H.L.; writing—review and editing, Y.-S.P.; visualization, H.-S.C.; supervision, Y.-S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This subject is supported by Korea Ministry of Environment as “The SS (Surface Soil conservation and management) projects; 2019 (2019002820001)”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented and used in the study are available on the request from the corresponding author.

Acknowledgments

The authors send special thanks to the Korea Environmental Industry & Technology Institute for continuous research project support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bieger, K.; Hörmann, G.; Fohrer, N. The impact of land use change in the Xiangxi Catchment (China) on water balance and sediment transport. Reg. Environ. Chang. 2013, 15, 485–498. [Google Scholar] [CrossRef]
  2. Arnold, J.; Srinivasan, R.; Muttiah, R.; Williams, J. Large area hydrologic modeling and assessment part I: Model development. J. Am. Water. Resour. Assoc. 1998, 34, 73–89. [Google Scholar] [CrossRef]
  3. Arnold, J.; Fohrer, N. SWAT2000: Current capabilities and research opportunities in applied watershed modelling. Hydrol. Process. 2005, 19, 563–572. [Google Scholar] [CrossRef]
  4. Guse, B.; Pfannerstill, M.; Fohrer, N. Dynamic modelling of land use change impacts on nitrate loads in rivers. Environ. Process. 2015, 2, 575–592. [Google Scholar] [CrossRef]
  5. Martin, K.L.; Hwang, T.; Vose, J.M.; Coulston, J.W.; Wear, D.N.; Miles, B.; Band, L.E. Watershed impacts of climate and land use changes depend on magnitude and land use context. Ecohydrology 2017, 10, e1870. [Google Scholar] [CrossRef]
  6. Tague, C.L.; Band, L.E. RHESSys: Regional Hydro-Ecologic Simulation System-An object-oriented approach to spatially distributed modeling of carbon, water, and nutrient cycling. Earth Interact. 2004, 8, 1–42. [Google Scholar] [CrossRef]
  7. Srivastava, A.; Kumari, N.; Maza, M. Hydrological response to agricultural land use heterogeneity using variable infiltration capacity model. Water Resour. Manag. 2020, 34, 3779–3794. [Google Scholar] [CrossRef]
  8. Aghsaei, H.; Dinan, N.M.; Moridi, A.; Asadolahi, Z.; Delavar, M.; Fohrer, N.; Wagner, P.D. Effect of dynamic land use/land cover change on water resources and sediment yield in the Anzali wetland catchment, Gilan, Iran. Sci. Total Environ. 2020, 712, 136449. [Google Scholar] [CrossRef] [PubMed]
  9. Harbor, J. A practical method for estimating the impact of land-use change on surface runoff, groundwater recharge, and wetland hydrology. J. Am. Plann. Assoc. 1994, 60, 95–108. [Google Scholar] [CrossRef]
  10. Bhaduri, B. A GIS-Based Model to Assess the Long-Term Impacts of Land Use Change on Hydrology and Nonpoint Source Pollution. Ph.D. Thesis, Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, IN, USA, 1998. [Google Scholar]
  11. Grove, M. Development and Application of a GIS-Based Model for Assessing the Long-Term Hydrologic Impacts of Land Use Change. Master’s Thesis, Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, IN, USA, 1997. [Google Scholar]
  12. Lim, K.J.; Engel, B.; Kim, Y.; Bhaduri, B.; Harbor, J. Development of the Long Term Hydrologic Impact Assessment (L-THIA) WWW systems. In Proceedings of the 10th International Soil Conservation Organization Meeting, Purdue University and USDA-ARS National Soil Erosion Research Laboratory, West Lafayette, IN, USA, 24–29 May 1999. [Google Scholar]
  13. Bhaduri, B.; Harbor, J.; Engel, B.; Grove, M. Assessing watershed-scale, long-term hydrologic impacts of land-use change using a GIS-NPS model. Environ. Manag. 2000, 26, 643–658. [Google Scholar] [CrossRef] [PubMed]
  14. Wilson, C.; Weng, Q. Assessing surface water quality and its relation with urban land cover changes in the Lake Calumet area, Greater Chicago. Environ. Manag. 2010, 45, 1096–1111. [Google Scholar] [CrossRef]
  15. Liu, Y.; Bralts, V.F.; Engel, B.A. Evaluating the effectiveness of management practices on hydrology and water quality at watershed scale with a rainfall-runoff model. Sci. Total Environ. 2015, 511, 298–309. [Google Scholar] [CrossRef] [PubMed]
  16. Ahiablame, L.M.; Engel, B.A.; Chaubey, I. Representation and evaluation of low impact development practices with L-THIA-LID: An example for site planning. Environ. Pollut. 2012, 1. [Google Scholar] [CrossRef]
  17. Eaton, T.T. Approach and case-study of green infrastructure screening analysis for urban stormwater control. J. Environ. Manag. 2018, 209, 495–504. [Google Scholar] [CrossRef] [PubMed]
  18. Li, F.; Chen, J.; Liu, Y.; Xu, P.; Sun, H.; Engel, B.A.; Wang, S. Assessment of the impact of land use/cover change and rainfall change on surface runoff in China. Sustainability 2019, 11, 3535. [Google Scholar] [CrossRef] [Green Version]
  19. ArcL-THIA Park, Y.S.; Lim, K.J.; Larry, T.; Engel, B.A.; Long-Term Hydrologic Impact Assessment. Department of Agricultural and Biological Engineering. 2013. Available online: http://npslab.kongju.ac.kr/models/LTHIA2013_Manual.pdf (accessed on 13 August 2018).
  20. Arnold, J.G.; Allen, P.M. Automated methods for estimating baseflow and ground water recharge from streamflow records. J. Am. Water Resour. Assoc. 1999, 35, 411–424. [Google Scholar] [CrossRef]
  21. Ahiablame, L.; Chaubey, I.; Engel, B.; Cherkauer, K.; Merwade, V. Estimation of annual baseflow at ungauged sites in Indiana USA. J. Hydrol. 2013, 476, 13–27. [Google Scholar] [CrossRef]
  22. Ryu, J.; Jang, W.S.; Kim, J.; Choi, J.D.; Engel, B.A.; Yang, J.E.; Lim, K.J. Development of a watershed-scale long-term hydrologic impact assessment model with the asymptotic curve number regression equation. Water 2016, 8, 153. [Google Scholar] [CrossRef] [Green Version]
  23. Environmental Geographic Information Service. Available online: https://egis.me.go.kr/main.do (accessed on 5 January 2021).
  24. Korea Meteorological Administration. Available online: https://www.weather.go.kr/w/index.do (accessed on 7 January 2021).
  25. Water Resources Management Information System. Available online: http://water.nier.go.kr/main/mainContent.do (accessed on 7 January 2021).
  26. Eckhardt, K. How to construct recursive digital filters for baseflow separation. Hydrol. Process. 2005, 19, 507–515. [Google Scholar] [CrossRef]
  27. Tetra Tech, Inc. User’s Guide Spreadsheet Tool for the Estimation of Pollutant Load (STEPL) Version 4.1; Tetra Tech, Inc.: Fairfax, VA, USA, 2011. [Google Scholar]
  28. Holland, J.H. Adaptation in Natural and Artificial Systems; University of Michigan Press: Ann Arbor, MI, USA, 1975. [Google Scholar]
  29. Lim, K.J.; Park, Y.S.; Kim, J.; Shin, Y.C.; Kim, N.W.; Kim, S.J.; Jeon, J.H.; Engel, B.A. Development of genetic algorithm-based optimization module in WHAT system for hydrograph analysis and model application. Comput. Geosci. 2010, 36, 936–944. [Google Scholar] [CrossRef]
  30. Toğan, V.; Daloğlu, T.A. An improved genetic algorithm with initial population strategy and self-adaptive member grouping. Comput. Struct. 2008, 86, 1204–1218. [Google Scholar] [CrossRef]
  31. Park, Y.S.; Engel, B.A.; Harbor, J. A web-based model to estimate the impact of best management practices. Water 2014, 6, 455–471. [Google Scholar] [CrossRef] [Green Version]
  32. Duda, P.B.; Hummel, P.R.; Donigian, A.S.; Imhoff, J.C. BASINS/HSPF: Model use, calibration, and validation. Trans. ASABE 2012, 55, 1523–1547. [Google Scholar] [CrossRef]
  33. Skaggs, R.W.; Youssef, M.A.; Chescheir, G.M. DRAINMOD: Model use, calibration, and validation. Trans. ASABE 2012, 55, 1509–1522. [Google Scholar] [CrossRef]
  34. Wang, X.; Williams, J.; Gassman, P.; Baffaut, C.; Izaurralde, R.; Jeong, J.; Kiniry, J. EPIC and APEX: Model use, calibration, and validation. Trans. ASABE 2012, 55, 1447–1462. [Google Scholar] [CrossRef]
  35. Moriasi, D.N.; Gitau, M.W.; Pai, N.; Daggupati, P. Hydrologic and water quality models: Performance measures and evaluation criteria. Am. Soc. Agric. Biol. Eng. 2015, 58, 1763–1785. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of study watershed.
Figure 1. Location of study watershed.
Water 13 02043 g001
Figure 2. Schematic flow of the study.
Figure 2. Schematic flow of the study.
Water 13 02043 g002
Figure 3. Scatter plot of separated monthly baseflow and estimated monthly baseflow: (a) Wsd-01; (b) Wsd-02; (c) Wsd-04; (d) Wsd-06; (e) Wsd-08; (f) Wsd-10; (g) Wsd-12; (h) Wsd-14; (i) Wsd-16; (j) Wsd-18; (k) Wsd-20.
Figure 3. Scatter plot of separated monthly baseflow and estimated monthly baseflow: (a) Wsd-01; (b) Wsd-02; (c) Wsd-04; (d) Wsd-06; (e) Wsd-08; (f) Wsd-10; (g) Wsd-12; (h) Wsd-14; (i) Wsd-16; (j) Wsd-18; (k) Wsd-20.
Water 13 02043 g003aWater 13 02043 g003bWater 13 02043 g003c
Figure 4. Flow duration curves of separated monthly baseflow and estimated monthly baseflow in calibration: (a) Wsd-01; (b) Wsd-02; (c) Wsd-04; (d) Wsd-06; (e) Wsd-08; (f) Wsd-10; (g) Wsd-12; (h) Wsd-14; (i) Wsd-16; (j) Wsd-18; (k) Wsd-20.
Figure 4. Flow duration curves of separated monthly baseflow and estimated monthly baseflow in calibration: (a) Wsd-01; (b) Wsd-02; (c) Wsd-04; (d) Wsd-06; (e) Wsd-08; (f) Wsd-10; (g) Wsd-12; (h) Wsd-14; (i) Wsd-16; (j) Wsd-18; (k) Wsd-20.
Water 13 02043 g004aWater 13 02043 g004b
Figure 5. Scatter plot of separated monthly baseflow and estimated monthly baseflow: (a) Wsd-03; (b) Wsd-05; (c) Wsd-07; (d) Wsd-09; (e) Wsd-11; (f) Wsd-13; (g) Wsd-15; (h) Wsd-17; (i) Wsd-19.
Figure 5. Scatter plot of separated monthly baseflow and estimated monthly baseflow: (a) Wsd-03; (b) Wsd-05; (c) Wsd-07; (d) Wsd-09; (e) Wsd-11; (f) Wsd-13; (g) Wsd-15; (h) Wsd-17; (i) Wsd-19.
Water 13 02043 g005aWater 13 02043 g005bWater 13 02043 g005c
Figure 6. Flow duration curves of separated monthly baseflow and estimated monthly baseflow in validation: (a) Wsd-03; (b) Wsd-05; (c) Wsd-07; (d) Wsd-09; (e) Wsd-11; (f) Wsd-13; (g) Wsd-15; (h) Wsd-17; (i) Wsd-19.
Figure 6. Flow duration curves of separated monthly baseflow and estimated monthly baseflow in validation: (a) Wsd-03; (b) Wsd-05; (c) Wsd-07; (d) Wsd-09; (e) Wsd-11; (f) Wsd-13; (g) Wsd-15; (h) Wsd-17; (i) Wsd-19.
Water 13 02043 g006aWater 13 02043 g006b
Table 1. Land uses in the study watersheds.
Table 1. Land uses in the study watersheds.
WatershedArea (ha)
UrbanAgricultureForestPastureWetlandBare landWaterTotal
Wsd-01590.3372.24637.727.126.141.30.15694.8
Wsd-022042.12741.01447.9199.036.1121.712.26600.0
Wsd-03196.31169.79256.982.79.1121.89.410,845.8
Wsd-04569.71536.99425.8142.763.4126.747.811,912.9
Wsd-05656.4479.511,134.4214.90.0142.82.312,630.2
Wsd-06435.73361.29002.8224.64.9155.436.013,220.6
Wsd-07131.41979.711,115.2281.013.172.2108.513,701.0
Wsd-08983.12266.711,494.5565.435.9149.045.515,540.0
Wsd-0920.21472.915,794.592.60.018.00.017,398.2
Wsd-10843.65017.511,771.9269.2112.4424.5131.618,570.7
Wsd-11165.22317.516,250.0239.00.355.876.119,103.8
Wsd-12343.85118.214,822.7200.685.8152.5111.820,835.4
Wsd-13770.95908.726,938.4934.1288.9289.4301.835,432.1
Wsd-14404.87297.326,403.61158.60.1384.5179.535,828.3
Wsd-15278.24779.935,251.6327.5120.3228.2342.641,328.3
Wsd-162454.45312.632,486.21044.5186.9655.1618.942,758.7
Wsd-171347.98703.731,043.61322.7259.1320.0380.643,377.7
Wsd-182098.621,737.323,612.52642.949.71319.6480.951,941.4
Wsd-19347.65614.773,684.2601.70.1371.194.080,713.3
Wsd-201763.317,163.2130,735.81857.5274.1644.03368.2155,805.9
Table 2. Statistics of monthly precipitation and daily stream flow.
Table 2. Statistics of monthly precipitation and daily stream flow.
WatershedMonthly Precipitation (mm)Daily Stream Flow (m3/s)
min.max.min.max.Mean
Wsd-010.5635.50.02118.651.31
Wsd-022.1738.10.16341.681.98
Wsd-034.1628.00.01241.933.20
Wsd-042.1738.10.01191.062.71
Wsd-051.5494.40.07133.193.18
Wsd-063.5822.00.011496.454.33
Wsd-078.1587.50.01247.784.03
Wsd-082.0469.50.13380.584.07
Wsd-091.0761.50.10189.603.13
Wsd-102.1738.10.13995.015.38
Wsd-117.6631.40.02618.545.49
Wsd-120.5492.60.05202.724.06
Wsd-130.9588.10.05312.416.78
Wsd-143.5822.00.01661.897.53
Wsd-150.8712.00.10894.3010.44
Wsd-165.7498.30.30697.6010.28
Wsd-171.8451.50.01702.767.73
Wsd-181.3616.70.131039.879.39
Wsd-190.7606.50.60736.2017.99
Wsd-204.1628.00.501160.0021.72
Table 3. Statistics describing the streamflow and separated baseflow.
Table 3. Statistics describing the streamflow and separated baseflow.
WatershedStreamflow (×106 m3)Baseflow (×106 m3)Mean Flow Percentage (%)
min.max.Meanmin.max.Mean
Wsd-010.13628.6003.5440.1048.8621.55243.8
Wsd-021.04544.8235.2110.76311.0332.24043.0
Wsd-030.00369.7178.1600.00220.8803.52043.1
Wsd-040.18557.8767.0110.03720.7243.01042.9
Wsd-050.71041.9938.3600.42815.6454.02648.2
Wsd-060.219214.74811.3750.17859.2784.75241.8
Wsd-070.53566.04710.2900.31424.5684.37842.6
Wsd-080.74398.04210.4320.49127.2054.25440.8
Wsd-090.16476.4218.1660.04626.3493.79246.4
Wsd-102.035172.18914.2311.06561.1506.88848.4
Wsd-110.328147.04114.4500.23654.0886.18742.8
Wsd-121.45775.85511.1201.16838.1425.95653.6
Wsd-131.456110.26518.2220.63261.6989.31251.1
Wsd-140.902151.18319.4630.17866.07810.20752.4
Wsd-150.130199.93026.5420.05278.2208.68132.7
Wsd-163.707117.33727.0341.97244.48413.59550.3
Wsd-170.21295.16819.1080.04133.3977.19437.6
Wsd-180.481280.72524.9980.421123.58413.50554.0
Wsd-197.379375.97848.0524.087194.16828.45659.2
Wsd-208.027524.90657.8265.825201.18530.41452.6
Table 4. Optimized coefficients for the Equation (2).
Table 4. Optimized coefficients for the Equation (2).
WatershedCURBNCAGRLCFRSTCPASTCWTLDCBARECWATR
Wsd-010.022410.100760.253730.092740.528900.092650.04994
Wsd-020.039810.480800.012560.184110.052220.219550.24430
Wsd-040.063140.194600.212890.240170.702070.034880.04919
Wsd-060.031370.490610.274210.079130.240280.191440.14365
Wsd-080.024370.359660.266960.082600.174130.093010.18265
Wsd-100.058510.424200.284130.231520.476020.218250.31505
Wsd-120.060500.232360.259960.153330.674170.116000.20427
Wsd-140.034000.390350.141210.282330.708070.175280.29423
Wsd-160.002890.335030.242790.121780.785990.236030.38384
Wsd-180.034940.342990.158100.171550.272990.232780.20131
Wsd-200.049690.478280.123890.314600.615440.041410.30796
Min.0.002890.100760.012560.079130.052220.034880.04919
Max.0.063140.490610.284130.314600.785990.236030.38384
Mean0.038330.348150.202770.177620.475480.150120.21603
Final value0.040.400.200.180.480.150.22
Table 5. R2 and NSE of separated and estimated monthly baseflow in the first watershed group.
Table 5. R2 and NSE of separated and estimated monthly baseflow in the first watershed group.
WatershedR2NSE
Wsd-010.7360.504
Wsd-020.6070.563
Wsd-040.6340.507
Wsd-060.6000.565
Wsd-080.7040.679
Wsd-100.7080.536
Wsd-120.6970.547
Wsd-140.6880.583
Wsd-160.8170.623
Wsd-180.7850.677
Wsd-200.6910.677
Table 6. R2 and NSE of separated and estimated monthly baseflow in the second watershed group.
Table 6. R2 and NSE of separated and estimated monthly baseflow in the second watershed group.
WatershedR2NSE
Wsd-030.6300.569
Wsd-050.7390.727
Wsd-070.6640.567
Wsd-090.6180.606
Wsd-110.6580.647
Wsd-130.6890.632
Wsd-150.7280.670
Wsd-170.7860.674
Wsd-190.6260.571
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, H.; Choi, H.-S.; Chae, M.-S.; Park, Y.-S. A Study to Suggest Monthly Baseflow Estimation Approach for the Long-Term Hydrologic Impact Analysis Models: A Case Study in South Korea. Water 2021, 13, 2043. https://doi.org/10.3390/w13152043

AMA Style

Lee H, Choi H-S, Chae M-S, Park Y-S. A Study to Suggest Monthly Baseflow Estimation Approach for the Long-Term Hydrologic Impact Analysis Models: A Case Study in South Korea. Water. 2021; 13(15):2043. https://doi.org/10.3390/w13152043

Chicago/Turabian Style

Lee, Hanyong, Hyun-Seok Choi, Min-Suh Chae, and Youn-Shik Park. 2021. "A Study to Suggest Monthly Baseflow Estimation Approach for the Long-Term Hydrologic Impact Analysis Models: A Case Study in South Korea" Water 13, no. 15: 2043. https://doi.org/10.3390/w13152043

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop