Evaluation of Seasonal Groundwater Quality Changes Associated with Groundwater Pumping and Level Fluctuations in an Agricultural Area, Korea

: This study was conducted to evaluate seasonal groundwater quality due to groundwater pumping and hydrochemical characteristics with groundwater level ﬂuctuations in an agricultural area in Korea. Groundwater levels were observed for about one year using automatic monitoring sensors, and groundwater uses were estimated based on the monitoring data. Groundwater use in the area is closely related to irrigation for rice farming, and rising groundwater levels occur during the pumping, which may be caused by the irrigation water of rice paddies. Hydrochemical analysis results for two separate times (17 July and 1 October 2019) show that the dissolved components in groundwater decreased overall due to dilution, especially at wells in the alluvial aquifer and shallow depth. More than 50% of the samples were classiﬁed as CaHCO 3 water type, and changes in water type occurred depending on the well location. Water quality changes were small at most wells, but changes at some wells were evident. In addition, the groundwater quality was conﬁrmed to have the effect of saltwater supplied during the 2018 drought by comparison with seawater. According to principal component analysis (PCA), the water quality from July to October was conﬁrmed to have changed due to dilution, and the effect was strong at shallow wells. In the study areas where rice paddy farming is active in summer, irrigation water may be one of the important factors changing the groundwater quality. These results provide a qualitative and quantitative basis for groundwater quality change in agricultural areas, particularly rice paddies areas, along with groundwater level and usage.


Introduction
In Korea, precipitation is concentrated in summer and varies greatly from season to season.In areas where there is a lot of rice paddy farming, irrigation water in summer is in high demand.If it rains less during that time, irrigation water should be drawn from other areas or secured from groundwater by developing wells in the area.When groundwater is pumped, its quality before and after pumping may change with the seasons.As groundwater pumping continues, the quality may change, which is caused by mixing from different aquifers, inflow of surface water from rivers or reservoirs [1][2][3][4], seawater intrusion [5][6][7], and so on.Also, groundwater quality may change temporally due to hydrologic processes, including seasonally varying rainfall and irrigation patterns [8][9][10][11].
Water 2021, 13 Hydrochemical facies in unconfined aquifers vary between young and old groundwater along the groundwater flow path, and the concentration of dissolved components, total dissolved solids (TDS), and electrical conductivity (EC) are related to the groundwater quality and residence time [1,12].In an agricultural region, nitrate in groundwater is decreased during dry/growing seasons of spring and fall due to the mixing of groundwater with deeper groundwater by heavy pumping, whereas it increases due to infiltration of water through the soil zone with highly concentrated nitrate [3].In addition, the chemical composition and groundwater flow patterns are changed by land use and groundwater use [2].In agricultural riparian zones, where intense groundwater pumping frequently occurs for the warming of greenhouses during the winter season, groundwater drawdown results in losing stream conditions and groundwater temperature decreases with considerable spatial variability, followed by hydrochemical and microbial community changes in groundwater [4,8].
In coastal aquifers, seawater intrusion by groundwater pumping has been reported in many studies [5], and the studies have been done to find the processes and factors associated with seawater intrusion [6].However, it is sometimes tough to distinguish the effect of pumping in these coastal aquifers.The aquifer in Fukuoka, Japan, was originally below the current seafloor, and deeper saline groundwater was isolated and not mixed with the shallower fresh groundwater, then the increased salinity by pumping did not result in seawater intrusion, but high-saline old groundwater was introduced [13].
Many cases in the literature show that groundwater quality changes with the seasons [14].Although many papers have been published on the change in groundwater quality, there are not many that have been dealt with the change in groundwater quality based on level fluctuations caused by groundwater pumping.This study was conducted to observe seasonal groundwater quality due to pumping and to evaluate hydrochemical characteristics with groundwater level variations in an agricultural area where groundwater is used intensively at specific times.Particularly, rice paddies in Korea require a lot of water in summer, and if it rains less in that time, groundwater must be used more for irrigation.Therefore, we wanted to evaluate the groundwater quality in different seasons based on the long-term monitoring for groundwater levels and hydrochemical analysis.

Study Area
The study area is located about 4 km from the west coast in the mid-western region of the Korean Peninsula at longitude 126 • 29 30"-126 • 31 10" E and latitude 36 • 32 43"-36 • 33 45" N; it is 2.83 km 2 (283.3 ha) of land that includes Yanggokri and a part of Sinri in Hongseong County, Chungnam Province.A regional water supply system has been established from the water source several kilometers away, so there is little problem with the supply for domestic use.However, agricultural water is not enough, particularly under drought conditions.Although a small reservoir is nearby, it is hard to get irrigation water during severe drought and groundwater has to be supplied or drawn from other regions.
The area consists of hilly mountains with an altitude of fewer than 100 m and relatively flat farming areas.The geology consists of Precambrian schist, gneiss, and Quaternary sediments [15,16].The upper layer is an alluvial layer mainly composed of porous media, and the lower layer is composed of bedrock.Therefore, the upper porous layer constitutes an alluvial aquifer, and the lower fractured bedrock constitutes a bedrock aquifer.
For agricultural land use, rice farming areas (paddy fields) account for 17.5% with 49.5 ha, and upland farming areas (upland fields) account for 13.3% with 37.7 ha.The area classified as forest or grassland is about 189.0 ha, accounting for 66.7% of the total area.In upland farming areas, various crops such as red pepper, sesame, soybean, pumpkin, and green onion have been cultivated (Figure 1).As of December 2018, the number of people living in the area was estimated to be about 120 based on the Annual Statistics of Hongseong County (http://www.hongseong.go.kr/).The total water supply rate of the whole county is 92.5%, the daily water supply is 255 L/day/person on average, and the amount of water supply of the area is a little more than the average.There are 107 wells in the area, 71 of which are for agricultural use and 36 for domestic use, based on data provided by the county administration.Annual groundwater use is reported to be 1846 m 3 /day for agricultural use and 274 m 3 /day for domestic use, indicating that most of the water is used for agriculture.However, it is difficult to trust these figures of groundwater use, because it is not actually monitored.Moreover, since seasonal use is unknown, the groundwater quality associated with pumping cannot be evaluated.
In 2018, due to the prolonged summer drought after the rainy season, irrigation water became scarce for rice farming.The local government urgently decided to supply water for irrigation from Hongseongho Lake.However, it is an artificial lake in the coastal area, storing water from a stream that flows into the sea, which is located about 4 km from the study area.At that time, Hongseongho Lake water contained a lot of salt, so the supplied water quality had to be checked regularly to make sure it did not exceed the salinity management criterion (<1280 ppm).From 8 August, Hongseongho Lake was supplied through a waterway.On 13 August, the salinity of the supplied water at the site was checked as 1320 ppm, and the water supply was stopped.This situation was repeated.Fortunately, at the end of August, about 200 mm of precipitation fell, and the drought situation was over.

Monitoring and Estimation of Seasonal Groundwater Use
To evaluate seasonal changes in groundwater quality, 55 wells were identified onsite in the study area, and 25 wells were finally selected for investigation.Table 1 shows the coordinates, elevation, usage, depth, and pump permit of the investigation wells, and automatic monitoring and hydrochemical analysis wells are marked.It is difficult to determine the exact structure of each well, but it is assumed that the screen of each well was installed according to main aquifer as shown in Table 1.As of December 2018, the number of people living in the area was estimated to be about 120 based on the Annual Statistics of Hongseong County (http://www.hongseong.go.kr/).The total water supply rate of the whole county is 92.5%, the daily water supply is 255 L/day/person on average, and the amount of water supply of the area is a little more than the average.There are 107 wells in the area, 71 of which are for agricultural use and 36 for domestic use, based on data provided by the county administration.Annual groundwater use is reported to be 1846 m 3 /day for agricultural use and 274 m 3 /day for domestic use, indicating that most of the water is used for agriculture.However, it is difficult to trust these figures of groundwater use, because it is not actually monitored.Moreover, since seasonal use is unknown, the groundwater quality associated with pumping cannot be evaluated.
In 2018, due to the prolonged summer drought after the rainy season, irrigation water became scarce for rice farming.The local government urgently decided to supply water for irrigation from Hongseongho Lake.However, it is an artificial lake in the coastal area, storing water from a stream that flows into the sea, which is located about 4 km from the study area.At that time, Hongseongho Lake water contained a lot of salt, so the supplied water quality had to be checked regularly to make sure it did not exceed the salinity management criterion (<1280 ppm).From 8 August, Hongseongho Lake was supplied through a waterway.On 13 August, the salinity of the supplied water at the site was checked as 1320 ppm, and the water supply was stopped.This situation was repeated.Fortunately, at the end of August, about 200 mm of precipitation fell, and the drought situation was over.

Monitoring and Estimation of Seasonal Groundwater Use
To evaluate seasonal changes in groundwater quality, 55 wells were identified on-site in the study area, and 25 wells were finally selected for investigation.Table 1 shows the coordinates, elevation, usage, depth, and pump permit of the investigation wells, and automatic monitoring and hydrochemical analysis wells are marked.It is difficult to determine the exact structure of each well, but it is assumed that the screen of each well was installed according to main aquifer as shown in Table 1.Groundwater level fluctuations are mainly caused by factors such as precipitation, groundwater pumping, tides, and so on.In the study area, long-term observations of groundwater levels are necessary because pumping for irrigation frequently happens.Automatic water level and quality monitoring sensors (TD-Diver, Van Essen Instruments, Delft, The Netherlands) were installed at 20 monitoring wells to assess the groundwater level fluctuations and pumping influence around the wells.Groundwater level monitoring was done automatically at 10 min intervals for approximately one year from September 2018 to September 2019.The automatic monitoring data were corrected by manually measured groundwater levels once a month.
The amount of groundwater use was estimated from the groundwater level data recorded by the automatic water level sensors.During a single pumping event, the pumping rate was checked, and it was multiplied by the pumping duration based on the groundwater level data.Estimated groundwater use is not exact because of the varying pumping rate during pumping and the monitoring interval, which is 10 min in this study.However, these errors could be ignored because groundwater uses in different seasons, and a long-term context must be evaluated rather than individual pumping events.

Hydrochemical Analysis
In order to evaluate the hydrochemical characteristics and seasonal changes of the groundwater in the study area, groundwater sampling from 21 wells in two different seasons, summer (17 July 2019) and fall (1 October 2019) was conducted, and measurement of groundwater levels was conducted at the same time to understand their spatial distribution in each time.The sampled groundwater was transported to the laboratory in a box with ice.Field parameters for groundwater samples, including temperature (T), pH, dissolved oxygen (DO), electrical conductivity (EC), and redox potential (Eh), were measured in the field using a digital meter for Intelligent Digital Sensors (IDS) (Multi3620 IDS, Xylem Analytics GmbH, Weilheim, Germany).The collected samples were filtered through a 0.2 µm membrane filter for instrumental analysis.

Multivariate Statistical Analysis
Principal component analysis (PCA) was used to identify major geochemical principal components (PCs) controlling groundwater chemistry.The PCs were extracted from 14 hydrochemical variables, including field parameters and dissolved components in the groundwater sample (T, pH, DO, EC, Ca 2+ , Mg 2+ , Na + , K + , SiO 2 (aq), HCO  2. The graphs in Figure 2 show groundwater level fluctuations of the main pumping wells and their surrounding observation wells by regions.Figure 2a shows the pumping wells are D-08 and D-11, and D-09 and D-11 in Figure 2b.In Figure 2c, there is no pumping well, and in Figure 2d, D-03 and D-05 are pumping wells.Except for the pumping wells, the wells shown in Figure 2 are observation wells.Precipitation is also shown in Figure 2c.Although the recovery of groundwater levels by precipitation was observed locally at some wells, fluctuations of levels showed a significant effect due to pumping.
Figure 2a shows fluctuations in the groundwater levels of D-08 and D-11 pumping wells and D-01, D-02, and C-09 observation wells.D-08, D-11, D-01, and D-02 wells are installed in the bedrock aquifer, and C-09 is installed in the alluvial aquifer.D-08 is the main supply well for irrigation in the area.Based on the groundwater level fluctuation data, pumping increased from mid-April and continued for about two months until mid-July.From 22 May to 16 July, the groundwater level drawdown was lower than the depth of the sensor installed in the well, and groundwater levels were recorded at the sensor location, about -6.0 m.Therefore, it is believed that the water levels actually dropped to a lower level than this.The effect of pumping from D-08 also led to a drop in groundwater level at nearby observation wells D-01 and D-02.At D-01, the level continuously lowered from 12.62 m on 12 May to 0.15 m on 7 June, showing a 12.47 m drop.At D-02, the level went from 12.69 m to 4.20 m at the same time as D-01, showing an 8.49 m drop.Based on the groundwater level data, D-01 and D-02 were affected by pumping from D-08, and they seemed to be closely related to each other.Groundwater levels from another pumping well, D-11, showed rapid drop and recovery due to pump and stop.Here, groundwater level data below the sensor installation depth, about −5.0 m, were missing, but a correlation of levels at D-11, D-01, and D-02 were not evident.C-09 is an observation well installed in the alluvial aquifer, unlike the other wells in the bedrock aquifer.The groundwater levels at C-09 fluctuated between 13.68 and 15.60 m in mean sea level elevation, showing a 1.2 m drop from 14.9 m on 12 May to 13.7 m on 7 June.C-09 is the closest observation well to D-08, the main pumping well in Figure 3a, but showed a different fluctuation pattern compared to D-01 and D-02, indicating poor connectivity between bedrock aquifer and alluvial aquifer.It is observed that groundwater levels rose as of 8 June despite continuous pumping from mid-April to early August, as shown in Figure 2a.Even though there was only a little rain during that time, as seen by the precipitation data shown in Figure 2c, the rise in groundwater levels must be attributed to other factors.The most likely factor is the recharge of the irrigation water in rice paddies.Water supplied to rice paddies from the pumping goes back into the ground and recharges the groundwater in the aquifer with a   Table 2 shows the estimation results of monthly groundwater use.The total amount of groundwater of these five wells was 54625 m 3 /year.D-08 pumped the most groundwater for one year, 26417 m 3 .The maximum monthly groundwater use was estimated at D-03, 8648 m 3 in May 2019.In summer 2019 (June to August), the amount of precipitation was 172 mm, which was only about 27% of the normal amount of the season.Therefore, groundwater use during that time could have been large compared with other years.It is observed that groundwater levels rose as of 8 June despite continuous pumping from mid-April to early August, as shown in Figure 2a.Even though there was only a little rain during that time, as seen by the precipitation data shown in Figure 2c, the rise in groundwater levels must be attributed to other factors.The most likely factor is the recharge of the irrigation water in rice paddies.Water supplied to rice paddies from the pumping goes back into the ground and recharges the groundwater in the aquifer with a time delay.The groundwater level increase as of mid-June can be seen in Figure 2a and the other graphs in Figure 2.Moreover, the groundwater level rise is observed in wells installed in both the alluvial aquifer and bedrock aquifer.The effect of recharge on paddy fields into aquifers can be found in various studies [18][19][20].Groundwater pumping frequently occurred from mid-April to late June, and the levels at D-10 were lowered to about −22.0 m from 9 May to 19 June.Due to pumping from D-09 and D-10, the groundwater levels at D-12 went down to about 1.5 m from May to mid-June.D-well is installed in bedrock and D-12 in the alluvial aquifer, but the groundwater level fluctuations between the two are correlated.Similarly, the levels at C-07 in the alluvial aquifer were from 5.64 to 7.98 m, and the effect of pumping from D-10 was evident at the well.Groundwater levels at D-10 were lowered by about 1.68 m from early May to mid-June.
Figure 2c shows the groundwater level fluctuations at C-02, C-03, C-04, C-05, C-06, and C-08 observation wells.The effects of both precipitation and pumping from nearby wells on these wells are illustrated.C-08 among the wells in Figure 2c is located at the highest altitude and the groundwater levels are the highest.The levels at C-05 reflected the effect of nearby groundwater pumping.The distance between C-05 and the nearby pumping well is about 25 m, and they are thought to be connected to each other.The groundwater level fluctuations at C-04 were similar to those at C-05, reflecting the pumping effect.
Figure 2d shows the groundwater level fluctuations at D-03, D-05, and D-13 pumping wells and D-14 and C-01 observation wells.Pumping from D-03 for irrigation of paddy fields was very often from May to July, and groundwater was also pumped from November to March, which was mainly for upland field irrigation.Groundwater levels were not recorded below the sensor installation depth, and the maximum drawdown of levels at D-03 is estimated to be over 54.2 m.The groundwater level at D-05 was affected by pumping from 22 May to 16 June, and was constant at about 21.3 m.At C-01, the groundwater level dropped by about 2 m due to pumping from D-05, and was at about 18.3 m.The level at D-14 was the highest at about 36.2 to 36.5 m with little fluctuation.The levels at D-13 fluctuated due to the individual pumping.

Spatial Distribution and Factors for Groundwater Level Changes
Figure 3 shows the spatial distribution of groundwater levels in two seasons, the first survey in July 2019 and the second in October 2019.The levels were only described in plain areas because there are no observation wells in forest areas.The depths of groundwater levels were 0.69 to 16.27 m (average 5.07 m) below the land surface in July, and 0.10 to 7.91 m (average 2.90 m) in October.The groundwater levels in mean sea level elevation were 2.38 to 32.75 m (average 15.18 m) in July and 7.02 to 32.80 m (average 15.71 m) in October.The groundwater levels in October were about 53 cm higher, on average, than in July.It is certain that the pumping influence in the eastern part of the area was reduced in October.In the southern part, the groundwater levels in October were higher than in July.Although the levels around D-09 and D-10 were still down in October, there was a slight rise by about 1.2 m compared with those in July.
Groundwater use in the area is closely related to irrigation for rice farming.Due to the large supply of irrigation water for rice paddies in summer, groundwater levels decrease by pumping in July.However, irrigation water in October, when rice is harvested, is rarely required.Precipitation is another factor that affects changes in groundwater levels in the area.Precipitation, which was 37.5 mm in July and 112 mm in September (the month prior to the second survey on 1 October 2019), may also have caused groundwater levels to increase.Additionally, irrigation water for rice paddies can flow back into the ground, raising the groundwater levels.Monitoring data also suggest that the irrigation water in rice paddies entered the aquifer to increase groundwater levels.

Monthly Groundwater Use
Monthly groundwater use was estimated based on the fluctuation data at D-03, D-05, D-08, D-09, and D-10 wells, which are major agricultural wells in the area.The groundwater use for a certain period was calculated by multiplying the pumping rate by pumping duration.The pumping rate of D-03, D-05, D-08, D-09, and D-10 wells for a single pumping event is 242, 90, 248, 95, and 174 m 3 /day, respectively.
Table 2 shows the estimation results of monthly groundwater use.The total amount of groundwater of these five wells was 54625 m 3 /year.D-08 pumped the most groundwater for one year, 26417 m 3 .The maximum monthly groundwater use was estimated at D-03, 8648 m 3 in May 2019.In summer 2019 (June to August), the amount of precipitation was 172 mm, which was only about 27% of the normal amount of the season.Therefore, groundwater use during that time could have been large compared with other years.There are a total of 25 wells for agricultural use in the area and the pump permit of those wells is 1918 m 3 /day, which is 700,070 m 3 /year.However, the actual amount used for one year is reported as 28,414 m 3 , so this amount is only about 4% of the pump permit.In this study, 54,625 m 3 from five wells (D-03, D-05, D-08, D-09, and D-10) was used over one year, from October 2018 to September 2019, equivalent to 192% (almost twice as much) of the total reported groundwater use.The estimated amount of groundwater use of the five wells is around 21.2% of the pump permit, or 705 m 3 /day (257,325 m 3 /year).It can be deduced that actual use was quite different from the reported amount, and only a small portion of the pump permit was used.

Hydrochemical Analysis and Groundwater Quality Evaluation 4.3.1. Groundwater Quality Changes and Water Types
The main causes of water quality change in the study area are various such as precipitation, groundwater pumping, and pollution from agricultural activities.In particular, when the drought in the summer of 2018 was severe, the temporary supply of salty water from Hongseongho Lake as irrigation water is also thought to affect water quality changes.
The hydrochemical analysis was properly conducted with a charge balance error (CBE) less than 5.5%, and the results are shown in Table 3 for the first survey on 17 July 2019 and the second on 1 October 2019.Total dissolved solids (TDS) amounted to 81-1409 mg/L (average 339 mg/L) in July and 78-709 mg/L (average 259 mg/L) in October, showing slightly lower TDS in groundwater in October than in July.Concentrations of major ions, both cations and anions, were also lower in October than in July.In addition, electrical conductivity (EC) was measured between 160-2630 µS/cm (average 598 µS/cm) in July and 153-1361 µS/cm (average 430 µS/cm) in October.4 show groundwater quality changes in the different seasons based on piper diagrams.More than 50% of the samples were classified as CaHCO 3 (II) water type, and water type changes occurred depending on the well location.Water quality changes were small at most wells, but changes at some wells were evident.Comparing the groundwater quality in July and October, the water type at B-11 and D-04 was changed from CaHCO 3 (II) to mixed CaMgCl (III); at C-07, D-09, and D-10 from mixed CaMgCl (III) to CaHCO 3 (II); and at D-01 from mixed CaNaHCO 3 (IV) to NaCl (IV).It seems certain that the groundwater in October was more diluted than in July.Such changes in groundwater quality may have been caused by precipitation and recharge of surface water, such as irrigation water of rice paddies.When the groundwater level is lowered by pumping, irrigation water may flow more into aquifers.Finally, intensive pumping in the area may have induced the irrigation water into aquifers, changing the groundwater quality.The flow of irrigation water into aquifers was reflected in the groundwater level rise despite continuous pumping from mid-June to early August.

Evaluation of Characteristics of Groundwater Quality
In order to see whether there were any seasonal changes in groundwater quality, the It seems certain that the groundwater in October was more diluted than in July.Such changes in groundwater quality may have been caused by precipitation and recharge of surface water, such as irrigation water of rice paddies.When the groundwater level is lowered by pumping, irrigation water may flow more into aquifers.Finally, intensive pumping in the area may have induced the irrigation water into aquifers, changing the groundwater quality.The flow of irrigation water into aquifers was reflected in the groundwater level rise despite continuous pumping from mid-June to early August.

Evaluation of Characteristics of Groundwater Quality
In order to see whether there were any seasonal changes in groundwater quality, the wells were categorized into two groups, group C and group B/D, according to the first letter of the well label, and box plots were constructed, shown in Figure 5. Group C wells are all installed at shallow depths in the alluvial aquifer, except for C-06 and C-07, and group B/D wells are installed deep in the bedrock.Although C-06 and C-07 are installed at about 50 m depth, it is reasonable to assume that they are affected by the alluvial layer because well screens penetrate both the upper alluvial aquifer and the lower bedrock.Most of the hydrochemical parameter values of group C are higher (temperature, TDS, EC, Ca 2+ , Mg 2+ , Na + , K + , HCO 3 − , Cl − , Br − , and SO 4 2− ) than those of the relatively deep wells of group D/B.pH, Eh, and SiO 2 (aq) values are relatively higher in group D/B than group C. DO and NO 3 − in group B/D show large variations.NO 3 − concentration, as the indicator of contamination from agricultural activity or livestock waste, is lower than the drinking water quality standards of Korea (NO 3 -N: 10 mg/L; NO 3 − : 44.3 mg/L), but the groundwater samples of C-08, B-20, and B-11 in both the first and second survey exceeded the standards.The NO 3 − concentration at B-20 and B-11 in the first survey was 48.1 mg/L and 41.1 mg/L, and in the second survey was 55.3 mg/L and 55.2 mg/L, respectively.They showed a higher concentration in October than in July.On the other hand, the groundwater sample at C-08 showed 60.5 mg/L in the first survey and 40.7 mg/L in the second survey, and the NO 3 − concentration was reduced from July to October.The seasonal variation of groundwater quality in group C was very clearly compared with group B/D.In particular, the parameters indicating dissolved components in groundwater (TDS, EC, Ca 2+ , Mg 2+ , Na + , K + , HCO 3 − , Cl − , Br − , and SO 4 2− ) were lower in October than in July.It may have been caused by precipitation or surface water infiltration, and the irrigated water used for rice paddies may be an important factor in reducing dissolved components in groundwater in the area.However, group B/D does not have large seasonal variation compared with group C. The main difference between the groups is whether the main aquifer for well development is in the upper alluvium (soil) or the lower bedrock.Group C wells are installed in the alluvial layer and are highly affected by agricultural activity, and it seems that a lot of matter in the soil might be dissolved in the groundwater.In addition, group C wells are developed at the upper layer, where the inflow of rainfall and surface water can occur more easily than at group B/D wells, and then the dilution effect can be more clearly seen.
Figure 6 shows the relationships between Cl − and Na + , and Cl − and Br − in each well to evaluate the characteristics and sources of groundwater.The seawater data were extracted from the literature [21].The groundwater samples are clearly lined up with groups C and B/D to the seawater as the endpoint in the graph of Cl − versus Na + and Br − concentration.Group C wells are lined up between group B/D and seawater.This suggests that the groundwater quality is affected by saltwater more in group C than in group B/D.As such, the possible seawater effects on groundwater occur due to the fact that the saline water was supplied from Hongseongho Lake during the summer drought in 2018.The effects are shown particularly in alluvial groundwater at shallow depths.The group B/D samples were farther from the seawater than group C samples, and the effect of seawater in deep wells seems to be relatively small.Comparing the groundwater quality in July and October, the groundwaters tend to move farther away from seawater, especially in Cl − and Na + in October.Over time, the influence of seawater seems to have decreased.
Figure 7 shows the analysis results of oxygen and hydrogen isotopes to identify the source of groundwater.The water samples were from the first survey, and no analysis was conducted for the second survey samples.The groundwater samples were close to the local meteoric water line (LMWL) [22] during the summer, the rainy season, and the groundwater seems to have been mostly recharged in summer.The samples of group C and some group B/D wells reflect the evaporation effect away from the LMWL.Particularly, the evaporation effect was the largest at C-03.Samples showing the evaporation effect were from shallow depths, 0.74-4.93m from the surface, and had relatively high temperatures (15.9-21.7 • C).The presence of evaporation effects at shallow depths means that external factors may easily cause groundwater quality to change.The wells in the elliptical region shown in Figure 7 have a large seasonal change in groundwater quality.Figure 7 shows the analysis results of oxygen and hydrogen isotopes to identify the source of groundwater.The water samples were from the first survey, and no analysis was conducted for the second survey samples.The groundwater samples were close to the local meteoric water line (LMWL) [22] during the summer, the rainy season, and the groundwater seems to have been mostly recharged in summer.The samples of group C and some group B/D wells reflect the evaporation effect away from the LMWL.Particularly, the evaporation effect was the largest at C-03.Samples showing the evaporation effect were from shallow depths, 0.74-4.93m from the surface, and had relatively high temperatures (15.9-21.7 °C).The presence of evaporation effects at shallow depths means that external factors may easily cause groundwater quality to change.The wells in the elliptical region shown in Figure 7 have a large seasonal change in groundwater quality.Table 5 shows a correlation matrix of the 14 hydrochemical parameters in groundwater samples.EC shows a high correlation with major dissolved ions such as Ca 2+ , Mg 2+ , Na + , K + , Cl − , SO 4 2− , etc., and DO is highly correlated to EC, Na + , K + , Cl − , and Br − .In particular, Br − , Cl − , Mg 2+ , and Na + , which indicate an association with seawater, are highly correlated.Table 5 shows a correlation matrix of the 14 hydrochemical parameters in groundwater samples.EC shows a high correlation with major dissolved ions such as Ca 2+ , Mg 2+ , Na + , K + , Cl − , SO4 2− , etc., and DO is highly correlated to EC, Na + , K + , Cl − , and Br − .In particular, Br − , Cl − , Mg 2+ , and Na + , which indicate an association with seawater, are highly correlated.8 show the PCA results of groundwater samples from both surveys.PCA results for all samples, including the first and second surveys, show that the three major eigenvalues are 56.34% for PC1, 13.98% for PC2, and 11.19% for PC3, explaining 81.51% of the total variance.Among the 14 hydrochemical parameters, PC1 confirms that DO, EC, Ca 2+ , Mg 2+ , Na + , K + , Cl − , SO4 2− , and Br − are the main contributors, while PC2    8 show the PCA results of groundwater samples from both surveys.PCA results for all samples, including the first and second surveys, show that the three major eigenvalues are 56.34% for PC1, 13.98% for PC2, and 11.19% for PC3, explaining 81.51% of the total variance.Among the 14 hydrochemical parameters, PC1 confirms that DO, EC, Ca 2+ , Mg 2+ , Na + , K + , Cl − , SO 4 2− , and Br − are the main contributors, while PC2 is affected by pH, SiO 2 (aq), and HCO 3 -, and PC3 is affected by NO 3 − .PC1 may be considered as the amount of main dissolved components in groundwater (Table 6).Figure 8a shows PC1 versus PC2 of the samples.It clearly distinguished from groups C and B/D, and the groundwater quality of C-03 shows the dramatic change from July to October.Table 7 shows the variation of principal components (PCs) during seasonal changes.The water quality of most samples changed in the direction that the value of PC1 was reduced from July to October.A reduced PC1 value is thought to indicate dilution.In addition, the water quality of most samples, except the C-03 sample, changed upward, as shown in Figure 8a.The upward direction for the increased value of PC2 means that pH, HCO 3 − , and SiO 2 (aq) concentrations increased from July to October, and this might be indicated by the degree of water-rock or water movement and mixing in the groundwater flow system.Figure 8b shows PC1 versus PC3 of the samples.PC3 values also tended to decrease as PC1 values decreased from July to October.The parameter representing PC3 is NO 3 , and it decreased together with the reduced PC1.
According to PCA, the water quality from July to October can change due to the dilution, and the effect is strong in group C. The period between July and October is when there is much precipitation.Therefore, precipitation may recharge much into the aquifer, and the recharged water may dilute the groundwater.Additionally, a lot of groundwater is pumped for irrigating rice paddies, and the irrigation water can flow back into the aquifer and influence the groundwater to be diluted.In the study areas where rice farming is active in summer, groundwater quality changes seem to be largely caused by both precipitation and the irrigated water.
is affected by pH, SiO2 (aq), and HCO3 -, and PC3 is affected by NO3 − .PC1 may be considered as the amount of main dissolved components in groundwater (Table 6).Figure 9 is presented as a conceptual diagram in order to easily understand the factors affecting the groundwater quality in the study area.From the PCA results, dilution (PC1), water-rock interaction (PC2), and contamination (PC3) were derived as the factors that cause water quality changes in July and October.Dilution can occur by precipitation, surface water infiltration, and the irrigated water to rice paddies which is induced by pumping or naturally flowing back into the aquifer.In addition, water-rock interaction occurs in the aquifer under groundwater flow control.Contaminations by fertilizer and livestock waste also affect groundwater quality in the study area.High saline groundwaters at some wells are attributed to the irrigated water from the outer water source during the past severe drought.From the study, we confirm that the changes in groundwater quality have occurred, and they are closely related to irrigation during the rice-growing season.The groundwater in the study area is influenced by ① dilution by precipitation, surface water infiltration and the irrigated water to the rice paddy field, and ② water-rock interaction under groundwater flow control, and ③ contamination by fertilizer and livestock waste.Additionally, high saline groundwaters at some wells are attributed to seawater influence by water supplied from outer water source during the past severe drought in 2018.

Conclusions
This study was conducted in priority to find out whether changes in water quality due to intensive pumping were observed in rice farming areas.Furthermore, if such a change in water quality was observed, it was attempted to determine the cause of the occurrence.Through long-term monitoring of the groundwater levels, it was confirmed that intensive pumping occurred for a certain period (from May to August) to supply irrigation water for rice farming.Groundwater levels fluctuated due to frequent pumping The groundwater in the study area is influenced by 1 dilution by precipitation, surface water infiltration and the irrigated water to the rice paddy field, and 2 water-rock interaction under groundwater flow control, and 3 contamination by fertilizer and livestock waste.Additionally, high saline groundwaters at some wells are attributed to seawater influence by water supplied from outer water source during the past severe drought in 2018.

Conclusions
This study was conducted in priority to find out whether changes in water quality due to intensive pumping were observed in rice farming areas.Furthermore, if such a change in water quality was observed, it was attempted to determine the cause of the occurrence.Through long-term monitoring of the groundwater levels, it was confirmed that intensive pumping occurred for a certain period (from May to August) to supply irrigation water for rice farming.Groundwater levels fluctuated due to frequent pumping for irrigation of rice paddies in summer, and groundwater uses were estimated from major pumping wells in the area based on automatic groundwater level monitoring data from September 2018 to September 2019, about one year.Groundwater use in the area is closely related to irrigation for rice farming.Besides precipitation, irrigation water for rice paddies is another factor that increases groundwater levels, which is clearly shown in the rise in levels during pumping.
The groundwater quality analysis results showed a large change in water quality in some wells, but in some wells, the change was insignificant.Water quality types are classified using the major cations and anions dissolved in groundwater.It seems certain that the groundwater in October is more diluted than in July.The seasonal variation of groundwater quality in group C of alluvial wells at shallow depth was very clearly compared with those in group B/D deep wells.The parameter values indicating dissolved components in groundwater (TDS, EC, Ca 2+ , Mg 2+ , Na + , K + , HCO 3 − , Cl − , Br − , SO 4 2− ) are lower in October than in July.Group B/D does not have large seasonal variation.In the isotope analysis results, it was found that the evaporation effect appeared in shallow wells and thus was greatly influenced by external factors.
The PCA results show that the three major eigenvalues were 56.34% for PC1, 13.98% for PC2, and 11.19% for PC3, explaining 81.51% of the total variance.PC1 can be related to the amount of main dissolved components in groundwater.The water quality of most samples changed in the direction of the PC1 value reduction.The PC2 may indicate the degree of water-rock or water-aquifer interaction, and the PC3, which is referred to NO 3 − concentration, may be related to contamination.Another factor influencing groundwater quality in the study area is the saline water irrigation from other regions, which was supplied temporarily during the summer drought in 2018.It was evaluated from the comparisons with the Na + , Cl − , and Br -components of seawater and groundwater.
This study is meaningful because it quantitatively evaluated the change of groundwater quality in relation with irrigation water by groundwater pumping, precipitation, and other possible factors (saline water supply during the summer drought in the previous year in this study) in rice farming areas for the first time.

Figure 1 .
Figure 1.Study area and well locations.

Figure 1 .
Figure 1.Study area and well locations.

Figure 2a shows fluctuations
Figure 2a shows fluctuations in the groundwater levels of D-08 and D-11 pumping wells and D-01, D-02, and C-09 observation wells.D-08, D-11, D-01, and D-02 wells are installed in the bedrock aquifer, and C-09 is installed in the alluvial aquifer.D-08 is the main supply well for irrigation in the area.Based on the groundwater level fluctuation data, pumping increased from mid-April and continued for about two months until mid-July.From 22 May to 16 July, the groundwater level drawdown was lower than the depth of the sensor installed in the well, and groundwater levels were recorded at the sensor location, about -6.0 m.Therefore, it is believed that the water levels actually dropped to a lower level than this.The effect of pumping from D-08 also led to a drop in groundwater level at nearby observation wells D-01 and D-02.At D-01, the level continuously lowered from 12.62 m on 12 May to 0.15 m on 7 June, showing a 12.47 m drop.At D-02, the level went from 12.69 m to 4.20 m at the same time as D-01, showing an 8.49 m drop.Based on the groundwater level data, D-01 and D-02 were affected by pumping from D-08, and they seemed to be closely related to each other.Groundwater levels from another pumping well, D-11, showed rapid drop and recovery due to pump and stop.Here, groundwater level data below the sensor installation depth, about −5.0 m, were missing, but a correlation of levels at D-11, D-01, and D-02 were not evident.C-09 is an observation well installed in the alluvial aquifer, unlike the other wells in the bedrock aquifer.The groundwater levels at C-09 fluctuated between 13.68 and 15.60 m in mean sea level elevation, showing a 1.2 m drop from 14.9 m on 12 May to 13.7 m on 7 June.C-09 is the closest observation well to D-08, the main pumping well in Figure 3a, but showed a different fluctuation pattern compared to D-01 and D-02, indicating poor connectivity between bedrock aquifer and alluvial aquifer.It is observed that groundwater levels rose as of 8 June despite continuous pumping from mid-April to early August, as shown in Figure2a.Even though there was only a little rain during that time, as seen by the precipitation data shown in Figure2c, the rise in groundwater levels must be attributed to other factors.The most likely factor is the recharge of the irrigation water in rice paddies.Water supplied to rice paddies from the pumping goes back into the ground and recharges the groundwater in the aquifer with a

Figure 3 .
Figure 3. Spatial distribution of groundwater levels based on (a) the first survey in July 2019 and (b) the second survey in October 2019.Groundwater levels are expressed as the altitude above mean sea level.

4. 2 .
Monthly Groundwater Use Monthly groundwater use was estimated based on the fluctuation data at D-03, D-05, D-08, D-09, and D-10 wells, which are major agricultural wells in the area.The groundwater use for a certain period was calculated by multiplying the pumping rate by pumping duration.The pumping rate of D-03, D-05, D-08, D-09, and D-10 wells for a single pumping event is 242, 90, 248, 95, and 174 m 3 /day, respectively.

Figure 3 .
Figure 3. Spatial distribution of groundwater levels based on (a) the first survey in July 2019 and (b) the second survey in October 2019.Groundwater levels are expressed as the altitude above mean sea level.

Figure
Figure 2b shows the groundwater level fluctuations at D-09 and D-10 pumping wells and D-12 and C-07 observation wells.Rapid groundwater drawdown and recovery at D-09 and D-10 are shown, which was caused by the pump and stop.Groundwater levels below the sensor installation depth of −21.75 m at D-09 and −12.33 m at D-10 were not observed.A groundwater level drop of about 33.2 m occurred at D-09, and 22.64 m at D-10.Groundwater pumping frequently occurred from mid-April to late June, and the levels at D-10 were lowered to about −22.0 m from 9 May to 19 June.Due to pumping from D-09 and D-10, the groundwater levels at D-12 went down to about 1.5 m from May to mid-June.D-well is installed in bedrock and D-12 in the alluvial aquifer, but the groundwater level fluctuations between the two are correlated.Similarly, the levels at C-07 in the alluvial aquifer were from 5.64 to 7.98 m, and the effect of pumping from D-10 was evident at the well.Groundwater levels at D-10 were lowered by about 1.68 m from early May to mid-June.Figure2cshows the groundwater level fluctuations at C-02, C-03, C-04, C-05, C-06, and C-08 observation wells.The effects of both precipitation and pumping from nearby wells on these wells are illustrated.C-08 among the wells in Figure2cis located at the highest altitude and the groundwater levels are the highest.The levels at C-05 reflected the effect of nearby groundwater pumping.The distance between C-05 and the nearby pumping well is about 25 m, and they are thought to be connected to each other.The groundwater level fluctuations at C-04 were similar to those at C-05, reflecting the pumping effect.Figure2dshows the groundwater level fluctuations at D-03, D-05, and D-13 pumping wells and D-14 and C-01 observation wells.Pumping from D-03 for irrigation of paddy fields was very often from May to July, and groundwater was also pumped from November to March, which was mainly for upland field irrigation.Groundwater levels were not recorded below the sensor installation depth, and the maximum drawdown of levels at D-03 is estimated to be over 54.2 m.The groundwater level at D-05 was affected by pumping from 22 May to 16 June, and was constant at about 21.3 m.At C-01, the groundwater level dropped by about 2 m due to pumping from D-05, and was at about 18.3 m.The level at D-14 was the highest at about 36.2 to 36.5 m with little fluctuation.The levels at D-13 fluctuated due to the individual pumping.

Figure 4 .
Figure 4. Piper diagrams indicating water type changes from the first survey in July 2019 to the second survey in October 2019.

Figure 4 .
Figure 4. Piper diagrams indicating water type changes from the first survey in July 2019 to the second survey in October 2019.

Figure 5 .
Figure 5. Box plots for groundwater quality changes indicating seasonal variations.X axes indicate group C and B/D of wells.Group C wells are installed mainly at shallow depths in the alluvial aquifer, and group B/D wells are installed deep in the bedrock.The first survey results in July 2019 of group C wells are marked as C1, and the second survey results in October 2019 as C2.Similarly, the first survey results of group B/D wells are marked as B/D1, and the second survey results as B/D2.Most of the hydrochemical parameter values of group C are higher (temperature, TDS, EC, Ca 2+ , Mg 2+ , Na + , K + , HCO3 − , Cl − , Br − , and SO4 2− ) than those of the relatively deep wells of group D/B.pH, Eh, and SiO2 (aq) values are relatively higher in group D/B than group C. DO and NO3 − in group B/D show large variations.NO3 − concentration, as the indicator of contamination from agricultural activity or livestock waste, is lower than the

Figure 5 .
Figure 5. Box plots for groundwater quality changes indicating seasonal variations.X axes indicate group C and B/D of wells.Group C wells are installed mainly at shallow depths in the alluvial aquifer, and group B/D wells are installed deep in the bedrock.The first survey results in July 2019 of group C wells are marked as C1, and the second survey results in October 2019 as C2.Similarly, the first survey results of group B/D wells are marked as B/D1, and the second survey results as B/D2.

Figure 6 .
Figure 6.Evaluation of characteristics and sources of groundwater quality by irrigation water influenced by seawater (SW): (a) in July 2019 and (b) in October 2019.Seawater data were modified from the literature [21].

Figure 6 .
Figure 6.Evaluation of characteristics and sources of groundwater quality by irrigation water influenced by seawater (SW): (a) in July 2019 and (b) in October 2019.Seawater data were modified from the literature [21].

Figure 7 .
Figure 7. Isotopic analysis results for groundwater samples from the first survey in July 2019.Dotted red line indicates evaporation line, and wells in elliptical region generally have large seasonal change in groundwater quality.

Figure 7 .
Figure 7. Isotopic analysis results for groundwater samples from the first survey in July 2019.Dotted red line indicates evaporation line, and wells in elliptical region generally have large seasonal change in groundwater quality.

Figure 8 .
Figure 8. Results of principal component analysis (PCA): (a) PC1 vs. PC2 and (b) PC1 vs. PC3.Arrow direction indicates water quality change from the first survey in July 2019 to the second survey in October 2019.Each well label represents well ID and sampled survey; for example, C1-01 means C-01 well sampled in the first survey, C2-01 sampled in the second survey.

Figure 8 .
Figure 8. Results of principal component analysis (PCA): (a) PC1 vs. PC2 and (b) PC1 vs. PC3.Arrow direction indicates water quality change from the first survey in July 2019 to the second survey in October 2019.Each well label represents well ID and sampled survey; for example, C1-01 means C-01 well sampled in the first survey, C2-01 sampled in the second survey.

Figure 9 .
Figure 9.The conceptual diagram to represent the groundwater quality changes during the pumping in the study area.The groundwater in the study area is influenced by ① dilution by precipitation, surface water infiltration and the irrigated water to the rice paddy field, and ② water-rock interaction under groundwater flow control, and ③ contamination by fertilizer and livestock waste.Additionally, high saline groundwaters at some wells are attributed to seawater influence by water supplied from outer water source during the past severe drought in 2018.

Figure 9 .
Figure 9.The conceptual diagram to represent the groundwater quality changes during the pumping in the study area.The groundwater in the study area is influenced by 1 dilution by precipitation, surface water infiltration and the irrigated water to the rice paddy field, and 2 water-rock interaction under groundwater flow control, and 3 contamination by fertilizer and livestock waste.Additionally, high saline groundwaters at some wells are attributed to seawater influence by water supplied from outer water source during the past severe drought in 2018.

Table 2 .
Monthly groundwater use of main pumping wells in the study area.(unit: m 3 ).

Table 2 .
Monthly groundwater use of main pumping wells in the study area.(unit: m 3 ).
-indicates no pumping recorded.

Table 3 .
Hydrochemical analysis results of groundwater.

Table 4 and
Figure

Table 5 .
Correlation matrix of hydrochemical parameters of groundwater samples.

Table 5 .
Correlation matrix of hydrochemical parameters of groundwater samples.

Table 6 and
Figure

Table 6 .
Principal component loadings, eigenvalues, and explained variance for five extracted components of groundwater samples.

Table 7 .
Variation of principal components (PCs) during seasonal change.Symbols and indicate increase and decrease of PCs, respectively.
Water 2021, 13, x FOR PEER REVIEW 19 of 21