==== Front Int J Environ Res Public Health Int J Environ Res Public Health ijerph International Journal of Environmental Research and Public Health 1661-7827 1660-4601 MDPI 10.3390/ijerph17239126 ijerph-17-09126 Article Groundwater Chemical Characteristics and Controlling Factors in a Region of Northern China with Intensive Human Activity Ren Chaobin 12 Zhang Qianqian 3* 1 School of Chemical Engineering, Zhengzhou University, Zhengzhou 450001, China; robin_ren@gs.zzu.edu.cn 2 School of Civil Engineering, Nanyang Institute of Technology, Nanyang 473004, China 3 Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang 050061, China * Correspondence: zhangqianqian@mail.cgs.gov.cn; Tel.: +86-18348837711 07 12 2020 12 2020 17 23 912614 10 2020 04 12 2020 © 2020 by the authors.2020Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).The identification of groundwater chemical characteristics and the controlling factors is of major concern in water environment studies. In this study, we identified the groundwater chemical characteristics, evolution laws and main control factors in a region severely affected by human activities using hydrochemical and multivariate statistical techniques. The results showed that the concentrations of NO3− and TH were the primary pollution factors in the region with intensive human activity because of high concentration and over the standard rates. The major types of groundwater chemistry were HCO3·SO4–Ca and HCO3·SO4–Ca·Mg. The sulfate-type water was as high as 75.0%, 69.2% and 41.2% in the three hydrogeological units. In addition, there were Cl-type and Na-type waters, indicating that the groundwater in this area has been significantly affected by human activities. A principal component analysis (PCA) indicated that the three factors affecting groundwater hydrochemistry in the study area are domestic sewage and fertilizer, water–rock interactions and industrial wastewater. Therefore, we suggest that the government and water environment management departments should prevent the discharge of domestic and industrial wastewater without standardized treatment first in order to effectively prevent the further deterioration of groundwater quality in this area. groundwaterhydrochemistryhuman activitycontrol factors ==== Body 1. Introduction Groundwater is one of the most vital resources for human living and production, especially in arid and semi-arid regions [1]. In recent years, with the rapid growth of China’s economy and population, processing capacity of sewage and garbage has not kept pace, and has led to unregulated discharge of sewage (especially in rural areas), the excessive application of pesticide and fertilizer, overexploitation of groundwater, and so on [2,3,4]. This has caused the deterioration of groundwater quality and changed groundwater chemical characteristics. It is well known that changes in the groundwater chemistry not only disrupt many ecological processes [5], but also affect groundwater use. Therefore, fully understanding the evolution law of groundwater chemistry is vital for the sustainable use of groundwater [6,7]. Previous studies found that the evolution of groundwater chemistry is mainly controlled by natural factors (including hydrogeological conditions, aquifer lithology, interactions of groundwater with rock, etc.) and human activities (such as domestic sewage, agricultural fertilizer and overexploitation of groundwater) [8,9]. In recent years, researchers have devoted to study the effect factors on groundwater chemistry [10,11,12]. Huang et al. [8] found that the main factors that control groundwater quality were seawater intrusion, water–rock interactions, lateral flow of river water, sewage discharge and agricultural pollution. Gu et al. [13] identified that the factors influencing groundwater quality were natural mineral dissolution and mine water mixing, domestic sewage and septic tank contamination and agricultural fertilizer contamination. Selecting the correct analysis method is important for identifying the main factors that control groundwater hydrochemistry. Principal component analysis (PCA) is a mathematical procedure that uses an orthogonal transformation to convert a set of observations of possibly correlated variables into a set of linearly uncorrelated variables called principal components [3]. Therefore, we could identify the unobserved latent controlling factors of groundwater chemistry. Nowadays, PCA has been certified as a powerful tool for determining the controlling factors of groundwater hydrochemistry [14,15] and identifying pollution sources of groundwater quality [16,17]. The Hutuo River alluvial-pluvial fan is located in the western portion of the North China Plain. It is an important source of water in Shijiazhuang, and is one of the important receiving areas of the South–North Water Transfer Project. Thus, the groundwater quality directly affects the drinking water security of the residents in the area. However, with rapid urbanization, industrialization and population growth, the groundwater chemistry has been seriously affected by human activities (such as unregulated discharge of domestic sewage and industrial wastewater and overexploitation of groundwater) [18]. Therefore, it is significant for deepening the understanding of groundwater hydrochemical evolution to study the controlling factors of groundwater hydrochemical evolution in this region. The main aims of this study were to: (1) evaluate the characteristics of groundwater hydrochemistry in a region with intensive human activity; (2) analyze the major ion sources using hydrochemistry and multivariate statistical techniques; (3) identify main factors of controlling groundwater chemistry using PCA. The results could be helpful in developing effective water quality protection strategies and could deepen the understanding of groundwater hydrochemical evolution law in a region with intensive human activity. 2. Materials and Methods 2.1. Study Area The study area is located in the Piedmont inclined plain at the eastern foot of Taihang Mountain, starting from Gangnan Reservoir in the west (the upstream recharge area of Hutuo River alluvial-proluvial fan) to Gaocheng in the east, and extending 5 km from Hutuohe River to the north and south sides. The geographical coordinates are 113°51′05″–114°56′00″ and 37°51′00″–38°24′40″. The administrative area includes Shijiazhuang City, Gaocheng District, Luquan District, Lingshou County, Zhengding County, Wuji County and Pingshan County, with a total area of 2442 km2 (Figure 1) and a total population of about 5.5 million. It has a temperate semi-humid and semi-arid monsoon climate, with average temperature and rainfall of 13.3–15.0 °C and 400–750 mm, respectively [1]. The topography is high in the west and low in the east. The landform is an alluvial-pluvial fan group formed by alluvial-pluvial action, and is a sloping plain formed by river and lake facies deposition. The western part is a hilly area with low and medium mountains, and the eastern part is flat North China Plain. The terrain slope is steeper in the west (5–6‰) and other areas are gentler (less than 2‰). The land use types in the study area include farmland (40%), city (11.60%), surface water (2.35%), town (13.09%), and wasteland (mainly including mountains and wasteland) (34.34%) (Figure 1). The aquifer in the study area belongs to the Quaternary thick aquifer system of Hebei Plain. The stratigraphic lithology is mainly composed of mild clay, sandy loam, sandy mild clay and pebbles, gravel and sand with different particle sizes [19]. The western part of the aquifer is made of a thick layer of sandy gravel; the sandy gravel in the eastern part becomes finer and eventually transitions to mainly coarse sand (Figure 2). The aquifer in the study area is rich in water, and the flow direction of groundwater is the same as that of Hutuo River, which flows from the southwest to northeast. The groundwater mainly consists of bedrock fissure water and loose rock pore water, and the main recharge sources are rainfall vertical infiltration, surface water lateral seepage and flood lateral seepage in rainy season [20]. The study area can be divided into three hydrogeological units, namely the fissure pore water unit in the plain of Hutuo River valley between Gangnan and Huangbizhuang reservoir (HRVP), the hydrogeological unit at the top of Hutuo River alluvial-pluvial fan (THRVPF) and the hydrogeological unit in the central of Hutuo River alluvial-pluvial fan (CHRVPF). The HRVP is located in the transition zone from the mountain area to the plain area. The groundwater depth is shallow (between 2 and 21 m), and it has good conductivity and water richness. In the THRVPF area, the groundwater depth is about 12–40 m and the water conductivity is good. Currently, the main exploitation aquifer is the second water group due to the first aquifer group has been drained. In the CHRVPF area, the groundwater depth is about 40–50 m and the lithology of the aquifer is mainly coarse sand. The third and fourth water groups are the current exploitation aquifer. In recent years, due to the needs of local social and economic development, the groundwater has been intensively exploited, and two large reservoirs have been built in the upstream of the study area, resulting in the recharge source of groundwater being blocked. The depth of groundwater dropped rapidly. At present, one of the largest descending funnels in China has been formed [20]. 2.2. Groundwater Sampling and Laboratory Analyses In October 2015, 46 groups of groundwater samples were collected in the Hutuo River alluvial-pluvial fan area (Figure 1). The groundwater was taken from civil wells and agricultural irrigation wells, with a buried depth of 4.0–50.0 m. One 1.5 L and one 500 mL polyethylene plastic bottle were used for the analysis of anions and cations, respectively. The water samples for cation analysis were added to HNO3 to adjust the pH < 2. Before sampling, rinse the sampling bottles with purified water and original water sample for three times, respectively, and put the samples into the refrigerator for test. The following 16 water chemical parameters were analyzed: pH, total dissolved solids (TDSs), potassium (K+), sodium (Na+), calcium (Ca2+), magnesium (Mg2+), nitrate (NO3−), nitrite (NO2−), ammonia (NH3+), bicarbonate (HCO3−), carbonate (CO32−), chloride (Cl−), sulfate (SO42−), iron (Fe), manganese (Mn) and total hardness (TH). The values of pH were measured in the field using a HQ40D multiparameter instrument (United States of America). Concentrations of the cations and trace elements (Fe and Mn) in the water samples were measured using an inductively coupled plasma atomic emission spectrometer (Agilent 7500ce ICP-MS, Tokyo, Japan), while anion analyses were carried out using spectrophotometry (Perkin-Elmer Lambda 35, Waltham, MA, USA). TDSs were measured using gravimetric methods, HCO3− was measured using acid–base titration and TH was determined using the ethylene diamine tetraacetic acid titration method. All analyses were carried out at the Groundwater Mineral Water and Environmental Monitoring Center at the Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences. 2.3. Data Analysis Since the data obtained in this study did not satisfy the normal distribution, we used the non-parametric Kruskall–Wallis and Mann–Whitney U test methods to analyze the significant differences of the concentration of groundwater hydrochemical parameters in different hydrogeological units and land use types. The Spearman correlation analysis was used to identify the correlation between hydrochemical indicators. The principal component analysis was used to identify the main factors of controlling groundwater hydrochemistry. R (3.6.2) (Core Team, Vienna, Austria) and origin (9.0)(OriginLab Corporation, Northampton, MA, USA) were used for data analysis. 2.4. Multivariate Data Analysis Contrast Coefficient Variance Variance is used to describe the dispersion degree of a random variable to its mathematical expectation (mean). Since the absolute value ranges of the chemical components of groundwater may be very different, it is not possible to compare the variances. Therefore, the ratio of each variable (Xi) of each sample to the mean value of the variable (Xm) is calculated first, which is defined as the contrast coefficient v, Vi = Xi/(Xm) Since the average value of the contrast coefficient is 1, each variable is converted into another variable with the same mean (its contrast coefficient), and then the variance of the contrast coefficient of each variable is calculated to compare the variance. The calculation formula of contrast coefficient variance (Vσi2) is as follows: Vσi2 = 1n∑i=1n(Vi−Vmi)2 where: Vi is the contrast coefficient of the i-th chemical component of groundwater; Vmi is the average value of the contrast coefficient of the i-th chemical component of groundwater. 3. Results 3.1. Main Ionic Characteristics The groundwater hydrochemical characteristics of the Hutuo River alluvial-pluvial fan are shown in Table 1. The pH of the groundwater in the study area is 7.38–8.40, and the average values in three hydrogeological units (valley plain area, top of alluvial-pluvial fan and central of alluvial-pluvial fan) are 7.73, 7.80 and 7.63, respectively. There is no excessive water sample, and the groundwater is weak alkaline water. The TDSs value is 152.2–1913 μs/cm. The mean values and over-standard rates of the three hydrogeological units are 950.4 μs/cm and 25%, 720.4 μs/cm and 11.5%, 412.8 μs/cm and 0%, respectively. The main cation in groundwater is Ca2+, with a concentration between 20.6 and 345.4 mg/L. The average values of Ca2+ in three hydrogeological units are 211.2, 143.8 and 69.8 mg/L, respectively. The order of cation concentration is Ca2+ > Na+ > Mg2+ >K+. The concentration of NH4+ is below the detection limit (0.04 mg/L) in three hydrogeological units in the study area. The main anion of groundwater is HCO3−. Its concentration range is between 141.0 and 429.3 mg/L. The mean concentration of HCO3− in three hydrogeological units is 315.7, 311.7 and 249.3 mg/L, respectively. The order of anion concentration is: HCO3− > SO42− > NO3− > Cl− > NO2− in the valley plain area, and HCO3− > SO42− > Cl− > NO3− > NO2− in the top and middle of alluvial-pluvial fan. CO32− was not detected in the three hydrogeological units in the study area. It is worth noting that the concentration of NO3− is 1.76–465.0 mg/L, and the average concentrations in the three hydrogeological units are 176.2, 75.8 and 9.10 mg/L, with over-the-standard (the Grade III standard for groundwater quality in China [21]) rates of 62.5%, 34.6% and 0%, respectively. It is obvious that groundwater in the valley plain area and the top of the alluvial-pluvial fan has been seriously polluted by NO3−. The concentration range of TH is 76.1–1058 mg/L. The average concentrations of TH in the three hydrogeological units are 650.6, 507.3 and 278.0 mg/L, respectively. The over-the-standard (the Grade III standard for groundwater quality in China [21]) rates are 87.5%, 73.1% and 0%, which is the most important pollution factor in the groundwater in this area. In addition, the concentration range of Fe is 0.01–1.313 mg/L. The average concentrations of Fe in the three hydrogeological units are 0.148, 0.067 and 0.214 mg/L, respectively, and the over-the-standard rates are 12.5%, 3.8% and 8.3%, respectively. The spatial variation of hydrochemical indicators is relatively obvious: Hutuo River valley plain area > top of alluvial-pluvial fan > central of alluvial-pluvial fan, which may be related to the groundwater level depth, the lithology of aeration zone and the intensity of human activities [1]. From the valley plain area to the middle of the alluvial-pluvial fan, the buried depth of the groundwater level increases gradually, the particles in the aeration zone become thinner and the permeability becomes poor, which makes it difficult for pollutants to enter the groundwater and the concentrations of the hydrochemical indexes are relatively low. 3.2. Hydrochemical Type It can be seen from Figure 3 that the main hydrochemical types of groundwater in the study area are HCO3·SO4-Ca and HCO3·SO4-Ca·Mg, accounting for 21.7% and 28.3%, respectively. The proportions of SO4-type water in the three hydrogeological units (valley plain area, top and middle of alluvial-pluvial fan) are 75.0%, 69.2% and 41.2%, respectively. In addition, Cl-type water also accounts for certain proportions in the three hydrogeological units, which are 50.0%, 26.9% and 8.3%, respectively. In this study, the proportions of HCO3-Ca(Mg)-type in the three hydrogeological units are 0%, 11.5% and 41.7%, respectively. 3.3. Analysis of Hydrochemistry A Gibbs diagram is often used to reflect the main control effects of groundwater hydrochemical types (rock weathering, evaporation concentration and precipitation) [22]. The water samples that are greatly affected by atmospheric precipitation are located in the lower right corner of Gibbs figure; the samples that are mainly controlled by rock weathering are located in the middle left corner and the samples; those that are controlled by evaporation crystallization are located in the upper right corner [23]. It can be seen from Figure 4 that most of the water samples from different hydrogeological units fall in the rock weathering area, indicating that rock weathering is the main control factor of groundwater hydrochemistry in this area, and it is less affected by evaporative concentration and atmospheric precipitation. 4. Discussion 4.1. Groundwater Chemical Characteristics in a Region with Intensive Human Activity Groundwater chemical characteristics are influenced by long-term geological history, seepage migration process and human activities [24]. Nowadays, China has experienced rapid urbanization and industrialization, and the influence of human activities on groundwater chemistry is increasing year by year. In the Hutuo River alluvial-pluvial fan region, the groundwater has been intensively exploited and has formed one of the largest descending funnels in China [20]. In addition, the main land use type is farmland in this area—thus, agricultural fertilization and pesticide spraying had serious impact on groundwater chemistry [25]. According to our research, TH and NO3− of the groundwater had significantly exceeded the Grade III standard for groundwater quality in China. Thus, TH and NO3− of groundwater in the study area have been significantly affected by human activities, such as domestic sewage and agricultural fertilizer [25,26]. This result is consistent with previous research. Zhang et al. (2019) found that NO3− was the main impact indicators for poor-quality groundwater in urbanized areas [27]. González Pérez et al. (2020) also found that intensive human activities have serious affected the level of NO3− in groundwater [28]. Egbi et al. (2020) also confirmed that the sources of NO3− in groundwater primarily originate from human contributions in the Lower Volta River basin of Ghana [29]. Therefore, TH and NO3− were the most important parameters, indicating that groundwater is affected by human activities. SO42− is also a vital parameter that reflected the effect of human activities on groundwater chemistry [30]. Indeed, in this study region, especially in the HRVP area, the mean concentration and over-the-standard (the Grade III standard for groundwater quality in China) rate of groundwater SO42− were 176.7% and 12.5%, respectively. Previous research has confirmed that SO42− was an important indicator in highly anthropogenic regions, which can evaluate groundwater vulnerability and pollution risk [31]. Zhang et al. (2020) also found that the contribution proportion of sewage to SO42− in groundwater was 35.5%–42.7% in the Hutuo River basin [18]. Therefore, SO42− is also a vital indicator representing the effects of human activity. In addition, the main hydrochemical type of groundwater in this study area was HCO3·SO4-Ca(Mg). However, previous studies have found that, in 1950, the hydrochemical type of groundwater in the Shijiazhuang area was mainly HCO3-Ca(Mg) [24]. In recent years, the hydrochemical type of groundwater has changed to HCO3·SO4-Ca(Mg), and there is Cl-type water and Na-type water. It can be seen that the groundwater in this area has been significantly affected by human activities, which could influence the safety of the drinking water. 4.2. Analysis of Major Ion Sources 4.2.1. Correlation Analysis If there is a strong correlation between the chemical components in the water environment, it means that they have a common material source. Based on the correlation analysis of groundwater quality parameters in the study area (Table 2), we found that pH had a significant correlation with other indicators. However, TH and TDS showed significant positive correlations with Na+, Ca2+, Mg2+, NO3−, Cl−, HCO3− and SO42−. It indicated that these components may come from common material sources, such as domestic sewage, agricultural fertilizer, industrial wastewater [25,29]. Because it was in the rainy season, during the process of rainwater recharge, a large number of pollutants from the surface were carried and infiltrated into the aquifer, which promoted water–rock interactions, and increased the concentration of groundwater ions [1]. K+ was only positively correlated with Na+. Na+, Ca2+, Mg2+, Cl−, NO3−, HCO3− and SO42− all showed significant positive correlations, indicating that they came from the same source, such as agriculture runoff, urban runoff and water–rock interactions [8,25]. The correlation between Fe and Mn was very significant, but they did not show correlation with other indexes, which proved that they have a common source. 4.2.2. Source Judgment by Using Contrast Coefficient Variance The factors that control the chemical composition of groundwater include natural conditions (such as water–rock interactions) and human activities. Under natural conditions, the chemical composition of groundwater will not change much in a complete groundwater system. However, the hydrochemical composition of a groundwater system affected by human activities often shows a great degree of change [32]. Due to the different impact characteristics on the chemical components of groundwater, the contrast coefficient variance is used in this study to distinguish whether the formation of the chemical components of groundwater is mainly dominated by the evolution of natural conditions or by human activities, and to compare the degree of influence of each component by human activities. On the whole, the variance values of the contrast coefficients of pH, TDS, TH, Ca2+, Mg2+ and HCO3− in the three hydrogeological units are small (Table 3). They are mainly controlled by the evolution of natural hydrochemistry. Because of the shallow groundwater depth, the Hutuo River valley plain unit is greatly affected by human activities and the variance value of the contrast coefficient of each component in groundwater is greater than that of the top and middle hydrogeological units of the Hutuo River alluvial-pluvial fan. It showed that: in the valley plain area, σ2 (K+), σ2(Na+), σ2(Cl−), σ2(NO3−), σ2(Fe) and σ2(Mn) have large values (σ2 > 0.5), which indicated that these components are greatly affected by human activities. In the top of the alluvial-pluvial fan, the values of σ2(Na+), σ2(Cl−), σ2(Mn) and σ2(Fe) are higher (σ2 > 0.5); in the middle of the alluvial-pluvial fan, the values of σ2(Na+), σ2(Cl−), σ2(NO3−) and σ2(Fe) exceeded 0.5, demonstrating that these components are significantly affected by human activities. In the three hydrogeological units, the variance value of the SO42− contrast coefficient is moderately low—0.199, 0.245 and 0.310, respectively—illustrating that it is less affected by human activities and mainly controlled by the evolution of natural hydrochemistry. 4.3. Controlling Factors of Groundwater Chemical In order to prevent groundwater pollution, the main influencing factors must be identified accurately. In this study, we selected 13 water quality parameters (pH, TDS, TH, K+, Na+, Ca2+, Mg2+, NO3−, Cl−, SO42−, HCO3−, Fe and Mn), and identified the main control factors of groundwater pollution in the Hutuo River alluvial-pluvial fan area by PCA. Before principal component analysis, a Kaiser–Meyer–Olkin (KMO) test and Barlett spherical test were carried out on all the data. The results showed that the KMO value was 0.713 and Barlett spherical test value was 963.07 (p < 0.001), indicating that the original data were suitable for PCA. Based on the eigenvalue being greater than 1, three main control factors (Figure 5) that caused groundwater pollution in the Hutuo River alluvial-proluvial fan area were identified. These factors can explain 74.01% of all variables, indicating that they basically contain all the information of the above 13 indicators. The contribution rate of the first main factor (FC1) was TH, Ca2+, NO3−, TDS, Na+, K+ and Cl−, showing a strong positive correlation with FC1, and SO42− showed a moderate positive correlation with PC1. It can be seen that PC1 reflects that groundwater is mainly affected by human activities. The possible sources of Cl− in nature were animal manure, domestic sewage, rising saltwater and road snow melting salt [33]. Since there was no rising saltwater that intruded into the aquifer in the study area, and it was in autumn without snowfall, it can be surmised that neither rising saltwater nor snowmelt salt were the main sources of Cl− in the groundwater in the area. Therefore, the higher concentration of Cl− is most likely to come from animal manure and domestic sewage. The main sources of NO3− included domestic sewage, agricultural fertilizers, manure and atmospheric deposition [34]. Previous studies have found that the nitrate pollution of groundwater is mainly related to the use of agricultural fertilizers, discharge of domestic sewage and sewage irrigation [35,36]. During the investigation and sampling, it was found that in the villages at the top of the Hutuo River alluvial-pluvial fan, a large amount of domestic sewage was directly discharged into the Hutuo River channel. In addition, the lithology of the stratum at the top of the fan was relatively coarse, and the sewage was easily infiltrated into the aquifer. As a result, domestic sewage was the main influencing factor of the groundwater in the study area. Furthermore, the main land use type was agricultural land (38.62%), so agricultural fertilizers were also an important source. In conclusion, we can conclude that the first factor (PC1) controlling groundwater chemistry in this area is human activity (domestic sewage and agricultural fertilizer). The contribution rate of the second main factor (FC2) was 13.78%. HCO3− and Mg2+ showed a strong positive correlation with PC2, pH showed a medium negative correlation with PC2 and SO42− showed a medium positive correlation with PC2. In general, HCO3−, Mg2+ and SO42− in groundwater were mainly from natural sources (dissolution of dolomite and limestone) [33]. Based on the previous analysis of the contrast coefficient, HCO3−, Ph, Mg2+ and SO42− in the groundwater of the study area were mainly controlled by the chemical evolution of natural water. In recent years, with the development of industry and agriculture, the groundwater in the Hutuo River alluvial-pluvial fan area has been intensively exploited. In addition, two large reservoirs have been built in the upstream of the study area, which cut off the recharge source of groundwater in the area, and continuously reduced the groundwater level. A depression funnel was formed in some parts, breaking the original balance of the groundwater dynamic field, thus triggering a series of physicochemical reactions. The equilibrium is conducive to the dissolution of carbonate minerals (calcite and dolomite), and then the concentration of HCO3−, Ca2+ and Mg2+ increased [1]. Therefore, the second factor controlling groundwater in this area is a natural factor (water–rock interaction). The contribution rate of the third main factor (FC3) is 8.77%. Fe and Mn have strong positive correlations with PC3. In the Hutuo River alluvial-pluvial fan area, the average concentrations of Fe in the groundwater of three hydrogeological units are 0.148, 0.067 and 0.214 mg/L, respectively, and there were four samples that exceeded the national groundwater quality standard of 0.3 mg/L. Previous studies have found that high concentrations of Fe in the water environment were mainly from industrial pollution [37,38]. There was a steel factory upstream of the study area, and the sampling point near it (5) has the highest Fe concentration (0.552 mg/L). The sampling point near Gaocheng industrial park (35) downstream of the study area had an Fe concentration of 1.313 mg/L, which further confirms our hypothesis. Therefore, PC3 shows that the groundwater in the study area is polluted by industrial sewage. Thus, the third factor controlling the groundwater in the study area is the human activity (industrial sewage). 5. Managerial Implications and Future Research The results indicate that high concentrations of NO3− and TH in groundwater are present in the region with intensive human activity. Local governments and environmental managers should pay due attention to NO3−, which can cause not only eutrophication and toxic algal blooms [39], but also increase the risk of diseases, such as methemoglobinemia, blue baby syndrome and stomach cancer [40]. Therefore, groundwater in the region must be treated before drinking. In addition, the main factors controlling the groundwater hydrochemistry in this area are domestic sewage, industrial sewage and agricultural fertilizer. Thus, it is suggested that the government and water environment management department should first stop the substandard discharge of domestic and industrial sewage, and promote optimization of fertilization strategies (such as formula fertilization based on soil testing) in order to effectively prevent the further deterioration of groundwater quality in this area. In this study, we have identified the main pollutants and factors controlling groundwater hydrochemistry in the region with intensive human activity. However, we have not quantitatively identified the source of the various pollutants. Nowadays, a quantitative traceability technique combined with isotope and source apportionment models (SIAR model) has been widely used in water environments. For example, previous researchers combined the isotope NO3− (δ15N and δ18O) and SIAR model to identify the NO3− sources [41], and used the isotope SO42− (δ34S and δ18O) and SIAR model to identify the SO42− sources [42]. Therefore, in order to formulate effective controlling measures for controlling groundwater pollution, it is necessary to quantitatively identify the sources of pollutants in groundwater in this region. 6. Conclusions In this research, we applied hydrochemistry (Piper and Gibbs diagrams) and multivariate statistical techniques (variance of contrast coefficient, correlation analysis and PCA) to study the hydrochemistry characteristics and the main control factors of groundwater in the region with intensive human activity—the Hutuo River alluvial-pluvial fan of northern China. The results showed that the groundwater quality was relatively poor, among which the concentration range of NO3− and TH are 1.76–465.0 mg/L and 76.1–1058 mg/L, respectively. The over-the-standard rates in three hydrogeological units are 62.5%, 34.6% and 0%, and 87.5%, 73.1% and 0%, respectively, which were the most important parameters indicating that groundwater is affected by human activities. Rock weathering is the main control function of groundwater hydrochemistry in this area. The main types of groundwater hydrochemistry are HCO3·SO4–Ca and HCO3·SO4–Ca·Mg, accounting for 21.7% and 28.3%, respectively. SO4-type water in three hydrogeological units is as high as 75.0%, 69.2% and 41.2%, respectively. In addition, there are Cl- and Na-type waters, indicating that the groundwater in this area has been significantly affected by human activities. The variance analysis of contrast coefficient shows that K+, Na+, NO3−, Cl−, Fe and Mn are greatly influenced by human activities, and the other indexes are mainly influenced by natural causes. The results of PCA show that the three factors controlling groundwater hydrochemistry in the study area are domestic sewage and chemical fertilizer (human activity), water–rock interactions (natural factor) and industrial sewage (human activity). In summary, the main factors controlling the groundwater hydrochemistry in this area are domestic sewage, industrial sewage and agricultural fertilizer. Therefore, in order to effectively prevent the further deterioration of groundwater quality in this area, it is suggested that the government and the water environment management department should first stop the substandard discharge of domestic and industrial sewage and promote optimization of fertilization strategies. In addition, it is necessary to quantitatively identify the sources of pollutants in the future in order to formulate effective controlling measures for controlling groundwater pollution. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Author Contributions Conceptualization, Q.Z. and C.R.; methodology, Q.Z.; software, C.R.; validation, Q.Z.; formal analysis, C.R.; investigation, C.R., and Q.Z.; resources, C.R.; data curation, C.R.; writing—original draft preparation, C.R., and Q.Z.; writing—review and editing, Q.Z.; visualization, C.R.; supervision, Q.Z.; project administration, Q.Z.; funding acquisition, Q.Z. All authors have read and agreed to the published version of the manuscript. Funding This research was funded by the Projects of China Geological Survey, China, grant number DD20190331, the National Key R&D Program of China, grant number 2018YFC1803001, the Belt and Road Fund on Water and Sustainability, China, grant number U2019nkms01, and Open Foundation of the State Key Laboratory of Urban and Regional Ecology of China, grant number SKLURE2019-2-3. Conflicts of Interest The authors declare no conflict of interest. Figure 1 Monitoring sites of groundwater in Hutuo River alluvial-pluvial fan. Figure 2 Hydrogeological cross-section of Hutuo River pluvial fan. Figure 3 Piper diagram of groundwater in Hutuo River alluvial-pluvial fan. Figure 4 Gibbs diagram for the groundwater in Hutuo River alluvial-pluvial fan. Figure 5 Loadings of 13 selected variables on varimax rotated factors in the Hutuo River alluvial-pluvial fan. ijerph-17-09126-t001_Table 1Table 1 Descriptive statistics of groundwater chemical parameters in three hydrogeological units in Hutuo River alluvial-pluvial fan. Parameters pH TDS K+ Na+ Ca2+ Mg2+ HCO3− Cl− SO42− NO3− NO2− Fe Mn TH (I) Hutuo River valley plain unit (n = 8) Min 7.64 500.7 0.6 10.1 135.2 10.23 187.0 51.09 62.4 68.8 0.002 0.017 0.001 379.8 Max 7.88 1913 22.0 180.8 345.4 53.1 381.4 322.4 322.3 465.0 0.007 0.552 0.012 1058 Mean 7.73 950.4 5.80 50.9 211.7 29.7 315.7 121.5 176.7 176.2 0.003 0.148 0.004 650.6 SD 0.08 455.6 7.66 56.1 68.3 15.8 70.0 99.3 78.7 139.7 0.002 0.180 0.004 219.7 Exceed standard rate (%) 0 25 — 0 — — — 12.5 12.5 62.5 0 12.5 0 87.5 (II) Top of the Hutuo River pluvial fan (n = 26) Min 7.42 241 0.27 9.08 53.5 11.7 171.7 15.9 21.8 5.04 0.002 0.010 0.001 182.2 Max 8.40 1508 7.05 307.2 224.1 79.1 429.3 343.5 318.2 216.0 0.041 0.324 0.013 814.2 Mean 7.80 720.4 2.00 47.3 143.8 36.0 311.7 89.9 151.4 75.8 0.005 0.067 0.003 507.3 SD 0.22 280.7 1.26 58.0 41.6 14.6 68.8 72.6 74.9 45.0 0.009 0.084 0.004 157.7 Exceed standard rate (%) 0 11.5 — 3.8 — — — 3.8 3.8 34.6 0 3.8 0 73.1 (III) Center of the Hutuo River pluvial fan (n = 12) Min 7.38 152.2 1.26 17.7 20.6 5.95 141.0 1.06 5.27 1.76 0.002 0.010 0.001 76.1 Max 7.97 626.5 2.56 139 110.1 41.5 334.2 135.3 168.6 21.8 0.004 1.313 0.040 445.9 Mean 7.63 412.8 2.13 36.5 69.8 25.2 249.3 42.6 81.4 9.10 0.002 0.214 0.005 278.0 SD 0.18 124.9 0.50 37.1 22.3 10.8 56.13 35.0 45.3 7.35 0.001 0.399 0.011 96.0 Exceed standard rate (%) 0 0 — 0 — — — 0 0 0 0 8.3 0 0 Standard 6.5–8.5 1000 — 200 — — — 250 250 88.6 3.29 0.3 0.1 450 Note: Mean: average value; SD: standard deviation; Min: minimum value; Max: maximum value; Standard is grade III standard for groundwater quality in China (GB/T14848–2017). ijerph-17-09126-t002_Table 2Table 2 Spearman’s correlation matrix of water quality of groundwater in the Hutuo River alluvial-pluvial fan. Parameters pH TDS K+ Na+ Ca2+ Mg2+ Cl− SO42− HCO3− NO3− Fe Mn TH pH 1.000 TDS −0.079 1.000 K+ −0.066 0.215 1.000 Na+ −0.088 0.687 ** 0.320 * 1.000 Ca2+ −0.067 0.917 ** 0.060 0.439 ** 1.000 Mg2+ −0.177 0.752 ** 0.196 0.540 ** 0.648 ** 1.000 Cl− −0.108 0.863 ** 0.170 0.706 ** 0.741 ** 0.598 ** 1.000 SO42− 0.015 0.826 ** 0.231 0.709 ** 0.746 ** 0.634 ** 0.681 ** 1.000 HCO3− −0.165 0.823 ** 0.257 0.656 ** 0.742 ** 0.808 ** 0.631 ** 0.670 ** 1.000 NO3− 0.043 0.868 ** 0.222 0.370 * 0.909 ** 0.585 ** 0.683 ** 0.649 ** 0.645 ** 1.000 Fe 0.148 0.080 0.111 0.051 0.076 0.067 0.063 −0.006 0.127 0.030 1.000 Mn 0.091 0.032 0.205 −0.018 0.048 −0.012 −0.018 −0.007 0.157 0.029 0.737 ** 1.000 TH −0.098 0.956 ** 0.130 0.519 ** 0.964 ** 0.803 ** 0.770 ** 0.756 ** 0.820 ** 0.889 ** 0.083 0.040 1.000 Note: * p < 0.05; ** p < 0.01. ijerph-17-09126-t003_Table 3Table 3 Variance analysis of contrast coefficient of groundwater hydrochemical compositions in the Hutuo River alluvial-pluvial fan. Parameters HRVP THRAPF CHRAPF pH 0.000 0.001 0.001 TH 0.114 0.097 0.119 TDS 0.230 0.152 0.092 K+ 1.744 0.397 0.056 Na+ 1.216 1.503 1.033 Ca2+ 0.104 0.084 0.102 Mg2+ 0.285 0.165 0.183 Cl− 0.668 0.653 0.675 SO42− 0.199 0.245 0.310 HCO3− 0.049 0.049 0.051 NO3− 0.629 0.353 0.652 Fe 1.474 1.562 3.480 Mn 1.011 1.541 0.292 Note: HRVP: Hutuo Rvier Valley Plain; THRAPF: Top of Hutuo River alluvial-pluvial fan; CHRAPF: Center of Hutuo River alluvial-pluvial fan. ==== Refs References 1. Zhang Q. Wang H. Wang Y. Yang M. Zhu L. Groundwater quality assessment and pollution source apportionment in an intensely exploited region of northern China Environ. Sci. Pollut. Res. 2017 24 16639 10.1007/s11356-017-9114-2 2. Jiang Y. Wu Y. Groves C. Yuan D. Kambesis P. Natural and anthropogenic factors affecting the groundwater quality in the Nandong karst underground river system in Yunan, China J. Contam. Hydrol. 2009 109 49 10.1016/j.jconhyd.2009.08.001 19717207 3. Qin R. Wu Y. Xu Z. Xie D. Zhang C. Assessing the impact of natural and anthropogenic activities on groundwater quality in coastal alluvial aquifers of the lower Liaohe River Plain, NE China Appl. Geochem. 2013 31 142 10.1016/j.apgeochem.2013.01.001 4. Xue D. Botte J. Baets B.D. Accoe F. Nestler A. Taylor P. Cleemput O.V. Berglund M. Boeckx P. Present limitations and future prospects of stable isotope methods for nitrate source identification in surface- and groundwater Water Res. 2009 43 1159 10.1016/j.watres.2008.12.048 19157489 5. Moore W.S. The subterranean estuary: A reaction zone of ground water and sea water Mar. Chem. 1999 65 111 125 10.1016/S0304-4203(99)00014-6 6. Tziritis E. Skordas K. Kelepertsis A. The use of hydrogeochemical analyses and multivariate statistics for the characterization of groundwater resources in a complex aquifer system. A case study in Amyros River basin, Thessaly, central Greece Environ. Earth Sci. 2016 75 339 10.1007/s12665-015-5204-y 7. Pazand K. Khosravi D. Ghaderi M.R. Rezvanianzadeh M.R. Identification of the hydrogeochemical processes and assessment of groundwater in a semi-arid region using major ion chemistry: A case study of Ardestan basin in Central Iran Groundw. Sustain. Dev. 2018 6 245 254 10.1016/j.gsd.2018.01.008 8. Huang G. Sun J. Zhang Y. Chen Z. Liu F. Impact of anthropogenic and natural processes on the evolution of groundwater chemistry in a rapidly urbanized coastal area, South China Sci. Total Environ. 2013 463 209 10.1016/j.scitotenv.2013.05.078 23810860 9. Lin C.Y. Abdullah M.H. Praveena S.M. Yahaya A.H.B. Musta B. Delineation of temporal variability and governing factors influencing the spatial variability of shallow groundwater chemistry in a tropical sedimentary island J. Hydrol. 2012 432 26 10.1016/j.jhydrol.2012.02.015 10. Li X. Tang C. Cao Y. Li D. A multiple isotope (H, O, N, C and S) approach to elucidate the hydrochemical evolution of shallow groundwater in a rapidly urbanized area of the Pearl River Delta, China Sci. Total Environ. 2020 724 137930 10.1016/j.scitotenv.2020.137930 32268282 11. Güler C. Kurt M.A. Alpaslan M. Akbulut C. Assessment of the impact of anthropogenic activities on the groundwater hydrology and chemistry in Tarsus coastal plain (Mersin, SE Turkey) using fuzzy clustering, multivariate statistics and GIS techniques J. Hydrol. 2012 414–415 435 451 12. Li Z. Yang Q. Yang Y. Ma H. Wang H. Luo J. Bian J. Martin J.D. Isotopic and geochemical interpretation of groundwater under the influences of anthropogenic activities J. Hydrol. 2019 576 685 697 10.1016/j.jhydrol.2019.06.037 13. Gu H. Chi B. Li H. Jiang J. Qin W. He W. Assessment of groundwater quality and identification of contaminant sources of Liujiang basin in Qinhuangdao, North China Environ. Earth Sci. 2015 73 6477 10.1007/s12665-014-3870-9 14. Boateng T.K. Opoku F. Acquaah S.O. Akoto O. Groundwater quality assessment using statistical approach and water quality index in Ejisu-Juaben Municipality, Ghana Environ. Earth Sci. 2016 75 489 10.1007/s12665-015-5105-0 15. Herojeet R. Rishi M.S. Lata R. Sharma R. Application of environmetrics statistical models and water quality index for groundwater quality characterization of alluvial aquifer of Nalagarh Valley, Himachal Pradesh, India Sustain. Water Resour. Manag. 2016 2 39 10.1007/s40899-015-0039-y 16. Gholizadeh M.H. Melesse A.M. Reddi L. Water quality assessment and apportionment of pollution sources using APCS-MLR and PMF receptor modeling techniques in three major rivers of South Florida Sci. Total Environ. 2016 566–567 1552 10.1016/j.scitotenv.2016.06.046 27317134 17. Zhang Q. Wang X. Wan W. Hou P. Li R. Ouyang Z. The Spatial-Temporal Pattern and Source Apportionment of Water Pollution in a Trans-Urban River Pol. J. Environ. Stud. 2015 24 841 18. Zhang Q.Q. Wang H.W. Lu C. Tracing sulfate origin and transformation in an area with multiple sources of pollution in northern China by using environmental isotopes and Bayesian isotope mixing model Environ. Pollut. 2020 265 115105 10.1016/j.envpol.2020.115105 32806467 19. Foster S. Garduno H. Evans R. Olson D. Tian Y. Zhang W. Han Z. Quaternary aquifer of the north china plain—Assessing and achieving groundwater resource sustainability Hydrogeol. J. 2004 12 81 93 10.1007/s10040-003-0300-6 20. Li Y. Zhang Z. Fei Y. Chen H. Qian Y. Yu D. Investigation of quality and pollution characteristics of groundwater in the Hutuo River Alluvial Plain, North China Plain Environ. Earth Sci. 2016 75 581 10.1007/s12665-016-5366-2 21. Ministry of Natural Resources of the People’s Republic of China Standard for Groundwater Quality, (GB/T14848-2017) General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China Beijing, China 2017 (In Chinese) 22. Li P. Wu J. Qian H. Assessment of groundwater quality for irrigation purposes and identification of hydrogeochemical evolution mechanisms in Pengyang County, China Environ. Earth Sci. 2013 69 2211 10.1007/s12665-012-2049-5 23. Tang J.P. Zhang Q. Hu Y. Groundwater chemical characteristics and analysis of their controlling factors in an alluvial fan of Jianjiang River Environ. Sci. 2019 40 3089 (In Chinese) 24. Zhang Q. Miao L. Wang H. Hou J. Li Y. How Rapid Urbanization Drives Deteriorating Groundwater Quality in a Provincial Capital of China Pol. J. Environ. Stud. 2020 29 441 10.15244/pjoes/103359 25. Zhang Q. Wang H. Assessment of sources and transformation of nitrate in the alluvial-pluvial fan region of north China using a multi-isotope approach J. Environ. Sci. 2020 89 9 10.1016/j.jes.2019.09.021 31892404 26. Meghdadi A. Javar N. Quantification of spatial and seasonal variations in the proportional contribution of nitrate sources using a multi-isotope approach and Bayesian isotope mixing model Environ. Pollut. 2018 235 207 10.1016/j.envpol.2017.12.078 29289831 27. Zhang F.E. Huang G.X. Hou Q.X. Liu C.Y. Zhang Y. Zhang Q. Groundwater quality in the Pearl River Delta after the rapid expansion of industrialization and urbanization: Distributions, main impact indicators, and driving forces J. Hydrol. 2019 577 124004 10.1016/j.jhydrol.2019.124004 28. González Pérez D.M. Martín Martín J.M. Guaita Martínez J.M. Sáez-Fernández F.J. An Analysis of the Cost of Water Supply Linked to the Tourism Industry. An Application to the Case of the Island of Ibiza in Spain Water 2020 12 2006 10.3390/w12072006 29. Egbi C.D. Anornu G.A. Ganyaglo S.Y. Appiah-Adjei E.K. Li S.L. Dampare S.B. Nitrate contamination of groundwater in the Lower Volta River Basin of Ghana: Sources and related human health risks Ecotoxicol. Environ. Saf. 2020 191 110227 10.1016/j.ecoenv.2020.110227 31981956 30. Ma M.D.C. Esteller M.V. Exposito J.L. Hirata R. Impacts of urbanization on groundwater hydrodynamics and hydrochemistry of the Toluca Valley aquifer (Mexico) Environ. Monit. Assess. 2014 186 2979 24590229 31. Busico G. Kazakis N. Colombani N. Mastrocicco M. Voudouris K. Tedesco D. A modified SINTACS method for groundwater vulnerability and pollution risk assessment in highly anthropized regions based on NO3 − and SO4 2− concentrations Sci. Total Environ. 2017 609 1512 1523 10.1016/j.scitotenv.2017.07.257 28800693 32. Meng F.A. Liang X.J. Hao Y. Wang Y.H. Lou Y. Li H.W. Dynamic Characteristics of Groundwater in Taoer River Fan Water Sav. Irrig. 2016 4 65 (In Chinese) 33. Liu C.Q. Li S.L. Lang Y.C. Xiao H.Y. Using delta15N- and delta18O-values to identify nitrate sources in karst ground water, Guiyang, southwest China Environ. Sci. Technol. 2006 40 6928 10.1021/es0610129 17153996 34. Jin Z. Xue Q. Chen L. Jin M. Li F. Using dual isotopes to evaluate sources and transformations of nitrate in the West Lake watershed, eastern China J. Contam. Hydrol. 2015 177 64 10.1016/j.jconhyd.2015.02.008 25835546 35. Bohlke J.K. Groundwater recharge and agricultural contamination Hydrogeol. J. 2002 10 153 10.1007/s10040-001-0183-3 36. Chen J. Tang C. Sakura Y. Yu J. Fukushima Y. Nitrate pollution from agriculture in different hydrogeological zones of the regional groundwater flow system in the North China Plain Hydrogeol. J. 2005 13 481 10.1007/s10040-004-0321-9 37. Bu H. Tan X. Li S. Zhang Q. Temporal and spatial variations of water quality in the Jinshui River of the South Qinling Mts., China Ecotoxicol. Environ. Saf. 2010 73 907 10.1016/j.ecoenv.2009.11.007 20047760 38. Juahir H. Zain S.M. Yusoff M.K. Hanidza T.T. Armi A.M. Toriman M.E. Mokhtar M. Spatial water quality assessment of Langat River Basin (Malaysia) using environmetric techniques Environ. Monit. Assess. 2011 173 625 10.1007/s10661-010-1411-x 20339961 39. Nestler A. Berglund M. Accoe F. Duta S. Xue D. Boeckx P. Taylor P. Isotopes for improved management of nitrate pollution in aqueous resources: Review of surface water field studies Environ. Sci. Pollut. Res. 2011 18 519 533 10.1007/s11356-010-0422-z 40. Past en-Zapata E. Ledesma-Ruiz R. Harter T. Ramı´rez A.I. Mahlknecht J. Assessment of sources and fate of nitrate in shallow groundwater of an agricultural area by using a multi-tracer approach Sci. Total Environ. 2014 470–471 855 864 10.1016/j.scitotenv.2013.10.043 41. Guo Z.F. Yan C.Z. Wang Z.S. Xu F.F. Yang F. Quantitative identification of nitrate sources in a coastal peri-urban watershed using hydrogeochemical indicators and dual isotopes together with the statistical approaches Chemosphere 2020 243 125364 10.1016/j.chemosphere.2019.125364 31760285 42. Torres-Martínez J.A. Mora A. Knappett P.S.K. Ornelas-Soto N. Mahlknecht J. Tracking nitrate and sulfate sources in groundwater of an urbanized valley using a multi-tracer approach combined with a Bayesian isotope mixing model Water Res. 2020 182 115962 10.1016/j.watres.2020.115962 32629319