Next Article in Journal
Numerical Investigation of the Sediment Load Exchange between a Coastal Mud Bank and Its Neighbouring Estuary
Previous Article in Journal
Reductive Sequestration of Chromate with Pyrite-Loaded nZVI@biochar Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogeochemical Characteristics and Sulfate Source of Groundwater in Sangu Spring Basin, China

1
Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang 050000, China
2
Key Laboratory of Groundwater Science and Engineering, Ministry of Natural Resources, Shijiazhuang 050000, China
*
Authors to whom correspondence should be addressed.
Submission received: 7 July 2024 / Revised: 5 October 2024 / Accepted: 8 October 2024 / Published: 11 October 2024
(This article belongs to the Section Hydrogeology)

Abstract

:
The Sangu Spring Basin is located in an important economic area, and groundwater is the main source of water for local life and industry. Understanding the sources of chemical components in groundwater is important for the development and utilization of groundwater. In this paper, we analyzed the origin of the chemical components of groundwater and their evolution in the Sangu Spring Basin using statistical analysis, Piper diagrams, Gibbs diagrams, ion ratios, and combined hydrochemistry–isotope analyses. The results show that the groundwater in the Sangu Spring Basin is mainly derived from atmospheric precipitation, that the groundwater in stagnant and confined environment zones was formed under colder climatic conditions, and that the surface water (SW) has a close hydraulic relation with the groundwater. Water–rock interaction is the main factor controlling the composition of groundwater. The compositions of groundwater are mainly derived from carbonate weathering, silicate weathering, and dissolution of gypsum. Na+ and K+ in groundwater mainly come from the dissolution of albite and potassium feldspar, rather than rock salt. Ion exchange occurs in karst groundwater (KGW) and fissure groundwater (FGW), and ion exchange is dominated by the exchange of Mg2+ and Ca2+ in the groundwater with Na+ and K+ in the rock or soil. Sulfate in groundwater is derived from dissolution of gypsum, infiltration of atmospheric precipitation, and leakage of SW. Groundwaters with the highest sulfate content are located in the vicinity of SW, as a result of receiving recharge from SW seepage. Groundwaters with higher sulfate contents are located in the stagnant and deeply buried zones, where sulfate is mainly derived from the dissolution of gypsum. SW seepage recharges groundwater, resulting in increased levels of Cl, NO3 and SO42− in groundwater. These insights can provide assistance in the protection and effective management of groundwater.

1. Introduction

The Sangu Spring Basin is one of the important coal-producing areas in China and has an important economic status. Coal mining and chemical production are the main economic activities in the study area. The Sangu Spring Basin is located in arid and semi-arid areas, and, coupled with the lack of SW resources and poor SW quality, this means that groundwater has become the main source of local industrial water and domestic water [1,2]. Groundwater is an essential element of the hydrological cycle and plays a key role in the base flow of rivers and the maintenance of terrestrial aquatic ecosystems [3,4]. There is a close hydrological connection between groundwater and SW in the Sangu Spring Basin, so it is of great significance to study the evolution of groundwater hydrochemistry for the management of groundwater.
Hydrochemical components and isotopes of groundwater are often used to analyze problems in hydrogeology [5,6,7,8]. Gibbs diagrams, Piper diagrams, deuterium–oxygen relationship diagrams, ion ratios, and combined hydrochemistry–isotope analyses are commonly used to analyze the origin and evolution of the chemical composition of groundwater [9,10]. Stable isotopes are the tools for studying the origin of groundwater, recharge elevations, recharge mechanisms, and hydraulic relations, and for calculating the mixing ratio of different water sources [11,12,13]. Combined use of δ34S-SO42−, δ18O-SO42−, δ18O-H2O, δD-H2O, and hydrochemistry–isotopes can be used to study the sources of sulfate in groundwater or SW and calculate the mixing ratios between different water sources [14,15,16]. Therefore, the method of combining hydrochemistry with isotopes can determine the geochemical factors and mechanisms that control the chemical composition of groundwater.
Previous studies on groundwater in the Sangu Spring Basin have mainly focused on the determination of spring boundaries [17,18,19], evaluation of groundwater quality [20,21], evaluation of water resources [17,22], decline of the groundwater level [23,24], protection of aquifers [2,25,26], and the impact of coal mining on the groundwater environment [27,28]. Only Zhang et al. [29] have studied the hydrochemical characteristics and evolution mechanism of KGW in the Sangu Spring Basin. At present, there is no systematic and comprehensive investigation on the hydrogeochemical behavior of groundwater in the Sangu Spring Basin. This paper studies the sources of hydrochemical components and their evolution of pore groundwater (PGW), FGW, and KGW in the Sangu Spring Basin by using statistical analysis, Gibbs diagrams, Piper diagrams, ion ratios, and joint hydrochemical–isotope analysis. The research results provide a scientific basis for the protection, pollution prevention, exploitation, and management of groundwater in the spring area.

2. Study Area

2.1. General Setting

The Sangu Spring Basin is located in the southeast of Shanxi province (Figure 1), China, with an area of 3214 km2. The Sangu Spring is exposed in the valley of the Dan River, which has now been flooded by the Qingtianhe Reservoir. The Sangu Spring Basin is surrounded by high mountains, with undulating intermountain basins and hills in the middle. The main rivers in Sangu Spring Basin are the Dan River, Baishui River, and Dongdan River, which belong to the Yellow River system. There are two medium-sized reservoirs, Renzhuang Reservoir and Qingtianhe Reservoir, as well as numerous small reservoirs in the Sangu Spring Basin. The Sangu Spring Basin belongs to the warm, temperate, semi-humid continental monsoon climate with four distinct seasons. The average annual temperature is 13.7 °C. The average annual precipitation is 593.66 mm and the average annual evaporation is 1633.95 mm [1,30].

2.2. Geological and Hydrogeological Setting

The outcrop strata in the Sangu Spring Basin include Cambrian, Ordovician, Carboniferous, Permian, and Quaternary strata. The Cambrian strata are distributed in a small area and only appear in deep gullies in the southeast of the spring region. The Ordovician and Carboniferous strata are widely distributed in the area, with large areas exposed in the western, central, eastern, and southwestern parts of the study area. The Quaternary strata are mainly exposed in the hills and basins between Gaoping City and Jincheng City. The Taiyuan Formation of the Carboniferous and the Shanxi Formation of the Permian systems are coal-bearing strata, which are widely distributed in the central and northern parts of the spring region.
The main aquifers (generalized) from bottom to top in the Sangu Spring Basin are Cambrian limestone aquifers, Ordovician limestone aquifers, Carboniferous thin limestone and sandstone aquifers, Permian sandstone aquifers, and Quaternary loose rock aquifers (Figure 2). KGW is mainly contained in the Cambrian and Ordovician strata, with the Zhangxia Formation of the Middle Cambrian and the Majiagou Formation of the Middle Ordovician as the main aquifers. The upper Cambrian and lower Ordovician Sanshanzi formations are a water-resisting layer in the region. FGW is stored in the Carboniferous and Permian strata, with the Taiyuan Formation limestone of the Carboniferous, the sandstone of the Shanxi Formation, and the Sandstone of the Shihezi Formation of the Permian as the main water-bearing layers. The Carboniferous Benxi Formation bauxite mudstone and Permian mudstone are water-resisting layers. PGW is stored in the Quaternary strata, and the aquifers are sand and gravel layers.
KGW is the most widely distributed and best water-rich groundwater in the Sangu Spring Basin, and it is also the most important source of water supply. FGW is the groundwater that is most affected by coal mining, and most of the springs’ flow rates become smaller or dry, and the water supply function becomes weaker gradually. PGW is mainly located in the mountain basins and river terrace in the middle of the spring area, with the smallest distribution area, thin aquifer thickness, and uneven water-richness.
Characteristics of the supplement, runoff, and draining of KGW are affected by various factors such as geological structure, stratigraphic lithology, distribution of modern water networks, and human activities [1,17,30].
The recharge sources of KGW include infiltration of atmospheric precipitation, leakage of surface water, and leaking recharge of the overlying aquifer. It is recharged by direct infiltration of atmospheric precipitation in the exposed areas of limestone, by leakage of surface water in the valleys and reservoirs of the limestone area, and by leaking recharge from the overlying aquifers in the overlying and buried areas.
There are three strong-flow zones in the KGW of the Sangu Spring Basin (Figure 1): one along the middle and lower valleys of the Dan River, one along the Gaoping–Bagong–Beishidian–Jincheng City line, and another along the valley of the Dongdan River [1,17,18].
The drainage methods of KGW include springs and artificially extracted groundwater. Drainage of spring water is mainly concentrated in the Guobi Spring (flow rate 0.4–0.66 m3/s) and the Sangu Spring (flow rate 3.5–7.0 m3/s). KGW is the main source of water for industrial production and domestic use in the Sangu Spring area, and groundwater exploitation is mainly distributed in the north–central part of the spring area, with an extraction volume of about 1.1568 × 108 m3/a.

3. Materials and Methods

3.1. Sample Collection and Treatment

From 2014 to 2015, 234 groundwater samples, 27 SW samples, and 12 atmospheric precipitation samples were collected, with the sampling point locations shown in Figure 1. Prior to the collection of well-water samples, the wells were first pumped until the electrical conductivity (EC) of the well water was stabilized. The samples were filtered through a 0.45 μm membrane filter and then placed into 1500 mL and 550 mL polyethylene bottles for analysis of major ionic constituents and trace elements. Nitric acid was added to the 550 mL samples until the pH was below 2, and the 1500 L samples did not need to be treated. The samples were stored below 4 °C until they were analyzed. Samples for stable isotope analysis (δ18O-H2O and δD-H2O) were collected into 50 mL glass bottles with a sealed mouth.
From 234 groundwater samples, 152 samples were selected for δ18O-H2O and δD-H2O analysis. A total of 12 atmospheric precipitation samples were all analyzed for δ18O-H2O and δD-H2O. Sixteen of the 27 SW samples were selected for δ18O-H2O and δD-H2O analysis.
A total of 33 samples were selected from 234 groundwater samples for δ18O-SO42− and δ34S-SO42− analysis. Two samples were selected from 27 SW samples for δ18O-SO42− and δ34S-SO42− analysis. The δ18O-SO42− and δ34S-SO42− analyses were performed on 5 samples from 12 atmospheric precipitation samples. Samples of δ18O-SO42− and δ34S-SO42− isotopes were collected in 5 L polyethylene containers and acidified with HCl to pH 2–3, and then saturated BaCl2 solution was added to precipitate BaSO4.

3.2. Analytical Methods

The hydrochemical components were determined by the Groundwater-Mineral Water and Environmental Monitoring Center, Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences. Major cations and trace elements were analyzed using a PerkinElmer Inductively Coupled Plasma Emission Spectrometer Optima Model 8300 (accuracy ±1%). The reliability of the water chemistry data was assessed by examining the ionic balance, with the ion charge imbalances within ±5%.
The concentration of HCO3 was determined by phenolphthalein titration. Except for HCO3, major anions were analyzed using a Thermo Scientific Dionex ICS-4000 (Dionex, Sunnyvale, CA, USA; accuracy: ±1%). The reliability of the water chemistry data was assessed by checking that the ionic balances were within ±5%.
The δ18O-H2O and δD-H2O isotopes were determined by the Picarro L 2130-i analyzer (Picarro, Santa Clara, CA, USA) of the Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences. The data are expressed in δ-symbols (δ) and have a precision of ±0.1‰ for δ18O-H2O and ±1‰ for δD-H2O relative to Vienna Standard Mean Ocean Water (V-SMOW).
δ18O-SO42− and δ34S-SO42− were analyzed in the Analytical Laboratory Beijing Research Institute of Uranium Geology (ALBRIUG). δ18O-SO42− was determined by coupling CO production with a Finnigan MAT 253 on a Thermal Conversion Elemental Analyzer (TCEA, Thermo Fisher Scientific, Bremen, Germany) at a temperature of 1450 °C. δ34S-SO42− was determined by coupling SO2 production with a Finnigan MAT 253 Isotope Ratio Mass Spectrometer (IRMS) on a TCEA (Thermo Fisher Scientific, Bremen, Germany). The working standards for δ18O-SO42− were GBW 04409, GBW 04410, and NBS 28. The working standards for δ34S-SO42− measurements were GBW 04414 and GBW 04415. δ34S-SO42− values corresponded to the Vienna Canyon Diablo Troilite (VCDT) standard, and δ18O-SO42− values correspond to the V-SMOW standard. The accuracy of the δ18O-SO42− and δ34S-SO42− analyses was better than ±0.5‰ and ±0.2‰, respectively, based on repeated analyses of internal standards.
Maximum, minimum, mean, and standard deviation of water chemistry components and deuterium–oxygen isotopes were calculated using IBM SPSS Statistics 26 statistical analysis software.

4. Results

4.1. Hydrochemical Characterization of Groundwater

The basic characteristics of the groundwater are shown in Table 1 and Figure 3.
A total of 137 samples of KGW were collected in the Sangu Spring Basin. The pH of the KGW ranged from 7.20 to 8.57, and the TDS ranged from 218.7 to 2411.0 mg/L. The hydrochemical types were dominated by HCO3-Ca, HCO3-Ca·Mg, HCO3·SO4-Ca·Mg, SO4·HCO3-Ca·Mg, SO4·HCO3-Na·Ca, and HCO3·SO4-CaHCO3-Ca. HCO3 and SO42− were the dominant anions, and Ca2+, Mg2+, and Na+ were the dominant cations.
A total of 74 samples of FGW were collected in the Sangu Spring Basin. The pH of the FGW ranged from 6.71 to 8.51, and the TDS ranged from 245.80 to 2424.50 mg/L. The hydrochemical types were dominated by HCO3-Ca, HCO3·SO4-Ca, SO4·HCO3-Ca, and SO4·HCO3-Ca·Na. HCO3 and SO42− were the dominant anions, and Ca2+ and Na+ were the dominant cations.
A total of 23 PGW samples were collected from the Sangu Spring Basin. The pH of PGW ranged from 7.35 to 8.28, and the TDS was 330.5–1226.0 mg/L. The hydrochemical types were dominated by HCO3·SO4-Ca, SO4·HCO3-Ca (Na), SO4·HCO3-Na (Ca), HCO3-SO4-Ca, SO4-HCO3-Ca (Na), and SO4-HCO3-Na (Ca). HCO3 and SO42− were the dominant anions, and Ca2+ and Na+ were the dominant cations.

4.2. Isotope Characterization of Groundwater

The δD-H2O and δ18O-H2O values of the KGW ranged from −83.0‰ to −42.0‰ and from −11.2‰ to −4.5‰, respectively. The δD-H2O and δ18O-H2O values of the FGW ranged from −83.0‰ to −57.0‰ and from −11.5‰ to −8.1‰, respectively. The δD-H2O and δ18O-H2O values of the PGW ranged from −69.0‰ to −48.0‰ and from −9.6‰ to −5.8‰, respectively. The δD-H2O and δ18O-H2O values of the SW ranged from −70.0‰ to −42.0‰ and from −9.1‰ to −3.9‰, respectively (Table 2). The δD-H2O and δ18O-H2O contents of SW and PGW were higher, while the δD-H2O and δ18O-H2O contents of some KGW and FGW were lower. The δD-H2O and δ18O-H2O contents of KGW and FGW had a wide range and were closely related to their hydrogeological conditions, and groundwaters that directly receive infiltration recharge from atmospheric precipitation or seepage recharge from SW had higher δD-H2O and δ18O-H2O contents.

5. Discussion

5.1. Groundwater Origins Indicated by δD-H2O and δ18O-H2O

The isotopes of δD-H2O and δ18O-H2O are commonly used to analyze the origin of groundwater [13]. The relationship between δD-H2O and δ18O-H2O in the samples is shown in Figure 4. The local meteoric water line (LMWL) is fitted to δD = 7.78δ18O-H2O + 10.33 [1], and its slope is slightly lower than the global meteoric water line (GMWL:δD = 8δ18O-H2O + 10) [31,32]. All groundwater samples were distributed along the atmospheric precipitation line (Figure 4), showing the origin of groundwater atmospheric precipitation.
SW is mainly distributed in Group 4, and some SW showed a drift of δ18O-H2O (Figure 4). SW close to the LMWL is mainly formed by the recent atmospheric precipitation recharge or received groundwater recharge, such as the middle and lower reaches of the Baishui River and the lower reaches of the Dan River. SW with δ18O-H2O drift is mainly affected by evaporation. For example, J53 is a sample of reservoir water from Shan’erdong Reservoir, and the water level of the reservoir was at the stagnant level when the sample was taken, with less recharge and strong evaporation, so there was a strong δ18O-H2O drift.
There are one PGW sample and one KGW sample in Group 4 (Figure 4). P11 is shallow PGW, located in the first terrace of Dan River. Its recharge sources are atmospheric precipitation and river-water infiltration (the river water comes from the reservoir). It is assumed that the higher deuterium–oxygen isotope content is due to the combined effect of recharge from the river water and evaporation. P10-f is KGW, and it is assumed that its higher deuterium–oxygen isotope content is due to the recharge of SW from Renzhuang Reservoir. The staggered distribution of SW and groundwater indicates a close hydraulic relationship between SW and groundwater.
The mixed distribution of SW, PGW, FGW, and KGW in Group 3 (Figure 4) indicates that there is a close hydraulic relationship between them. For example, the higher deuterium and oxygen isotope content of Je35, D112, etc. in KGW is due to the fact that they receive recharge from overburdened PGW, which in turn receives seepage recharge from the Dan River, and the river water is affected by evaporation. For instance, the shallow KGW (XB90, XA165, E184) directly receives infiltration recharge from atmospheric precipitation, has a short cycle path, a shallow water level, and is affected by a certain degree of evapotranspiration, resulting in a high deuterium–oxygen isotope content.
Group 2 is dominated by KGW (Figure 4), located in a zone of weak runoff, with slow water turnover and tritium levels of less than 7 TU. The groundwater may contain groundwater recharged during colder climates.
The groundwater in Group 1, which includes two KGW and two FGW samples, has the lowest deuterium–oxygen isotope content (Figure 4). The two KGW samples in Group 1 are located on the boundary of the Sangu Spring Basin, which is in the stagnant-flow area of KGW, the water alternation is extremely slow, and the groundwaters may have formed in a period of colder climate [6,7,32]. The two FGW samples in Group 1 are water discharged during the extraction of coalbed methane (ZC38, ZC112), and the groundwaters are in a confined environment, speculated to have been formed during a period of colder climate [33].

5.2. Causes of Groundwaters

The Gibbs diagram can be used to study the major diagenetic mechanisms of groundwater evolution. At present, this method is widely used to identify the sources of groundwater chemical composition [34,35]. The Gibbs diagram is generally divided into rock weathering dominant zones, evaporation dominant zones, and precipitation dominant zones, and the water–rock action dominant zones can be further divided into carbonate dissolution dominant zones and silicate dissolution dominant zones [36]. Groundwater samples from the Sangu Spring Basin were mainly distributed in areas where rock weathering is dominant (Figure 5), indicating that groundwater is mainly controlled by water–rock interaction. The Na+/(Na+ + Ca2+) of groundwater varies in a larger range, while TDS varies in a smaller range, indicating that cation exchange may influence the variation of Na+ and Ca2+ content in groundwater [37]. Rock weathering is the main factor affecting the hydrochemical components of groundwater and does not prove that groundwater formation mechanisms are not affected by human activities. The influence of human activities includes NO3, SO42−, and Cl directly or indirectly caused by acid rain, sewage, and chemical fertilizer running into groundwater.
The relationship diagrams for Mg2+/Na+ versus Ca2+/Na+ and HCO3/Na+ versus Ca2+/Na+ were plotted (Figure 6), through which the sources of dissolved substances in the groundwater can be further determined [38,39,40,41]. KGW, FGW, and PGW are located on the weathering line of silicates–carbonates and near the end element (Figure 6), indicating that carbonate weathering and silicate weathering are the main controlling factors of groundwater components in the Sangu Spring Basin. FGW is located on the weathering line of silicate–carbonate, which is consistent with the nature of the surrounding rocks of the aquifer. The aquifers of FGW are sandstone and limestone, especially in the Carboniferous strata, where limestone is the main aquifer. The PGW is also located in the weathering line of silicate–carbonate, which is consistent with the hydrogeological conditions. PGW is mainly located in intermountain basins and river terraces, and groundwater from the surrounding mountains is the main source of recharge.

5.3. Major Ion Ratios of Groundwater and Their Implications for Mineral Dissolution

Ion ratios are a common method for analyzing the source of groundwater components and have been widely used [7,35,42,43]. The chemical reaction equations for carbonate and silicate weathering are shown in Equation (1) through Equation (3). As a result of weathering, the anions increase HCO3, SO42−, and NO3 and the cations increase Ca2+, Mg2+, Na+, and K+ in the chemical fractions of groundwater.
Carbonate minerals weathering:
4 ( C a 1 x M g x ) C O 3 ( s ) + H 2 C O 3 + H 2 S O 4 + H N O 3 4 ( 1 x ) C a 2 + ( aq ) + 4 xM g 2 + ( a q ) + 5 H C O 3 ( aq ) + N O 3 ( aq ) + S O 4 2 ( a q )
Silicate minerals weathering:
C a A l 2 S i 2 O 8 ( s ) + C O 2 ( g ) + 2 H 2 O ( l ) + H 2 S O 4 + H N O 3 2 C a 2 + ( a q ) + H C O 3 ( a q ) + N O 3 ( a q ) + S O 4 2 ( a q ) + 4 S i O 2 + A l ( O H ) 3 + 3 A l O O H
4 N a ( K ) A l S i 3 O 8 ( s ) + C O 2 ( g ) + 2 H 2 O ( l ) + H 2 S O 4 + H N O 3 4 N a + ( K + ) ( a q ) + H C O 3 ( a q ) + N O 3 ( a q ) + S O 4 2 ( a q ) + 12 S i O 2 + A l ( O H ) 3 + 2 A l O O H
The Ordovician strata in the Sangu Spring Basin contain gypsum intercalations, and the dissolution of gypsum can lead to the entry of Ca2+ and SO42− into the groundwater (Equation (4)). If Ca2+ and SO42− in KGW are mainly derived from gypsum dissolution, then the γCa2+/γSO42− (γCa2+ represents the milligram equivalent of calcium ions in groundwater, and γSO42− represents the milligram equivalent of sulfate ions in groundwater, in meq/L.) ratio is 1. KGW samples are not distributed along the line of gypsum dissolution (Figure 7a), which indicates that gypsum dissolution is not one of the main sources of Ca2+ and SO42− in KGW.
C a 2 S O 4 ( s ) C a 2 + ( aq ) + S O 4 2 ( a q )
If only carbonate weathering occurs, then γ(Ca2+ + Mg2+): γ(HCO3 + CO32− + SO42− + NO3) = 1:1, and the groundwater should be distributed along the carbonate weathering line (Figure 7c). If only silicate weathering occurs, then γ(Ca2+ + Na+ + K+): γ(HCO3 + CO32− + SO42− + NO3) = 1:1. If carbonate weathering, silicate weathering, and sulfate dissolution occur at the same time, then γ(Ca2+ + Mg2+ + Na+ + K+ − Cl): γ(HCO3 + CO32− + SO42− + NO3) = 1:1. γ(Ca2+ + Mg2+ + Na+ + K+ − Cl) represents the cationic component of carbonate weathering, silicate weathering, and sulfate dissolution into the groundwater. In the plot of γ(Ca2+ + Mg2+ + Na+ + K+ − Cl) versus γ(HCO3 + CO32− + SO42− + NO3) (Figure 7b), the PGW, FGW, and KGW are all distributed along the line of 1:1, which indicates that the chemical fractions of the groundwater in the Sangu Spring Basin watershed mainly originated from carbonate weathering, silicate weatherig, and dissolution of sulfate.
Weathering of salt rocks causes increases in Na+, K+, and Cl content in groundwater (Equations (5) and (6)). The groundwater is not distributed along the dissolution line of the salt rocks (Figure 7d), which indicates that the dissolution of the salt rocks is not the main source of Na+ and K+ in the groundwater in the Sangu Spring Basin. Dissolution of carbonic, sulfuric, and nitric acids with sodium feldspar (NaAlSi3O8) and potassium feldspar (KAlSi3O8) (Equation (3)) is the source of Na+ and K+ content in groundwater.
N a C l ( s ) N a + ( aq ) + C l ( a q )
                K C l ( s ) K + ( aq ) + C l ( a q )
Precipitation contains SO2 and NOx, and their participation promotes the weathering of carbonate and silicate minerals (Equations (1) and (2)). The correlation coefficients between (Ca2+ + Mg2+) and (HCO3 + SO42− + CO32− + NO3) in KGW, FGW, and PGW are 0.85, 0.78, and 0.76, respectively, showing a significant positive correlation of level 0.01, indicating that the dissolution of minerals in groundwater was mainly controlled by carbonic acid, sulfuric acid, and nitric acid.

5.4. Ion Exchange

The exchange of Ca2+ and Mg2+ in groundwater with Na+ and K+ in soil or rock is a common cause of changes in ionic composition in groundwater. The equations for cation exchange occurring in groundwater are shown in Equations (7) and (8). In the Gibbs plots, the Na+/(Na+ + Ca2+) ratios of some groundwater samples are greater than 0.5, indicating that cation exchange has some influence on the change in cation content in groundwater [37,44]. Scherer (1977) [45] proposed a mathematical relationship that can better represent the ion exchange that occurs between groundwater and aquifer rocks, called the chlor–alkali index (CAI). The chlor–alkali index method provides a quantitative analysis of ion exchange and reverse ion exchange processes (Equations (7) and (8)) and is widely used [31,35,46,47].
C a 2 + ( M g 2 + ) ( aq ) + 2 N a ( K ) X ( s ) 2 N a + ( K + ) ( aq ) + C a ( M g ) X 2 ( s )
C a ( M g ) X 2 ( s ) + 2 N a + ( K + ) ( aq ) C a 2 + ( M g 2 + ) ( aq ) + 2 N a ( K ) X ( s )
C A I 1 = ( C l ( N a + + K + ) ) / C l
C A I 2 = ( C l ( N a + + K + ) ) / ( H C O 3 + C O 3 2 + S O 4 2 + N O 3 )
If the chemical reaction in Equation (7) occurs, CAI1 < 0, and CAI2 < 0, ion exchange occurs, and Mg2+ and Ca2+ in groundwater are exchanged with Na+ and K+ in rock or soil. If the chemical reaction in Equation (8) occurs, CAI1 > 0, and CAI2 > 0, the Na+ and K+ in groundwater are exchanged with Mg2+ and Ca2+ in rock or soil. The ranges of CAI1 and CAI2 in KGW were −30.12–5.86 and −0.58–0.08, respectively (Figure 8a). A total of 69.34% of KGW experienced ion exchange, which means that Mg2+and Ca2+ in KGW were exchanged with Na+ and K+ in rock or soil. A total of 8.03% of KGW experienced reverse ion exchange, which means that Na+ and K+ in KGW were exchanged with Mg2+and Ca2+ in rock or soil. The ranges of CAI1 and CAI2 in FGW were −40.41–0.46 and −0.66–0.10, respectively. A total of 85.14% of the FGW experienced ion exchange, which means that Na+ and K+ were exchanged with Mg2+ and Ca2+ in rocks or soils. A total of 12.16% of the FGW experienced reverse ion exchange, which means that Na+ and K+ were exchanged with Mg2+ and Ca2+ in rocks or soils. The ranges of CAI1 and CAI2 in PGW were −11.74–0.52 and −0.48–0.10, respectively. A total of 69.57% of the PGW had ion exchange, which means that Mg2+ and Ca2+ in the PGW were exchanged with Na+ and K+ in the rock or soil. A total of 26.09% of the PGW had reverse ion exchange, which means that Na+ and K+ were exchanged with Mg2+and Ca2+ in the rock or soil.
In recent years, the relationship plot of (K+ + Na+ − Cl) versus ((Ca2+ + Mg2+) − (HCO3 + SO42− + CO32− + NO3) has often been used to determine whether cation exchange occurs in groundwater [31,48,49]. KGW, FGW, and PGW are basically distributed along the −1:1 line (Figure 8b), indicating that ion exchange occurs in groundwater. Ion exchange is dominated by the exchange of Mg2+ and Ca2+ in groundwater with Na+ and K+ in rock or soil. In a small portion of groundwater, Na+ and K+ are exchanged with Mg2+ and Ca2+ in rocks or soils. Although the PGW samples are distributed along the −1:1 line, the CAI1 and CAI2 of most of the samples are close to zero and the ion exchange is weak.

5.5. Sources of Sulfate

Sources of sulfate in groundwater include natural and anthropogenic sources [50,51,52]. Natural sources include atmospheric deposition, dissolution of sulfate minerals (e.g., gypsum), sulfide oxidation, seawater intrusion, and volcanic eruption [53,54]. Anthropogenic sources include domestic wastewater, industrial wastewater, mining activities, and the application of chemical pesticides, etc. [14,55,56,57,58,59,60]. The sulfur and oxygen isotopes of sulfates from various sources in groundwater have significant differences. Combined with the hydrogeological conditions, the sources of sulfate in groundwater can be identified by the comprehensive use of water chemistry and multi-isotope analysis [61,62,63,64].
In Figure 9a, no water samples fall within the experimental zone of sulfate generation from sulfide oxidation [65], indicating that sulfate in groundwater, SW, and atmospheric precipitation does not originate from sulfide oxidation. This is also illustrated by the absence of samples within the sulfate oxidation zone in Figure 9d.
KGW in Zone I in Figure 9b has higher SO42−content and the heaviest δ34S-SO42− [66]. These KGWs are located in the deeply buried groundwater zone in the northern part of the Sangu Spring Basin and in the stagnant-flow zone near the spring boundary. Due to the large circulation depth of the KGW, poor groundwater runoff conditions, and the long water–rock interaction time, the sulfate mainly comes from the dissolution of gypsum. The KGW located in Zone I in Figure 9b falls in the gypsum dissolution zone in Figure 9d, which also indicates that the sulfate in the KGW comes from the dissolution of gypsum. According to geological data, gypsum deposits are widespread in the Middle Ordovician strata in the Sangu Spring Basin. Gypsum has been dissolved in areas with good water circulation conditions, but is still present in deeply buried areas. Gypsum was found in the borehole where the ZK05 water sample was collected, in the Ordovician strata in the form of fissure filling, and the δ34S-SO42− content of the gypsum was 27.5. Figure 9c also proves that the KGW in the deeply buried zone and stagnant-flow zone has the heaviest δ34S-SO42− and the oldest age.
Groundwater in Zone II in Figure 9b has lower SO42−content and heavier δ34S-SO42−. The two FGW samples are the groundwater discharged from the extraction of coalbed methane. The groundwater is in a reducing environment, and the reducing effect of sulfur bacteria leads to an increase in δ34S-SO42− and a simultaneous decrease in SO42− content in the groundwater.
The groundwater with the lowest SO42− content and the lightest δ34S-SO42− is located in zone III of Figure 9b, and the atmospheric precipitation is also located in this zone. The sulfate in KGW and FGW in this area mainly comes from atmospheric precipitation.
The groundwater with the highest SO42−content and lighter δ34S-SO42− is located in zone IV of Figure 9b, and sample E141 from the Baishui River is also located in this zone. Due to the low SO42− content in precipitation, atmospheric precipitation is unlikely to be the main source of SO42− in these groundwater samples [67]. KGW in this area is affected by SW leakage, resulting in higher SO42− content and lighter δ34S-SO42− in groundwater. Groundwater in Zone IV has a higher deuterium–oxygen isotope content and is in Group 3 with E141 in the δD-H2O versus δ18O-H2O relationship diagram.
Groundwater in the other areas in Figure 9b has higher SO42− content and lighter δ34S- SO42−; it is speculated that sulfate in groundwater comes from atmospheric precipitation, SW bodies, and gypsum dissolution. These samples in Figure 9d are also located between the atmospheric precipitation, SW, and gypsum end members, and are subject to the combined influence of atmospheric precipitation infiltration, SW body seepage, and gypsum dissolution.
In conclusion, the sulfate in the groundwater of the Sangu Spring Basin is mainly influenced by the dissolution of gypsum, the infiltration of atmospheric precipitation, and the seepage of SW bodies. Sulfate in the deeply buried and stagnant-flow areas of KGW is mainly derived from the dissolution of gypsum, and groundwater recharged by SW seepage tends to have a high sulfate content. FGW in reducing environments may have experienced the reducing effect of sulfur bacteria.
The hydrochemical components of groundwater in the Sangu Spring Basin are mainly affected by water–rock interaction and SW leakage. Further research on the water–rock interaction of groundwater and the hydraulic connection between surface water and groundwater is recommended. The protection and management of surface water quality should be strengthened, and already polluted surface water should be treated and monitored to avoid groundwater pollution caused by surface water seepage.

6. Conclusions

Using a combination of Piper diagrams, Gibbs diagrams, ion ratios, and hydrochemistry–isotope analyses is an effective means of determining the source of groundwater components and their influencing factors. The hydrochemical types of groundwater were mainly HCO3-Ca, HCO3·SO4-Ca, HCO3-Ca·Mg, HCO3·SO4-Ca·Mg, SO4·HCO3-Ca·Mg, SO4·HCO3-Ca (Na), and SO4·HCO3-Na·Ca. HCO3 and SO42− were the dominant anions, and Ca2+, Mg2+ and Na+ were the dominant cations. All groundwater samples were distributed along the atmospheric precipitation line, suggesting an atmospheric precipitation origin for the groundwater. Groundwater in the stagnant zone and closed environment had the lightest δ18O-H2O and δD-H2O, indicating that the groundwater in this area was formed during a period of cooler climate. Groundwater with heavier δ18O-H2O and δD-H2O mainly receives recharge from SW.
The Gibbs diagram shows that the water–rock interaction is the main factor influencing the hydrochemical components of the groundwater, and the relationship diagrams for Mg2+/Na+ vs. Ca2+/Na+ and HCO3/Na+ vs. Ca2+/Na+ show that carbonate weathering and silicate weathering are the main sources of the chemical components of the groundwater. The ion ratio analysis also shows that the chemical components of groundwater mainly come from carbonate weathering and silicate weathering. The dissolution of minerals in groundwater is mainly controlled by carbolic acid, sulfuric acid, and nitric acid. Dissolution of salt rocks is not the main source of Na+ and K+ in groundwater in the Sangu Spring Basin, but dissolution of sodium feldspar (NaAlSi3O8) and potassium feldspar (KAlSi3O8) is the main source of Na+ and K+ in groundwater.
The ion exchange of groundwater is dominated by the exchange of Mg2+ and Ca2+ with Na+ and K+ in rocks or soil, and a small part of groundwater Na+ and K+ is exchanged with Mg2+ and Ca2+ in rocks or soil. Ion exchange mainly occurs in KGW and FGW, while ion exchange in PGW is weak.
Sulfate in the groundwater of the Sangu Spring Basin is affected by the combined effects of gypsum dissolution, infiltration of larger air precipitation, and leakage of SW. The sulfate in the KGW of the deeply buried and stagnant-flow zones is mainly derived from the dissolution of gypsum and has a higher sulfate content. Groundwater recharged by leakage from SW tends to have the highest sulfate content. The reducing effect of sulfur bacteria may have occurred in the FGW which is in a reducing environment.
In view of the environmental problems linked to groundwater in the Sangu Spring Basin, future work should strengthen the study of the hydraulic connection between KGW and SW, strengthen the protection of the quality of SW, and carry out monitoring of the hydrochemical anomaly zone of KGW, which is of great significance to the management of KGW resources in the Sangu Spring Basin.

Author Contributions

Methodology, Z.B. and X.H.; software, X.H. and Z.W.; validation, Z.B. and X.L.; formal analysis, X.L. and Z.B.; investigation, X.H., Z.B., Z.W., X.Z., C.G. and C.Z.; data curation, Z.W.; writing—original draft preparation, Z.B. and X.H.; writing—review and editing, Z.B. and X.H.; visualization, Z.B.; supervision, X.L.; project administration, X.H. and Z.W.; funding acquisition, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the China Geological Survey project (No. 12120114010701, DD20160296, 1212011140028, DD20190252, DD20221812, and DD20230539).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

We thank Jianfei Ma, Changchang Fu, Ming Gao and Juanjuan Li in our team for their help in the investigation of groundwater and the environment.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hou, X.W.; Li, X.Q.; Wang, Z.X.; Bai, Z.X.; Gui, C.L.; Zuo, X.F.; Zhang, C.C. Hydrogeological and Environment Geological Survey Report of Jindong Energy Base; Institute of Hydrogeology and Environmental Geology: Shijiazhuang, China, 2016; pp. 45–189. (In Chinese) [Google Scholar]
  2. Wang, Z.X.; Li, X.Q.; Hou, X.W.; Zhang, C.C.; Gui, C.L.; Zuo, X.F. Karst aquifer protection evaluation to the Sangu Spring Basin under the condition of coal mining. Carsologica Sin. 2019, 38, 28–29, (In Chinese with English abstract). [Google Scholar]
  3. Winter, T.C. Relation of streams, lakes, and wetlands to groundwater flow systems. Hydrogeol. J. 1999, 7, 28–45. [Google Scholar] [CrossRef]
  4. Glok-Galli, M.; Vadillo-Pérez, I.; Jiménez-Gavilán, P.; Ojeda-Rodriguez, I.; Urresti-Estala, B.; Martinez, D.E. Application of hydrochemical and multi-isotopic (87Sr/86Sr, δ13C-Dic, δ2H-H2O, δ18O-H2O) tools to determine contamination sources and processes in the Guadalhorce River Basin, southern Spain. Sci. Total Environ. 2022, 828, 154424. [Google Scholar] [CrossRef] [PubMed]
  5. Clark, I.D.; Fritz, P. Environmental Isotopes in Hydrogeology Lewis; Springer: New York, NY, USA, 1997. [Google Scholar]
  6. Liu, P.; Hoth, N.; Drebenstedt, C.; Sun, Y.J.; Xu, Z.M. Hydro-geochemical paths of multi-layer groundwater system in coal mining regions—Using multivariate statistics and geochemical modeling approaches. Sci. Total Environ. 2017, 600–601, 1–14. [Google Scholar] [CrossRef]
  7. Fu, C.; Li, X.; Ma, J.; Liu, L.; Gao, M.; Bai, Z. A hydrochemistry and multi-isotopic study of groundwater origin and hydrochemical evolution in the middle reaches of the Kuye River basin. Appl. Geochem. 2018, 98, 82–93. [Google Scholar] [CrossRef]
  8. Gao, M.; Li, X.; Qian, J.; Wang, Z.; Hou, X.; Fu, C.; Ma, J.; Zhang, C.; Li, J. Hydrogeochemical characteristics and evolution of karst groundwater in Heilongdong Spring basin, Northern China. Water 2023, 15, 726. [Google Scholar] [CrossRef]
  9. Olea-Olea, S.; Silva-Aguilera, R.A.; Alcocer, J.; Escolero, O.; Morales-Casique, E.; Florez-Peñaloza, J.R.; Almora-Fonseca, K.A.; Oseguera, L.A. Water-rock interaction processes in groundwater and flows in a Maar Lake in central Mexico. Water 2024, 16, 715. [Google Scholar] [CrossRef]
  10. Alogayell, H.M.; EL-Bana, E.M.M.; Abdelfattah, M. Groundwater quality and suitability assessment for irrigation using hydrogeochemical characteristics and pollution indices: A case study of North Al-Quwayiyah Governorate, Central Saudi Arabia. Water 2023, 15, 3321. [Google Scholar] [CrossRef]
  11. Keesari, T.; Roy, A.; Mohokar, H.; Pant, D.; Sinha, U.K. Characterization of mechanisms and processes controlling groundwater recharge and its quality in drought-prone region of central India (Buldhana, Maharashtra) using isotope hydrochemical and end-member mixing modeling. Nat. Resour. Res. 2020, 29, 1951–1973. [Google Scholar] [CrossRef]
  12. Qu, S.; Wang, G.C.; Shi, Z.M.; Xu, Q.Y.; Guo, Y.Y.; Ma, L.; Sheng, Y.Z. Using stable isotopes (D, 18O, 34S and 87Sr/86Sr) to identify sources of water in abandoned mines in the Fengfeng coal mining district, northern China. Hydrogeol. J. 2018, 26, 1443–1453. [Google Scholar] [CrossRef]
  13. Mohamed, O.A.; Tiziano, B.; Abdillahi, E.A.; Mahamoud, A.C.; Moussa, M.A.; Omar, A.D.; Youssouf, D.S.; Nima, M.E.; Ali, D.K.; Ibrahim, H.K.; et al. Origin of nitrate and sulfate sources in volcano-sedimentary aquifers of the East Africa Rift System: An example of the Ali-Sabieh groundwater. Sci. Total Environ. 2022, 804, 150072. [Google Scholar]
  14. Jin, M.G.; Zhang, J.; Zhang, Z.X.; Cao, M.D.; Huang, X. A review on source identification of dissolved sulfate in groundwater: Advances, problems and development trends. Bull. Geol. Sci. Technol. 2022, 41, 160–171. [Google Scholar]
  15. Hosono, T.; Siringan, F.; Yamanaka, T.; Umezawa, Y.; Onodera, S.i.; Nakano, T.; Taniguchi, M. Application of multiisotope ratios to study the source and quality of urban groundwater in Metro Manila, Philippines. Appl. Geochem. 2010, 25, 900–909. [Google Scholar] [CrossRef]
  16. Puig, R.; Soler, A.; Widory, D. Characterizing sources and natural attenuation of nitrate contamination in the Baix Ter Aquifer System (NE Spain) using a multi-isotope approach. Science. Total Environ. 2017, 580, 518–532. [Google Scholar] [CrossRef] [PubMed]
  17. Han, X.R.; Shi, J.; Tang, J.S.; Li, Q.S. Method of study on karst groundwater systems in north China—Case hsistory of Sangu Spring Basin. Carsologica Sinica 1990, 9, 197–210, (In Chinese with English abstract). [Google Scholar]
  18. Wang, Z.X.; Li, X.Q.; Hou, X.W.; Zhang, C.C.; Gui, C.L.; Zuo, X.F. Distribution Characteristics and Subsystem Identification of Karst Groundwater in Sangu Spring Basin. Earth Environ. 2020, 48, 228–239, (In Chinese with English abstract). [Google Scholar]
  19. Li, X.L.; Wu, Y.Z. Study on Karstic Groundwater Boundary Dividing in Sangu Spring Area. Ground Water 2003, 25, 224–240, (In Chinese with English abstract). [Google Scholar]
  20. Xu, Z.F.; Zhang, Z.X.; Liu, X.X. Water environmental quality assessment and strategy for water pollution control in Sangu spring basin. Ground Water 2012, 34, 87–90, (In Chinese with English abstract). [Google Scholar]
  21. Liu, X.H.; Zhang, Y.B. On the karst underwater water quality assessment of San-guquan in Jincheng of Shanxi. Shanxi Archit. 2009, 35, 173–174, (In Chinese with English abstract). [Google Scholar]
  22. Liu, X.H.; Zhang, Y.B. Groundwater resources evaluation on Sangu spring region of Jincheng of Shanxi Province. J. Taiyuan Univ. Sci. Technol. 2009, 30, 261–263, (In Chinese with English abstract). [Google Scholar]
  23. Zhang, Z.X.; Liu, X.H.; Zhang, Y.B. Research on dynamic characteristics and attenuation causes of flow Rate of Sangu Spring. Sci—Tech Inf. Dev. Econ. 2010, 20, 168–170, (In Chinese with English abstract). [Google Scholar]
  24. Tang, Q.F. Analysis of the decline factors in karst groundwater level of Sangu Spring Basin. Shanxi Water Resour. 2002, 50–51. (In Chinese) [Google Scholar]
  25. Li, H.Q. Analysis on the groundwater ecological status in the Sanguquan spring-feeding area of Jincheng City and restoration countermeasure. Shanxi Hydrotech. 2018, 2, 61–67. [Google Scholar]
  26. Qin, H.Z. Prediction of influence of coal mining on groundwater with numerical method—A case study of the Zhaozhuang Coal Mine in Jincheng area. Environ. Prot. Sci. 2013, 39, 53–56, (In Chinese with English abstract). [Google Scholar]
  27. Chen, F. Impact from Coal Winning on Spring Area Environment. Coal Geol. China 2020, 32, 34–37, (In Chinese with English abstract). [Google Scholar]
  28. Wang, Z.X.; Hou, X.W.; Li, X.Q.; Zhang, C.C.; Gui, C.L.; Zuo, X.F. Study on the influence of coal mining to north karst area. Yellow River 2019, 41, 75–82, (In Chinese with English abstract). [Google Scholar]
  29. Zhang, C.C.; Hou, X.W.; Li, X.Q.; Wang, Z.X.; Gui, C.L.; Zuo, X.F. Hydrogeochemical characteristics and evolution mechanism of karst groundwater in the catchment area of the Sangu Spring. Hydrogeol. Eng. Geol. 2021, 48, 62–71. [Google Scholar]
  30. Zhang, C.C.; Hou, X.W.; Li, X.Q.; Wang, Z.X.; Gui, C.L.; Zuo, X.F. Numerical simulation and environmental impact prediction of karst groundwater in Sangu Spring Basin, China. J. Groundw. Sci. Eng. 2020, 8, 210–222. [Google Scholar]
  31. Dansgaard, W. Stable isotopes in precipitation. Tellus 1961, 16, 436–468. [Google Scholar] [CrossRef]
  32. Craig, H. Isotopic variation in meteoric waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef]
  33. Yin, L.H.; Hou, G.C.; Dou, Y.; Tao, Z.P.; Li, Y. Hydrogeochemical and isotopic study of groundwater in the habor lake basin of the Ordos plateau, NW China. Environ. Earth Sci. 2011, 64, 1575–1584. [Google Scholar] [CrossRef]
  34. Umarani, P.; Ramu, A.; Kumar, V. Hydrochemical and statistical evaluation of groundwater quality in coastal aquifers in Tamil Nadu, India. Environ. Earth Sci. 2019, 78, 452. [Google Scholar] [CrossRef]
  35. Wang, X.D.; Zheng, W.D.; Tian, W.; Gao, Y.M.; Wang, X.Z.; Tian, Y.Q.; Li, J.S.; Zhang, X.Y. Groundwater hydrogeochemical characterization and quality assessment based on integrated weight matter-element extension analysis in Ningxia, upper Yellow River, northwest China. Ecol. Indic. 2022, 135, 108525. [Google Scholar] [CrossRef]
  36. Gibbs, R.J. Mechanisms controlling world water chemistry. Science 1970, 170, 1088–1090. [Google Scholar] [CrossRef]
  37. Liu, F.; Song, X.F.; Yang, L.H.; Han, D.M.; Zhang, Y.H.; Ma, Y.; Bu, H.M. The role of anthropogenic and natural factors in shaping the geochemical evolution of groundwater in the Subei Lake basin, Ordos energy base, northwestern China. Sci. Total Environ. 2015, 538, 327–340. [Google Scholar] [CrossRef]
  38. Gaillardet, J.; Dupr’e, B.; Louvat, P.; All`egre, C.J. Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers. Chem. Geol. 1999, 159, 3–30. [Google Scholar] [CrossRef]
  39. Xiao, J.; Jin, Z.D.; Wang, J.; Zhang, F. Hydrochemical characterisrics, controlling factors and solute sources of groundwater within the Tarim River Basin in the extreme arid region, NW Tibetan Tibertan Platiau. Quat. Int. 2015, 380–381, 237–246. [Google Scholar] [CrossRef]
  40. Ángel, R.P.; Álvaro, P.C.; Richard, S.A.; Diana, M.C.; Jorge, L.P.; Johonathan, S.C. Hydrochemical characterization of groundwater in the Loja Basin (Ecuador). Appl. Geochem. 2019, 104, 1–9. [Google Scholar]
  41. Irina, T.; Aleksei, K.; Alexander, Z.K.; Igor, C. Chemical composition of groundwater in abandoned coal mines: Evidence of hydrogeochemical evolution. Appl. Geochem. 2022, 137, 105210. [Google Scholar]
  42. Tanushree; Kumari, R. A hydrochemical and remote sensing approach to decrypt the groundwater salinization in the coastal district of Sabarmati basin, Gujarat. Groundw. Sustain. Dev. 2021, 15, 100673. [Google Scholar] [CrossRef]
  43. Owen, D.D.R.; Cox, M.E. Hydrochemical evolution within a large alluvial groundwater resource overlying a shallow coal seam gas reservoir. Sci. Total Environ. 2015, 523, 233–252. [Google Scholar]
  44. Xin, Y.; Hou, S.C.; Xiang, L.; Yu, Y.-X. Adsorption and substitution effects of Mg on the growth of calcium sulfate hemihydrate: An ab initio DFT study. Appl. Surf. Sci. 2015, 357, 1552–1557. [Google Scholar] [CrossRef]
  45. Schoeller, H.J. Geochemistry of Groundwater. Groundwater Studies—An International Guide for Research and Practice; UNESCO: Paris, France, 1977; pp. 1–18. [Google Scholar]
  46. Chaouki, K.; Meriem, T.; Lassâad, D. An appraisal of Chott El Jerid brine encroachment in the Tozeur-south shallow aquifer Geoelectrical and hydrochemical approach. J. Appl. Geophys. 2021, 190, 1–16. [Google Scholar]
  47. Ismail, A.H.; Hassan, G.; Sarhan, A. Hydrochemistry of shallow groundwater and its assessment for drinking and irrigation purposes in Tarmiah district, Baghdad governorate, Iraq. Groundw. Sustain. Dev. 2020, 10, 100300. [Google Scholar] [CrossRef]
  48. Li, Z.J.; Yang, Q.C.; Yang, Y.S.; Ma, H.Y.; Wang, H.; Luo, J.N.; Bian, J.M.; Martin, J.D. Isotopic and geochemical interpretation of groundwater under the influences of anthropogenic activities. J. Hydrol. 2019, 576, 685–697. [Google Scholar] [CrossRef]
  49. Awaleh, M.O.; Boschetti, T.; Soubaneh, Y.D.; Kim, Y.; Baudron, P.; Kawalieh, A.D.; Ahmed, M.M.; Daoud, M.A.; Dabar, O.A.; Kadieh, I.H.; et al. Geochemical, multi-isotopic studies and geothermal potential evaluation of the complex Djibouti volcanic aquifer (republic of Djibouti). Appl. Geochem. 2018, 97, 301–321. [Google Scholar] [CrossRef]
  50. Krouse, H.R.; Grinenko, V.A. Stable Isotopes: Natural and Anthropogenic Sulphur in the Environment; John Wiley & Sons, Ltd.: Chichester, UK; New York, NY, USA, 1991. [Google Scholar]
  51. Zak, D.; Hupfer, M.; Cabezas, A.; Jurasinski, G.; Audet, J.; Kleeberg, A.; McInnes, R.; Kristiansen, S.M.; Petersen, R.J.; Liu, H.J.; et al. Sulphate in freshwater ecosystems: A review of sources, biogeochemical cycles, ecotoxicological effects and bioremediation. Earth-Sci. Rev. 2020, 212, 103446. [Google Scholar] [CrossRef]
  52. Benadela, L.; Bekkoussa, B.; Gaidi, L. Multivariate analysis and geochemical investigations of groundwater in a semi-arid region, case of superficial aquifer in Ghriss Basin, Northwest Algeria. J. Groundw. Sci. Eng. 2022, 10, 233–249. [Google Scholar]
  53. Fernando, W.A.M.; Ilankoon, I.M.S.K.; Syed, T.H.; Yellishetty, M. Challenges and opportunities in the removal of sulphate ions in contaminated mine water: A review. Miner. Eng. 2018, 117, 74–90. [Google Scholar] [CrossRef]
  54. Galhardi, J.A.; Bonotto, D.M. Hydrogeochemical features of SW and groundwater contaminated with acid mine drainage (AMD) in coal mining areas: A case study in southern Brazil. Environ. Sci. Pollut. Res. 2016, 23, 18911–18927. [Google Scholar] [CrossRef]
  55. Kinnunen, P.; Kyllonen, H.; Kaartinen, T.; Mäkinen, J.; Heikkinen, J.; Miettinen, V. Sulphate removal from mine water with chemical, biological and membrane technologies. Water Sci. Technol. 2018, 194–205. [Google Scholar] [CrossRef] [PubMed]
  56. Zhang, J.; Jin, M.G.; Cao, M.D.; Huang, X.; Zhang, Z.X.; Zhang, L. Sources and behaviors of dissolved sulfate in the Jinan karst spring catchment in northern China identified by using environmental stable isotopes and a Bayesian isotope-mixing model. Appl. Geochem. 2021, 134, 105109. [Google Scholar] [CrossRef]
  57. Vitòria, L.; Otero, N.; Soler, A.; Canals, A. Fertilizer characterization: Isotopic data (N, S, O, C, and Sr). Environ. Sci. Technol. 2004, 38, 3254–3262. [Google Scholar] [CrossRef] [PubMed]
  58. Jurado, A.; Vàzquez-Suñé, E.; Soler, A.; Tubau, I.; Carrera, J.; Pujades, E.; Anson, I. Application of multi-isotope data (O, D, C and S) to quantify redox processes in urban groundwater. Appl. Geochem. 2013, 34, 114–125. [Google Scholar] [CrossRef]
  59. Otero, N.; Soler, A.; Canals, Á. Controls of δ34S and δ18O in dissolved sulphate: Learning from a detailed survey in the Llobregat River (Spain). Appl. Geochem. 2008, 23, 1166–1185. [Google Scholar] [CrossRef]
  60. Strawn, D.G.; Bohn, H.L.; O’Connor, G.A. Soil Chemistry, 4th ed.; John Wiley & Sons: New York, NY, USA, 2015; p. 381. [Google Scholar]
  61. Man, K.; Ma, Z.M.; Xu, X.J. Research on the mechanism of sulfate pollution of roundwater in Jiaozuo area. Appl. Mech. Mater. 2014, 665, 436–439. [Google Scholar] [CrossRef]
  62. Strebel, O.; Böttcher, J.; Fritz, P. Use of isotope fractionation of sulfate-sulfur and sulfate- oxygen to assess bacterial desulfurication in a sandy aquifer. J. Hydrol. 1990, 121, 155–172. [Google Scholar] [CrossRef]
  63. 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. [Google Scholar] [CrossRef]
  64. Vengosh, A.; Lindberg, T.T.; Merola, B.R.; Ruhl, L.; Warner, N.R.; White, A.; Di Giulio, R.T. Isotopic imprints of mountaintop mining contaminants. Environ. Sci. Technol. 2013, 47, 10041–10048. [Google Scholar] [CrossRef]
  65. Stempvoort, D.R.V.; Krouse, H.R. Controls of δ18O in sulphate: Review of experimental data and application to specific environments. Environ. Asp. Mine Wastes Short Course Ser. 1994, 31, 447–479. [Google Scholar]
  66. Huang, Q.B. The Carbon Sequestration Effect in Semi-Arid Karst Area: A Case Study of Liulin Spring Catchment, Shanxi Province. Ph.D. Thesis, China University of Geosciences, Wuhan, China, 2019. (In Chinese with English abstract). [Google Scholar]
  67. Zhou, J.W.; Zhang, Y.P.; Zhou, A.G.; Liu, C.F.; Cai, H.S.; Liu, Y.D. Application of hydrochemistry and stable isotopes (δ34S, δ18O and δ37Cl) to trace natural and anthropogenic influences on the quality of groundwater in the piedmont region, Shijiazhuang, China. Appl. Geochem. 2016, 71, 63–72. [Google Scholar] [CrossRef]
Figure 1. Exposed strata and distribution of sampling points in the Sangu Spring Basin. Note: In the legend, Karst groundwater, Surface water, Pore groundwater and Fissure groundwater all represent sampling locations.
Figure 1. Exposed strata and distribution of sampling points in the Sangu Spring Basin. Note: In the legend, Karst groundwater, Surface water, Pore groundwater and Fissure groundwater all represent sampling locations.
Water 16 02884 g001
Figure 2. Hydrogeological cross-section (A,B in Figure 1) of the study area.
Figure 2. Hydrogeological cross-section (A,B in Figure 1) of the study area.
Water 16 02884 g002
Figure 3. Piper diagram of groundwater in the Sangu Spring Basin.
Figure 3. Piper diagram of groundwater in the Sangu Spring Basin.
Water 16 02884 g003
Figure 4. δD −H2O versus δ18O − H2O plot of groundwater in the Sangu Spring Basin.
Figure 4. δD −H2O versus δ18O − H2O plot of groundwater in the Sangu Spring Basin.
Water 16 02884 g004
Figure 5. Gibbs diagrams representing controlling factors of groundwater quality, expressed in mg⋅L−1 (a) TDS vs. (Na/(Na + Ca)) and (b) TDS vs. (Cl/(Cl + HCO3)).
Figure 5. Gibbs diagrams representing controlling factors of groundwater quality, expressed in mg⋅L−1 (a) TDS vs. (Na/(Na + Ca)) and (b) TDS vs. (Cl/(Cl + HCO3)).
Water 16 02884 g005
Figure 6. Normalized bivariate diagrams for (a) Mg2+/Na+ vs. Ca2+/Na+ and (b) HCO3/Na+ vs. Ca2+/Na+.
Figure 6. Normalized bivariate diagrams for (a) Mg2+/Na+ vs. Ca2+/Na+ and (b) HCO3/Na+ vs. Ca2+/Na+.
Water 16 02884 g006
Figure 7. Plots of (Ca2+ + Mg2+) vs. (HCO3 + SO42− + CO32− + NO3) (a), (Ca2+ + Mg2+ + Na+ + K+ − Cl) vs. (HCO3 + SO42− + CO32− + NO3) (b), Ca2+ vs. SO42− (c), (Na+ + K+) vs. Cl (d), Mg2+ vs. SO42− (e) and Ca2+ vs. Mg2+ (f).
Figure 7. Plots of (Ca2+ + Mg2+) vs. (HCO3 + SO42− + CO32− + NO3) (a), (Ca2+ + Mg2+ + Na+ + K+ − Cl) vs. (HCO3 + SO42− + CO32− + NO3) (b), Ca2+ vs. SO42− (c), (Na+ + K+) vs. Cl (d), Mg2+ vs. SO42− (e) and Ca2+ vs. Mg2+ (f).
Water 16 02884 g007
Figure 8. Relationship plots of CAI1 vs. CAI2 (a) and ((Ca2++Mg2+) − (HCO3 + SO42− + CO32− + NO3)) vs. (K+ + Na+ − Cl) (b), expressed in meq L−1.
Figure 8. Relationship plots of CAI1 vs. CAI2 (a) and ((Ca2++Mg2+) − (HCO3 + SO42− + CO32− + NO3)) vs. (K+ + Na+ − Cl) (b), expressed in meq L−1.
Water 16 02884 g008
Figure 9. Plots of δ18O − SO42− (‰) vs. δ18O − H2O (‰) (a), δ34S − SO42− (‰) vs. SO42− (mg/L) (b), δ34SSO4 (‰) vs. 14C (PMC) (c), and δ34S − SO42− (‰) vs. δ18O − SO42− (‰) (d).
Figure 9. Plots of δ18O − SO42− (‰) vs. δ18O − H2O (‰) (a), δ34S − SO42− (‰) vs. SO42− (mg/L) (b), δ34SSO4 (‰) vs. 14C (PMC) (c), and δ34S − SO42− (‰) vs. δ18O − SO42− (‰) (d).
Water 16 02884 g009
Table 1. Hydrochemical components of groundwater in the Sangu Spring Basin.
Table 1. Hydrochemical components of groundwater in the Sangu Spring Basin.
K+Na+Ca2+Mg2+ClSO42−HCO3NO3TDSpH
KGW (n = 137)
Min0.263.5449.3712.590.3621.86159.400.00218.707.20
Max20.40198.00487.70130.00232.201413.00441.10378.502411.008.57
Mean1.5043.61126.2333.4224.18230.89301.2123.72653.187.73
Std2.4552.8860.6415.7126.07207.2539.6438.04337.310.25
FGW (n = 74)
Min0.074.6254.128.434.5823.610.801.98245.806.71
Max7.78282.60307.0051.40230.10846.00465286.002424.508.51
Mean0.9667.49136.6023.3230.85253.30293.5335.75776.237.79
Std1.0974.4456.5710.3335.11213.3485.121.56408.260.33
PGW (n = 23)
Min0.2810.8058.3012.048.8140.38250.001.63330.507.35
Max9.45196.00286.9060.38101.00521.00478.60227.501226.008.28
Mean1.4757.00145.7326.3937.49233.21324.1045.52796.107.78
Std1.9257.8154.3211.2625.24135.7650.8947.37231.400.26
Note on units: concentrations are expressed in milligram per liter (mg/L), except for pH.
Table 2. δD-H2O and δ18O-H2O isotope composition of all the water samples.
Table 2. δD-H2O and δ18O-H2O isotope composition of all the water samples.
KGW (n = 112)FGW (n = 32)
MinMaxMeanStdMinMaxMeanStd
δD-H2O (‰)−83.0−42.0−70.54.03−83.0−57.0−67.74.68
δ18O-H2O (‰)−11.2−4.5−9.60.57−11.57−8.1−9.20.67
PGW (n = 8)SW (n = 16)
MinMaxMeanStdMinMaxMeanStd
δD-H2O (‰)−69.0−48.0−63.06.74−70.0−42.0−56.38.16
δ18O-H2O (‰)−9.6−5.8−8.51.17−9.1−3.9−7.41.38
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bai, Z.; Hou, X.; Li, X.; Wang, Z.; Zhang, C.; Gui, C.; Zuo, X. Hydrogeochemical Characteristics and Sulfate Source of Groundwater in Sangu Spring Basin, China. Water 2024, 16, 2884. https://doi.org/10.3390/w16202884

AMA Style

Bai Z, Hou X, Li X, Wang Z, Zhang C, Gui C, Zuo X. Hydrogeochemical Characteristics and Sulfate Source of Groundwater in Sangu Spring Basin, China. Water. 2024; 16(20):2884. https://doi.org/10.3390/w16202884

Chicago/Turabian Style

Bai, Zhanxue, Xinwei Hou, Xiangquan Li, Zhenxing Wang, Chunchao Zhang, Chunlei Gui, and Xuefeng Zuo. 2024. "Hydrogeochemical Characteristics and Sulfate Source of Groundwater in Sangu Spring Basin, China" Water 16, no. 20: 2884. https://doi.org/10.3390/w16202884

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