Next Article in Journal
Evolution of Water–Sediment Situation and Attribution Analysis in the Upper Yangtze River, China
Next Article in Special Issue
An Assessment of Geospatial Analysis Combined with AHP Techniques to Identify Groundwater Potential Zones in the Pudukkottai District, Tamil Nadu, India
Previous Article in Journal
Investigating Climate Change Effects on Evapotranspiration and Groundwater Recharge of the Nile Delta Aquifer, Egypt
Previous Article in Special Issue
Appraisal of Groundwater Vulnerability Pollution Mapping Using GIS Based GOD Index in Tiruchendur, Thoothukudi District, India
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Analysis of Net Groundwater Recharge Using Water Budget and Climate Change Scenarios

1
Department of Geological Sciences, BK21 School of Earth and Environmental Systems, Pusan National University, Busan 46241, Republic of Korea
2
Geoscience Co., Busan 48058, Republic of Korea
3
Radwaste Technology & Research Institute, Korea Radioactive Waste Agency, Gyeongju 38062, Republic of Korea
4
Department of Spatial Information Engineering, Pukyong National University, Busan 48513, Republic of Korea
5
Institute of Environmental Studies, Pusan National University, Busan 46241, Republic of Korea
*
Author to whom correspondence should be addressed.
Water 2023, 15(3), 571; https://doi.org/10.3390/w15030571
Submission received: 14 November 2022 / Revised: 20 January 2023 / Accepted: 28 January 2023 / Published: 1 February 2023
(This article belongs to the Special Issue Drought and Groundwater Development)

Abstract

:
Estimating the groundwater recharge rate is essential in all groundwater-related fields, including groundwater development, use, management, modeling, and contamination analysis. In this study, we proposed a combined method of water budget and climate change scenario for estimating the net groundwater recharge rate in the Nakdong River watershed (NRW), South Korea. For the climate change scenario method, the representative concentration pathway (RCP) 4.5 and 8.5 climate scenarios were adopted. First, using the water budget method from 2009 to 2018, the net groundwater recharge rate (NGRR) of 12.15–18.10% relative to annual precipitation (AP) was obtained, subtracting direct runoff (DR) of 21.18–25.32% relative to AP, evapotranspiration (EP) of 40.53–52.29% relative to AP, and baseflow of 12.42–17.84% relative to AP, from the AP (865–1494 mm). The average annual NGRR of the NRW was 200 mm (15.59%). Second, the mean NGRRs from 2009 to 2100 under the RCP 4.5 and RCP 8.5 scenarios were anticipated as 8.73% and 7.63%, respectively. The similarity between the water budget and climate change scenarios was confirmed using data from 2009 and 2018. According to the simple climate change scenario, it is predicted that annual precipitation will increase over the years while the groundwater level and net groundwater recharge rate will decrease. Nonetheless, the estimated NGRR by the water budget method in this study possesses uncertainty due to using potential ET instead of actual ET which should be estimated by considering soil water content.

1. Introduction

Groundwater originates mainly from rainfall and is recharged in aquifers. Watershed groundwater recharge is important for its development, use, conservation, and management. Moreover, it is critical for accurately estimating the groundwater recharge rate, which local and federal governments use for relevant policymaking. The groundwater recharge rate is one of the most challenging factors in assessing the quantitative evaluation of groundwater resources [1]. The groundwater recharge rate is calculated using methods, such as baseflow separation, water budget, water table fluctuation (WTF), groundwater modeling, climate change scenarios, and artificial intelligence. The net groundwater recharge rate (NGRR) using the water budget method is the amount of precipitation and groundwater inflow from other basins, excluding direct runoff, evapotranspiration, baseflow, pumping quantity, groundwater outflow into other basins, and coastal groundwater discharge [2,3]. The WTF method, a physical estimation method, estimates the water table rise in the phreatic aquifer and the groundwater recharge rate due to rainfall [4]. The hybrid water table fluctuation method (hybrid-WTF) combines the WTF method and the model of water content in unsaturated zones by considering the delay time of water reaching the water table from the land surface, as well as baseflow and groundwater level drops simultaneously in the condition of a permeable unsaturated zone and the underlying aquifer [5].
Moon et al. [6] estimated the groundwater recharge rate of four major river basins in South Korea using an improved WTF method based on groundwater level. In the rainy season (June–September) of 2005 in the uppermost watershed, Choi et al. [7] estimated average groundwater recharge rates of 24.75%, 13.4%, and 14.0%, respectively, using the WTF method, chloride mass balance method, and baseflow separation methods with chloride tracers. Chung et al. [8] calculated the groundwater recharge in the Mihocheon stream basin in South Korea using the SWAT-MODFLOW model [9], which is a surface water-groundwater integration model. In the Gangcheon Stream basin in Gyeonggi Province, South Korea, Noh et al. [10] calculated average groundwater recharge rates of 17.81% and 18.13%, respectively, using the hybrid-WTF method [5] and groundwater modeling. Bae [11] calculated the groundwater recharge rate using the natural resources conservation service curve number (NRCS-CN) and baseflow separation methods in five sub-watersheds and obtained a 30 mm/yr higher groundwater recharge rate using the baseflow separation and NRCS-CN methods. Lee and Bae [12] examined the variability in groundwater recharge rates due to urbanization and land cover changes. Viji et al. [13] calculated CN, including the antecedent moisture condition (AMC), for the Kundapallam watershed in India, considering land use and hydrological soil grouping (HSG). Ling et al. [14] proposed an improved soil conservation service (SCS)-CN method using nonparametric statistical analysis of rainfall and runoff data from the Wangjiaqiao watershed in China. Shi et al. [15] proposed the SCS-CN method by considering topographic slope, soil moisture, and storm duration factors, and applied the method to the Loess Plateau watershed in China.
Groundwater recharge rates are anticipated to change due to future climate change. Dripps and Bradbury [16] evaluated groundwater resources according to climate change by applying an improved soil water balance (SWB) model by considering soil moisture content, land cover, topography, climate data, and geographic information. Holman [17] analyzed the climate and socioeconomic impacts, uncertainty, and vulnerability of the groundwater recharge rate. Schindler et al. [18] investigated groundwater level variability according to soil type. They proposed a management plan for agricultural groundwater use in northeastern and central Germany, utilizing soil data at 368 monitoring points from 1951 to 2000 as well as using climate change scenarios from 2001 to 2055. Patricia et al. [19] proposed a standardized anomaly index for evaluating the impact of drought on a watershed scale in the Philippines, based on the spatiotemporal characteristics of groundwater level change. Jang et al. [20] estimated the groundwater level change from 2000 to 2012 in four major river basins in South Korea. Using climate change scenarios, they then predicted the groundwater level change from 2000 to 2100. Lee et al. [21,22] proposed coupled model development between groundwater recharge rate and climate change in river watersheds in Korea and calculated groundwater rate by using Visual HELP3. By using climate change scenarios for the Korean Peninsula and a geographic information system (GIS), Lee and Lee [23] estimated the groundwater recharge rates of 27.37% (1970 to 2000), 27.43% (2001 to 2030), 26.06% (2031 to 2050), and 27.88% (2051 to 2100) in the NRW. Lee and Lee [24] revealed the groundwater recharge mechanism in the Wonju and Uiwang areas through a time-series analysis of groundwater levels and rainfall. Kim et al. [25] conducted a time series analysis of the groundwater levels and precipitation on Geumjeongsan Mountain in Busan City. They suggested that the groundwater recharge rate of the bedrock aquifer is governed by fracture zones more than the direct rainfall effect.
The objective of this study is to estimate the net groundwater recharge rate (NGRR) from 2009 until 2100 by using a combined analysis of the water budget method and the climate change scenarios of representative concentration pathways (RCP) 4.5 and 8.5. Using data of 2009–2018, baseflow was calculated using streamflow data at the estuary of the NRW, and the groundwater recharge rate was estimated based on the water balance. The long-term groundwater recharge until 2100 was predicted using precipitation from the RCP 4.5 and RCP 8.5.

2. Study Area

The Nakdong River watershed (NRW) consists of eight sub-basins in the upstream area (Andong Dam, Imha Dam, Mainstream I, Naeseongcheon, Yeonggang, Byeongseongcheon, Ssanggyecheon, and Wicheon Basin), seven sub-basins in the midstream region (Mainstream II, Gamcheon, Mainstream III, Geumho River upstream, Geumho River downstream, Main stream IV, Hoecheon Basin), six subbasins in the downstream western region (Hwanggang, Mainstream V, Namgang upstream, Yangcheon, Namgang Dam, and Namgang lower basin), six sub-basins in the downstream eastern region (Main stream VI, Miryanggang, Mulgeum, Yangsancheon, and Hagueon, Seobun) Nakdong River basin). The Nakdong River, which originates in Taebaek City, Gangwon Province, passes through Yeongnam Province and, flows into the South Sea through Busan Metropolitan City, with the longest length (510 km) and the second largest area (23,647 km2) in South Korea.
A total of 122 wells under the National Groundwater Monitoring Network (NGMN) are distributed in the NRW (Figure 1). The surface elevation of the wells was the highest at 735.6 m (ND-3-SS well) and the lowest at 1.4 m (ND-3YN well). The highest groundwater level was 735.1 m, average mean sea level (amsl) (ND-3-SS well), and the lowest was −6.8 m amsl (ND-4-DC well), with an average groundwater level of 98.6 m amsl (Table 1). The highest groundwater temperature was 20.6 °C (ND-3-DGB well) and the lowest was 11.0 °C (ND-1-NA well), with the average groundwater temperature of all the wells being 15.0 °C. In addition, the average groundwater electrical conductivity (EC) of all the wells was 536 μS/cm, with the highest reading of 11,694.6 μS/cm (ND-4-DC well near the coast of the South Sea) and the lowest reading of 83.9 μS/cm (ND-4-USDC well).
Based on the precipitation distribution map of the NRW, prepared from precipitation data from 30 meteorological stations (Figure 2) from 1996 to 2019, the mean annual precipitation in the NRW was 1277 mm. The precipitation data of the 30 meteorological stations corresponded with the nearest monitoring well. The highest precipitation of 1989 mm/yr was observed in the ND-4-SH well, and the lowest precipitation of 917 mm/yr was observed in the ND-2-PC well. The precipitation in each well was used to calculate the net groundwater recharge.

3. Methods

3.1. Method of Water Budget

In the hydrological cycle, the net groundwater recharge rate (NGRR) is the water that reaches the saturation zone through the unsaturated zone, and is calculated as the amount excluding direct runoff (DR), baseflow (BF), the outflow (OU) from the aquifer to the sea, the inflow (RN) of other geological layers, evapotranspiration (ET), as well as pumping (DO), if it occurs, from precipitation (P) [2,3]:
NGRR = PDRETBFDOOU + RN
Direct runoff (DR) is the amount that flows through the land surface and enters streams. If there is no RN and DO with no consideration of OU, then Equation (1) is transformed into:
NGRR = PDRETBF
Direct runoff (DR) is inversely proportional to the amount of infiltration and directly proportional to the water content in the soil layers. When there is no DR data, the DR is generally calculated using the SCS-CN method [27,28]. The SCS-CN method, DR should be greater than 0, and P should be greater than 0.2:
D R = ( P 0.2 S ) 2 P + 0.8 S
where S is governed by the antecedent soil moisture condition (AMC) and depends on the runoff curve number (CN), which is inversely proportional to DR. The AMC values are determined by 5 or 3-day accumulated antecedent precipitation.
Evapotranspiration (ET) has been estimated in several ways [29,30,31,32,33,34]. Actual evapotranspiration (AET) is a fraction of potential evapotranspiration (PET) that is assumed as ET in the saturated condition of the soil. In this study, the PET (ET0, mm/d) as an approximation of the AET was estimated by using the Penman–Monteith equation that is widely used in the world [35] and that considers atmospheric temperature, humidity, wind speed, and radiation. Additionally, Yang et al. [36] estimated the ratio of AET over PET to be 0.96 in the coastal area of Busan City belonging to the NRW.
Baseflow (BF), the water provided by riverside aquifers, is estimated by hydrograph separation [37,38] and the recession of BF is influenced by topography, watershed shape, soil, geological media, etc. The volume of BF discharged from riverside aquifers is given by Meyboom [37]:
V t p = Q 0 t 1 2.3026
Finally, the baseflow (BF) during seasons adjacent to A and B is expressed by:
B F = 2 ( Q B Q A ) t 1 2.3026
where QA is the baseflow at the critical time of recession after the hydrograph peak in season A and QB is the baseflow at the critical time of recession after the hydrograph peak in season B. Time t1 is the time taken for the baseflow recession to decline by one log cycle.

3.2. Method of Climate Change Scenarios

Climate change scenarios calculate future climate factors, such as temperature, precipitation, wind, and humidity, by applying changes in radiative forcing caused by anthropogenic causes (greenhouse gases, aerosols, and land use changes) to the Earth system model. There are six climate change scenarios, A1 (A1FI, A1T, and A1BI), A2, B1, and B2 depending on the projected carbon dioxide emissions [39]. The IPCC (2007) determined the greenhouse gas concentration based on the amount of radiation exerted on the atmosphere by human activities according to the RCP. The Korea Meteorological Administration (KMA) predicted future climate change up to 2100 according to the RCP scenario after simulating past climates for natural and anthropogenic forcing from 1860 to 2005. The KMA developed a climate change scenario for the Korean Peninsula (www.climate.go.kr, accessed on 13 November 2022) using meteorological data (averaged from 2001 to 2010) obtained from the Automated Surface Observing System (ASOS) and the Automatic Weather System (AWS). The 1-km climate change scenario for South Korea was obtained by applying a 12.5-km resolution climate change scenario to the PRISM based downscaling estimation model (PRIDE). The calculated 1-km grid data is converted into administrative district-level data using GIS techniques. This 1-km climate change scenario can reflect the effects of complex topography in South Korea, which cannot be realized in the global climate model.
In this study, using the climate change scenarios of RCP 4.5 and RCP 8.5, NGRRi is calculated by Jang et al. [20]:
N G R R i = P i · f
where recharge ratio f is determined by:
f = Δ H i · S y / P i
where Sy is the specific yield. ∆Hi is the groundwater level difference between Himax (the highest monthly groundwater level in the ith year) and Himin (the lowest groundwater level in the ith year) for each year between 2009 and 2018. The values of f are computed by using Equation (7) with ∆Hi and Pi at each point in the historical period of 2009–2018. Finally, the mean f value of all the wells for the period of 2009–2018 is applied to each well with corresponding precipitation for the period of 2019–2100. Pi is annual precipitation of the ith year.

4. Results

4.1. Net Groundwater Recharge Rate by the Water Budget

The net groundwater recharge rate (NGRR) by the water budget was calculated for 264 standard sub-watersheds (Figure 1). For the determination of direct runoff (DR) by the SCS-CN method, the SCS-CN values were obtained using the land cover map and soil map of the NRW. The 21 land cover types were discovered in the study area, including forest 67.17% (broadleaf forest 30.58%, coniferous forest 25.39%, mixed forest 11.20%), residential areas 1.01%, industrial areas 1.32%, commercial areas 1.17%, cultural and sports areas 0.02%, transportation areas 0.95%, public facilities areas 0.01%, rice paddy fields 7.57%, farm fields 11.04%, greenhouse cultivation areas 2.39%, orchard 0.43%, other cultivation areas 0.21%, natural grasslands 1.51%, artificial grasslands 0.58%, inland wetlands 0.01%, coastal wetlands 0.01%, natural bare fields 1.05%, and artificial bare fields 0.44%. Drainage grades were divided into A (55.38%), B (20.58%), C (3.03%), and D (21.01%). The S values were calculated using the SCS-CN values, considering the area corresponding to each drainage class. Finally, annual DR values from 2009 to 2018 were calculated as 183.21 mm (21.18%) to 372.25 mm (25.32%) by using Equation (3). The average annual DR for the ten years is 308.28 mm, which is 24.14% of the average annual precipitation of 1277 mm (Table 2).
Annual ET0 was calculated using the Penman-Monteith method, by using the daily average temperature (14.5–15.7 °C), daily maximum temperature (32.5–37.3 °C), daily minimum temperature (−12.8–−6.0 °C), solar radiation (5.99–7.70 h), atmospheric pressure (100.7 to 100.9 kPa), daily average water vapor pressure (1.23 to 1.27 kPa), and daily average wind speed (3.1–3.4 m/s) from 2009 to 2018 (Table 3). Precipitation, temperature, and solar radiation amounts were obtained from 30 meteorological stations in the NRW. Due to the difficulty in directly measuring or calculating longwave radiation [40], 30% of the radiation energy was applied as the net radiation value according to the latitude effect provided by FAO. The annual PET from 2009 to 2018 ranged from 452.27 mm (40.53%) to 605.59 mm (52.29%), and the annual average ET for the ten years was 567.33 mm, which corresponded to 44.43% of the 10-year average annual precipitation (1277 mm).
The streamflow rates of peak and re-peak (2763–21,657 m3/d), numbered from 1 to 10, with a time interval of 1.6–3.2 months between two peaks at the estuary of the NRW (Figure 3), were used for the calculation of the BF values. Using Equation (5), from 2009 to 2018, the BF values of the NRW were between 122.64 mm (12.42%) and 254.81 mm (17.84%), and the 10-year average BF value was 202.32 mm, which was 15.84% of the 10-year annual average precipitation (Table 2).
From 2009 to 2018 in the NRW, the NGRR values were estimated to range from 12.15% to 18.10%. The average of the NGRR values in the NRW from 2009 to 2018 is 200 mm (15.66%) when DR of 308 mm (24.12%), PET of 567 mm (44.40%), and BF of 202 mm (15.82%) are subtracted from the average AP (1277 mm).

4.2. Net Groundwater Recharge Rate under the Climate Change Scenarios

Before the RCP 4.5 and RCP 8.5 climate change scenarios, the recharge ratio (f), using the observed 2009–2018 data, was calculated by multiplying ∆Hi and specific yield (Sy) in Equation (7). The Sy value was determined to be 0.12 as a result of matching NGRR values with the NGRR values by the water budget method. In contrast, the mean Sy value of the NRW was reported to be 0.0134 by Moon et al. [41]. Next, the NGRR values from 2019 to 2100 were computed using the predicted precipitation data and f values in Equation (6), based on RCP 4.5 and RCP 8.5 climate change scenarios with a resolution of 1 km. Based on the RCP 4.5 and RCP 8.5 climate change scenarios, Figure 4 and Table 4 show mean AP (mm/yr) for the periods of 2009–2018, 2019–2039, 2040–2059, 2060–2079, and 2080–2100. The AP of the RCP 4.5 displays an increasing tendency from 2040 to 2079 and then a decreasing tendency from 2080 to 2100, while the AP of the RCP 8.5 displays an increasing tendency from 2019 until 2100.
According to the climate change scenarios of RCP 4.5 and RCP 8.5, it is predicted that the decline rate of the groundwater level by RCP 4.5 and RCP 8.5 will be 0.023 m/yr and 0.024 m/yr, respectively, reaching a decline of 2.3 m and 2.4 m in 2100, respectively. In general, the decrease in groundwater level over a long period can be attributed to climate change and various man-made effects of groundwater pumping, irrigation, land-use change, etc. [4]. In this study, it is judged that the groundwater level decline with the annual precipitation decrease in the long term is caused by a decrease in infiltration rate into the subsurface as well as an increase in precipitation and an increase in direct runoff during the wet season relative to smaller precipitation in the dry season [20]. Moreover, a temperature increase of 3.5–4.0 °C until 2100 on the Korean peninsula [42] will increase the ET rate [43]. On the other side, urbanization can decrease ET rates by increasing runoff [44]. The NGRR values will be decreased to 7.05% based on the RCP 4.5 and to 5.90% based on the RCP 8.5 climate change scenarios in the period of 2080–2100 (Table 4), with a slightly decreasing tendency (Figure 5). The mean NGRR from 2009 to 2100 was estimated to be 8.73% in the RCP 4.5 scenario and 7.63% in the RCP 8.5 scenario, being similar to each other. Precipitation in summer for the periods of 2009–2019, 2020–2039, 2040–2059, 2060–2079, and 2080–2100 by using the RCP 4.5 and 8.5 climate change scenarios is intensified when going to 2100 (Table 5; Figure 6).

5. Discussion

Thus far, there has been no case of calculating the NGRR based on the long-term water budget method and simple climate change scenarios for the total area of the NRW, except in this study. Moon et al. [6] estimated the groundwater recharge rate in the NRW using the WTF method and obtained 6.1%, which is much smaller than the 15.54% of the NGRR according to the water budget method in this study. In this study, the observed streamflow rate was compared with the baseflow and calculated direct runoff by the SCS-CN method, showing good matching such that the observed streamflow rate of 323.59–622.42 mm/yr and the baseflow plus the calculated direct runoff of 321.89–617.99 mm/yr from 2009 to 2018 resulted in slightly greater values of the observed streamflow rate than those of the calculated direct runoff and baseflow (Figure 7).
The NGRR (200 mm, 15.66% of the AP) of the NRW in this study was compared with the average annual groundwater recharge rate (196 mm, 14.9% of the AP) in the Yanggok-ri subwatershed in the Geum River watershed from 2001 to 2018 [45]. Hence, the NGRR of the NRW is slightly greater than that of the Yanggok-ri subwatershed in the Geum River watershed, owing to the differences in the hydrogeological and climatic conditions. In addition, the NGRR of the NRW of 12.42–17.84% during 2009–2018 showed a much smaller variation than that of 3.6–28.2% in the Yanggok-ri subwatershed of the Geum River watershed during 2001–2018. In addition, the larger the watershed, the smaller the change in the NGRR. In this study, the mean of the NGRRs from 2009 to 2100 was 8.73% and 7.63%, respectively, using the RCP 4.5 and RCP 8.5 scenarios. In contrast, the NGRR of 26.7 mm to 432.5 mm (average 174.6 mm) in the Yanggok-ri subwatershed of the Geum River watershed from 2019 to 2100 corresponds to 2.8–45.1% of the precipitation in 2100, based on the climate change scenario RCP 8.5 [38].
Looking at the change in the NGRR in the period of 2080–2100, the NGRR predicted by the climate change scenarios of the RCP 4.5 and RCP 8.5 is predicted to decrease from 13.73% to 7.05% and 14.43% to 5.90%, respectively. The NGRRs by the RCP 4.5 and RCP 8.5 scenarios are almost the same during 2009–2018, with a slight difference from the NGRR by the water budget analysis (Figure 8). Therefore, it is judged that the estimated values of NGRR in the climate change scenarios RCP 4.5 and RCP 8.5 are reasonable. Figure 9A shows the spatial distribution of the average NGRRs in the water budget analysis during 2009–2018. The NGRR showed high values in the eastern and southern regions but exhibited relatively low values in the inland areas, closely resembling precipitation distribution. In addition, the NGRR in the upstream area of the NRW was relatively smaller than that in the downstream area. On the other side, in Figure 9B, the NGRR distribution in 2080–2100 using the RCP 4.5 scenario shows overall higher values in the eastern region than those in the inland and southern areas, displaying a different spatial distribution of the NGRR from that using the water budget analysis in the period of 2009–2018 and showing smaller NGRRs than those by the water budget analysis.
The estimated NGRR by the water budget analysis in this study possesses uncertainty due to using PET instead of AET, which can be estimated by considering soil water content [36,46,47]. A future study needs to reveal a more accurate NGRR in the NRW which was not possible due to insufficient data on soil water content. For the climate change scenarios, there is uncertainty about the specific yield (Sy), which is an average value (0.12) of the NRW since Sy is spatially different depending on geological formations in the NRW. Finally, the water budget analysis and the climate change scenarios also imply uncertainty due to not considering the effect of land cover change.

6. Conclusions

In this study, the change in the NGRR of the NRW was estimated using the water budget method and the climate change scenarios of RCP 4.5 and RCP 8.5. By using the water budget method and precipitation data from 2009 to 2018, direct runoff, evapotranspiration, and baseflow were estimated prior to determining the NGRR. The DR in the NRW by the SCS-CN method was calculated to be 183.21 mm (21.18%) to 372.25 mm (25.32%). The PET obtained by the Penman–Monteith method ranged from 452.27 mm (40.53%) to 605.59 mm (52.29%). The baseflow in the NRW using the baseflow separation method was 122.64 mm (12.42%) to 254.81 mm (17.84%).
The mean annual NGRR was determined to be 200 mm (15.66%) by subtracting the mean annual DR of 308 mm (24.12%), the mean annual PET of 567 mm (44.40%), and the mean annual BF of 202 mm (15.82%) from the mean annual precipitation of 1277 mm in the NRW for ten years (2009–2018), with an increasing trend in the NGRR over the ten years. The regional distribution of the NGRR showed high values in the eastern and southern regions but exhibited relatively low values in the inland areas, closely resembling precipitation distribution. In addition, the NGRR in the upstream area of the NRW was relatively smaller than that in the downstream area.
Under a decreased tendency of precipitation and the constant recharge ratio of the NRW, the mean NGRR for the entire NRW over the long term until 2100 is predicted to decrease from 13.73 to 7.05% based on the RCP 4.5 and from 14.43 to 5.90% based on the RCP 8.5. The NGRR distribution in 2080 – 2100 using the RCP 4.5 scenario shows overall high values in the eastern region and lower values in the inland and southern areas, similar to the spatial distribution of the NGRR, which displays a different spatial distribution of the NGRR from that using the water budget analysis in the period of 2009–2018, with smaller NGRR values than those by the water budget analysis.
The accuracy of the prediction of the NGRR according to the climate change scenarios RCP 4.5 and RCP 8.5 was verified by comparing the 2009 to 2018 water budget analysis and climate change scenarios. However, the robustness of the estimation of the NGRR is constrained by the uncertainty of the PET, the specific yield, and the land cover change. Therefore, it is judged that the estimated NGRR based on the climate change scenarios RCP 4.5 and RCP 8.5 was reasonable. The NGRR of the NRW presented in this study will be useful in the long-term development and management of groundwater resources in watersheds worldwide as well as in the NRW. Future research will focus on factors such as change in runoff, baseflow, and evapotranspiration over a long period.

Author Contributions

Conceptualization, S.-M.Y., H.-T.J. and S.-Y.H.; methodology, H.-T.J., J.-Y.C., S.-Y.H. and J.K.; software, S.-M.Y., H.-T.J. and J.-Y.C.; validation J.K. and S.-Y.H.; writing, original draft preparation, S.-M.Y., H.-T.J. and S.-Y.H.; writing, review and editing, J.-Y.C. and S.-Y.H.; visualization, S.-M.Y.; supervision, S.-Y.H.; project administration, S.-Y.H.; funding acquisition, S.-Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the National Research Foundation of Korea (NRF) grant funded by the Ministry of Science and ICT (2020R1A2B5B02002198).

Data Availability Statement

Data is unavailable due to privacy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sophocleous, M.A. Combining the soilwater balance and water-level fluctuation methods to estimate natural groundwater recharge: Practical aspects. J. Hydrol. 1991, 124, 229–241. [Google Scholar] [CrossRef]
  2. Fetter, C.W. Applied Hydrology, 4th ed.; Prentice Hall: Englewood Cliffs, NJ, USA, 2003; 676p. [Google Scholar]
  3. Schwartz, F.W.; Zhang, H. Fundamentals of Groundwater; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2003; 583p. [Google Scholar]
  4. Healy, R.W.; Cook, P.G. Using groundwater levels to estimate recharge. Hydrogeol. J. 2002, 10, 91–109. [Google Scholar] [CrossRef]
  5. Park, E.; Parker, J.C. A simple model for water table fluctuations in response to precipitation. J. Hydrol. 2008, 356, 344–349. [Google Scholar] [CrossRef]
  6. Moon, S.-K.; Woo, N.C.; Lee, K.S. Statistical analysis of hydrographs and water-table fluctuation to estimate groundwater recharge. J. Hydrol. 2004, 292, 198–209. [Google Scholar] [CrossRef]
  7. Choi, I.H.; Woo, N.C.; Kim, S.-J.; Moon, S.-K.; Kim, J. Estimation of the groundwater recharge rate during a rainy season at a headwater catchment in Gwangneung, Korea. Korean J. Agric. For. Meteorol. 2007, 9, 75–87. [Google Scholar] [CrossRef]
  8. Chung, I.-M.; Kim, N.W.; Lee, J. Estimation of groundwater recharge by considering runoff process and groundwater level variation in watershed. J. Soil Groundw. Environ. 2007, 12, 19–32. [Google Scholar]
  9. Kim, N.W.; Chung, I.-M.; Won, Y.S. An integrated surface water-groundwater modeling by using fully combined SWAT-MODFLOW model. J. Korean Soc. Civil Eng. 2006, 26, 481–488. [Google Scholar]
  10. Noh, D.N.; Park, H.J.; Cheong, J.-Y.; Hamm, S.-Y. Groundwater recharge analysis and comparison using hybrid water-table fluctuation method and groundwater modeling: A case of Gangcheon basin in Yeoju City. J. Geol. Soc. Korea 2018, 54, 169–181. [Google Scholar] [CrossRef]
  11. Bae, S.K. Applicability of NRCS-CN method for the estimation of groundwater recharge. J. Korean Soc. Civ. Eng. B 2005, 25, 425–430. [Google Scholar]
  12. Lee, S.H.; Bae, S.K. Groundwater balance in urban area. J. Environ. Sci. 2011, 20, 1553–1560. [Google Scholar]
  13. Viji, R.; Prasanna, P.R.; Ilangovan, R. Gis Based SCS-CN Method for estimating runoff In Kundahpalam watershed, Nilgries district, Tamilnadu. Earth Sci. Res. J. 2015, 19, 59–64. [Google Scholar]
  14. Ling, L.; Yusop, Z.; Yap, W.-S.; Tan, W.L.; Chow, M.F.; Ling, J.L. A calibrated, watershed-specific SCS-CN Method: Application to Wangjiaqiao watershed in the three gorges area, China. Water 2019, 12, 60. [Google Scholar] [CrossRef]
  15. Shi, W.; Wang, N. An improved SCS-CN method incorporating slope, soil moisture, and storm duration factors for runoff prediction. Water 2020, 12, 1335. [Google Scholar] [CrossRef]
  16. Dripps, W.R.; Bradbury, K.R. A simple daily soil-water balance model for estimating the spatial and temporal distribution of groundwater recharge in temperate humid areas. Hydrogeol. J. 2007, 15, 433–444. [Google Scholar] [CrossRef]
  17. Holman, I.P. Climate change impacts on groundwater recharge- uncertainty, shortcomings, and the way forward? Hydrogeol. J. 2005, 14, 637–647. [Google Scholar] [CrossRef] [Green Version]
  18. Schindler, U.; Steidl, J.; Müller, L.; Eulenstein, F.; Thiere, J. Drought risk to agricultural land in Northeast and Central Germany. J. Plant Nutr. Soil Sci. 2007, 170, 357–362. [Google Scholar] [CrossRef]
  19. Patricia, A.J.S.; Wang, L.; Koike, T. Modeling the hydrologic responses of the Pampanga River basin, Philippines: A quantitative approach for identifying droughts. Water Resour. Res. 2011, 47, 1–21. [Google Scholar]
  20. Jang, S.; Hamm, S.-Y.; Yoon, H.S.; Kim, G.B.; Park, J.H.; Kim, M.S. Predicting long-term change of groundwater level with regional climate model in South Korea. Geosci. J. 2015, 19, 503–513. [Google Scholar] [CrossRef]
  21. Lee, J.H.; Jeon, S.W.; Lee, M.J.; Hong, H.J. Coupled Model Development between Groundwater Recharge Rate Quantity and Climate Change in River Watershed; Korea Environment Institute: Sejong City, Republic of Korea, 2009; 09-06-52(4); 142p. [Google Scholar]
  22. Lee, J.H.; Jeon, S.W.; Lee, M.J.; Hong, H.J. Coupled Model Development between Groundwater Recharge Rate Quantity and Climate Change in River Watershed II; Korea Environment Institute: Sejong City, Republic of Korea, 2010; 10-02-97(5); 149p. [Google Scholar]
  23. Lee, M.J.; Lee, J.H. Coupled model development between groundwater recharge quantity and climate change using GIS, Green Growth Research Report. J. Korean Assoc. Geogr. Inf. Stud. 2011, 14, 36–51. [Google Scholar] [CrossRef]
  24. Lee, J.Y.; Lee, K.K. A Comparative study on characteristics of waterlevel responses to rainfall in the two aquifer systems. J. Soil Groundw. Environ. 2002, 7, 3–14. [Google Scholar]
  25. Kim, T.W.; Hamm, S.-Y.; Cheong, J.Y.; Ryu, S.M.; Lee, J.H.; Son, K.T.; Kim, N.H. Time series and groundwater recharge analyses using water fluctuation data in mountain Geumjeong area. J. Environ. Sci. 2008, 17, 257–267. [Google Scholar]
  26. Jeon, H.-T. A Study of Long-Term Net Recharge Rate of Groundwater in Korea. Ph.D. Dissertation, Pusan National University, Busan, Republic of Korea, 2021; 167p. [Google Scholar]
  27. Morel-Seytoux, H.J.; Verdin, J.P. Extension of the SCS Rainfall-Runoff Methodology for Ungauged Watersheds; Colorado State University: Fort Collins, CO, USA, 1981. [Google Scholar]
  28. Aller, L.; Bennett, T.; Lehr, J.H.; Petty, R.H.; Hackett, G. DRASTIC: A Standardized System for Evaluating Ground Water Pollution Potential Using Hydrogeologic Settings; Prepared for U.S. EPA Office of Research and Development, Ada, OK; National Water Well Association: Dublin, OH, USA, 1987. [Google Scholar]
  29. Thornthwaite, C.W. Report of the Committee on Transpiration and Evaporation, 1943–1944. Eos Trans. Am. Geophys. Union 1944, 25, 683–693. [Google Scholar]
  30. Penman, H.L. Natural evapotranspiration from open water, bare soil, and grass. Proc. R. Soc. Lond. Ser. A 1948, 193, 120–145. [Google Scholar]
  31. Blaney, H.F.; Criddle, W.D. Determining Water Requirements in Irrigated Areas from Climatological and Irrigation Data; USDA (SCS) TP-96; US Department of Agriculture: Washington, DC, USA, 1950.
  32. Turc, L. Evaluation des besoins en eau d’irrigation. evapotranspiration potentielle, formulation simplifie et mise a jour. Ann. Agron 1963, 12, 13–49. [Google Scholar]
  33. Doorenbos, J.; Pruitt, W.O.; Aboukhaled, A. Crop water requirements. In Food and Agriculture Organization Irrigation and Drainage Paper 24, Technical Report; FAO: Rome, Italy, 1977; 156p. [Google Scholar]
  34. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapo-Transpiration-Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; Volume 300, p. D05109. [Google Scholar]
  35. Liu, Z. Estimating land evapotranspiration from potential evapotranspiration constrained by soil water at daily scale. Sci. Total Environ. 2022, 834, 155327. [Google Scholar] [CrossRef]
  36. Yang, S.I.; Kang, D.H.; Kwon, B.H.; Kim, B.W. Influence of land use and meteorological factors for evapotranspiration estimation in the coastal urban area. J. Environ. Sci. Int. 2010, 19, 295–304. [Google Scholar] [CrossRef]
  37. Meyboom, P. Estimating groundwater recharge from stream hydrographs. J. Geophys. Res. 1961, 66, 1203–1214. [Google Scholar] [CrossRef]
  38. Rorabaugh, M.I. Estimating changes in bank storage and ground-water contribution to stream flow. Inter. Assoc. Sci. Hydrol. 1964, 63, 432–441. [Google Scholar]
  39. IPCC. Climate Change 2007: The Physical Science Basis. In The Fourth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2007. [Google Scholar]
  40. Choi, D.-H.; Lee, B.-Y.; OH, H.-Y. Long and short wave radiation and correlation analysis between downtown and suburban area(Ⅱ). J. Korean Sol. Energy Soc. 2013, 33, 101–110. [Google Scholar] [CrossRef]
  41. Moon, S.-K.; Woo, N.C.; Lee, K.S. A study on the relation between types and recharges of groundwater: Analysis on national groundwater monitoring network data. J. Soil Groundw. Environ. 2002, 7, 45–59. [Google Scholar]
  42. Keem, M.; Ko, I.; Kim, S. An analysis of the effect of climate change on Nakdong River flow condition using CGCM’s future climate information. J. Korean Soc. Water Qual. 2009, 25, 863–871. [Google Scholar]
  43. Snyder, R.; Moratiel, R.; Song, Z.; Swelam, A.; Jomaa, I.; Shapland, T. Evapotranspiration response to climate change. Acta Hortic. 2011, 922, 91–98. [Google Scholar] [CrossRef]
  44. Mazrooei, A.; Reitz, M.; Wang, D.; Sankarasubramanian, A. Urbanization impacts on evapotranspiration across various spatio-temporal scales. Earth’s Future 2021, 9, e2021EF002045. [Google Scholar] [CrossRef]
  45. Ha, K.; Park, C.; Kim, S.; Shin, E.; Lee, E. Groundwater recharge evaluation on Yangok-ri area of Hongseong using a distributed hydrologic model (VELAS). Econ. Environ. Geol. 2021, 54, 161–176. [Google Scholar] [CrossRef]
  46. Filgueiras, R.; Almeida, T.S.; Mantovani, E.C.; Dias, S.H.B.; Fernandes-Filho, E.I.; Cunha, F.F.; Venancio, L.P. Soil water content and actual evapotranspiration predictions using regression algorithms and remote sensing data. Agric. Water Manag. 2020, 241, 106346. [Google Scholar] [CrossRef]
  47. Senay, G.B.; Kagone, S.; Velpuri, N.M. Operational global actual evapotranspiration: Development, evaluation, and dissemination. Sensors 2020, 20, 1915. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Location of the wells of the national groundwater monitoring network and the standard subwatersheds in the Nakdong River watershed.
Figure 1. Location of the wells of the national groundwater monitoring network and the standard subwatersheds in the Nakdong River watershed.
Water 15 00571 g001
Figure 2. Location of the weather monitoring stations in the Nakdong River watershed [26].
Figure 2. Location of the weather monitoring stations in the Nakdong River watershed [26].
Water 15 00571 g002
Figure 3. Streamflow in the Nakdong River watershed from 2009 until 2018. The numbers indicate the peaks for estimating baseflow.
Figure 3. Streamflow in the Nakdong River watershed from 2009 until 2018. The numbers indicate the peaks for estimating baseflow.
Water 15 00571 g003
Figure 4. Predicted precipitation of 2009 to 2100 by the RCP 4.5 and 8.5 climate scenarios in the Nakdong River watershed.
Figure 4. Predicted precipitation of 2009 to 2100 by the RCP 4.5 and 8.5 climate scenarios in the Nakdong River watershed.
Water 15 00571 g004aWater 15 00571 g004b
Figure 5. Predicted net groundwater recharge rate in the Nakdong River watershed for 2019–2100.
Figure 5. Predicted net groundwater recharge rate in the Nakdong River watershed for 2019–2100.
Water 15 00571 g005aWater 15 00571 g005b
Figure 6. Seasonal trend of precipitation for the periods of 2009–2019, 2020–2039, 2040–2059, 2060–2079, and 2080–2100 using the RCP 4.5 and 8.5 climate change scenarios.
Figure 6. Seasonal trend of precipitation for the periods of 2009–2019, 2020–2039, 2040–2059, 2060–2079, and 2080–2100 using the RCP 4.5 and 8.5 climate change scenarios.
Water 15 00571 g006
Figure 7. The observed streamflow and calculated direct runoff during 2009–2018.
Figure 7. The observed streamflow and calculated direct runoff during 2009–2018.
Water 15 00571 g007
Figure 8. The NGRR (%) relative to the corresponding year from the water budget analysis and RCP 4.5 and 8.5 scenarios from 2009 until 2018.
Figure 8. The NGRR (%) relative to the corresponding year from the water budget analysis and RCP 4.5 and 8.5 scenarios from 2009 until 2018.
Water 15 00571 g008
Figure 9. Distribution map of the average NGRRs in the Nakdong River watershed by (A) the water budget analysis of the period 2009–2018 and (B) the RCP 4.5 scenario in 2080–2100.
Figure 9. Distribution map of the average NGRRs in the Nakdong River watershed by (A) the water budget analysis of the period 2009–2018 and (B) the RCP 4.5 scenario in 2080–2100.
Water 15 00571 g009
Table 1. Well information statistics in the Nakdong River watershed for 2009–2018.
Table 1. Well information statistics in the Nakdong River watershed for 2009–2018.
Surface
Elevation
(m)
Well
Depth
(m)
Groundwater
Level
(m, amsl)
Temp.
(°C)
EC
(μS/cm)
Max735.6168.0735.120.611,695
Min1.440.0−6.811.083.9
Arithmetic mean105.679.798.615.0536
Median78.970.072.415.0320
Standard deviation105.419.8104.61.21083
Range734.2128.0741.99.611,611
Table 2. Estimated NGRR, DR, ET0, and BF of the Nakdong River watershed, with the ratios (%) to the annual average precipitation (P) during 2009–2018.
Table 2. Estimated NGRR, DR, ET0, and BF of the Nakdong River watershed, with the ratios (%) to the annual average precipitation (P) during 2009–2018.
YearP, mmDR, mm (%)ET0, mm (%)BF, mm (%)NGRR, mm (%)
20091191 274.91 (23.08)567.52 (47.65)196.95 (16.54)151.62 (12.73)
20101345 315.14 (23.43)590.48 (43.90)239.94 (17.84)199.44 (14.83)
20111494 372.25 (24.92)605.59 (40.53)254.81 (17.06)261.35 (17.49)
20121464 367.67 (25.11)602.13 (41.13)254.37 (17.38)239.83 (16.38)
20131072 254.72 (23.76)547.28 (51.05)133.09 (12.42)136.91 (12.77)
20141363 345.11 (25.32)592.25 (43.45)189.08 (13.87)236.56 (17.36)
20151025 238.26 (23.24)506.74 (49.44)155.47 (15.17)124.53 (12.15)
20161465 363.33 (24.80)604.17 (41.24)236.25 (16.13)261.25 (17.83)
2017865 183.21 (21.18)452.27 (52.29)122.64 (14.18)106.88 (12.36)
20181482 368.18 (24.84)604.89 (40.82)240.64 (16.24)268.29 (18.10)
Table 3. Factors for the estimation of ET0.
Table 3. Factors for the estimation of ET0.
Average Daily Temp.
(°C)
Daily Max. Temp.
(°C)
Daily Min. Temp.
(°C)
Solar
Radiation (h)
Atmospheric Pressure
(kPa)
Average Daily Water Vapor Pressure
(kPa)
Average Daily
Wind Speed
(m/s)
200915.232.5−7.65.99100.71.253.4
201014.934.1−8.16.25100.81.263.3
201114.633−12.86.37100.91.233.3
201214.534.5−9.97.12100.71.273.3
201315.335−10.77.7100.71.243.4
201415.132.9−66.89100.81.253.2
201515.433.5−7.87.04100.81.243.1
201615.737.3−10.26.99100.81.263.1
201715.236.2−7.77.55100.81.253.2
201815.136.4−9.97.41100.81.243.2
Table 4. Estimated NGRR (%), annual precipitation (P, mm), recharge ratio (f), and specific yield (Sy) in the Nakdong River watershed for the periods of 2009–2018, 2019–2039, 2040–2059, 2060–2079, and 2080–2100 by the RCP 4.5 and RCP 8.5 climate change scenarios.
Table 4. Estimated NGRR (%), annual precipitation (P, mm), recharge ratio (f), and specific yield (Sy) in the Nakdong River watershed for the periods of 2009–2018, 2019–2039, 2040–2059, 2060–2079, and 2080–2100 by the RCP 4.5 and RCP 8.5 climate change scenarios.
PeriodRCP 4.5RCP 8.5
(Starting Year–Ending Year)PSyNGRRPSyNGRR
2009–20181455.10.1213.731445.20.1214.43
2019–20391312.60.127.741387.40.126.31
2040–20591493.10.127.811424.70.125.78
2060–20791538.50.127.301474.80.125.72
2080–21001414.40.127.051585.00.125.90
Table 5. Maximum, minimum, and mean precipitation (mm/mon) based on the RCP 4.5 and RCP 8.5 climate change scenarios for the periods of 2009–2018, 2019–2039, 2040–2059, 2060–2079, and 2080–2100.
Table 5. Maximum, minimum, and mean precipitation (mm/mon) based on the RCP 4.5 and RCP 8.5 climate change scenarios for the periods of 2009–2018, 2019–2039, 2040–2059, 2060–2079, and 2080–2100.
PeriodRCP 4.5RCP 8.5
(Starting Year–Ending Year)MaxMinMeanMaxMinMean
2009–2018520.2August 20141.9November 2017107.6520.2August 20141.9November 2017107.6
2019–2039502.3July 20205.7December 2020109.9580.8July 20314.9November 2035115.9
2040–2059702.7July 20524.8February 2042124.4657.0July 20405.8December 2058118.7
2060–2079824.1July 20776.6February 2079128.2549.1July 20613.5April 2072122.9
2080–2100980.3July 21004.8February 2098117.9790.2July 20924.5December 2089132.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yun, S.-M.; Jeon, H.-T.; Cheong, J.-Y.; Kim, J.; Hamm, S.-Y. Combined Analysis of Net Groundwater Recharge Using Water Budget and Climate Change Scenarios. Water 2023, 15, 571. https://doi.org/10.3390/w15030571

AMA Style

Yun S-M, Jeon H-T, Cheong J-Y, Kim J, Hamm S-Y. Combined Analysis of Net Groundwater Recharge Using Water Budget and Climate Change Scenarios. Water. 2023; 15(3):571. https://doi.org/10.3390/w15030571

Chicago/Turabian Style

Yun, Sul-Min, Hang-Tak Jeon, Jae-Yeol Cheong, Jinsoo Kim, and Se-Yeong Hamm. 2023. "Combined Analysis of Net Groundwater Recharge Using Water Budget and Climate Change Scenarios" Water 15, no. 3: 571. https://doi.org/10.3390/w15030571

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