Next Article in Journal
Dam Renovation to Prolong Reservoir Life and Mitigate Dam Impacts
Next Article in Special Issue
Uncertainty Analysis of Numerical Simulation of Seawater Intrusion Using Deep Learning-Based Surrogate Model
Previous Article in Journal
Short-Range Movement Pattern of Amphidromous Lagoon Fish Schools: Ecological Applications
Previous Article in Special Issue
Ecological and Health Risk Assessment of Heavy Metals in Soils from Recycled Lead Smelting Sites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrochemical Characteristics and Hydrogeochemical Simulation Research of Groundwater in the Guohe River Basin (Henan Section)

1
College of Geosciences and Engineering, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
2
Collaborative Innovation Center for Efficient Utilization of Water Resources, Zhengzhou 450046, China
3
Sinopec North China Oil and Gas Branch, Yulin 719000, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(9), 1461; https://doi.org/10.3390/w14091461
Submission received: 29 March 2022 / Revised: 28 April 2022 / Accepted: 28 April 2022 / Published: 3 May 2022

Abstract

:
With the implementation of the policy of ecological protection and high-quality development of the Yellow River Basin, the Guohe River Basin, which is close to the middle reaches of the Yellow River Basin, is also an important part of future development. Mathematical statistics, the Piper diagram, ion proportion coefficient method, Gibbs diagram and reverse hydrogeochemical simulation are used to analyze the chemical characteristics and evolution of groundwater in the Guohe River Basin (Henan Section). The dominant ions in the study area are HCO3 and Na+, and the three-layer aquifer has obvious zoning characteristics. The results show that the chemical types of shallow groundwater in this area are HCO3−Na • Mg • Ca, intermediate HCO3-Na and deep HCO3−Na. Using the ion proportion coefficient method, it is found that Na+, Ca2+, and Mg2+ in the groundwater aquifer undergo cation exchange in the aquifer. According to the reverse hydrogeochemical simulation, gypsum in the three aquifers is in a dissolved state, carbonate and sulfide ores in the shallow layer are dissolved, dolomite and halite in the intermediate layer are dissolved, calcite and sulfide ores are precipitated and carbonate, halite and sulfide ores in the deep layer are precipitated; the hydrogeochemical evolution model is established to find that Ca2+ in groundwater displaces Na+ in the aqueous medium. This research can provide a scientific basis for the rational development and utilization of groundwater and ecological protection in the Yellow River Basin.

1. Introduction

The chemical composition of groundwater is mainly affected by geological conditions, the extent of chemical weathering of rocks, characteristics of aquifers and human activities [1,2]. Its formation and evolution are the long-term result of groundwater diagenesis. With the circulation of groundwater, the process is more complex and changeable [3]. The quality of groundwater in an area depends on the chemical composition of rock stratum and make-up water, the interaction between soil and water, the rock present in unsaturated areas and its interaction with make-up water, the recharge time of the aquifer and the action of the aquifer itself [4,5]. The definition of hydrogeochemistry was proposed in 1979 as a science with which to study various chemical interactions and physical phenomena between the hydrosphere and its surrounding natural environment [6]. In 1993, Shen Zhaoli et al. [7] compiled the basis of hydrogeochemistry and established the discipline system of hydrogeochemistry.
The relevant information of groundwater circulation characteristics from different angles can be acquired by studying the chemical composition, and distribution characteristics of groundwater can be obtained [8,9]. Lloyd, S. Adams [5,10] and others used a variety of conventional statistical methods (such as mathematical statistical analysis, factor analysis, etc.) to conduct mathematical statistical analysis and interpretation. Hitchon et al. [11], Gupta et al. [12] and S. Kumar [13] added multivariate analysis to the interpretation of groundwater chemistry. After decades of development, hydrogeochemical theory has been continuously improved, and research methods such as ion proportion coefficient method and water rock interaction type have been continuously enriched [14,15]. Rabelani M. [16] studied the hydrogeochemical characteristics of groundwater in Grand Giani City, Limpopo Province, South Africa using a Piper diagram, Gibbs diagram and ion ratio. Kaur L. et al. [17] investigated the chemical composition of groundwater in the Panipa special district (semi-arid alluvial area) of Haryana, India by using a Piper diagram and ion proportion relationship. N.M. Refat Nasher [18] used factor analysis and ion ratio to analyze the hydrogeochemical characteristics of groundwater in three basins in Bangladesh. The depth of ion ratio analysis is adopted to describe the interpretation of calcite and dolomite to the traditional hydrogeochemical concept [19]. Remy R. [20] used groundwater and surface water samples collected in different seasons to explain the hydrogeochemical process of Kattumannarkoil Taluk. Aher S. et al. [21] explored the hydrogeochemical characteristics of the Pravara River Basin using groundwater samples collected on the cross-section of the basin. Accurate assessment of groundwater hydrogeochemical characteristics and groundwater quality is helpful and necessary for sustainable groundwater protection and management [22,23].
In recent years, scholars have used isotope analysis and hydrogeochemical simulation methods to study hydrochemical characteristics to good effect [24,25,26]. Ren et al. [27] found that groundwater in any area has a unique chemical composition, which is caused by many processes such as the interaction between water and rock during recharge and groundwater flow through the aquifer and the long-term water storage of the aquifer. Wei et al. [28], Subba Rao et al. [29] and Wu et al. [30,31] aimed to delineate the process responsible for the geochemical changes of aquifer water environment. Prasanna et al. [32] used the method of digital system analysis to calculate the saturation index of groundwater. Parkhurst et al. [33] began the classification of mineral saturation index in groundwater samples using PHREEQC in 2013.
Through hydrochemical analysis of the shallow groundwater in the Guohe River Basin (Anhui section), Ma Tao et al. [34] found that the shallow groundwater in the study area is susceptible to human activities and atmospheric rainfall. Zheng Tao et al. [35] investigated the chemical characteristics and genetic mechanism of groundwater in the central area of the Guohe River Basin. They found that there are significant vertical differences in the chemical composition of groundwater, and the formation of groundwater hydrochemical characteristics is affected by water–rock interaction, cation exchange and human activities. At present, there are few studies on hydrochemistry that have focused on the upper reaches of the Guohe River Basin. In the context of ecological protection and high-quality development in the Yellow River Basin, the research aims to determine the chemical characteristics of shallow, mid-deep and deep groundwater in the Guohe River Basin (Henan Section), and to identify its salient hydrogeochemical formation, assess the transport of groundwater components through reverse hydrogeochemical simulation, quantitatively study the chemical components of groundwater using the hydrogeochemical evolution model and reveal the chemical formation process of groundwater and the spatial transport trends in its groundwater hydrochemistry. The results may provide a scientific basis for rational development and use of groundwater and ecological protection therein.

2. Materials and Methods

2.1. Regional General Situation

The study area is located in the east of Henan Province, in the upper reaches of the Guohe River Basin, and forms an important part of the Huaihe River Basin, covering an area of about 14,000 km2, with geographic coordinates of 113°22′–115°39′ E and 33°22′–34°50′ N. The Guohe River Basin is located in the warm temperate zone and has obvious monsoonal climate characteristics. According to the precipitation data from meteorological stations in Kaifeng, Shangqiu and other places, the annual average precipitation over the basin is 638 to 690 mm, and the trend of precipitation increases slightly toward the southeast, consistent with the direction of the Guohe River. The humidity is the highest from June to September, and the annual average temperature is 13 to 15 ℃. The evaporation intensity decreases from west to east, and the evaporation is the strongest from May to June. The lithology of the Pre-Cenozoic Middle Jixian System is mainly metamorphic sandstone and schist. The Cambrian lithology is mainly limestone and dolomite, with a thickness of 21.2 to 218.29 m. The Ordovician-Cambrian is mainly composed of limestone and dolomite with a total thickness of more than 400 m. The main lithology of the Cenozoic Paleogene is claystone, sandy claystone, etc., and sand layers are distributed alternately therein. In some areas, there are sand-argillaceous limestones, and the strata are in angular unconformity contact with the Cambrian and Ordovician strata.

2.2. Sample Collection and Research Methods

The project team collected 180 sets of groundwater samples in the study area from September to October 2020, including 108 sets of shallow groundwater, 34 sets of mid-deep groundwater and 38 sets of deep groundwater (Figure 1). The water samples were sent to the Experimental Testing Center of the Natural Resources Monitoring Institute of Henan Province for testing. The sampling data analysis included seven constant inorganic components (K+, Na+, Ca2+, Mg2+, Cl, SO42−, and HCO3), as well as pH, water temperature, COD, total hardness and TDS. The assay for Na+ relied upon atomic absorption spectrophotometry. The tests for Ca2+, Mg2+, Cl and SO42− relied upon ion chromatography; HCO3 was determined by on-site titration; WTW portable multi parameter water quality analyzers were used to measure the TDS, water temperature and pH value of spring water. The measurement accuracy is 1 mg/L, 0.1 ℃ and 0.01 pH units, respectively; COD was measured by the three-position fluorescence method.
The geological data of the study area show that the boundary between shallow and middle deep groundwater in the Guohe River Basin is set at about 50 m, at 50 to 300 m the middle aquifer (the Kaifeng control depth is about 400 m) is found, and the deep control depth is between 300 and 600 m. Shallow groundwater samples come from rural decentralized water supply wells and farmland irrigation wells, with a depth of 8 to 50 m; the intermediate groundwater comes from rural centralized water supply wells with a depth of 50 to 300 m; the deep groundwater comes from the water supply well serving the urban waterworks, with a depth of more than 300 m.
Mathematical statistics could be used to study the variational characteristics of the coefficient of variation and salinity of the three-layer aquifers in the study area. The main source of ions was explored by using a Gibbs diagram to determine the hydrolithic effect of ions in groundwater. Depending on the water chemistry, TDS partition and hydrogeological conditions, using PHREEQC software (Parkhurst and others, USDS), the research calculated mineral saturation index (SI) using the reverse hydrogeochemical simulation, simulated the shallow, intermediate, and deep groundwater migration paths, traced dissolved mineral phases as they underwent precipitation and mapped the Guohe River Basin along the flow profile on the hydrogeochemical evolution model, thus revealing the hydrogeochemical formation process of groundwater in the South of the Guohe River Basin.

2.3. Statistical Analysis

Descriptive statistics provide a general view of the hydrochemical characteristics [36]. To prevent subjective and objective errors in the test, the cation anion balance test method is adopted here, and the test formula is determined as follows:
E = Σ m c Σ m a Σ m c + Σ m a × 100 %
where E denotes the relative error, %; mc and ma are the milligram equivalent concentrations of cations and anions, meq/L.
If the relative error is ± 5%, it indicates that the data are valid and reliable. Substituting the measured data into Formula (1), the collected 108 groups of shallow groundwater, 34 groups of intermediate groundwater and 38 groups of deep groundwater meet the relative error criterion of ± 5%; that is, they satisfy the prerequisite conditions.

3. Results and Discussion

3.1. Analysis of Groundwater Chemical Types

3.1.1. Statistical Characteristics of Main Indicators: Analysis of Hydrochemical Types and Characteristics

A mathematical statistical analysis was conducted on the groundwater sample data in was performed (Table 1), and the average content distribution diagram of major ions and indicators was drawn (Figure 2). The coefficient of variation (C.V.) can represent the dispersion degree of various indexes in groundwater and the complexity of influencing factors of groundwater chemical composition and its evolution [37]. The coefficient of variation of K+, Cl and SO42− in the shallow layer is large, indicating that the shallow layer is strongly affected by human activities. It can be seen that the coefficient of variation of Na+, Ca2+ and Mg2+ in the middle and deep layer is large. This phenomenon may be affected by the recharge of shallow groundwater or cation exchange adsorption. From shallow to intermediate to deep, the groundwater salinity in the study area presents a trend of high to low to low. This zonation indicates that the vertical component difference of the three aquifers may be caused by the surrounding rock environment and the intensity of water–rock interaction. Meanwhile, the salinity of the shallow layer is generally higher than that of the intermediate layer, which is possibly influenced by human activities and farmland irrigation. The results of this study differ from those of Zheng Tao et al. [35], wherein an analysis of the low → high → low trend of groundwater stratification mineralization in The Guohe River Basin (Anhui Section) is revealed.

3.1.2. Analysis of Hydrochemical Types and Characteristics

At present, many scholars use Piper three-line diagrams to analyze the ion composition and variation characteristics of surface water and groundwater [38,39]. D. Gamvroula and D et al. [40] used this method to explain in detail the distribution and characteristics of macronutrient ions in the groundwater of Megara Basin. Piper tri-maps of shallow, intermediate, and deep layers in the study area were drawn (Figure 3). The results showed that Na+, Ca2+ and Mg2+ ions are abundant and evenly distributed in shallow groundwater, and the hydrochemical types are mainly HCO3-Na•Mg•Ca, accounting for 35.19% overall. The hydrochemical types of some sampling sites are Cl•HCO3•SO4-Na•Mg and HCO3-Ca•Mg, accounting for 14.82% overall. HCO3 is abundant in the intermediate groundwater, and the main cations are Na+. The hydrochemical type of HCO3-Na is 44.12% overall. Some are of the HCO3•Cl•SO4-Na type, accounting for 17.65% overall, and the proportion of HCO3•Cl-Na type is 14.71%. The type of deep groundwater is similar to that of middle deep groundwater, mainly HCO3-Na type, accounting for 60.53% overall with some HCO3•Cl-Na type present, accounting for 23.68% overall.
Along the direction of groundwater flow, the order of cation contents in shallow groundwater in the upper reaches of the study area is Ca2+ > Mg2+ > Na+, the order of cation contents in the lower reaches is Na+ > Mg2+ > Ca2+, and the order of anions is HCO3 > SO42− > Cl. The order of anionic contents in the intermediate groundwater is HCO3 > Cl> SO42−, and the order of cations is Na+ > Ca2+ > Mg2+. The order of cations in the deep groundwater is different from that in the other two layers, so the exchange effect of the deep groundwater is weak. Meanwhile, HCO3 and Na+ are found to play a dominant role in groundwater with increasing depth, and their basic hydrochemical type is characterized by the transformation from HCO3-Na•Mg•Ca type in shallow groundwater to HCO3-Na type at intermediate depth, and then to HCO3-Na type in deep groundwater, with zoning characteristics.

3.2. Analysis of Water-Rock Interaction

3.2.1. Ion Proportion Relation and Component Origin Analysis

γNa+/γCl is known as groundwater genetic coefficient, which can be used to reflect the sources of Na+ and Cl in groundwater [7]. The leachate value of halite-bearing strata is about 1 [41]. According to the proportion coefficient of water ions in the study area, Figure 4a shows that the ratio of γNa+/γCl in groundwater is mostly greater than 1, indicating that Na+ is not completely derived from the dissolution of halite, but also the dissolution of other sodium-bearing ores. Meanwhile, it is found that γNa+/γCl is present as follows: deep > intermediate > shallow layers. This indicates that the deep and intermediate layers are relatively closed. It can be seen in Figure 4b that HCO3 is not completely derived from the dissolution of calcite. The shallow, intermediate and deep underground water samples in the study area are distributed on the side of the isoratio line that is biased towards HCO3, indicating that the alternating adsorption of Na+ and Ca2+ cations or desulfuration may occur in the study area, which reduces Ca2+ and increases HCO3. Figure 4c indicates that the shallow groundwater and intermediate groundwater are distributed on both sides of the contour line, indicating that gypsum leaching may have occurred. The deep groundwater is all below the contour line, biased to the side of SO42−, suggesting that cation alternating adsorption may have occurred to decrease Ca2+ and increase Na+ content. The deep groundwater is all below the isoratio line, which is on the side of SO42−. In addition to the cation alternation adsorption, which can reduce Ca2+, the dissolution of other sulfate minerals such as mirabilite, or even the contribution of sulfide oxidation to SO42− content cannot be ruled out. Figure 4d shows a positive correlation between Ca2++Mg2+ and HCO3 in the shallow layer, indicating calcite and dolomite solubilization. Some of the shallow groundwater is distributed on one side of Ca2++Mg2+, suggesting that the solubilization of other calcium and magnesium minerals, such as gypsum and tremolite, may have occurred. The groundwater in the intermediate layer is inclined to HCO3, implying that the alternate adsorption of cations may occur in the study area. According to Figure 4e, water samples from intermediate layers in the study area have γ (Ca2++Mg2+)/γ (HCO3 + SO42−) ratios of less than 1, indicating that the main sources of Ca2+ and Mg2+ are the dissolution of evaporite and silicate minerals, suggesting that the negative charges of HCO3 and SO42− are balanced by another cation such as Na+. HCO3 and SO42− may be controlled by a series of hydrochemical actions such as silicate hydrolysis, sulfide oxidation and cation-alternating adsorption, especially in intermediate and deep groundwater. The obvious enrichment characteristics of HCO3 and SO42− warrant further discussion. Figure 4d and e show that there is significant enrichment of Ca2+ and Mg2+ in the shallow layer. This enrichment is probably affected by human activities, but the enrichment of Ca2+ and Mg2+ in the shallow layer caused by the evaporation concentration cannot be ruled out, which also warrants further discussion.
Through [γ(Ca2++Mg2+)-γ(SO42−+HCO3)]/γ(Na+-Cl) ratio analysis of the cation exchange, adsorption is found to occur; γ(Na+-Cl) represents the change of Na+ content caused by other sodium-bearing minerals except halite dissolution. [γ(Ca2++Mg2+)-γ(SO42− +HCO3)] represents the changes of Ca2+ and Mg2+ contents caused by other substances other than gypsum, dolomite, and calcite. As shown in Figure 4f, the slope in shallow groundwater is −0.925, the slope in intermediate groundwater is −1.032, and the slope in deep groundwater is −0.984 (all close to −1), indicating that the change is linear, and that cation exchange adsorption has occurred at the spatial scale of groundwater.

3.2.2. Gibbs Diagrams and Their Implications

The chemical natural formation mechanism of groundwater mainly includes atmospheric precipitation, water-rock interaction and evaporation concentration [42,43]. The Gibbs diagrams of shallow, intermediate and deep groundwater in the study area demonstrate that the shallow groundwater body in the study area is jointly affected by water–rock interaction and evaporation (albeit mainly by water–rock interaction); the sampling points in the intermediate and deep groundwater are situated in the middle area. The chemistry of intermediate groundwater is controlled by water–rock interaction; that is, the formation mechanism of groundwater runs from shallow to deep, and weathering of the rock is continuously intensified (Figure 5).

3.3. Reverse Hydrogeochemical Simulation

In the reverse hydrogeochemical simulation, the selection of appropriate simulation paths and possible mineral phases is the key to the establishment of the model, so as to determine the main SI and the transformation amount of each mineral phase on the simulation path, so as to study the hydrogeochemical formation [44]. According to the hydrogeological conditions prevailing across the study area [45,46,47], calcite (CaCO3), dolomite (CaMg (CO3)2), gypsum (CaSO4•2H2O), halite (NaCl), CO2, O2, H2S, NaX and CaX2 were selected as the reaction mineral phases.

3.3.1. Mineral Saturation Index

The content of CO2 is used as a proxy for the degree of dissolution of carbonate; H2S and O2 can distinguish the degree of dissolution of pyrite, tremolite and other sulfur-bearing minerals; according to the Gibbs diagram analysis results, calcite and dolomite were added as the main reactive mineral phases. NaX and CaX2 are added according to results of the ion proportion coefficient method. The mineral SI indicates the saturation state of minerals relative to groundwater. The SI values at some representative points in the study area are listed in Table 2.
When SI > 0, it indicates that the mineral facies is in a supersaturated state and tends to be precipitated in the groundwater aquifer (otherwise, it will be further dissolved). The results show that the groundwater aquifers in the Guohe River Basin (Henan Section) are divided into distinct zones, calcite and dolomite in shallow and deep groundwater is probably precipitated, and calcite and dolomite in intermediate aquifers may be further dissolved in groundwater. Gypsum and halite always maintain their ability to dissolve further at the vertical spatial scale of the groundwater profile.

3.3.2. Mineral Saturation Index

According to the TDS, isoline is the main factor in the selection of reverse hydrogeochemical simulation paths; the shallow, intermediate and deep simulation paths in the study area are shown in Figure 6.
The measured pH, water temperature and major ion components (Table 3) are applied to PHREEQC for reverse hydrogeochemical simulation, and the results are summarized in Table 4.
The shallow simulation path (W-Q06→W-Q51) mainly involves the dissolution of calcite, dolomite, gypsum and sulfur-bearing minerals. In the simulation path for the intermediate layers, gypsum, salt and dolomite dissolution and precipitation of calcite and sulfide ores are the main factors influencing the hydrogeochemistry. From the simulation starting point (W-Z02) to its end (W-Z36), the contents of Na+, Cl and SO42− increase significantly; the contents of Ca2+ and Mg2+ rise slightly, and the contents of other ions change little. The dissolution of gypsum and the precipitation of calcite, dolomite, halite and sulfide ores occur in the deep simulation path. The contents of Na+, HCO3 and SO42− in the water grow from the simulated starting point (W-S08) to its end (W-S30). The dissolution of NaX and precipitation of CaX2 in the three-layer groundwater simulation path reveal that the cation exchange reaction between Ca2+ and Na+ is common in the aquifer, which also verifies the conclusion obtained by the ion proportional coefficient method. From the vertical space, the hydrochemical type changes from shallow HCO3-Na•Mg•Ca type to intermediate HCO3-Na type, and then to the deep groundwater dominated by almost all HCO3-Na type; Na+ and HCO3 gradually become the dominant ions in groundwater. The proportions of Ca2+, Mg2+, Cl and SO42− in the water decrease with increasing depth. The hydrogeochemical formation of groundwater in the Guohe River Basin (south section of the River) is quantitatively revealed by the amount of mineral dissolved deposits calculated by reverse hydrogeochemical simulation.

3.3.3. Analysis of the Hydrogeochemical Evolution Model

The hydrogeochemical evolution model on spatial scale was developed by combining geological and hydrogeological conditions, hydrogeochemical processes and main anions in groundwater at the time of groundwater hydration. To reveal the horizontal transport of anions and cations along the groundwater flow path, and through the changes in their relative concentrations with depth, the changes of main anions and cations in the vertical spatial scale are analyzed, and the spatial transport associated with the prevailing groundwater hydrochemistry is studied (Figure 7).
The results show that shallow groundwater is susceptible to atmospheric precipitation, human activities and evaporation and concentration, resulting in strong fluctuation of cation and anion concentrations therein. The concentrations of anions and cations in the intermediate layer are generally increasing, which is mainly caused by leaching and cation alternation adsorption, in a manner consistent with the reverse hydrogeochemical simulation results. Na+ and HCO3 in the deep layer augment significantly along the flow direction, indicating that the deep layer is affected by leaching and cation alternating adsorption.
According to the analysis in Section 3.1, the Guohe River Basin (Henan Section) has obvious zoning characteristics. The hydrochemical types of the shallow layer along the direction of groundwater flow vary, indicating that the shallow layer is affected by evaporation and concentration, and the intermediate and deep layers gradually change to HCO3-Na type along the direction of groundwater flow (Figure 7). Leaching and cation alternating adsorption are the main hydrogeochemical processes affecting the three-layer aquifer in the basin.

4. Conclusions

The hydrochemical types of groundwater in this river basin (Henan Section) are as follows: an HCO3-Na•Mg•Ca type predominates in the shallow layer, an HCO3-Na type predominates in intermediate and deep layers and the groundwater has zoning characteristics. According to Gibbs diagrams, the lithogenesis of the three aquifers becomes more intense with increasing depth in the groundwater, and the weathering of the rock is gradually intensified. The shallow layer is affected by water–rock interaction, evaporation and concentration, while the intermediate and deep layers are mainly affected by the action of water on the rock.
The ion proportion coefficient method shows that Na+ in groundwater is not completely derived from the dissolution of halite, and there may be dissolution of other sodium-bearing ores; Ca2+ and Mg2+ are mainly caused by the dissolution of evaporite and silicate minerals. The alternation adsorption of cations occurred in each aquifer, which reduced the content of Ca2+ and Mg2+ and increased the content of Na+ in the solution.
PHREEQC shows that the SI index shows a general trend, such that SI > 0 and the groundwater is in a supersaturated state. Its shallow path is characterized by the dissolution of calcite, dolomite, gypsum and sulfur-bearing minerals. In the deep path, dolomite, gypsum and halite were dissolved, while calcite and sulfur minerals (pyrite, etc.) are precipitated. The deep path is characterized by the dissolution of gypsum and the precipitation of calcite, dolomite, halite, sulfur minerals and other minerals. The results are consistent with the calculated mineral SI data, indicating that water–rock interaction plays a major role in the geochemistry of these groundwater aquifers. Meanwhile, the NaX and CaX2 dissolved precipitates of the three aquifers confirm quantitatively the cation-alternation adsorption between Ca2+ and Na+.
According to the hydrogeochemical evolution model, the dissolution and precipitation of mineral facies exert a significant influence on the evolution of groundwater chemistry. Along the direction of groundwater flow, the aquifers show an increasing trend of HCO3 and Na+ concentration, and in the vertical scale, the hydrochemical type is stable as HCO3-Na along the direction of groundwater flow. The influence of cation alternation adsorption is significant, which verifies the trend in transformation of hydrochemistry type with increasing depth and qualitatively reveals the formation of the prevailing hydrogeochemistry in the Guohe River Basin (Henan Section). The ion concentrations in shallow groundwater fluctuate, and there is a small deviation from the results predicted by mineral simulation, which may be caused by human activities, evaporation and concentration.

Author Contributions

Conceptualization, F.Y.; methodology, F.Y.; investigation, F.Y. and D.Z.; data curation, X.L.; writing—original draft preparation, D.Z.; writing—review and editing, F.Y. and D.Z.; supervision, Z.L.; funding acquisition, F.Y. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos 41402225 and 41972261), and the Key Scientific Research in Universities in Henan Province (Grant No. 21A170014).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Domenico, P.A. Concepts and Models in Groundwater Hydrology; McGraw-Hill: New York, NY, USA, 1972; pp. 1–105. [Google Scholar]
  2. Schuh, W.M.; Klinkebiel, D.L.; Gardner, J.C. Tracer and nitrate movement to groundwater in the northern great plains. J. Environ. Qual. 1997, 26, 1335–1347. [Google Scholar] [CrossRef]
  3. Freeze, A.R.; Cherry, J.A. Groundwater; Prentice-Hall, Inc.: Englewood Cliffs, NJ, USA, 1979; p. 604. [Google Scholar]
  4. Klovan, J.E. R- and Q-Mode Factor Analysis. In Concepts in Geostatistics; McCammon, R.B., Ed.; Springer: Berlin/Heidelberg, Germany, 1975; pp. 21–69. [Google Scholar]
  5. Lloyd, J.W.; Heathcote, J.A. Natural Inorganic Hydrochemistry in Relation to Groundwater; Clarendon Press: Oxford, UK, 1985; p. 15. [Google Scholar]
  6. Zhang, Z. Developing hydrogeology. Hydrogeol. Eng. Geol. 1979, 1, 12–14. [Google Scholar]
  7. Shen, Z.; Zhu, W.; Zhong, Z. Fundamentals of Hydrogeochemistry; Geological Press: Beijing, China, 1993; pp. 1–25. [Google Scholar]
  8. Gu, H.; Chi, B.; Wang, H.; Zhang, Y.; Wang, M. Relationship between surface water and groundwater in the Liujiang basin-Hydrochenmical constrains. Adv. Earth Sci. 2017, 32, 789–799. [Google Scholar]
  9. Fuchu, Z.; Bin, W.; Fan, G.; Mingliang, D.; Liantong, X. Hydrochemical characteristics of groundwater and evaluation of water quality in arid area of Northwest China: A case study in the plain area of Kuitun River Basin. Arab. J. Geosci. 2021, 14, 2099. [Google Scholar] [CrossRef]
  10. Adams, S.; Titus, R.; Pietersen, K. Hydrochemical characteristics of aquifers near Sutherland in the Western Karoo, South Africa. J. Hydrol. 2001, 241, 91–103. [Google Scholar] [CrossRef] [Green Version]
  11. Hitchon, B.; Billings, G.k.; Klovan, J.E. Geochemistry and origin of formation waters in the western Canada sedimentary basin—III. Factors controlling chemical composition. Geochim. Cosmochim. Acta 1971, 35, 567–598. [Google Scholar] [CrossRef]
  12. Gupta, L.P.; Subramanian, V. Geochemical factors controlling the chemical nature of water and sediments in the Gomti River, India. Environ. Geol. 1998, 36, 102–108. [Google Scholar] [CrossRef]
  13. Kumar, S.; Venkatesh A, S.; Singh, R. Geochemical signatures and isotopic systematics constraining dynamics of fluoride contamination in groundwater across Jamui district, Indo-Gangetic alluvial plains, India. Chemosphere 2018, 205, 493. [Google Scholar] [CrossRef]
  14. Mukherjee, I.; Singh U, K. Hydrogeochemical characterizations and quality evaluation of groundwater in the major river basins of a geologically and anthropogenically driven semi-arid tract of India. Sci. Total Environ. 2022, 805, 150323. [Google Scholar] [CrossRef]
  15. Kaur, L.; Rishi M, S.; Sharma, S. Hydrogeochemical characterization of groundwater in alluvial plains of River Yamuna in Northern India: An insight of controlling processes. J. King Saud Univ. Sci. 2019, 31, 1245–1253. [Google Scholar] [CrossRef]
  16. Mudzielwana, R.; Gitari, M.W.; Akinyemi, S.A. Hydrogeochemical characteristics of arsenic rich groundwater in Greater Giyani Municipality, Limpopo Province, South Africa. Groundw. Sustain. Dev. 2020, 10, 100336. [Google Scholar] [CrossRef]
  17. Nasher, N.M.R.; Ahmed, M.H. Groundwater geochemistry and hydrogeochemical processes in the Lower Ganges-Brahmaputra-Meghna River Basin areas, Bangladesh. J. Asian Earth Sci. X 2021, 6, 100062. [Google Scholar]
  18. Vinnarasi, F.; Srinivasamoorthy, K.; Saravanan, K. Hydrogeochemical characteristics and risk evaluation of potential toxic elements in groundwater from Shanmughanadhi, Tamilnadu, India. Environ. Res. 2022, 204, 112199. [Google Scholar] [CrossRef]
  19. Gomo, M. Conceptual hydrogeochemical characteristics of a calcite and dolomite acid mine drainage neutralised circumneutral groundwater system. Water Sci. 2018, 32, 355–361. [Google Scholar] [CrossRef] [Green Version]
  20. Rumuri, R. Dataset on hydrogeochemical characteristics of groundwater and surface water in Kattumannarkoil Taluk, India. Data Brief 2020, 32, 106058. [Google Scholar] [CrossRef]
  21. Aher, S.; Deshmukh, K.; Gawali, P. Hydrogeochemical characteristics and groundwater quality investigation along the basinal cross-section of Pravara River, Maharashtra, India. J. Asian Earth Sci. X 2022, 30, 100082. [Google Scholar] [CrossRef]
  22. Chu, W.; Chen, F.; Xiong, W.; Ge, Z.; Yuzhe, Z. Hydrogeochemical characterization and quality assessment of groundwater using self-organizing maps in the Hangjinqi gasfield area, Ordos Basin. NW China. Geosci. Front. 2021, 12, 781–790. [Google Scholar]
  23. Zhaoshi, W.; Xijun, L.; Kuanyi, L. Water quality assessment of rivers in Lake Chaohu Basin (China) using water quality index. Ecol. Indic. 2021, 20, 121. [Google Scholar]
  24. Peng, K.; Liu, W.; Wei, S.; Liu, C.; Chen, Y.; Dong, H.; Su, D.; Yuan, W.; Han, L. Study on the recharge source of karst groundwater in Jinan city based on hydrogeochenmical and isotopic characteristics. Carsologica Sin. 2020, 39, 650–657. [Google Scholar]
  25. Gao, Z.; Liu, J.; Li, Y.; Li, N.; Wang, M.; Wang, Z.; Wang, Z.; Liu, M. Hydrochenmical characteristics and hydrogeochenmical simulation of pore groundwater in Lhasa valley area. J. Shandong Univ. Sci. Technol. (Nat. Sci.) 2020, 39, 1–10. [Google Scholar]
  26. Jia, Z.; Zang, H.; Hobbs, P.; Zheng, X.; Xu, Y.; Wang, K. Application of inverse modeling in a study of the hydrogeochemical evolution of karst groundwater in the Jinci Spring region, northern China. Environ. Earth Sci. 2017, 76, 312. [Google Scholar] [CrossRef]
  27. Ren, X.; Li, P.; He, X.; Su, F.; Elumalai, V. Hydrogeochemical processes affecting groundwater chemistry in the central part of the Guanzhong Basin, China. Arch. Environ. Contam. Toxicol. 2021, 80, 74–91. [Google Scholar] [CrossRef]
  28. Wei, M.; Wu, J.; Li, W.; Zhang, Q.; Su, F.; Wang, Y. Groundwater Geochemistry and its Impacts on Groundwater Arsenic Enrichment, Variation, and Health Risks in Yongning County, Yinchuan Plain of Northwest China. Expo. Health. 2021, 1–20. [Google Scholar] [CrossRef]
  29. Subba Rao, N.; Ravindra, B.; Wu, J. Geochemical and health risk evaluation of fluoride rich groundwater in Sattenapalle Region, Guntur district, Andhra Pradesh, India. Human Ecol. Risk Assess. 2020, 26, 2316–2348. [Google Scholar] [CrossRef]
  30. Wu, J.; Li, P.; Wang, D.; Ren, X.; Wei, M. Statistical and multivariate statistical techniques to trace the sources and affecting factors of groundwater pollution in a rapidly growing city on the Chinese Loess Plateau. Human Ecol. Risk Assess. 2020, 26, 1603–1621. [Google Scholar] [CrossRef]
  31. Wu, J.; Zhang, Y.; Zhou, H. Groundwater chemistry and groundwater quality index incorporating health risk weighting in Dingbian County, Ordos basin of northwest China. Geochemistry 2020, 80, 125607. [Google Scholar] [CrossRef]
  32. Prasanna, M.V.; Chidambaram, S.; Shahul Hameed, A. Srinivasamoorthy, K. Study of evaluation of groundwater in Gadilam basin using hydrogeochemical and isotope data. Environ. Monit. Assess. 2010, 168, 63–90. [Google Scholar] [CrossRef] [Green Version]
  33. Parkhurst, D.L.; Appelo, C.A.J. Description of input and examples for phreeqc version 3—A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. In U.S. Geological Survey Techniques and Methods; Book 6, Chaper A43; U.S. Geological Survey: Denver, CO, USA, 2013; p. 497. [Google Scholar]
  34. Ma, T.; Zhang, C. Chemical characteristics of shallow groundwater in Guo River area. Ground Water 2021, 43, 22–24. [Google Scholar]
  35. Zheng, T.; Jiao, T.; Hu, B.; Gong, J.; Hou, X.; Wang, H. Hydrochemical characteristics and origin of groundwater in the central Guohe river Basin. Environ. Sci. 2021, 42, 766–775. [Google Scholar]
  36. He, X.-L.; Wu, Y.-H.; Zhou, J.; Bing, H.-J. Hydro-chemical Characteristics and Quality Assessment of Surface Water in Gongga Mountain Region. Environ. Sci. 2016, 37, 3798–3805. [Google Scholar]
  37. Wang, X.; Xu, Y.J.; Zhang, L. Watershed scale spatiotemporal nitrogen transport and source tracing using dual isotopes among surface water, sediments and groundwater in the Yiluo River Watershed, Middle of China. Sci. Total Environ. 2022, 833, 155180. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, X.; Zhang, G.; Xu, Y.J. Impacts of the 2013 extreme flood in Northeast China on regional groundwater depth and quality. Water 2015, 7, 4575–4592. [Google Scholar] [CrossRef] [Green Version]
  39. Wang, X.H.; Liu, C.L.; Hou, H.B.; Wang, X.Y. Identifying watershed-scale spatiotemporal groundwater and surface water mixing function in the Yiluo River, middle of China. Environ. Sci. Pollut. Res. 2021, 28, 11053–11065. [Google Scholar] [CrossRef]
  40. Gamvroula, D.; Alexakis, D.; Stamatis, G. Diagnosis of groundwater quality and assessment of contamination sources in the Megara basin (Attica, Greece). Arab. J. Geosci. 2013, 6, 2367–2381. [Google Scholar] [CrossRef]
  41. Li, H.; Xia, L.; Jiang, G. Analysis on geothermal fluid ion ratios in Tianjin area. Ground Water 2015, 37, 35–37. [Google Scholar]
  42. Tang, C.; Zheng, X.; Liang, Y. Hydrochemical Characteristics and Formation Causes of Ground Karst Water Systems in the Longzici Spring Catchment. Environ. Sci. 2020, 41, 2087–2095. [Google Scholar]
  43. Sun, P.; Yu, S.; Mo, F.; He, S.; Lu, J.; Yuan, Y. Hydrochemical characteristics and influencing factors in different geological background: A case study in Darongjiang and Lingqu Basin, Guangxi, China. Environ. Sci. 2016, 37, 123–131. [Google Scholar]
  44. Xu, D. Study on Chemical Evolution Characteristics and Control Factors of Groundwater in Suxian Mining Area. Master’s Degree, Hefei University of Technology, Hefei, China, 2017. [Google Scholar]
  45. LI, X. Study on Groundwater Circulation Characteristics in Guohe River Basin (Henan Section). Master’s Degree, North China University of Water Resources and Electric Power, Zhengzhou, China, 2020. [Google Scholar]
  46. Zhang, L. Study on Evaluation and Sustainable Development and Utilization of Geothermal Water Resources in Kaifeng City. Master’s Degree, Henan Polytechnic University, Jiaozuo, China, 2011. [Google Scholar]
  47. Jiang, X.; Li, J.; Guo, L.; Zhang, T.; Zhu, J. Chemical characteristics and formation mechanism of shallow groundwater in the northern Henan plain. Bull. Geol. Sci. Technol. 2021, 40, 290–300. [Google Scholar]
Figure 1. Distribution of hydrochemical sampling points in the Guohe River Basin.
Figure 1. Distribution of hydrochemical sampling points in the Guohe River Basin.
Water 14 01461 g001
Figure 2. Average content distribution of main ions and indicators in the three-layer aquifer.
Figure 2. Average content distribution of main ions and indicators in the three-layer aquifer.
Water 14 01461 g002
Figure 3. Piper tri-maps: shallow, intermediate, and deep groundwater in the Guohe River Basin (Henan Section).
Figure 3. Piper tri-maps: shallow, intermediate, and deep groundwater in the Guohe River Basin (Henan Section).
Water 14 01461 g003
Figure 4. Ion proportion coefficients of shallow, intermediate, and deep groundwater in the study area.
Figure 4. Ion proportion coefficients of shallow, intermediate, and deep groundwater in the study area.
Water 14 01461 g004
Figure 5. Gibbs diagram of shallow, intermediate and deep groundwater in the study area.
Figure 5. Gibbs diagram of shallow, intermediate and deep groundwater in the study area.
Water 14 01461 g005
Figure 6. Shallow, intermediate and deep simulation paths of Guohe River Basin (Henan Section).
Figure 6. Shallow, intermediate and deep simulation paths of Guohe River Basin (Henan Section).
Water 14 01461 g006
Figure 7. Hydrogeochemical evolution model of the three-layer groundwater aquifer along the chosen profile in the Guohe River Basin (unit: mg/L).
Figure 7. Hydrogeochemical evolution model of the three-layer groundwater aquifer along the chosen profile in the Guohe River Basin (unit: mg/L).
Water 14 01461 g007
Table 1. Concentration and content statistics of main hydrochemical indicators (unit: mg/L).
Table 1. Concentration and content statistics of main hydrochemical indicators (unit: mg/L).
TypeProjectK+Na+Ca2+Mg2+ClSO42−HCO3HardnessTDSTDS [35]
ShallowMean2.55147.6689.7773.84117.25126.50686.91528.12917823
S.D.4.79158.9449.3132.68122.55159.28191.52225.39507321
Min.0.3015.9015.3024.9012.303.82359.00141400289
Max.41.801420.00384.00205.00696.001287.001640.00178043721620
C.V. 11.881.080.550.441.051.260.280.430.550.39
IntermediateMean1.58203.0118.1816.2677.8676.28398.59112.356221322
S.D.0.3771.4521.0818.6665.9952.01114.69126.74178264
Min.0.6037.003.301.7021.5021.50251.0015366869
Max.2.40353.0087.8088.90320.00211.00807.0058510801710
C.V.10.240.351.161.150.850.680.291.130.290.20
DeepMean11.168.1960.8064.73372.6161.63560.001.56200792
S.D.15.1010.1342.3540.3960.5178.82136.330.3358139
Min.1.0055.103.001.4018.8022.10244.0013347561
Max.2.40330.0094.6057.90188.00193.00490.00475950970
C.V. 11.351.240.700.620.161.280.240.210.290.18
1 (The coefficient of variation is dimensionless).
Table 2. Mineral SI of representative points of shallow, intermediate and deep layers in Guohe River Basin (Henan Section).
Table 2. Mineral SI of representative points of shallow, intermediate and deep layers in Guohe River Basin (Henan Section).
PointSICalciteSIDolomiteSIGypsumSIHalitePointSICalciteSIDolomiteSIGypsumSIHalite
W-Q060.380.60−2.44−7.62W-S08−0.08−0.07−3.30−6.91
W-Q510.411.60−1.40−4.70W-S160.110.40−3.12−6.71
W-Q880.540.88−2.01−8.01W-S220.451.14−2.96−6.73
W-Z02−0.34−0.51−2.91−6.91W-S240.200.61−3.22−6.78
W-Z35−1.13−2.18−3.39−7.01W-S300.190.54−3.25−6.81
W-Z36−0.56−0.85−2.04−5.57
Table 3. Ion table of constant components at reverse hydrogeochemical simulation point (unit: mg/L).
Table 3. Ion table of constant components at reverse hydrogeochemical simulation point (unit: mg/L).
PointK+Na+Ca2+Mg2+ClSO42−HCO3pHTemp.TDS
W-Q062.3032.2092.3033.4027.9013.10472.007.2719.0453
W-Q512.401420.0034.00114.006961287.001640.007.5718.74372
W-Z021.80130.0012.209.6034.0024.30294.007.5719.5403
W-Z351.10164.004.403.1021.5021.50359.007.1816.6420
W-Z361.75336.0022.0021.70320.00154.00282.007.2120.11080
W-S081.60176.004.402.5025.5028.90344.008.1823.5438
W-S161.30190.005.003.0038.6041.20335.008.2828.6477
W-S301.20230.003.702.2025.4045.50461.008.4326.4552
Table 4. Reverse hydrogeochemical simulation results of Guohe River Basin (Henan Section) (unit: mmol/L).
Table 4. Reverse hydrogeochemical simulation results of Guohe River Basin (Henan Section) (unit: mmol/L).
Mineral PhaseShallow Simulation PathIntermediate Simulation PathDeep Simulation Path
W-Q06→W-Q51W-Z02→W-Z35W-Z35→W-Z36W-S08→W-S16W-S16→W-S30
CaCO312.5501.098−4.343−1.398−1.264
CaMg (CO3)23.408−0.2680.7670.021−0.033
CaSO4·2H2O3.246−0.1093.6931.4762.284
NaCl18.940−0.3538.4300.370−0.452
CO2−0.0891.4951.5521.2103.398
O220.1700.159−4.624−2.696−4.478
H2S10.0800.079−2.312−1.348−2.239
NaX41.0801.833−0.6480.1672.038
CaX2−20.540−0.9160.324−0.083−1.019
(Positive values indicate dissolution and negative values denote precipitation).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, F.; Zhou, D.; Li, Z.; Li, X. Hydrochemical Characteristics and Hydrogeochemical Simulation Research of Groundwater in the Guohe River Basin (Henan Section). Water 2022, 14, 1461. https://doi.org/10.3390/w14091461

AMA Style

Yu F, Zhou D, Li Z, Li X. Hydrochemical Characteristics and Hydrogeochemical Simulation Research of Groundwater in the Guohe River Basin (Henan Section). Water. 2022; 14(9):1461. https://doi.org/10.3390/w14091461

Chicago/Turabian Style

Yu, Furong, Dongxu Zhou, Zhiping Li, and Xiao Li. 2022. "Hydrochemical Characteristics and Hydrogeochemical Simulation Research of Groundwater in the Guohe River Basin (Henan Section)" Water 14, no. 9: 1461. https://doi.org/10.3390/w14091461

APA Style

Yu, F., Zhou, D., Li, Z., & Li, X. (2022). Hydrochemical Characteristics and Hydrogeochemical Simulation Research of Groundwater in the Guohe River Basin (Henan Section). Water, 14(9), 1461. https://doi.org/10.3390/w14091461

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop