Next Article in Journal
A Numerical Validation of 3D Experimental Dam-Break Wave Interaction with a Sharp Obstacle Using DualSPHysics
Next Article in Special Issue
Automatic Calibration for CE-QUAL-W2 Model Using Improved Global-Best Harmony Search Algorithm
Previous Article in Journal
Current Challenges and Advancements on the Management of Water Retreatment in Different Production Operations of Shale Reservoirs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Groundwater Abstraction on Hydrological Responses during Extreme Drought Periods in the Boryeong Dam Catchment, Korea

1
Department of Civil and Environmental Engineering, University of Science & Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 34113, Korea
2
Department of Land, Water and Environment Research, Korea Institute of Civil Engineering & Building Technology, 283 Goyang-daero, Ilsanseo-gu, Goyang-si 10223, Korea
*
Author to whom correspondence should be addressed.
Water 2021, 13(15), 2132; https://doi.org/10.3390/w13152132
Submission received: 22 June 2021 / Revised: 26 July 2021 / Accepted: 27 July 2021 / Published: 2 August 2021
(This article belongs to the Special Issue Advances and Challenges in Hydrological Modeling and Engineering)

Abstract

:
Groundwater withdrawal results in a significant depletion of groundwater storage due to the frequency and intensity of droughts and increasing irrigation demands. To ensure the sustainable use of groundwater resources, it is necessary to accurately simulate the groundwater behavior of catchments using a surface–groundwater integrated runoff model. Most of the existing catchment runoff models have been applied to surface water management; thus, integrated runoff analysis studies that consider the interaction between surface water and groundwater are required. Due to the intensive agricultural sector in Korea and the position of rice as the staple in the Korean diet, more than 50% of groundwater abstraction is used for irrigation. Therefore, it is very important to understand the hydrological interrelationships between agricultural areas and the entire watershed. This study aimed to compare and analyze the groundwater levels in the mountainous areas and paddy field areas in the Boryeong Dam catchment through a surface–groundwater integrated runoff simulation using the Catchment Hydrologic Cycle Assessment Tool model, and to compare the hydrological responses in wet years (2010–2012) and dry years (2014–2016). The maximum difference in the monthly groundwater level in the dry years compared to the wet years was 1.07 m at the forest catchment and 0.37 m at the paddy catchment. These results indicate that the impact of drought on the groundwater level of paddy catchments is not significant compared to the forest catchments; however, drought slows the recovery of the groundwater level before the rainy season, thereby limiting the agricultural groundwater use in the catchment.

1. Introduction

Climate change has been shown to be a principal factor changing the hydrological cycle on a global scale [1]. The drought cycle pattern in South Korea has been shortened recently, and droughts of varying intensity have occurred every year since 2013 [2]. The Boryeong Dam, which serves as a major water source in the middle-western region of Korea, recorded its lowest storage volume in 2015 due to the continuous shortage of precipitation since 2014. The dam levels were stabilized by limiting domestic water consumption in the region and constructing a water supply tunnel of 21.9 km, connecting the Geum River to the Boryeong Dam catchment [3]. The Ministry of Agriculture, Food and Rural Affairs of Korea [4] has classified the drought warning steps according to the storage of the dam into four stages: the concern stage (less than 70% storage in a normal year), warning stage (less than 60%), caution stage (less than 50%), and serious stage (less than 40%). The water storage volume of the Boryeong Dam in 2015 decreased to 18.87%, reaching the serious stage due to the extreme drought that occurred in the central region of Korea [5]. The drought in the central region continued, and the water storage rate reached 8.29% in 2017—the lowest water storage level since the construction of the Boryeong Dam [5]. Jung et al. [6] analyzed the drought index using unstructured data of the climate, streamflow, and rainfall of Boryeong city in the dry years of 2014–2016. Yu et al. [7] analyzed the drought transitions in meteorological and agricultural droughts, including the 2015 drought, and Kim et al. [8] assessed the water shortage of the Boryeong Dam by applying future climate change scenarios.
As the frequency and intensity of droughts have increased, the importance of groundwater as a sustainable water resource has been recognized [9]. According to the Third Comprehensive Long-Term Water Resources Plan (2016) of the Ministry of Land, Infrastructure, and Transport of Korea, the groundwater level is decreasing every year due to groundwater consumption and urbanization. Groundwater use, which was 2.57 billion m3 in 1994, increased by ~60% to 4.1 billion m3 in 2014, indicating a very high dependence on groundwater use in Korea [10]. Groundwater is an important water resource that dominates streamflow, especially in the dry season, and its variability directly contributes to the health of a stream [11]. Brunner et al. [12] analyzed MODFLOW simulations related to surface water–groundwater interactions. Mukherjee et al. [13] showed that the depletion of streams is related to the reduction in baseflow caused by groundwater storage depletion in adjacent aquifers.
Groundwater storage is not only affected by meteorological factors but can also be significantly affected by groundwater extraction for agricultural purposes [14], and groundwater withdrawals for agricultural purposes have a significant impact on the groundwater system of the watershed itself [15]. In addition, the paddy field area accounts for approximately 8.21% of the Korean territory [16], and about 51% of the total amount of groundwater consumption is used only for agriculture [17]. Therefore, the impact of groundwater withdrawals for agricultural purposes on the hydrological cycle of the catchment must be closely examined to accurately estimate the impact of drought on groundwater levels in paddy field areas in Korea. Since the aquifer is dynamically connected to the nearby sub-catchments, it is necessary to quantitatively analyze not only the groundwater movement in paddy fields but also the hydrological interactive responses between the paddy field areas and the whole watershed. Leng et al. [18] performed and analyzed Community Land Model 4.0 (CLM4) simulations driven by downscaled/bias-corrected historical simulations and future projections from general circulation models (GCMs) to investigate the effects of irrigation on global water resources. Ferguson and Maxwell [19] evaluated and compared the impacts of large-scale climate change and local-scale water management practices in an agricultural watershed using a fully integrated model of groundwater, surface water, and land surface processes. Sorooshian et al. [20] investigated the impacts of agricultural irrigation on hydrological processes in California using a regional climate model (RCM), an offline land surface model output, and available in situ observations and remote sensing data. Kannan [21] developed a modeling approach for canal irrigation systems and irrigation best management practices (BMPs) to simulate the hydrological balance of irrigated watersheds based on the water requirements of crops, the number and frequency of irrigation systems, and critical crop water requirement stages. Ozdogan et al. [22] described an original technique for applying satellite-derived irrigation data within a land surface model (LSM) and evaluated the effects on irrigation in the United States. Essaid et al. [23] examined changes in streamflow, groundwater discharge to the stream, and stream temperature resulting from irrigation practices using a watershed-scale surface water and groundwater flow integrated model.
The Boryeong Dam catchment is representative of most Korean catchments, with a steep and mountainous catchment including paddy field areas. Even though this catchment is not a regional-scale watershed, with an area of only 163 km2, it has complex interactive hydrological cycle processes due to human activities such as groundwater abstraction, agricultural reservoirs, and export and import water supply systems in the catchment. However, most studies on hydrological runoff modeling in Korea have focused on the surface runoff analysis and have not suggested human activities including groundwater pumping for irrigation demand and water intake systems in the watershed.
In this study, a runoff simulation of the Boryeong Dam catchment was conducted using a physical parameter-based and surface–groundwater integrated hydrological model for paddy fields in Korea, including a paddy module that considers the ponding depth and levee height of the irrigation area of the watershed. The groundwater levels of the paddy and forest areas of the watershed were analyzed and compared with two consecutive periods of wet years (2010–2013) and drought years (2014–2016).

2. Materials and Methods

2.1. Study Area Description

The study area is the Boryeong Dam catchment located in Chungnam Province (the middle-western region of South Korea) (Figure 1), and the area of the catchment is 163 km2, with an average slope of 40.19% and an average elevation of 232.56 m. The average annual rainfall of the watershed is 1.244 mm, and the average monthly temperature of the catchment during the study period was −4.5 to 27.3 °C, indicating that the seasonal temperature fluctuation was very large. The Boryeong Dam, located in the Ungcheon stream of the Geum River, is a multi-purpose dam that was constructed in 2000. The Boryeong Dam has a water storage capacity of 116.9 million m3 and supplies seven adjacent cities and counties [24].

2.2. Data Collection

The runoff simulation was conducted for the years 2000 to 2019, and 20 years of daily rainfall, streamflow, and meteorological data were collected to establish the Catchment Hydrologic Cycle Assessment Tool (CAT) model. Rainfall data were collected from the Boryeong Dam station, Dohwadam station, and Samsan station operated by the Korea Water Resources Corporation (K-water) [5], and the meteorological data used for estimating the evapotranspiration using the Penman–Monteith method, such as minimum and maximum temperature (°C), humidity (%), sunshine hours, and wind speed (m/s), were collected from the Boryeong Meteorological Station operated by the Korea Meteorological Administration [25].
Soil parameters required for model simulation were extracted from a 1:25,000 scale soil map of the Rural Development Administration of Korea [26], and groundwater aquifer parameters were extracted from the Boryeong Region Groundwater Research Report 2007 [5,10]. The land cover map of the Ministry of Environment [27] was used for classification into forest, urban, and paddy fields, and the input data for the model simulation were derived by dividing the watershed into pervious and impervious areas.

2.3. CAT Model

The CAT is a physical parameter-based distributed hydrological model that allows the quantitative evaluation of the long- and short-term water cycle of catchments [28]. It is a link–node connecting model designed to simulate runoff, including infiltration, evapotranspiration, and groundwater flow for each spatial unit, and it divides the hydrological cycle process into pervious and impervious areas. The major physical parameters required for the initial simulation are slope, soil type, land use, aquifer, and river information. The concept and development of the basic model were based on the unconfined aquifer and single soil layer assumptions. The model categorizes the incoming rainfall that falls on pervious, impervious, and paddy field areas, providing the surface flow, infiltration, or evapotranspiration (Figure 2).
The CAT model has been applied in various studies for different catchments. Jang et al. [29,30] assessed the long-term hydrological behavior of agricultural reservoirs in the Idong catchment and analyzed the impact of future climate change on the hydrological responses and runoff in the Gyeongan–Cheon River Basin. Birhanu et al. [31] assessed the results of five hydrological models, including the CAT, applied in 10 catchments in Korea. Choi et al. [32] analyzed a short-term CAT runoff simulation and carried out a sensitivity analysis of soil parameters using three infiltration methods. Lee and Cho [33] analyzed the hydrological cycle of four catchments in Ulsan city, Korea, using the CAT, and Miller et al. [34] investigated the changes in stormwater runoff caused by the transformation of rural landscapes into peri-urban areas in the UK using the CAT.
To simulate the runoff process in the paddy field (Figure 3), the soil and groundwater layers were divided in the same way as the pervious area of the watershed. Artificial drainage facilities can be included in the soil layer to accommodate underground culvert drainage and pipe drainage; however, only surface drainage by levee height was considered in this study according to the field situation of the catchment. Surface drainage occurs when the ponding depth is greater than the height of the surface drain water threshold. The surface drain equations are
Q s =   α H s H p ( H s > H p )
Q s = 0 ( H s H p )
where Q s is the discharge from the surface (m3/s), α is the drainage coefficient of the surface drain levee in the paddy (mm0.5/h), H s is the ponding depth of the paddy (m), and H p is the height of the surface drain levee of the paddy (m).
The hydrological interrelation between surface water and groundwater was calculated according to Darcy’s law, with flows based on the hydraulic conductivities, river stage, and groundwater level. The interrelation between the river level and the groundwater level was calculated using Equations (3) and (4). Equation (4) was applied when the river level was higher than the groundwater level in the vicinity; otherwise, Equation (3) was applied.
Q r = K s r A r
Q r = K s r ( h H r b r ) A r
Q g = K s r h x · l · T
Q i n Q o u t = A · S d h d t
where Q r is the inflow into the river or recharge from the river (m3/s), Q g is the groundwater flow (m3/s), Q i n and Q o u t are the inflow and outflow of the aquifer (m3/s), respectively, K s r is the saturated hydraulic conductivity of the riverbed (m/s), A r is the area of the riverbed (m2), b r is the riverbed thickness (m), h is the groundwater level (m), H r is the riverbed elevation (m), h / x is the slope of the groundwater level,   l is the connected length between catchments (m), T is the average aquifer thickness (m), A is the catchment area (m2), S is the storage coefficient, and d h / d t is the rate of level change.

2.4. Model Calibration by Infiltration Methods and Comparison Criteria

The Boryeong Dam catchment was assumed to be a single catchment to calibrate and compare the model performances by three infiltration methods. Since the type of infiltration method applied to the model could significantly affect the results of the runoff simulation, it was necessary to simulate the runoff model for each infiltration method in advance and select the most suitable infiltration method for the study area. The CAT provides three infiltration analysis methods: the rainfall excess method, the Horton method, and the Green–Ampt method, with each infiltration method having unique parameters to be calibrated. Rainfall excess is that part of rainfall which is not lost to infiltration, interception, and depression storage. In the case of the rainfall excess infiltration method, it is assumed that infiltration occurs unconditionally until the groundwater zone is saturated using unsaturated vertical and horizontal hydraulic conductivities, and the saturated soil moisture increases since it is assumed that the excess precipitation is surface runoff. The Green–Ampt infiltration method is calculated using the soil parameters, and the Horton infiltration method is calculated using the infiltration capacity curve equation as a function of time using empirical values. In the case of the rainfall excess method, the simulated runoff hydrograph tends to be the closest to the observed hydrograph, and the simulated peak flow tends to be the most underestimated among the three infiltration methods. In the case of the Green–Ampt method, the runoff ratio increases sharply according to the occurrence of rainfall, and it starts to decrease when the rainfall ends. The Green–Ampt infiltration method is based on the physical soil parameters; however, the simulated results tend to be highly sensitive to rainfall characteristics.
The SCEUA-P global optimization algorithm [35], which is a module of PEST (Model-Independent Parameter Estimation) and a parameter optimization and uncertainty analysis package [36] associated with the CAT model, was applied for parameter calibration while considering its applications, robustness, and efficiency, and the sensitive parameters for runoff simulation were selected as the calibration parameters.
In the rainfall excess infiltration method, the saturated soil moisture (θs), saturated vertical hydraulic conductivity (ks), and saturated horizontal hydraulic conductivity (ksi) were selected for parameter calibration. The saturated soil moisture (θs), saturated hydraulic conductivity (ks), saturated horizontal hydraulic conductivity (ksi), and wetting front soil suction head (PSI) were calibrated using the Green–Ampt infiltration method, and the maximum and minimum infiltration capacities (fo, fc) were calibrated using the Horton infiltration method. The saturated hydraulic conductivity of the riverbed (ku_riv) and the storage coefficient of the aquifer (aqf_S) were manually calibrated for each of the three infiltration methods (Table 1).
The model was calibrated for the period 2000–2009 and validated for the period 2010–2019. The Nash–Sutcliffe efficiency (NSE) index is commonly used as an indicator of hydrological model performance; however, its over-sensitivity to extreme values has been recognized [37,38]. The multi-objective functions used for evaluating model performances in this study were the Kling–Gupta efficiency (KGE), the determination coefficient (R2), NSE, and the root mean square error (RMSE) [37]:
R 2 = 1 ( Q o b s Q s i m ) 2 ( Q o b s Q o b s ¯ ) 2
N S E = 1 ( Q o b s Q s i m ) 2 ( Q o b s Q m e a n ) 2
K G E = 1 ( R 2 [ Q o b s , Q s i m ] 1 ) 2 + ( S D [ Q s i m ] S D [ Q o b s ] 1 ) 2 + ( M [ Q s i m ] M [ Q o b s ] 1 ) 2
M S E = ( Q o b s Q s i m ) 2 n
where Q o b s , Q s i m , and Q m e a n are the observed, simulated, and mean observed streamflow, respectively, and S D , M , and n are the standard deviation, mean, and number of data points, respectively.

2.5. CAT Model Setup for 10 Sub-Catchments

In this study, the Boryeong Dam catchment was divided into five forest catchments and five paddy field catchments to analyze the groundwater hydrological response in each sub-catchment and the interrelationship between the paddy fields and the watershed. The watershed was divided into three regions (R1–R3) based on the daily groundwater pumping data [5]. The regions were then divided into 10 sub-catchments to analyze the groundwater level in detail. The watershed was assumed to consist of only paddy fields and forest areas. The sub-catchments were divided by the area ratio of paddy fields and forests of three regions according to the land use map [27]. The R1 (northwestern) region and the R3 (southern) region of the catchment were divided into one forest area and one paddy field area, while the R3 (northeastern) region of the watershed was divided into three forest catchments and three paddy catchments considering the amount of water intake from the agricultural reservoirs to the paddy fields (Table 2).
In addition, we considered additional potential factors that could affect the hydrological cycle of the watershed. As shown in Figure 4, we considered the water export to the Boryeong downtown through the Gaehwa weir, which is located outside the catchment in the northwestern region of the catchment. The Samsan and Hwasung agricultural reservoirs in the northeastern region were considered as the water supplied from these reservoirs directly affects the paddy fields in the catchment. In addition, a water supply from the Geum River into the watershed was also considered.
Groundwater movement connections between adjacent sub-catchments were considered by linking the upslope forest catchment and downslope paddy field catchment (Figure 5). The groundwater movement connection parameters were the aquifer slope (%), the average distance between sub-catchments (m), the length of the boundary line (m), and the hydraulic conductivity of the aquifer (mm/s). The F 1 catchment was connected to the P 1 catchment, the F 2, F 2-1, and F 2-2 catchments were connected to the P 1, P 2, and P 2-2 catchments, respectively, and the F 3 catchment was connected to the P 3 catchment.

2.6. Human Impact on Hydrological Cycle in Boryeong Dam Catchment

2.6.1. Water Supply Tunnel from Geum River

Due to the severe drought that occurred in 2015, the middle-western part of Korea was severely affected by water shortages. The water storage rate of the Boryeong Dam continued to decrease despite the implementation of restrictions on water use in the catchment. The average inflows of the Boryeong Dam in 2014 and 2015 were 2.32 m3/s and 1.93 m3/s, respectively, and remained at about 50% of the average inflow (4.33 m3/s) for 18 years from 1998 to 2015. In response, the Korean government constructed an emergency water supply tunnel facility that connects the lower reaches of the Geum River and the upper Boryeong Dam, which has been operational since February 2016. The length of the water supply tunnel is approximately 21 km, and it can supply up to 115,000 m3/day of water to the Boryeong Dam (Figure 6). The Geum River water supply was calculated using the daily inflow time series data provided by the Korea Water Resources Corporation [5] using the module for external water supply in the CAT.

2.6.2. Groundwater Pumping Rate

The groundwater use data for the 10-year period from 2008 to 2017 were provided by the National Groundwater Information Center [16], and the average daily groundwater use data were applied to the study period. The daily groundwater use for agricultural, industrial, and domestic water use was classified, and the volume of groundwater for each purpose was divided by the area ratio in each sub-catchment. Agricultural groundwater was applied to the paddy fields, and domestic and industrial groundwater use was applied to the forest areas. In the case of the R1 region, the agricultural groundwater use rate was 96% of the total groundwater consumption over 2000 to 2019, and the rates were 93.8% and 44.3% for the R2 and R3 regions, respectively. Table 3 shows the daily groundwater use per sub-catchment according to the area ratio. The average daily volume of groundwater pumping was applied as a constant value in the model. Groundwater pumping is only applied in the period from April to September since irrigation begins in early April and continues until the end of September in Korea.

2.6.3. Reservoirs in the Catchment

In the Boryeong Dam catchment, there is a Gaewha diversion weir structure in the northwestern area of the catchment, which supplies water to the Boryeong downtown, located out of the catchment. There are two reservoirs upstream of the northeastern area of the catchment that are used to supply agricultural water to nearby paddy fields: the Samsan reservoir supplies agricultural water to P 2-1, and the Hwasung reservoir supplies water to P 2-2. The intake volume of the Gaewha diversion weir is 5240 m3/day, and this was applied as a constant value. The agricultural water supply data from the Samsan and Hwasung reservoirs were collected as a time series from 2012 to 2018, and the intake water supply for the paddy fields was considered only during the cultivation period from April 1 to September 31. Table 4 shows the average daily intake rate of the reservoirs.

2.7. Wet and Dry Periods of the Catchment

To compare and analyze the hydrological responses and groundwater level changes in the wet period and the drought period in the catchment, consecutive wet and dry years were selected based on the average annual precipitation and the average annual runoff ratio within the validated period of 2010–2019. The 30-year average annual precipitation in the Boryeong Dam catchment is approximately 1244 mm, and the 20-year average annual runoff ratio is approximately 53.4%. Three consecutive years of 2010, 2011, and 2012 were selected as the wet period, with annual precipitation rates of 1467 mm, 1971 mm, and 1588 mm, respectively, and with runoff ratios of 64.8%, 75.3%, and 54.3 %, respectively, indicating higher than average precipitation and runoff ratios. The three consecutive years of 2014, 2015, and 2016 were selected as the dry period, with annual precipitation rates of 1089 mm, 1023 mm, and 1086 mm, respectively, and with runoff ratios of 40.4%, 36.2%, and 47.4%, respectively, indicating lower than average precipitation and runoff ratios (Figure 7).

3. Results

3.1. Model Calibration for Single Watershed Simulation Using Three Different Infiltration Methods

The parameter values before and after calibration for the rainfall excess method, Green–Ampt method, and Horton infiltration method are shown in Table 5. The θs, ks, and ksi parameters were calibrated using the rainfall excess method. The θs, ks, ksi, and PSI parameters were calibrated using the Green–Ampt method, and the fo and fc parameters were calibrated using the Horton method.
Table 6 shows the statistical results for the observed and simulated discharge before and after parameter optimization for each infiltration analysis method, assuming the Boryeong Dam catchment as a single catchment. In all cases, before and after parameter calibration, the NSE results, when applying the rainfall excess infiltration analysis method, show relatively high accuracy. After parameter calibration, the NSE results are 0.802 for rainfall excess, 0.625 for Green–Ampt, and 0.762 for the Horton infiltration analysis method, and the results for KGE are 0.827, 0.775, and 0.832, respectively. In the case of KGE, the results of the Horton infiltration method are similar or show slightly higher accuracy than the results of rainfall excess; however, the model performance when applying the rainfall excess infiltration method was considered to be satisfactory considering the R2 and RMSE results together. Therefore, we assumed the rainfall excess infiltration method to be the most suitable infiltration method for analyzing the groundwater level of the watershed in the detailed sub-catchment model setup.
A comparison of the observed and simulated streamflow for the calibrated and validated periods on a logarithmic scale is shown in Figure 8a. The R2, NSE, and KGE values were 0.803, 0.802, and 0.827, respectively, during the calibration period of 2000–2009, and the values of R2, NSE, and KGE were 0.781, 0.722, and 0.825, respectively, during the validation period of 2010–2019 (Figure 8b). The statistical results of the calibration period show a slightly higher accuracy for the model simulation compared to the results of the validation period; however, both the calibration and validation periods had satisfactory model performances.

3.2. Comparison of Streamflow in Wet and Dry Periods

We simulated the 10 sub-catchment models and compared the streamflow results against the observed streamflow at the outlet of the watershed. The annual comparisons of the observed and simulated streamflow at the outlet of the watershed in each year of the wet period are shown in Figure 9a, and those of the dry period are shown in Figure 9b. In the case of the wet period, the annual observed and simulated streamflows were 1491 m3/s and 1825 m3/s in 2010, 2435 m3/s and 2812 m3/s in 2011, and 1597 m3/s and 1634 m3/s in 2012, respectively. In the dry period, the annual observed and simulated streamflows were 754 m3/s and 834 m3/s in 2014, 617 m3/s and 701 m3/s in 2015, and 835 m3/s and 976 m3/s in 2016, respectively. The simulated streamflow discharge tended to be overestimated in the wet period compared to the observed streamflow; however, it was slightly underestimated during the dry period.
The comparisons of the observed and simulated streamflow for the entire wet period and the entire dry period are shown as a logarithmic scale in Figure 10a,b, respectively. The statistical results for the wet period show values of R2, NSE, KGE, and RMSE of 0.74, 0.73, 0.80, and 10.25, respectively, and the statistical results for the dry period show values of R2, NSE, KGE, and RMSE of 0.76, 0.76, 0.80, and 2.54, respectively. The results of both the wet and dry period simulations show satisfactory model performances, and the simulation of the dry period showed a slightly better model performance compared to the wet period.

3.3. Groundwater Level of Wet and Dry Periods

To analyze the differences in the groundwater level between the drought and wet periods, it is necessary to closely compare the monthly groundwater level fluctuations in each period. Table 7 shows the difference between the monthly average groundwater level for the three years of the wet period and the dry period in forest sub-catchments, and Table 8 shows the difference between the monthly average groundwater levels for the three years of the wet and dry periods in the paddy sub-catchments. The differences indicate that the groundwater level decreases in the forest sub-catchments were larger than the groundwater level decreases in the paddy sub-catchments. In the case of the forest sub-catchments, the decrease in the annual average groundwater level for each sub-catchment shows that the minimum value was 0.01 m in the F 2-1 catchment, and the maximum value was 0.38 m in the F 3 catchment. In the case of the paddy sub-catchments, the minimum value was 0.09 m in the P 2 and P 2-1 catchments, and the maximum value was 0.12 m in the P 3 catchment. The monthly groundwater level decrease in the dry period compared to the wet period in the forest sub-catchments shows that the maximum difference was 1.07 m in the F 3 catchment in September, and the monthly groundwater level decrease in the paddy sub-catchments shows that the maximum difference was 0.37 m in the P 3 catchment in September, indicating that the maximum annual decrease in the groundwater level in the forest sub-catchments was about 2.9 times greater than that of the paddy sub-catchments.
The difference in the groundwater level of the paddy sub-catchments in the comparison of the wet and dry periods was less than the difference in the groundwater level of the forest sub-catchments due to the reduced agricultural groundwater use during the dry period (Table 9). In the case of P 1, there were no significant differences in the use of agricultural groundwater between the wet and dry periods as the area of the catchment is smaller than that of the other paddy sub-catchments. However, in the case of the paddy sub-catchments P 2, P 2-1, P 2-2, and P 3, groundwater use was greater during the wet period since the paddy sub-catchment areas are larger, indicating that agricultural groundwater use sharply decreased during the dry period (Table 9). In the case of P 2, the annual average groundwater use during the wet and dry periods was 14,773 m3/day and 3164 m3/day, respectively, indicating that only 21.4% of the groundwater was used in the dry period compared to the wet period. In the case of P 2-1 and P 2-2, groundwater use in the dry period was approximately 21.4% of that of the wet period. In the case of P 3, the annual average groundwater use during the wet and dry periods was 22,636 m3/day and 2634 m3/day, respectively, indicating that only 11.6% of the groundwater was used in the dry period compared to the wet period. The results show that the available groundwater for agriculture is greatly affected by the drought season, and the efficient management of groundwater is very important, especially in groundwater-fed irrigation catchments.
In Korea, approximately 30% of the annual precipitation is concentrated over a month, from late June to late July, during the rainy season [25], and irrigation begins in early April and continues until the end of September. Therefore, to analyze the fluctuations in the groundwater level, it is necessary to closely examine the fluctuations in the groundwater level during the rainy seasons but also the irrigation period when agricultural groundwater is used. The monthly groundwater level fluctuations during the three years of the wet and dry periods are shown in Figure 11 for the forest sub-catchments and in Figure 12 for the paddy sub-catchments. In the case of the forest sub-catchments in the wet period, the groundwater level gradually decreased from January to June, then started to rise from the end of June—the rainy season in Korea—and reached the highest groundwater level in July. During the dry period in forest sub-catchments, the resulting graphs show a marked difference, showing a lower groundwater level than that of the wet period; however, there was no significant difference in the monthly groundwater level because of the lack of baseflow due to insufficient precipitation in the rainy season. In the case of the paddy catchments in the wet period, the groundwater level started to decrease from April, since the irrigation period is from April to September in Korea. The groundwater level then recovered rapidly, rising sharply in August in response to the precipitation of the rainy season starting from the end of June. The changes in the groundwater level in response to precipitation were slower and had a longer time lag than those of the forest catchments due to the agricultural groundwater pumping in the paddy catchments. During the dry period in the paddy catchments, the groundwater level started to decrease from April and could not recover quickly in July, as it did in the wet period, due to insufficient precipitation during the rainy season. The groundwater level gradually recovered from September after the irrigation period ended.

4. Conclusions and Summary

The impact of drought on groundwater levels over the Boryeong Dam catchment was analyzed in the present study. We conducted a single watershed runoff simulation for 2000–2019 using three infiltration methods to calibrate the parameters and determine the most suitable infiltration method for the study area. Runoff simulations with three infiltration methods were calibrated for 2000–2009 and validated for 2010–2019 against the observed streamflow. The results of the single watershed runoff simulation show that the rainfall excess infiltration method showed the most satisfactory model performance compared to the Horton and Green–Ampt infiltration methods, with values of 0.802 for NSE and 0.827 for KGE in the calibrated period, and 0.722 for NSE and 0.825 for KGE in the validated period. Then, the Boryeong Dam catchment model with the rainfall excess infiltration method was divided into 10 sub-catchments to analyze the hydrological cycle processes in paddy field and forest sub-catchments and the groundwater level changes in drought periods. The wet period of 2010–2012 and the dry period of 2014–2016 were selected based on the average annual precipitation and the average annual runoff ratio. In addition, the model considered the human impacts on the hydrological cycle over the study area, such as reservoirs, the agricultural groundwater pumping rate, and the import or export of water in the catchment.
The annual comparisons of the observed and simulated streamflow for each year showed that streamflow was overestimated by the model in the wet period and slightly underestimated in the dry period. In the case of the wet period, the statistical results of R2, NSE, KGE, and RMSE were 0.74, 0.73, 0.80, and 10.25, respectively, while in the case of the dry period, these values were 0.76, 0.76, 0.80, and 2.54, respectively, indicating that both model performances for the wet and dry periods showed a satisfactory model performance.
The differences in the groundwater level in the forest catchments were larger than those in the paddy field catchments, as groundwater use in the dry period in the paddy sub-catchments was 11.6% compared to 21.4% in the wet period.
During the wet period in the forest sub-catchments, the groundwater level continued to decrease from January to May and started to increase from June—the rainy season—showing the highest groundwater level in July. During the dry period, the groundwater levels did not significantly change in the forest catchments. In the case of the paddy field catchments, the groundwater level affected by the irrigation period recovered rapidly, showing a sharp increase in August in response to the rainy season during the wet period; however, it recovered slowly during the dry period due to the lack of precipitation in the rainy season. The agricultural groundwater use was limited during the drought period to 11.6–21.4% compared to the normal period, and the recovery of the groundwater level in the paddy field catchments took a longer time.
Groundwater is the main source of fresh water required for socio-economic development in the region for both domestic and agricultural use. The sustainable use of agricultural groundwater depends on understanding the dynamic behavior and responses of groundwater in paddy catchments and their relationship to other hydrological components. Since groundwater and surface water flows have strong interrelationships, they need to be considered and managed in combination. The use of integrated surface–groundwater modeling analysis is necessary to model groundwater availability under present and future climate conditions and to ensure the sustainability of groundwater use.
The results of this study reveal that agricultural groundwater availability is dependent on the prevailing climatic conditions. It is necessary to secure sustainable groundwater resources through efficient water management at the catchment scale and to compare the hydrological cycle processes between the mountainous catchments and paddy field catchments in Korea. Subsequent studies are needed to accurately analyze the hydrological cycle process of paddy field catchments by analyzing the relationship between precipitation, groundwater levels, baseflow, and groundwater recharge.

Author Contributions

Conceptualization, S.P., H.K.; Methodology, S.P., H.K. and C.J.; Formal analysis, S.P.; Investigation, C.J.; Resources, S.P., C.J.; Data curation, S.P., H.K.; Writing—original draft preparation, S.P.; Writing—review and editing, H.K.; Supervision, H.K., C.J.; Project administration, H.K.; Funding acquisition, H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Korea Institute of Civil Engineering and Building Technology, grant number 20210194-001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This study did not include any publicly available datasets.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, X.; Tang, Q.; Zhang, X.; Lettenmaier, D.P. Runoff sensitivity to global mean temperature change in the CMIP5 models. J. Geophys. Res. Lett. 2014, 41, 5492–5498. [Google Scholar] [CrossRef]
  2. National Drought Information Portal. Available online: https://www.drought.go.kr (accessed on 6 May 2021).
  3. Park, B.S.; Yoon, H.J.; Hong, Y.S.; Kim, S.P. A study on the optimal operation and policy of the Boryeong Dam diversion pipe line using the SWAT model. J. Korean Soc. 2020, 36, 546–558. [Google Scholar] [CrossRef]
  4. Ministry of Agriculture, Food and Rural Affairs of Korea. Available online: https://www.mafra.go.kr (accessed on 6 May 2021).
  5. Korea Water Resources Corporation. Available online: https://www.water.or.kr (accessed on 6 May 2021).
  6. Jung, J.H.; Park, D.H.; Ahn, J.H. Drought evaluation using unstructured data: A case study for Boryeong area. J. Korea Water Resour. Assoc. 2020, 53, 1203–1210. [Google Scholar] [CrossRef]
  7. Yu, M.S.; Cho, Y.H.; Kim, T.W.; Chae, H.S. Analysis of drought propagation using hydrometeorological data: From meteorological drought to agricultural drought. J. Korea Water Resour. Assoc. 2018, 51, 195–205. [Google Scholar] [CrossRef]
  8. Kim, W.J.; Jung, C.G.; Kim, J.U.; Kim, S.J. Water shortage assessment by applying future climate change for Boryeong Dam using SWAT. J. Korea Water Resour. Assoc. 2018, 51, 1195–1205. [Google Scholar] [CrossRef]
  9. Lee, J.M.; Park, J.H.; Chung, E.; Woo, N.C. Assessment of Groundwater Drought in the Mangyeong River Basin, Korea. Sustainability 2018, 10, 831. [Google Scholar] [CrossRef] [Green Version]
  10. Ministry of Land, Infrastructure and Transport. Available online: http://www.molit.go.kr (accessed on 6 May 2021).
  11. Lee, D.J.; Lee, J.M.; Kim, J.G.; Lim, K.J.; Engel, B.A.; Yang, J.E.; Jung, Y.H. Effects of slope magnitude and length on SWAT baseflow estimation. J. Irrg. Drain Eng. 2019, 145, 04018037. [Google Scholar] [CrossRef]
  12. Brunner, P.; Simmons, C.T.; Cook, P.G.; Therrien, R. Modeling surface water-groundwater interaction with MODFLOW: Some considerations. Groundwater 2020, 48, 174–180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mukherjee, A.; Bhanja, S.N.; Wada, Y. Groundwater depletion causing reduction of baseflow triggering Ganges river summer drying. Sci. Rep. 2018, 8, 12049. [Google Scholar] [CrossRef]
  14. Laveti, N.V.S.; Banerjee, A.; Kartha, S.A.; Dutta, S. Anthropogenic influence on monthly groundwater utilization in an irrigation dominated Ganga river sub-basin. J. Hydrol. 2021, 593, 125800. [Google Scholar] [CrossRef]
  15. McCallum, A.M.; Anderson, M.S.; Kelly, B.F.J.; Giambastiani, B.; Acworth, R.I. Hydrological investigations of surface water-groundwater interactions in a sub-catchment in the Namoi Valley, NSW, Australia. IAHS Publ. 2009, 329, 157–166. [Google Scholar]
  16. Statistics Korea. Available online: http://kostat.go.kr (accessed on 6 May 2021).
  17. National Groundwater Information Center. Available online: http://www.gims.go.kr (accessed on 6 May 2021).
  18. Leng, G.; Huang, M.; Tang, Q.; Leung, R. A modeling study of irrigation effects on global surface water and groundwater resources under a changing climate. J. Adv. Model. Earth Syst. 2015, 7, 1285–1304. [Google Scholar] [CrossRef]
  19. Ferguson, I.M.; Maxwell, R.M. Human impacts on terrestrial hydrology: Climate change versus pumping and irrigation. Environ. Res. Lett. 2012, 7, 044022. Available online: https://iopscience.iop.org/article/10.1088/1748-9326/7/4/044022 (accessed on 6 May 2021). [CrossRef] [Green Version]
  20. Sorooshian, S.; AghaKouchak, A.; Li, J. Influence of irrigation on land hydrological processes over California. J. Geophys. Res. 2014, 119, 137–152. [Google Scholar] [CrossRef] [Green Version]
  21. Kannan, N.; Jeong, J.; Srinivasan, R. Hydrologic modeling of a canal-irrigated agricultural watershed with irrigation best management practices: Case study. J. Hydrol. Eng. 2011, 16, 746–757. [Google Scholar] [CrossRef] [Green Version]
  22. Ozdogan, M.; Rodell, M.; Beaudoing, H.K.; Toll, D.L. Simulating the effects of irrigation over the United States in a Land Surface Model based on satellite-derived agricultural data. J. Hydromeorol. 2010, 11, 171–184. [Google Scholar] [CrossRef]
  23. Essaid, H.I.; Caldwell, R.R. Evaluating the impact of irrigation on surface water-groundwater interaction and stream temperature in an agricultural watershed. Sci. Total Environ. 2017, 599–600, 581–596. [Google Scholar] [CrossRef]
  24. Rim, C.S. Hydrologic simulation considering watershed characteristics in Boryeong-Dam watershed. J. Civ. Eng. 2016, 10, 67–72. [Google Scholar]
  25. Korea Meteorological Administration. Available online: http://hydro.kma.go.kr (accessed on 6 May 2021).
  26. Rural Development Administration. Available online: http://www.rda.go.kr (accessed on 6 May 2021).
  27. Ministry of Environment. Available online: http://me.go.kr (accessed on 6 May 2021).
  28. Kim, H.J.; Jang, C.H.; Noh, S.J. Development and application of the catchment hydrologic cycle assessment tool considering urbanization (Ⅰ)—Model development-. J. Korea Water Resour. Assoc. 2012, 45, 203–215. [Google Scholar] [CrossRef] [Green Version]
  29. Jang, C.H.; Kim, H.J.; Kim, J.T. Prediction of reservoir water level using CAT. J. Korean Soc. Agric. Eng. 2012, 54, 27–38. [Google Scholar] [CrossRef] [Green Version]
  30. Jang, C.H.; Kim, H.J.; Ahn, S.R.; Kim, S.J. Assessment of hydrological changes in a river basin as affected by climate change and water management practices, by using the CAT model. Irrig. And Drain. 2016, 65, 26–35. [Google Scholar] [CrossRef]
  31. Birhanu, D.; Kim, H.J.; Jang, C.H.; Park, S.H. Does the complexity of evapotranspiration and hydrological models enhance robustness? Sustainability 2018, 10, 2837. [Google Scholar] [CrossRef] [Green Version]
  32. Choi, S.W.; Jang, C.H.; Kim, H.J. Analysis of short-term runoff characteristics of CAT-PEST connected model using different infiltration analysis methods. J. Korea Acad. Industr. Coop. Soc. 2016, 17, 26–41. [Google Scholar] [CrossRef]
  33. Lee, S.H.; Cho, H.J. Analysis of water cycle at main streams in Ulsan using CAT model. J. Korea Water Resour. Assoc. 2014, 47, 1–10. [Google Scholar] [CrossRef]
  34. Miller, J.D.; Kim, H.J.; Kjeldsen, T.R.; Packman, J.; Grebby, S.; Dearden, R. Assessing the impact of urbanization on storm runoff in a peri-catchment using historical change in impervious cover. J. Hydrol. 2014, 515, 59–70. [Google Scholar] [CrossRef] [Green Version]
  35. Duan, Q.; Sorooshian, S.; Gupta, V.K. Optimal use of the SCE-UA global optimization method for calibrating watershed models. J. Hydrol. 1994, 158, 265–284. [Google Scholar] [CrossRef]
  36. Roy, T.; Gupta, H.V.; Serrat-Capdevila, A.; Valdes, J.B. Using satellite-based evapotranspiration estimates to improve the structure of a simple conceptual rainfall-runoff model. Hydrol.Earth Syst. Sci. 2017, 21, 879–896. [Google Scholar] [CrossRef] [Green Version]
  37. Moussa, R.; Chahinian, N. Comparison of different multi-objective calibration criteria using a conceptual rainfall-runoff model of flood events. Hydrol.Earth Syst. Sci. 2009, 13, 519–535. [Google Scholar] [CrossRef] [Green Version]
  38. Beven, K.; Binley, A. Glue: 20 years on. Hydrol. Process. 2014, 28, 5897–5918. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of the Boryeong Dam catchment including the rainfall stations, reservoirs, water level, and meteorological station.
Figure 1. Location of the Boryeong Dam catchment including the rainfall stations, reservoirs, water level, and meteorological station.
Water 13 02132 g001
Figure 2. Schematic diagram of the CAT model (CAT 3.0 user’s manual, 2017).
Figure 2. Schematic diagram of the CAT model (CAT 3.0 user’s manual, 2017).
Water 13 02132 g002
Figure 3. Hydrological cycle in the paddy field in the CAT model.
Figure 3. Hydrological cycle in the paddy field in the CAT model.
Water 13 02132 g003
Figure 4. CAT model setup of the Boryeong Dam catchment including 10 sub-catchments and water export and import.
Figure 4. CAT model setup of the Boryeong Dam catchment including 10 sub-catchments and water export and import.
Water 13 02132 g004
Figure 5. Groundwater movement connections between paddy field and forest catchments.
Figure 5. Groundwater movement connections between paddy field and forest catchments.
Water 13 02132 g005
Figure 6. Water supply tunnel from the lower reaches of the Geum River to the upper Boryeong Dam catchment.
Figure 6. Water supply tunnel from the lower reaches of the Geum River to the upper Boryeong Dam catchment.
Water 13 02132 g006
Figure 7. Comparison of annual precipitation, discharge, and runoff ratio in the wet and dry periods.
Figure 7. Comparison of annual precipitation, discharge, and runoff ratio in the wet and dry periods.
Water 13 02132 g007
Figure 8. (a) Comparison of observed, calibrated (2000–2009), and validated (2010–2019) streamflow of the rainfall excess infiltration method-applied model. (b) Scattered plot of calibration and validation against observed and simulated streamflow.
Figure 8. (a) Comparison of observed, calibrated (2000–2009), and validated (2010–2019) streamflow of the rainfall excess infiltration method-applied model. (b) Scattered plot of calibration and validation against observed and simulated streamflow.
Water 13 02132 g008
Figure 9. (a) The comparison of observed and simulated streamflow and precipitation for each year in the wet period. (b) The comparison of observed and simulated streamflow and precipitation for each year in the dry period.
Figure 9. (a) The comparison of observed and simulated streamflow and precipitation for each year in the wet period. (b) The comparison of observed and simulated streamflow and precipitation for each year in the dry period.
Water 13 02132 g009
Figure 10. (a) The observed and simulated streamflow plotted on a logarithmic scale for the wet period. (b) The statistical model performance of the wet period. (c) The observed and simulated streamflow plotted on a logarithmic scale for the dry period. (d) The statistical model performance of the dry period.
Figure 10. (a) The observed and simulated streamflow plotted on a logarithmic scale for the wet period. (b) The statistical model performance of the wet period. (c) The observed and simulated streamflow plotted on a logarithmic scale for the dry period. (d) The statistical model performance of the dry period.
Water 13 02132 g010
Figure 11. The monthly average groundwater level changes in the forest sub-catchments in the wet and dry periods.
Figure 11. The monthly average groundwater level changes in the forest sub-catchments in the wet and dry periods.
Water 13 02132 g011
Figure 12. The monthly average groundwater level changes in the paddy sub-catchments in the wet and dry periods.
Figure 12. The monthly average groundwater level changes in the paddy sub-catchments in the wet and dry periods.
Water 13 02132 g012
Table 1. Soil parameters for SCEUA-P global optimization, and river and aquifer parameters for manual calibration.
Table 1. Soil parameters for SCEUA-P global optimization, and river and aquifer parameters for manual calibration.
DescriptionParameterRangeUnit
Soil
Saturated soil moistures_per0.3–0.6-
Saturated vertical hydraulic conductivityks0.00001–0.08mm/s
Saturated lateral hydraulic conductivityksi0.00001–0.8mm/s
Wetting front soil suction headPSI49.5–316.3mm
Maximum (initial) infiltration capacityfo25–260mm/h
Minimum (asymptotic) infiltration capacityfc0–12mm/h
River
Saturated hydraulic conductivity of reverbedku_riv0.0001–0.000001mm/s
Aquifer
Storage coefficientS0.1–0.01-
Table 2. Areas and slopes of each sub-catchment in the Boryeong Dam catchment.
Table 2. Areas and slopes of each sub-catchment in the Boryeong Dam catchment.
RegionsR1R2R3
Sub-CatchmentsF 1P 1F 2F 2-1F 2-2P 2P 2-1P 2-2F 3P 3
Area (km2)40.10.944.81.82.15.80.80.861.24.8
Slope (-)0.260.150.400.300.280.190.160.150.230.11
F: forest, P: paddy.
Table 3. Daily groundwater use distributed according to the area ratio of sub-catchments in the irrigation period.
Table 3. Daily groundwater use distributed according to the area ratio of sub-catchments in the irrigation period.
Sub-CatchmentF 1F 2F 2-1F 2-2F 3
Groundwater use (m3/day)13284912022952
P 1P 2P 2-1P 2-2P 3
Groundwater use (m3/day)105510,0651388139214,553
Table 4. Amount of average daily intake rate for the Samsan and Hwasung agricultural reservoirs during the cultivation period.
Table 4. Amount of average daily intake rate for the Samsan and Hwasung agricultural reservoirs during the cultivation period.
Intake2012201320142015201620172018
Samsan (m3/day)5118443142573041184137443394
Hwasung (m3/day)5165430849164771347722633151
Table 5. Parameter values before and after calibration for the rainfall excess method, Green–Ampt method, and Horton infiltration method.
Table 5. Parameter values before and after calibration for the rainfall excess method, Green–Ampt method, and Horton infiltration method.
Rainfall ExcessGreen–AmptHorton
θsksksiθsksksiPSIfofc
Default0.37480.00060.0020.37480.00060.0021505150
Calibrated0.44610.0040.0010.49880.02390.0423248.649.88199.62
Table 6. Comparison of statistical results against observed and simulated streamflow using three infiltration analysis methods.
Table 6. Comparison of statistical results against observed and simulated streamflow using three infiltration analysis methods.
Rainfall ExcessR2NSEKGERMSE
Calibration0.8030.8020.8277.046
Validation0.7810.7220.8256.667
Green–AmptR2NSEKGERMSE
Calibration0.7110.6250.7759.691
Validation0.7220.4750.6169.163
HortonR2NSEKGERMSE
Calibration0.7680.7620.8327.719
Validation0.7840.7200.8136.688
Table 7. The differences in monthly average groundwater elevation between the wet and dry periods in the forest sub-catchments (m).
Table 7. The differences in monthly average groundwater elevation between the wet and dry periods in the forest sub-catchments (m).
MonthF 1F 2F 2-1F 2-2F 3
Jan0.03−0.06−0.38−0.030.06
Feb0.01−0.06−0.41−0.070.04
Mar0.02−0.04−0.42−0.080.04
Apr0.02−0.02−0.41−0.070.03
May−0.03−0.06−0.49−0.15−0.02
Jun0.01−0.01−0.45−0.120.02
Jul0.270.22−0.190.140.29
Aug0.650.530.260.550.68
Sep1.010.800.740.941.07
Oct0.900.610.750.850.99
Nov0.670.380.630.670.78
Dec0.470.200.470.490.58
Average0.340.210.010.260.38
Table 8. The differences in monthly average groundwater elevation between the wet and dry periods in the paddy sub-catchments (m).
Table 8. The differences in monthly average groundwater elevation between the wet and dry periods in the paddy sub-catchments (m).
MonthP 1P 2P 2-1P 2-2P 3
Jan0.040.060.050.060.10
Feb0.020.030.020.030.07
Mar0.000.01−0.010.000.02
Apr−0.01−0.01−0.02−0.010.00
May−0.03−0.03−0.05−0.04−0.04
Jun−0.03−0.03−0.05−0.03−0.04
Jul0.010.010.000.010.03
Aug0.150.130.140.160.18
Sep0.320.270.310.330.37
Oct0.290.250.300.310.35
Nov0.220.190.230.230.26
Dec0.170.150.180.180.21
Average0.100.090.090.100.12
Table 9. Groundwater pumping rate for the wet and dry periods by sub-catchment (m3/day).
Table 9. Groundwater pumping rate for the wet and dry periods by sub-catchment (m3/day).
WetDry
201020112012Avg.201420152016Avg.
F 18731710170314291449147414871470
F 2371396501423503537796612
F 2-11516201720223225
F 2-21718231923253628
F 36139339418291259131713331303
P 15591377137711041154116611991173
P 213,51114,78116,02814,7733286332328823164
P 2-11864203922112038453458398436
P 2-21868204422162043454459399437
P 322,32022,77822,81022,6362480269827252634
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Park, S.; Kim, H.; Jang, C. Impact of Groundwater Abstraction on Hydrological Responses during Extreme Drought Periods in the Boryeong Dam Catchment, Korea. Water 2021, 13, 2132. https://doi.org/10.3390/w13152132

AMA Style

Park S, Kim H, Jang C. Impact of Groundwater Abstraction on Hydrological Responses during Extreme Drought Periods in the Boryeong Dam Catchment, Korea. Water. 2021; 13(15):2132. https://doi.org/10.3390/w13152132

Chicago/Turabian Style

Park, Sanghyun, Hyeonjun Kim, and Choelhee Jang. 2021. "Impact of Groundwater Abstraction on Hydrological Responses during Extreme Drought Periods in the Boryeong Dam Catchment, Korea" Water 13, no. 15: 2132. https://doi.org/10.3390/w13152132

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