Characterization and Cause Analysis of Shallow Groundwater Hydrochemistry in the Plains of Henan Province, China

: With the development of the human population and society, groundwater environmental problems have become an important factor limiting global socioeconomic development, and the study of groundwater hydrochemical characteristics and pollution is a current hot issue. In this study, data regarding shallow groundwater quality in 76 instances were used to evaluate the quality of shallow groundwater in the plains of Henan Province, China, by using a combination of subjective and objective assignments, mathematical statistics, Piper trilinear diagram, Gibbs diagram, ion ratio analysis, and other methods to study the hydrochemical characteristics of groundwater and its formation mechanism. The results showed that the groundwater quality in most areas of Henan Plain is in good condition, and the proportion of samples with excellent grades and good grades is as high as 43.42% and 35.53%. The range of poor and extremely poor water quality is small, and only ﬁve samples are judged as poor and extremely poor grades, mainly distributed in Jiaozuo City, Xinxiang City, Zhoukou City, and Puyang City. The groundwater anionic hydrochemistry is mainly of the HCO 3 type, accounting for 61.84% of the samples and locally transformed downstream to HCO 3 · SO 4 , HCO 3 · SO 4 · Cl, HCO 3 · Cl · SO 4 , and Cl · SO 4 · HCO 3 . Cations are predominantly of the Ca/Mg and Ca– Mg/Mg–Ca type, and gradually transformed to the Na–Ca/Ca–Na and Na–Mg/Mg–Na type along the runoff direction. Water–rock interactions and anthropogenic factors dominate the hydrochemistry evolution, with major geochemical processes involving the precipitation of calcite and dolomite as well as the weathering dissolution of rock salt and ﬂuorite. Human activity is an important factor affecting the distribution of NO 3 –N and Fe 3+ . It is recommended that groundwater be continuously monitored to provide scientiﬁc data for sustainable groundwater quality management and that appropriate measures be developed to prevent further degradation of the groundwater environment.


Introduction
Groundwater is an important source of fresh water for human survival, with about one-third of the world's population drinking from groundwater and about 40% of food production also relying on groundwater [1,2]. However, as modernization accelerates around the world, increasingly intense agricultural, industrial, and human activities pose a great threat to groundwater ecology [3]. According to statistics, globally, about 80% of the population do not have access to safe drinking water sources, and approximately 2 million groundwater resource renewal studies, sustainable mining assessments, or risk assessments of individual sources of pollution [35][36][37] and have not yet undertaken a comprehensive and systematic assessment of groundwater quality. In view of this, the present study was conducted based on sampling data of shallow groundwater in the plains of Henan Province, using the entropy weight method with combined subjective assignment to determine the groundwater quality index, as well as mathematical statistics, Piper trilinear diagram, Gibbs diagram, correlation analysis, and ion ratio analysis to analyze its hydrochemical characteristics and formation mechanism and explore its pollution sources, in order to provide some scientific basis for the development, management and protection of shallow groundwater resources.

Data Sampling
Henan Province is located in the central part of China, with coordinates betweeñ 31 • 23 -36 • 22 N and~110 • 21 -116 • 39 E. It has a total area of 167,000 km 2 and is an important comprehensive transportation hub and center of human, logistics, and information flow in China. It has a continental monsoon climate with a transition from the north subtropical zone to the warm temperate zone, with remarkable climatic features, an average temperature of~12-16°C, and an average annual rainfall of~532.5-138.6 mm, mostly concentrated in the months of June to August [38]. Groundwater has become the main source of water supply in Henan Province owing to year-on-year increasing exploitation [39]. In this study, a total of 76 sets of shallow groundwater samples were collected in the plains of Henan Province, with a sampling depth of less than 50 m and a relatively uniform distribution of sampling points ( Figure 1). The samples were tested for K + , Na + , Ca 2+ , Mg 2+ , Fe 3+ , Cl − , SO 4 2− , HCO 3 − , F − , NO 3 -N, Mn, As, total hardness (TH), total dissolved solid (TDS), and pH. The detection limits for Mn and As were 0.001 mg/L. During data processing, samples with values less than the detection values for these two indicators were treated as 0.001 mg/L.

Research Methodology
After the completion of sample testing and analysis, the entropy weight method combined with the subjective assignment method was used to determine the relative weights of the evaluation indexes and to continue the traditional water quality index (WQI) method to classify the groundwater classes in the study area. The characteristics of the hydrochemistry parameters were analyzed by mathematical statistics to explore the spatial variability of groundwater chemical ions in the study area. The formation pathways of groundwater hydrochemistry were analyzed using a Gibbs diagram to reveal the ionic origins of the water bodies. The Piper trilinear diagram was used to analyze the hydrochemical composition of the water bodies in the study area and their distribution characteristics. Correlation analysis and ion ratio relationship analysis were conducted to reveal the main hydrochemical processes of groundwater quality assessment. The relative weights of the evaluation indicators were determined as follows: The sample data were standardized according to (Equation (1)).
where x ij is the concentration value for sample i of the jth chemical indicator, x j,max is the maximum value of the concentration of the jth chemical indicator, and x j,max is the maximum value of the concentration of the jth chemical indicator.  The weight of the ith sample was calculated under the jth chemical indicator after standardization: where n is the number of samples. The information entropy (Equation (3)), information utility value (Equation (4)), and entropy weight (Equation (5)) were calculated for each chemical indicator [40]: The average of the calculated entropy weights and subjective assigned weights was used as the final indicator weight for this evaluation to calculate the quality rating of each chemical indicator in each sample by Equation (6) and the groundwater quality index (GWQI) by Equation (7), where the subjective assigned weights were calculated according to the relative importance of the indicators provided in the literature [26,27] as follows: where C i,j is the concentration of sample i of chemical indicator j, and S j is the WHOrecommended limit value for chemical indicator j.

Statistical Characteristics of Groups
To assess the data characteristics of hydrochemical indicators in the study area, the minimum, maximum, and mean values of each chemical indicator for the 76 groups of tested samples and the coefficient of variation were calculated and analyzed in comparison with World Health Organization's (WHO) guidelines for drinking-water quality (Table 1), which state that if an indicator exceeds the limit value, it is detrimental to human health [41]. Moreover, the groundwater chemical indicator concentrations in the study area were interpolated using the ordinary kriging method (Figure 2), and the concentrations of each indicator were classified into 10 classes using the geometric partitioning method to determine more significant differences.  Table 1 and Figure 2 show that the concentrations of Fe 3+ and Mn in the samples ranged from 0.010 to 6.130 mg/L and from 0.001 to 2.060 mg/L, with mean values of 0.493 and 0.182 mg/L, respectively, and both mean values were higher than the WHOrecommended limits. The coefficients of variation of Fe 3+ and Mn were 2.062 and 1.797, respectively, which were relatively larger than those of the other chemical indicators, and Fe 3+ and Mn had a large uneven distribution of concentrations, which also indicates that these two chemical indicators are sensitive to environmental changes. The pH in the samples varied between 7.000 and 8.300, with a mean value of 7.423 and a coefficient of variation of 0.040, and the dispersion of the concentration distribution was small, indicating that the shallow groundwater in the study area was overall weakly alkaline, which meets WHO drinking water standards. The TDS in the samples varied between 333.090 and 1764.690 mg/L, with a mean value of 779.659 mg/L and a variation coefficient of 0.422, indicating a moderate degree of dispersion, whereby the WHO standard is 1000 mg/L; 27.63% of the samples exceeded the recommended limit. Total groundwater hardness is the sum of all metal ions except alkali metals, and the TDS in the tested samples varied between 25.000 and 1057.500 mg/L, with a mean value of 436.901 mg/L and a coefficient of variation of 0.464. The degree of dispersion was similar to that of TDS, with 28.94% of the samples exceeding the WHO-recommended limit value. The maximum concentration of arsenic in the samples was up to 0.062 mg/L, whereas the WHO-recommended limit was 0.010 mg/L; therefore, the maximum value far exceeded the recommended limit. Statistically, approximately 6.48% of the samples exceeded the limit for As, and the coefficient of variation for this indicator was 2.534, with a concentration of areas in the study area where As levels exceeded the limit. Nitrate nitrogen varied between 0.005 and 69.320 mg/L, with a relatively small mean value of 7.689 and a coefficient of variation of 1.797. The WHOrecommended limit value for fluoride is 1.5 mg/L, and 17.11% of the samples collected in this study area exceeded the recommended value. The average anion concentrations were HCO 3 − > SO 4 2− > Cl − , with HCO 3 − concentrations being significantly higher than those of the other two indicators. The mean values of all three concentrations were lower than the WHO-recommended limits, and the coefficient of variation was not significantly different, varying between 0.333 and 0.784. Furthermore, 26.32%, 9.21%, and 5.26% of the sample data exceeded the recommended limits for these three indicators, respectively. The average concentration of cations was Na + > Ca 2+ > Mg 2+ > K + , similar to that of anions, and the average values of all three concentrations were also below the WHO-recommended limits. The coefficient of variation of K + was as high as 2.696, which was quite different from that of the other three cations. In terms of spatial distribution, Na + , Mg 2+ , Cl − , SO 4 2− , HCO 3 − , F − , TDS, and pH concentrations were similarly distributed, with a gradual increase from the southwest to the northeast. K + was more abnormally distributed relative to the other cations, with high concentrations of K + in almost every city except Xinyang City, Hebi City, and Anyang City, with a patchy distribution, especially in Pingdingshan City. The concentration of TH was higher in the northwestern part of the study area, with a decreasing trend from the northwest to the southeast, and the smallest concentration of 25.000 mg/L in the southeastern part of Xinyang City. The high concentration of Fe 3+ was mainly distributed in Jiaozuo City, Xinxiang City, Puyang City, and the northwestern part of Zhoukou City, with obvious variability in distribution. The distribution of Mn and Fe 3+ was basically similar, but the overall concentration was significantly lower than that of Fe 3+ . In contrast, the distribution of As was opposite to that of NO 3 -N, and only a certain amount of As was present in central Xinxiang City, southeastern Anyang City, and southwestern Kaifeng City.

Distribution of Hydrochemistry Types
From the Piper trilinear diagram of the groundwater chemistry in Figure 3, it can be seen that most of the samples had the hydrochemistry type of HCO 3 , accounting for approximately 61.84% of the total samples; 32.89% of the samples had the hydrochemistry types of HCO 3 -Cl and HCO 3 -SO 4 , and they are mainly located in the central and northern part of the study area. Moreover, a small number of samples had the hydrochemistry type of a mixture of three anions and were mainly located in the central part of Xinxiang City and southwest of Hebi City, showing a trend of conversion from the HCO 3 type to the HCO 3 -SO 4 type, HCO 3 -SO 4 -Cl type, HCO 3 -Cl-SO 4 , and Cl-SO 4 -HCO 3 type. From the cationic point of view, the hydrochemistry types were dominated by the Ca/Na, Ca-Mg/Mg-Ca, and Na-Mg/Mg-Na types. The sample distribution of these water chemistry types had a certain zonation, as shown by the predominance of the Ca/Mg-and Ca-Mg/Mg-Ca-type water in the pre-mountain sloping plain, where As toward the central plain area, the Na+ content gradually increased, and the water chemistry type began to change to the Na-Ca/Ca-Na type, Na-Mg/Mg-Na type, and finally to the Na type, which predominated in the eastern low plain area. As shown by the prismatic diagram in Figure 3, the groundwater types were overall predominated by alkaline earth metals and weak acids.

Groundwater Quality Assessment Results
The relative weights of the indicators calculated according to the entropy weighting method and the literature [26,27] are shown in Table 2, where the missing items are in analogous values, and their average values were taken to calculate the water quality index according to Equations (6) and (7). The calculated water quality index for each sample was interpolated according to the inverse distance weighting method and graded according to Table 3 to obtain the distribution of the water quality index for the whole study area (Figure 4).  According to the numerical calculation and interpolation results of each groundwater sample, there were 33 samples with excellent grade, accounting for 43.42%; 27 samples with good grade, accounting for 35.53%; 11 samples with moderate grade, accounting for 14.47%; and fewer samples with poor and very poor grade of 3 and 2, respectively. The areas with excellent water quality in the study area were mainly distributed in Xinyang City, Zhumadian City, most of Xuchang City, central and northeastern Pingdingshan City, central Anyang City, and southeastern Shangqiu City, with a small distribution in other cities. The areas with good water quality occupied most of the entire study area. The area with moderate water quality was mainly distributed in the central and northeastern part of the study area, namely Kaifeng City, Xinxiang City, and Puyang City, and had a certain area distribution in Jiaozuo City. The area of poor and very poor water quality in the study area was small, mainly distributed in Jiaozuo City, Xinxiang City, Zhoukou City, and Puyang City, in a patchy distribution.

Ionic Change Pattern
The evolutionary process of hydrochemistry is related to groundwater runoff conditions and to the burial depth of the groundwater table [42]. Along the direction of groundwater runoff, the concentration of major groundwater anions and cations will generally increase as the groundwater runoff path grows and is subject to more persistent water-rock interactions. Moreover, the shallower the groundwater burial depth is, the larger the proportion of Cl − and SO 4 2− that will be affected by evaporation concentration, and these will even become the dominant ions. In this study, two typical profiles were selected, and Stiff diagrams were drawn to investigate the water chemistry change pattern according to the groundwater flow direction and burial depth distribution. The locations of the two water chemistry profiles are shown in Figure 5, and the corresponding Stiff diagrams are shown in Figure 6. P1 is the profile along the burial depth of groundwater level from shallow to deep, from JC004→JC028→JC072→JC091; P2 is the profile along the direction of groundwater runoff, from JC012→JC046→JC069→JC029. From P1 in Figure 7, it can be seen that the anion content in groundwater increased with the burial depth from shallow to deep, and the HCO 3 − content increased significantly, indicating that the effect of evaporation and concentration on HCO 3 − gradually decreased. The content of Cl − and SO 4 2− did not decrease, but increased slightly, which is obviously not in accordance with the change in water chemistry with the burial depth of groundwater under natural conditions, and it is presumed that the content of Cl − and SO 4 2− is influenced by stronger human activities. As shown by P2, along the direction of groundwater runoff, Na + showed an upward trend in general, but at the middle and lower reaches, it initially decreased and then increased; Ca 2+ had the same overall trend as Na + , but at the middle and lower reaches, there was an opposite trend to Na + , which might be due to the cation exchange adsorption of Na + and Ca 2+ in this region; furthermore, the change inMg 2+ content was basically the same, with no obvious rise and fall characteristics, which may be influenced by the use of magnesium fertilizer in agricultural production activities. It is obvious that the changes in groundwater chemistry in the study area do not fully conform to the natural laws, and to some extent, may also be influenced by stronger human production activities.

Analysis of Ion Sources
Correlation analysis between ions can, to some extent, reveal the causal relationship between ions [43]. To explore the source of ions, correlation analysis was performed on the chemical indicators corresponding to the 76 samples (Table 4). It can be seen that Na + has a moderate positive correlation with Cl − , SO 4 2− , F − , and HCO 3 − , indicating that Na + may have the same material source as all these anions. Na + in the study area shows a positive correlation with Cl − , which may be related to the dissolution of rock salt (NaCl) (Equation (8)). The evolutionary relationship of the groundwater phase can be inferred from the proportional relationship between ions, and as shown in Figure 7a, the scatter of Na + and Cl − was basically distributed on a 1:1 straight line, but most of it lay above the straight line. The amount of scattered Na + located above the straight line was larger than that of Cl − , indicating that the source of Na + is not only from the dissolution of rock salt but also may be related to the weathering dissolution of albite (Equation (9)) and the presence of mannite (Na 2 SO 4 ) in the inclusion zone (Equation (10)), and may also be related to the occurrence of cation exchange adsorption of Mg 2+ and Ca 2+ . Na + , SO 4 2− , and HCO 3 − show a stronger positive correlation, corroborating part of the above speculation. A small portion of the scatter lies below the 1:1 straight line, and this extra Cl − may come from anthropogenic pollution. Na + is positively correlated with F − , but F − is negatively correlated with Ca 2+ , suggesting that the formation of F − is related to the dissolved precipitation of fluorescent stone (CaF 2 ) (Equation (11)). NaCl→Na + + Cl −   4 2− , which may be related to the use of Mg fertilizer in agricultural production in the study area. The chemical reaction equations of calcite and dolomite reflect the ratio between dissolved ions (Equations (12) and (13) The above chemical reaction equation shows that the ionic ratio of Ca 2+ , Mg 2+ , and Ca 2+ to HCO 3 − from the dissolution of calcite and dolomite is between the 1:1 and 1:2 straight lines, and its magnitude depends on the amount of CO 2 involved in the reaction. As displayed in Figure 7b, a small portion of the scatter lay above the 1:1 straight line, indicating that a small portion of the Ca 2+ in the groundwater may have come from the dissolution of dolomite or anorthite; there was also a portion of the scatter between the 1:1 and 1:2 straight lines, and this portion of the Ca 2+ may have come from the dissolution of calcite. Most of the scatter lay below the 1:2 straight line, and the decrease in Ca 2+ in the groundwater confirms that the increase in Na + in Figure 7a may be related to cation exchange sorption. As shown in Figure 7c, a small sector of the scatter lay above the 1:1 straight line, which is similar to the distribution in Figure 7b, indicating the presence of other sources of Ca 2+ and Mg 2+ in a small portion of groundwater in the study area; a significant portion of the scatter lay between the 1:1 and 1:2 straight lines, indicating that dissolution of dolomite is an important source of Ca2 + and Mg 2+ in groundwater in the study area. Figure 7d shows the ionic ratio between Ca 2+ + Mg 2+ , and HCO 3 − + SO 4 2− , with most of the scatter located along the 1:1 straight line, and a small portion of the scatter located above the straight line, presumably due to excess Mg 2+ , which confirms the above point about magnesium fertilizer use. A significant portion of the scatter lies below the straight line, further suggesting that the reduction in Ca 2+ may be caused by cation exchange.
Nitrate pollution generally results from the heavy use of nitrogen fertilizers, large discharge of domestic and industrial effluents, accumulation of solid wastes and animal manure, and the potent role of human activities in the distribution of nitrate [44]. For this reason, in this study, to investigate the main sources of nitrate pollution in the study area, the distribution of construction land in the study area was spatially overlaid with the distribution of nitrate nitrogen concentration ( Figure 8). As shown in Figure 8, nitrate nitrogen was basically distributed at higher concentrations in places where construction land was more concentrated, and nitrate nitrogen content was lower in areas where construction land was sparser, such as the southern part of Xinyang City. The migration of Fe and Mn ions in groundwater is generally controlled by the composition of the water-bearing medium, runoff conditions, nature of the overlying soils, acid and alkaline conditions, and reducing environment [45], and because the properties of Fe and Mn are very similar, they often coexist in groundwater. The distribution of Fe and Mn concentrations in the study area was generally similar, except in Puyang City, Pingdingshan City, Xinxiang City, and the border between Zhoukou and Kaifeng City. These regions all had a high concentration of Fe 3+ , presumably due to industrial wastewater discharge; therefore, the total number of industrial sources in each city of the study area was counted in this study (Figure 9), and the top three cities in the order of the number of industrial sources were Jiaozuo City > Zhengzhou City > Puyang City > (Kaifeng City + Zhoukou City). Thus, it is clear that the regions with high Fe 3+ concentration in the study area all had a considerable number of industrial pollution sources. Although the number of industrial pollution sources in Pingdingshan City was relatively small, the city of Wugang, which is known as the "capital of iron and steel," has a large amount of Fe 3+ in the wastewater produced by iron and steel making. Arsenic in groundwater generally originates from the release of sediments produced by the weathering of rocks and can also arise from the discharge of industrial wastewater. The strongest correlation of As with Fe 3+ and the most similar distribution with HCO 3 − in the study area indicate that the strong reducing environment limits the formation of iron oxides and indirectly attenuates the adsorption of iron oxide on As, leading to the elevated As content [46].

Mechanisms of Hydrochemical Formation
The Gibbs diagram allows the classification of factors affecting the chemical composition of natural water into evaporative concentration, water-rock interaction, and atmospheric precipitation [47]. According to the Gibbs diagram (Figure 10), most of the groundwater samples in the study area were located in the water-rock interaction region, and very few samples were located in the evaporative concentration region. The relationship between TDS and Na/(Na + Ca) shows that as the value of Na/(Na + Ca) increased, TDS did not change much, further indicating the presence of cation exchange adsorption during groundwater runoff. The trend in the distribution of the samples shows that the value of Cl/(Cl + HCO 3 ) tended to increase with increasing TDS but remained within the range of water-rock interaction, which also indicates the presence of evaporative concentration of the shallow groundwater chemistry components in the study area with less influence. Given that the main formation mechanism of groundwater chemistry in the study area is water-rock interaction, and the main expression of water-rock interaction is the dissolution and precipitation of minerals, this study determined the water-rock interaction state of groundwater in the study area by calculating the saturation index of each sample relative to calcite, dolomite, fluorite, and rock salt. Figure 11 shows the relationship between the saturation index of groundwater samples relative to each mineral and TDS.
As shown in Figure 11, the saturation indices of groundwater samples relative to calcite and dolomite in the study area were mostly greater than 0, except for very few that were less than 0. The saturation indices ranged from −0.34 to 1.42 and from −0.79 to 2.61, respectively, indicating that precipitation of calcite and dolomite occurs in most areas of the study area. The saturation index scatter of groundwater samples relative to fluorite and rock salt all fell below the SI = 0 axis between −0.25 and −3.71 and from −8.03 to −5.6, respectively, suggesting that weathering dissolution of rock salt and fluorite occurs in groundwater in the study area. With the increase in TDS, each saturation index scatter showed an increasing trend, and the trend of fluorescent stone change was especially obvious, implying that rock weathering dissolution is an important factor affecting TDS in the study area.     Figure 11. Scatter plot of saturation index with TDS.

Discussions
After assessing the groundwater quality in the study area and analyzing the water chemistry distribution characteristics, change patterns, ion sources, and formation mechanisms, it can be found that the overall water quality in the study area is good, and the areas with poor water quality are mainly distributed locally in Jiaozuo City, Xinxiang City, Zhoukou City, and Puyang City. Most regions with poor water quality have the characteristics of concentrated emission of pollution sources, strong human activities, or long runoff paths. From the perspective of anions, the hydrochemical type in the study area is mainly HCO 3 type, and the ion content increases gradually along the runoff path. This is because HCO 3 − is mainly derived from the dissolution of carbonate rocks, the material source is sufficient, and this process is further enhanced by the dissolution of CO 2 produced by microbial respiration or in the air with atmospheric rainfall or river runoff [48]. Cl − and SO 4 2− in local groundwater are dominant ions, which are closely related to the weathering and dissolution of rock salt and mirabilite and may also be related to the pollution caused by human activities. From the perspective of cations, the distribution of hydrochemical types has a certain degree of zonation such that the Ca/Mg type and Ca-Mg/Mg-Ca type are dominant in the piedmont sloping plain, and its formation mechanism is related to the weathering and dissolution of calcite and dolomite. In the eastern plain area, the hydrochemical type began to change into the Na-Ca/Ca-Na type, Na-Mg/Mg-Na type, and finally changed into the Na type, indicating that a strong ion exchange occurred in the process of groundwater runoff, which was verified by the distribution of Ca 2+ and Na + in concentration, and these conclusions are similar to the early research results in the region [49,50]. The long history of human activities in Henan Province has made great changes in the hydrochemistry of the region, which is mainly divided into intuitive and indirect effects. It is more intuitive that K + , NO 3 -N, and Fe 3+ have abnormally high concentrations in some cities due to the use of fertilizers and sewage discharge, or the extensive use of dry toilets due to climatic influences. The indirect effect does not directly change the component content, but affects the hydrodynamic field through irrigation, mining, and other ways, thus affecting the evaporation concentration or water-rock interaction and resulting in some component concentration anomalies. Therefore, the traditional Gibbs graph cannot meet the research on the formation mechanism of groundwater chemistry under the influence of human beings, and some studies have shown that the Gibbs graph oversimplifies the interpretation of the aquifer system and ignores important processes in defining groundwater geochemistry [51]. A variety of methods are needed to have a more comprehensive understanding of groundwater chemical field changes.

Conclusions
Shallow groundwater is an important source of freshwater in Henan Province, China. In this study, we assessed the quality of shallow groundwater in the plain area of Henan Province based on a combination of subjective and objective empowerment and analyzed the groundwater hydrochemical characteristics and formation mechanisms, the following conclusions can be drawn: (1) A combination of the entropy weighting and subjective assignment methods was used to determine the relative weights of each evaluation index not only reflecting the data characteristics of the evaluation index itself but also to some extent combining the subjective judgment of the importance of the evaluation index. According to the results of groundwater quality assessment, the water quality in most areas of the study area is above medium to good, and the water quality is poor in areas with long runoff paths or strong human activities, mainly distributed locally in Jiaozuo City, Xinxiang City, Zhoukou City, and Puyang City, in a patchy distribution. (2) The predominance of HCO 3 type anions in the groundwater types in the study area indicates that carbonate dissolution plays a dominant role in regulating hydrochemistry. The presence of ions dominated by Cl − and SO 4 2− in local areas indicates that the weathering dissolution of rock salt and mannite as well as human activities also influence the hydrochemistry components simultaneously. From the perspective of cations, the distribution of hydrochemistry types has a certain zonation, and the precipitation of calcite and dolomite, the weathering dissolution of rock salt and fluorite, and the cation exchange and adsorption are the main factors affecting the hydrochemistry distribution of cations.
(3) Human activities can directly or indirectly affect changes in the chemical components of groundwater. For example, inappropriate behaviors such as overuse of fertilizers, industrial and agricultural wastewater discharge, manure accumulation, groundwater over-exploitation, and unreasonable irrigation can lead to the abnormal distribution of certain ions and deteriorate the pristine groundwater environment. Therefore, it is necessary to continuously monitor groundwater to provide scientific data for sustainable groundwater quality management and to develop corresponding measures to prevent further degradation of the groundwater environment. Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare that they have no known competing financial interests or personal relationships that could influence the work reported herein.