Next Article in Journal
Assessing Groundwater Mineralization Process, Quality, and Isotopic Recharge Origin in the Sahel Region in Africa
Previous Article in Journal
Assessing the Changes of Ecosystem Services in the Nansi Lake Wetland, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Control of Groundwater Flow Systems and Geochemical Processes on Groundwater Chemistry: A Case Study in Wushenzhao Basin, NW China

1
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
3
Institutions of Earth Science, Chinese Academy of Sciences, Beijing 100029, China
4
Xi’an Center of Geological Survey, China Geological Survey, Xi’an 710054, China
5
Key Laboratory for Groundwater and Ecology in Arid and Semi-arid Areas, Xi’an Center of Geological Survey, China Geological Survey, Xi’an 710054, China
6
Mining Engineering and Geology College, Xinjiang Institute of Engineering, Wulumuqi 830023, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 790; https://doi.org/10.3390/w11040790
Submission received: 31 March 2019 / Revised: 5 April 2019 / Accepted: 8 April 2019 / Published: 16 April 2019
(This article belongs to the Section Hydrology)

Abstract

:
The lowest reaches of a large-scale basin could be the discharge areas of local, intermediate and regional groundwater flow systems with significantly distinct travel distances and travel times. This study aims to delineate the groundwater chemical characteristics and the mechanism controlling the chemical evolution in the lowest reaches of the Wushenzhao Cretaceous basin, NW China. A total of 38 groundwater samples were collected and were chemically classified into five distinct water types by means of a Piper Plot. According to the hydrogeological setting and groundwater age, the spatial distribution of these water types is found to be associated with hierarchically nested groundwater flow systems (local and regional system): Types 1, 2, 3 and 4 belong to the local groundwater flow system, while type 5 belongs to the regional flow system. Graphical plots, stable isotopes and geochemical modeling techniques were used to interpret the observed compositions. The results show the dominance of carbonate and gypsum dissolution in type 1 waters; ion exchange in types 2, 3 and 4; and evaporite dissolution in type 5. In addition, human activities in the form of extensive irrigation also affect the chemical compositions of type 1 water. These findings are important for the sustainable management of groundwater resources in the study area.

1. Introduction

Geochemical processes taking place within the aquifer have a profound effect on groundwater chemical evolution [1,2]. An increased knowledge of groundwater flow and geochemical processes could contribute to the improved understanding of groundwater system in a region, which is of critical importance for the effective management and sustainable development of groundwater resources [3,4,5,6,7,8,9]. During the last decades, the integrated application of hydrochemistry and environmental isotopes has been proven to be a powerful tool for solving these hydrological problems [10,11,12,13].
Ordos Plateau is an important energy base located in the semiarid and arid regions, NW China (Figure 1a), which conducts plenty of industrial activities using substantial water resources [14]. Groundwater is the only reliable source of water for drinking, irrigation, industry and ecosystem [15]. In recent years, due to the rapid development of the energy base, the demand for water resources has increased tremendously. As a result, the groundwater level is increasingly declining, and the ecosystem environment is worsening, especially the Wushenzhao basin in the northeast of Plateau (Figure 1b) [16]. This may lead to the deterioration of groundwater quality in the future [8]. To help the sustainable development of groundwater (i.e., maintaining both yield and quality), a thorough understanding of the groundwater system, particularly the hydrochemical system, is essential in those Lake Basins.
The available geochemical investigations of groundwater have been just carried out at a regional scale in the Ordos Plateau [15,17,18], while detailed investigations on the groundwater geochemistry are rather scattered so far, and analyses conducted are limited to the major ions chemistry of groundwater in the Lake Basins, NE Ordos Plateau, such as Dake Lake, Mukai Lake and Hutongchahan (HTCH) lake (Figure 1b) [14,16,19,20,21,22]. In addition, the HTCH Lake Basin is chosen as the key study area (Figure 1c) because it is located in the lowest reaches of the Wushenzhao Basin and might serve as the fingerprint of nested flow systems. The purpose of this study is to delineate the flow and hydrochemical evolution of groundwater in the HTCH Lake Basin. A multi-tracer approach including ion chemistry and environmental isotopes (e.g., δ2H, δ18O, 87Sr/86Sr, δ34S, 3H and 14C), together with geochemical modeling techniques, was employed. This work is the first detailed study on the groundwater chemistry in the study area and correlates the hydrochemical pattern with the groundwater flow systems. Moreover, groundwater age as well as improved potentiometric head contour is also first reported. These findings will significantly improve the sustainable development of groundwater in the study area.

2. Study Area

The Ordos Plateau is located in the north part of Ordos Basin, NW China (Figure 1a). It is one of the most important energy and industry bases due to an abundance of minerals resources, such as coal, natural gas, petroleum and halite [18,23]. The Plateau can be divided into three groundwater systems by the Sishi Ridge and Xinzhao Divide, namely the Wuding system, the Dosit system and the Molin system, as shown in Figure 1b [15]. The HTCH Lake Basin is located in the northeast part of Ordos Plateau and the north part of the Wuding groundwater system (Figure 1b). It belongs to the continental semiarid to arid climate: The annual precipitation averages about 340 mm/year; the mean annual potential evaporation is around 2300 mm/year, which far exceeds the precipitation.
Quaternary lacustrine sediments, aeolian sands and Cretaceous formations can be observed in the study area (Figure 2a). The first two are relatively thin, while the latter has a huge thickness of about 700–1000 m, which is further divided into the Huanhe and Luohe groups from the top to bottom (Figure 2a). The major aquifer is the poorly consolidated Cretaceous sandstone with a high permeability in the Plateau. Lens of mudstone are present in the Cretaceous formation, but they cannot constitute a regional aquitard due to their limited extent and thin thickness. Therefore, the Cretaceous aquifer is usually considered as a relatively homogenous and unconfined system [18]. The dominant minerals of aquifer rocks are quartz, Na- and K-feldspars, and calcite and dolomite [14].
The groundwater flow is driven by the wavy topography in the Ordos Plateau, which develops hierarchically nested groundwater flow systems. The local flow system generally develops between groundwater divides and neighboring depressions with a shallow circulation depth of around 200 m; intermediate flow systems occurs between the major divides and the relatively larger lakes with a relatively deep circulation depth; and the regional flow system has a much larger lateral extent and a higher circulation depth (up to 500 m) when compared with the local and intermediate flow systems [14,15,18,20,22].
There is a hypothesis of groundwater flow systems in the lake-concentrated regions (Figure 1c): It also develops hierarchically nested groundwater flow systems that include local, intermediate and regional flow systems under the control of the groundwater-fed lakes as discharge areas, as shown in Figure 2a [24,25]. The local flow system usually occurs between groundwater divides and neighboring lakes (e.g., Dake Lake and Mukai Lake) with a relatively shallow circulation depth, while, intermediate and regional flow systems are present between the major divides and the relatively larger lakes such as HTCH Lake, with significantly deeper circulation depths and longer flow distances relative to local flow systems. As shown in Figure 2a, the study area (HTCH Lake Basin) located in the lower reaches of the basin probably develops local, intermediate and regional flow systems (Figure 2). According to the potentiometric contours map of the shallow groundwater (Figure 1d), groundwater flows along the cross section A′–A″ from northwest to southeast, which agrees well with the topography in the study area (Figure 1d).

3. Materials and Methods

3.1. Sampling and Analysis

An important sampling campaign was conducted in the HTCH basin during September 2017. A total of 38 groundwater samples were collected along the cross section A′–A″ following the groundwater flow direction. Thirty-three samples were gathered from active water wells with depths from about 10 m to 200 m, and 5 samples were collected from deep boreholes with depths varying between about 350 m to 600 m. In addition, a total of 5 lake water samples were taken from the HTCH Lake. The sampling sites are shown in Figure 1c,d. The shallow groundwater samples are considered belonging to the same groundwater flow system, namely the local flow system, according to the hydrogeological context reported in previous study [14,15]. Meanwhile, in order to better understand the chemical evolution of groundwater, a total of 5 rock core samples with different depths were collected from a deep borehole, as shown in Figure 1c.
All groundwater samples were collected at a well head using pumps installed in the wells. In both cases, the wells were allowed to discharge for a period of time prior to sampling to ensure the removal of stored water from the well column. In situ cellulose membrane filters (0.45 μm, Jinteng, Tianjin Municipality, China) were used to filter the samples for analyses of water chemistry and Sr isotope; raw water samples without filtering were collected in sampling polyethylene bottles for hydrogen, oxygen, sulfur, carbon isotope and tritium determinations. The sampling bottles were prerinsed with water samples three times before the final water was gathered. Finally, all water sample bottles were sealed with parafilm tape in order to prevent them from evaporation.
Total dissolved solids (TDS), pH, electrical conductivity (EC), Eh and temperature parameters were measured in the field by a portable Multiparameter Digital Meter (HQ40D, Hach, Loveland, CO, USA) that was calibrated using standard solution before use. These field parameters were recorded in situ until they reached a steady state. The analyses for the groundwater performed in the laboratory included water chemistry and isotopes (δ2H, δ18O, 87Sr/86Sr, δ34S, 14C and 3H). All lake water samples were measured for major ion chemistry (K+, Na+, Ca2+, Mg2+, Cl, SO42−, HCO3, CO32−, and NO3). All groundwater samples were analyzed for major ion chemistry and stable water isotopes (δ2H and δ18O). Minor ion concentrations and Sr isotope ratios were determined for 15 representative samples. Ten samples were selected for the carbon isotopic analysis (δ13C and 14C), and tritium (3H) was measured for 7 representative samples. Major ions except HCO3 and CO32− were analyzed using Ion Chromatography (Dionex ICS-2100, Sunnyvale, CA, USA) in the Water Isotope and Water–Rock Interaction Laboratory, Institute of Geology and Geophysics, Chinese Academy of Sciences (IGG, CAS). The alkalinity was measured by the titration method using 0.0048 M HCl immediately after returning the same laboratory. The charge balance error of major ions for all groundwater samples was within ±5%. The δ2H and δ18O compositions were analyzed using a liquid water isotope laser analyzer (Picarro L2120-i, Picarro, Inc., Santa Clara, CA, USA) in the same laboratory. The isotopic data was expressed in the standard delta notation (δ) as per mil (‰) relative to the Vienna Standard Mean Ocean Water (VSMOW, 0‰), and the measurement precision was ±0.5‰ for δ2H and ±0.1‰ for δ18O. The minor elements were measured by inductively coupled plasma mass spectrometry (ICP-MS NexION300D, PerkinElmer, Inc., Shelton, CT, USA) at the Beijing Research Institute of Uranium Geology, and the analytical precision is ±3% of the concentration based on the reproducibility of samples and standards. The 87Sr/86Sr ratios were determined using a thermal ionization mass spectrometer (Phoenix, ThermoFisher Scientific, Waltham, MA, USA) at the Beijing Research Institute of Uranium Geology with the analytical standard error (2σ) less than 0.00003. The δ34S values for dissolved SO4 in groundwater were measured by MAT 253 stable isotope mass spectrometry (ThermoFisher Scientific, Waltham, MA, USA) in the stable isotopes laboratory of IGG-CAS, and the analytical precision was approximately 0.2‰. The tritium (3H) was also measured at the Beijing Research Institute of Uranium Geology with the detection limit of 1.3 TU. The δ13C and 14C was determined in the Beta AnalyticTM in the United States by accelerator mass spectrometry (AMS) with an analytical precision of ±0.1‰ for δ13C and 0.1 pMC (percent Modern Carbon) for 14C. In addition, the aquifer mineralogical compositions were determined using an X-ray Diffraction (XRD) method at Lanzhou Center for Oil and Gas Resource, Institute of Geology and Geophysics, Chinese Academy of Sciences.

3.2. Inverse Geochemical Modeling

Geochemical modeling provides a useful tool for the quantitative evaluation of the chemical evolution in a particular groundwater system [26,27]. Inverse geochemical modeling is widely used to quantitatively simulate the key geochemical processes that can explain the differences observed between the initial and final water compositions along the flow path in a groundwater system [28,29]. In this study, the inverse geochemical modeling was conducted by means of the geochemical modeling program PHREEQC (version 2) developed by the United States Geological Survey (USGS) [30]. Furthermore, the saturation indices (SI) for minerals were also calculated with the help of this program.

4. Results and Discussion

The well locations, field data and major ion chemistry of groundwater and lake water are summarized in Table S1 in the Supplementary Material (SM). Table S2 in the SM presents the minor ion chemistry, isotopic compositions (δ2H, δ18O, 87Sr/86Sr, δ34S, 3H, and 14C) and SI, along with the calculated 14C age of groundwater in the study area.

4.1. Hydrochemical Pattern Resulting from Groundwater Flow Systems

Based on the available chemical data, the hydrochemical characteristics of HTCH lake water and groundwater were discussed.
The chemical data for lake water show that TDS, a comprehensive indicator of hydrochemistry, is relatively high and ranges from 4464.5 mg/L to 5165.7 mg/L with the mean value of 4719.2 mg/L. The pH values lie between 9.97 and 10.03 with an average of 10.00, indicating an alkaline nature. The water type is Na–Cl(HCO3), as shown in Piper plot (Figure 3). Na+ and Cl are the dominant elements: Na+ accounts for, on average, 88.6% of total cations, and Cl accounts for, on average, 51.3% of total anions. This is in accordance with hydrochemical characteristics of inland salt lake in other regions [16].
The TDS of groundwater range between 171.8 and 724.0 mg/L with mean and median values of 300.7 mg/L and 234.9 mg/L, respectively. Groundwater is slightly alkaline and has pH values varying from 7.69 and 9.39. The measured Eh values of groundwater range from −137.50 to −55.60 mV with the mean value of −90.34 mV. Ca2+, and the Na+ concentrations in groundwater vary in a wide range from 2.6 to 64.3 mg/L and 14.9 to 225.3 mg/L, respectively. The Cl and NO3 concentrations also exhibit a wide range from 3.5 to 60.9 mg/L and from ND (not detected) to 67.8 mg/L, respectively.
The Piper plot is widely used to graphically display the bulk major ion compositions of groundwater, and it is, therefore, used to classify groundwater [31]. Figure 3 shows that they are found to be chemically classified into five distinct types: Ca(Mg)–HCO3 (type 1), mixed cations-HCO3 (type 2), Na–HCO3 (type 3), Na–HCO3(SO4) (type 4) and Na–SO4 (type 5). In this classification scheme, type 1 through type 5 represent a natural chemical evolutionary sequence (Figure 3). The TDS and major ion concentrations for the 5 water types are also show illustrated in the Box plots of Figure 4. The calculated average values of the geochemical data for the 5 types of groundwater, including TDS, temperature, pH, Eh, major ions and stable isotopes (δD, δ18O and 87Sr/86Sr) are summarized in Table S3. The hydrochemical characteristics of each type of groundwater are discussed in detail below.
Type 1 waters have a Ca(Mg)–HCO3 composition and average TDS of 234 mg/L (Table S3). These waters contain more than 70% (in meq/L) Ca2+ + Mg2+ of the total cations content. HCO3 is the major anion component, which accounts for more than 60% of total anions contents. This type predominantly occurs in the topographic highs of the basin, as shown in Figure 2b. Type 2 waters are mixed cations-HCO3, and no dominant cation exists, belonging to the transitional type. The anion relationship is rHCO3 > rSO42− > rCl. Type 3 is Na–HCO3 water. Na+ is the major cation, which constitutes more than 54% of the total cations; HCO3, the major anion, accounts for more than 67% of the total anions. These waters are low in SO42 content that makes up less than 23% of the total anions. Type 4 is Na–HCO3(SO4) water with the mean TDS of 305.2 mg/L. Na+ is the predominant cation, and it accounts for greater than 77% of the total cation contents; HCO3 is the predominant anion, and it occurs in amounts greater than 48% of the total anion contents. The SO42− constitutes more than 23% of the total anions, which is considerably higher relative to type 3 waters. These three types, including type 2, 3 and 4, are predominantly located in the topographic lows, as shown in Figure 2b. Type 5 is Na–SO4 water, with the mean TDS value of 692.8 mg/L. Na+ contents range from 213.8 to 258.0 mg/L with the average value of 225.03 mg/L and constitute, on average, about 91% of total cations; SO42− concentrations vary between 322.2 and 608.7 mg/L with the mean of 387.11 mg/L, which accounts for, on average, 68% of total anions. The concentrations of TDS, Na+ and SO42− in type 5 are significantly higher than types 1, 2, 3 and 4, as illustrated in Figure 4. In summary, type 1 is characterized by and dominated by Ca2+, while types 3, 4 and 5 are characterized by Na+ among cations; types 1, 2, 3 and 4 are characterized by a predominance of HCO3, whereas type 5 is characterized by SO42 among anions (Figure 3). In addition, it is should be noted that the mean Cl concentrations are much higher in type 1 (16.4 mg/L) than type 2 (9.0 mg/L) and type 3 (8.2 mg/L), as shown in Table S3 and Figure 4e, and is not consistent with the general pattern that Cl contents increase along the flow line. This is attributed to the fact that some samples in type 1 waters have been strongly affected by anthropogenic contamination, such as agricultural fertilizer, as indicated by higher mean NO3 contents in type 1 (range: 13.6–67.8 mg/L; mean: 36.1 mg/L) than type 2 (range: 3.1–31.1 mg/L; mean: 18.5 mg/L) and type 3 (range: 13.6–180.1 mg/L; mean: 2.4 mg/L) (Table S3 and Figure 5). A similar contamination phenomenon caused by human activities is also observed in the regions close to the study area in the Ordos Plateau [20,32].
Broad relationships between the hydrochemical pattern and groundwater flow systems have been demonstrated by various workers. According to theory and observations, the groundwater flow from the recharge to discharge zones may be attended generally by an increase in the TDS and a systematic variation in the dominant cations and anions [33,34,35]. On the other hand, in a complex basin with the development of three types of groundwater flow systems (i.e., local, intermediate and regional), as described by Toth [36], groundwater from local flow systems may exhibit significantly distinct hydrochemical characteristics from those from intermediate/regional flow systems. This is because of the significant difference in the length of the flow paths and, hence, the contact time between groundwater and rock, especially in the lower reaches of a basin where a regional flow system lies below a local system [37,38]. Generally, systems may be characterized as local by low-TDS bicarbonate waters, as intermediate by intermediate-TDS-sulfate waters and as regional by high-TDS chloride-rich waters.
When we correlate the water types with the geographic locations, types 1, 2, 3 and 4 characterized by a predominance of HCO3 were found to be located in the shallow aquifer (<200 m), whereas all samples of type 5 dominated by SO42− occurred in the deep parts (around 350–600 m), as shown in Figure 2b and Figure 5. Such a vertical pattern of groundwater chemistry is in accordance with the propositions of Freeze [37] and Toth [38]. In addition, groundwater is of a meteoric origin, and its chemical composition is predominantly controlled by a water–rock interaction within the aquifer in the study area that is relatively simple and macroscopically homogeneous sandstone (see Section 4.3). Therefore, the significant geochemical differences between shallow and deep groundwater seem to support the hypothesis of nested flow systems with different flow distances and flow times (Figure 2). That is, the groundwater chemistry suggests that shallow groundwater (types 1, 2, 3 and 4) belongs to a local flow system with a relatively short travel time and distance, while deep groundwater (type 5) belongs to a regional flow system with a longer travel time and distance. The TDS and Cl content, traditional indicators of groundwater residence time, are significantly higher in deep groundwater compared to shallow groundwater, as shown in Figure 5a,g, which also support this interpretation.

4.2. Groundwater Origin

The environmental isotopes of deuterium (2H) and oxygen-18 (18O) are excellent tracers for determining the groundwater origin [39]. A local meteoric water line (LMWL: δ2H = 6.45δ18O − 6.51) for the Ordos Plateau basin defined by Yin et al. [23] provides the basis for the interpretation of δ2H and δ18O values for groundwater in this study. The lower slope of LMWL relative to the global meteoric water line (GMWL) (δ2H = 8δ18O − 10) is attributed to the secondary evaporation process during rainfall as a result of an arid environment.
The δ2H and δ18O values in shallow groundwater exhibit a relatively wide range from −86.5‰ to −57.5‰ and from −11.3‰ to −7.4‰, respectively, which is within the general range previously observed for shallow groundwater in the Ordos Plateau [23,40]. The plot of δ2H versus δ18O (Figure 6) shows that all groundwater samples lie close to the LMWL and GMWL, indicating that they are of local meteoric origin. However, they can be divided into two groups. One group (group 1) is the type 1 groundwater located in the recharge area of local flow system; the average stable isotopic composition (−63.3‰ for δ2H and −8.3‰ for δ18O) in these groundwater corresponds to that of modern rainfall in the Ordos Plateau with a weighted mean of −58.5‰ for δ2H and −8.1‰ for δ18O (Figure 7) [23]. This suggests that they are recharged primarily by modern precipitation or precipitation in the late Holocene period. The other group (group 2) includes the other four types of groundwater (types 2, 3, 4 and 5). They are depleted in heavy isotopes (mean −80.9‰ for δ2H and −10.4‰ for δ18O) in comparison with type 1 groundwater and present meteoric water (Figure 7), suggesting that they are primarily recharged during the late Pleistocene and early Holocene periods when climate was relatively colder and wetter than the present, after having considered that the altitude effect was not significant in the whole Ordos Plateau [23,40]. In addition, in group 2, type 5 waters are found to have isotopic compositions very similar to type 4 and eight samples in type 3 because they are both old waters, while type 2 waters and four samples in type 3 display the isotopic compositions within those of type 1 and those 8 samples in type 3 (type 4, 5), which may be caused by a mixing between type 1 and type 3 based on the chemical characteristics and hydrogeological setting.

4.3. Mechanisms Controlling Groundwater Chemical Evolution

This section aims to identify the dominant mechanisms responsible for the chemical compositions of groundwater of various types. Understanding geochemical processes in groundwater is a prerequisite for the correction of 14C ages and for recording groundwater recharge and flow. Processes controlling the hydrochemical evolution were determined using hydrochemistry, stable isotope and inverse geochemical modeling.

4.3.1. Constraints from Water Chemistry

Generally, the chemical compositions of groundwater are primarily controlled by various factors, including the chemical composition of recharge water, water–rock interaction, dissolution/precipitation of mineral phases, evaporation, mixing processes and human activities [1,2]. The Gibbs diagram [41], which exhibits the ratios of Na2+/(Na2+ + Ca2+) and Cl/(Cl + HCO3) as a function of TDS, is an effective tool to assess the natural mechanisms controlling water chemistry, including the atmospheric precipitation, the water–rock interaction and the evaporation-crystallization process. Figure 8 shows that all groundwater samples fall into the zone of rock dominance, suggesting that water–rock interactions are the major mechanism controlling groundwater chemistry. The Cretaceous sandstones consist mainly of quartz; K-feldspar; plagioclase; calcite; and minor amounts of dolomite, ankerite, siderite, pyrite, hematite, analcite, Glauber’s salt and clay minerals, as presented in the Table S4. Therefore, potential water–rock interactions include the dissolution/weathering of these minerals forming aquifer and cation exchange occurring on the clay minerals. In addition, the dissolution of halite and gypsum can also take place in the aquifer, although they are not found in the study area. This is due to the fact that they have been observed in the Cretaceous aquifer in the Ordos Plateau [15].
In order to further reveal the origin of solutes in groundwater, the relationship between the dissolved chemical constituents is employed [42]. The relationship between Na+ and Cl is commonly used as a first order indicator of water–rock interactions [43]. The plot of Na+ versus Cl (Figure 9a) shows that only three samples in type 1 groundwater samples (about 25% of type 1 waters) lie close to the 1:1 relation line expected from the dissolution of halite, indicating the occurrence of halite dissolution in these waters, while other samples cluster above the 1:1 line and molar Na/Cl ratios vary between 1.5 and 30.2 with a mean value of 8.7. This indicates the existence of additional sources contributing to the dissolved Na+, which potentially include the weathering of Na-bearing silicates, a dissolution of Na-bearing evaporites and cation exchange between Ca2+/Mg2+ in groundwater and Na+ on clay minerals after taking into consideration the mineralogical compositions of aquifer studied [44,45,46].
The plot of (Ca2+ + Mg2+) versus (HCO3 + SO42–) (Figure 9b) indicates that type 1 groundwater is distributed around the 1:1 relation line, suggesting that Ca, Mg and HCO3 in these waters are derived mainly from simple dissolution of calcite, dolomite and gypsum. The mean SI values for calcite and dolomite are 0.71 and 0.93, respectively, indicative of a slight oversaturation of groundwater with carbonates minerals and the importance of carbonates dissolution in these waters. These reactions commonly occur in the recharge areas in many sedimentary aquifer environments, and the reaction equations are as follows [2]:
CaCO3 + CO2 + H2O = Ca2+ + 2HCO3
CaMg(CO3)2 + 2CO2 + 2H2O = Ca2+ + Mg2+ + 4HCO3
CaSO4 = Ca2+ + SO42−
In contrast, type 2, 3, 4 and 5 waters exhibit a significant deviation from the 1:1 relation line and plot above it, indicating the simple dissolution of carbonate and gypsum alone could not interpret their chemical compositions. This may be caused by the cation exchange, silicates weathering or evaporites dissolution [47]. Figure 9c shows that majority of these waters plot above or below the theoretical line of weathering of plagioclase (y = x), suggesting silicate weathering is not the major source of Na+. As shown in Figure 9d, type 2, 3, 4 and 5 waters have a much higher Na/Ca ratio relative to type 1 water and TDS slightly increases as groundwater evolves from type 1 to types 2, 3 and 4. This seems to imply that a cation exchange has taken place as groundwater flows in the aquifer. The plot of (Na+ − Cl) against (Ca2+ + Mg2+) could be employed to evaluate the cation exchange [46,48]. If the relationship between these two parameters is linear with a slope of −1.0, suggesting a cation exchange is significant in controlling geochemical compositions. As shown in Figure 9e, types 1, 2, 3 and 4 are located near the 1:1 relation line; in contrast, type 5 groundwater lie away from the 1:1 line. This mirror image of mono- (Na+) and bivalent cations (Ca2+, Mg2+), along with the almost constant HCO3 contents along the groundwater flow path (Figure 5e), indicate that a cation exchange controls the hydrochemical evolution of shallow groundwater (types 1, 2, 3 and 4) along the flow path; the reaction equations are as follows:
NaX + Ca2+/Mg2+ = CaX2/MgX2 + Na+
The use of two Chloro-Alkaline indices (CAI-1 and CAI-2) proposed by Schoeller is helpful for indicating the cation exchange process [49]. CAI-1 and CAI-2 are defined as the ratio of (Cl − (Na + K)) to Cl and (HCO3 + SO4 + CO3 + NO3), respectively. The positive indices suggest that the ion exchange occurs between Na+ in the groundwater and Ca2+/Mg2+ in the aquifer sediments, while the negative values suggest that the reverse ion exchange takes place. As shown in Figure 9f, both of the indices (CAI-1 and CAI-2) are found to be below zero in the type 2, 3 and 4 groundwater samples, implying that a high concentration of Na+ and a low concentration of Ca2+ are attributed to the ion-exchange process in these waters.
Figure 9g shows that deep type 5 water displays a significant deviation from the 1:1 line expected from gypsum dissolution and plots above the 1:1 line, indicating that gypsum is not the major process in this type of water and that SO42− has additional sources. The SI values for gypsum in these waters range from −1.85 to −1.67 with an average of −1.79, furthering supporting this interpretation. However, these waters have a relatively higher average SI value for gypsum than type 1 (−2.33), type 2 (−2.46), type 3 (−2.84) and type 4 (−2.84), as shown in Table S3, suggesting a stronger degree of gypsum dissolution in this type than in the other four types. Figure 9h shows that type 5 groundwater plots relatively close to the 1:1 relation line, suggesting Na+, SO42− and dominant ions in type 5 primarily derive from the dissolution of evaporite including halite and Glauber’s salt (Na2SO4). The reaction equation of Glauber’s salt dissolution is as follows:
Na2SO4 = 2Na+ + SO42−
Another significant characteristic of the groundwater is the extremely low K+ content between 0.6 and 4.4 mg/L (mean 2.3) and the very low K/Na molar ratios of near 0. The low levels of potassium may be attributed to its tendency to participate in the formation of secondary minerals and the lower reactivity of K-feldspars compared to plagioclase [8,43].
The dissolved silica concentrations (range: 1.66–11.90 mg/L; mean: 7.14 mg/L) are generally lower than the typical case where silicate weathering is a major process in groundwater (approx. 20 to >100 mg/L) [9], suggesting that the silicate weathering is a minor source of Na+ in groundwater.

4.3.2. Constraints from Strontium Isotope

The Sr isotopes are extensively applied as tracers to effectively delineated mineral weathering reactions in subsurface hydrology [50,51,52]. This is attributed to the following properties. The strontium can be removed from the water through mineral precipitation or cation exchange reaction, but the Sr isotopic composition of water remains largely unaffected because its mass-dependent fraction is usually negligible in nature. Moreover, Sr2+ is geochemically similar to Ca2+ and usually resides in carbonates, gypsums and K-bearing or Ca-bearing silicates [51,53]. These common rock-forming minerals generally exhibit a wide and predictable range of 87Sr/86Sr ratios that are controlled by initial values, Rb/Sr ratios and the age of minerals. Silicate minerals, particularly potassic silicates (e.g., biotite and K-feldspar), that exceed a few million years old have more radiogenic 87Sr/86Sr ratios due to the addition of radiogenic 87Sr through the decay of 87Rb (t1/2 = 4.8 × 1010 years). In contrast, Ca-rich minerals such as calcite or gypsum display relatively less radiogenic 87Sr/86Sr ratios and remain relatively constant over time.
Groundwater primarily originates from the local precipitation in the study area, as indicated by δD and δ18O. Therefore, the strontium geochemistry of groundwater is initially controlled by the precipitation. The Sr concentrations in precipitation samples collected from the Otak meteorological station (Figure 1b), nearest to the study area, range from 10 to 250 μg/L with the volume-weighted mean (VWM) value of 50 μg/L [54]. The Sr concentrations in groundwater vary in a wide range from 377 to 1000 μg/L, with an average value of around 630 μg/L. The much higher Sr2+ concentration in groundwater compared to those observed for local precipitation demonstrates that Sr2+ in groundwater have other potential sources besides precipitation, which may include the dissolution of carbonate and silicate minerals that are major Sr-bearing minerals in the aquifer matrix.
The Sr sources might control the 87Sr/86Sr ratios and can be well constrained by the relationship between the Sr isotopic compositions and concentration ratios for groundwater [55]. Figure 10a shows their relationship in the plot of 87Sr/86Sr ratios versus 1/Sr. The geochemical compositions of three potential end-members are also shown: (a) carbonate fraction of aquifer rocks, represented by an 87Sr/86Sr ratio and Sr concentrations of 0.710148 and 309 μg/g, respectively [56]; (b) silicate fractions of aquifer rocks, characterized by a significantly more radiogenic 87Sr/86Sr ratio (0.716552) and a relatively lower Sr content (284 μg/g) relative to carbonate minerals [57]; and (c) local precipitation, characterized by an 87Sr/86Sr ratio and Sr concentration of 0.710760 and 309 μg/g, respectively [54]. Figure 10a shows almost all the groundwater samples points are located within the triangle consisting of local precipitation, carbonate and silicate fractions, and the plot is obviously closer to the carbonate fraction end-member than the other two end-members, which indicates the carbonate dissolution exerts the significantly more important control on the Sr2+ in groundwater than the other two sources. Figure 10b shows that 87Sr/86Sr ratios tend to gradually decrease with increasing Sr2+ concentrations in groundwater, furthering demonstrating that the major source of the dissolved Sr2+ is not the weathering of silicates with more radiogenic 87Sr/86Sr ratios. This is due to the fact that 87Sr/86Sr ratios will increase with increasing Sr2+ concentrations if silicate weathering is an important process controlling the groundwater hydrochemical evolution [9]. Figure 10b displays that types 3, 4 and 5 have relatively lower 87Sr/86Sr ratios compared to type 1, which may result from the gypsum dissolution as groundwater moves from recharge areas within the aquifer. This is supported by the fact that the mean SO42− concentrations are higher for type 3 (28.6 mg/L), type 4 groundwater (67.1 mg/L) and type 5 (342.8 mg/L) than type 1 groundwater (18.9 mg/L) as shown in Table S3.
As the silicate weathering is not significant within the aquifer, excess Na+ relative to the Cl concentration in type 3 and 4 groundwater is believed to be the primary result of the cation exchange process and the dissolution of Glauber’s salt. Types 3 and 4 are characterized by a predominance of Na+ and HCO3; therefore, cation exchange is the major source of excess Na in these waters. That is, Sr geochemical signatures further confirm the conclusion that a cation exchange significantly contributes to dissolved Na+ in type 3 and type 4 groundwater, deduced from the groundwater chemistry discussed above.

4.3.3. Constraints from Sulfur Isotope

Groundwater sulfate can be derived from various sources, which include the accession of inorganic SO4 via the atmospheric deposition of marine SO4 aerosols, the dissolution of SO4 minerals such as gypsum, the oxidation of reduced sulfide minerals such as pyrite or the mineralization of organic S compounds within the soil zone [58,59]. The sulfur isotopes values for SO434SSO4) of these sources are markedly different. The δ34S values of sedimentary reduced sulfides (e.g., pyrite) typically range from approximately −50 to +10‰, while the evaporite minerals tend to be relatively enriched in δ34S values between around 10 and 20‰ [13]. Therefore, the sources of SO42− in groundwater can be better determined by sulfur isotope values of SO42− [60,61].
In this study, an isotopic (δ34SSO4) analysis was carried out for type 5 groundwater (Na–SO4) with relatively higher SO42− concentrations in order to identify the SO42− source in groundwater. Groundwater primarily originates from the local precipitation in the study area. Hence, the sulfur geochemistry of groundwater is initially controlled by the precipitation. Precipitation has SO42− concentrations between 9.56 and 28.67 mg/L, and the mean δ34S value for SO4 in precipitation is around 3.7‰ in the Ordos Plateau [62]. The SO42− concentrations range from 322.2 to 367.4 mg/L in type 5 groundwater, which is significantly higher compared to local precipitation. This indicates the presence of additional SO42− sources in groundwater except precipitation. The relatively higher δ34S values of these waters (12.4‰ to 14.1‰, mean 13.0‰) than those of precipitation also support this conclusion. These waters have δ34S values within the typical range of evaporite minerals, indicating that SO42− in type 5 is primarily from the dissolution of evaporite such as gypsum and Glauber’s salt in the aquifer. This interpretation is consistent with the hydrochemical data discussed above.

4.3.4. Constraints from Inverse Geochemical Modeling

An inverse geochemical modeling was performed in order to further constrain and quantify the main natural geochemical processes identified above: the dissolution of calcite, dolomite, anhydrite, halite and Glauber’s salt; the precipitation of calcite and dolomite; and a cation exchange with clays. The weathering of silicate minerals is thought to be unimportant for the chemical evolution and, hence, is not considered in the models.
During the inverse modeling, at least two chemical analyses of groundwater at different points of the flow path and a set of plausible phases (mineral and/or gases) that potentially react along this path are required. In the models, the average chemical parameter values from the statistical analyses of each type of groundwater, as presented in Table S3, were used to represent the “initial” and “final” waters along the groundwater flow paths. A total of five models (models 1, 2, 3, 4 and 5) have been run for an interpretation of the chemical compositions of the five types of groundwater present in the study area. The five simulated paths are represented by the transition between the following end-members: (1) pure water to type 1 water; (2) type 1 water to type 2 water; (3) type 2 to type 3 water; (4) type 3 to type 4 water; and (5) type 4 to type 5 water. The phases used were calcite, dolomite, gypsum and halite in all the models based on the lithological and mineralogical information. CO2 gas was only included in Model 1 from the precipitation to type 1 water to simulate carbonate dissolution under the action of CO2 because only this system is supposed to be open for CO2 in the recharge area. CaX, MgX and NaX phases were only included in Models 2, 3 and 4 to simulate the cation exchange; Glauber’s salt was only considered in Model 5 to simulate Glauber’s salt dissolution. The constraints used included Ca2+, Mg2+, Na+, HCO3, SO42− and Cl in all models.
Table S5 presents the selected results of the inverse modeling from all possible models based on the statistical parameters (the sum of residuals and maximum fractional error) calculated by PHREEQC, representing different possible combinations of reactants and products that can account for the changes in the groundwater chemistry. This result demonstrates that the amount of a given mineral phase that must dissolve or the precipitation to produce the observed variations in groundwater chemistry between two end points along the flow path.
Model 1 is the geochemical path from the local precipitation to a typical recharge water of type 1 (Ca·Mg–HCO3); the pure water is assumed to be representative of precipitation in this model. It is summarized as follows:
Precipitation (pure water) + CO2 gas + Calcite + Dolomite + Gypsum + Halite→Ca·Mg–HCO3 type water
Model 2 is the geochemical path from type 1 to type 2 (mixed cations-HCO3) and is summarized as follows:
Ca·Mg–HCO3 type water + Na from cation exchange + Gypsum→mixed cations-HCO3 type water + calcite precipitation (or Ca loss for ion exchange) + dolomite precipitation (or Mg loss for ion exchange)
Model 3 reproduces the geochemical path from type 2 groundwater to type 3 (Na–HCO3) and is summarized as follows:
Mixed cations-HCO3 type water + Na from cation exchange + Gypsum→Na–HCO3 type water + Ca loss for ion exchange + Mg loss for ion exchange
Model 4 reproduces the geochemical path from type 3 water to type 4 (Na–HCO3·SO4) and is summarized as follows:
Na–HCO3 type water + Na from cation exchange + Gypsum + Halite→Na–HCO3·SO4 type water + calcite precipitation (Ca loss for ion exchange) + dolomite precipitation (Mg loss for ion exchange)
Model 5 represents the geochemical path from type 4 water to type 5 (Na–SO4), and it is summarized as follows:
Na–HCO3·SO4 type water + Glauber’s salt + Gypsum + (Halite)→Na–SO4 type water + calcite dolomite precipitation
Summarizing, the geochemical modeling results support the hypothetical geochemical processes derived from the qualitative analyses for the groundwater chemical and isotopic compositions. From the modeling results, it can be confirmed that the dissolution of carbonate and evaporite minerals, as well as a cation exchange within the aquifer primarily control the chemical evolution of groundwater in the study area.

4.4. Groundwater Age

Groundwater age can be used to further constrain groundwater origin and flow systems. Radiocarbon (14C) is the most widely applied dating method by far because of the almost ubiquitous presence of dissolved inorganic carbon (DIC) in groundwater [13].
In spite of the widespread use of 14C, the interpretation of 14C model age still is limited by many uncertainties in determining the 14C content of dissolved carbon in the recharge areas and in accounting for various chemical and physical processes affecting the 14C content along groundwater flow paths in aquifers [63,64,65,66]. Using the measurements of 14C of the DIC in groundwater, 14CDIC, the basic equation for groundwater dating by using 14C is as follows [63]:
t = 5730 ln 2 ln ( 14 C DIC 14 C 0 )
where t is the groundwater age; 5730 is the modern half-life of 14C; 14CDIC is the measured 14C value of the DIC in groundwater sample; and 14C0 is the initial 14C content of DIC after correction for geochemical reactions without radioactivity decay.
The estimation of 14C0 is essential for groundwater 14C dating. Numerous single-sample-based traditional geochemical adjustment models have been proposed to calculate 14C0 so far, such as Vogel, Tamer, Pearson, Fontes and Garnier models and so on [67,68,69,70]. These different models may give significantly distinct values of 14C0. The incorrect use of the models for chemical and isotopic measurements in groundwater will lead to a wide range of estimations of 14C0 and greatly limits the usefulness of radiocarbon as a dating tool for groundwater [63]. Therefore, choosing the most appropriate models is of critical importance for improving the estimation of 14C age in a particular groundwater system, depending on which models most completely account for the geochemical processes affecting the chemical and isotopic compositions of the waters [65].
The Han and Plummer diagram of δ13C versus 14C [63] can aid in determining which adjustment model is the most appropriate and should be applied for the specific groundwater system. Figure 11 shows a plot of 14C versus δ13C (Han and Plummer diagram) for the groundwater samples studied, indicating that all samples can be divided into two groups. One group is composed of type 1 waters; they plot close to the IAEA (International Atomic Energy Agency) line, suggesting they are young waters with 14C age close to zero. Another group includes type 4 and 5 waters, which have δ13CDIC ranging between δ13Ci (approx. −11.05‰) and δ13CE (approx. −0.15‰) and are plotted in the Eichinger area. This indicates that the carbon isotopic compositions of DIC in these waters have been affected by a carbon exchange between DIC in groundwater and solid carbonate under closed system, and the Eichinger model [71], accounting for a carbon exchange between solid carbonate and DIC under a closed system, should be applied to estimate 14C0 for type 4 and 5 waters. In fact, the models of Wigley and Evans also take into consideration the carbon exchange between DIC and carbonate mineral, but they are less precise compared to the Eichinger model because they neglect the contribution of CO2 (aq) to DIC (the models assume that DIC = HCO3). Therefore, Eichinger’s model considering the contribution of CO2 (aq) to DIC is selected to obtain a relatively accurate age in this study.
According to Han and Plummer, the Han and Plummer’s model and IAEA model are also applicable for groundwater samples studied [63]. Han and Plummer’s model is a revised version of the Fontes and Garnier model [72], which considers not only carbon mixing but also carbon exchange occurring both in the saturated and unsaturated zone. This model is a combination of three models: Tamers’, Mook’s and Eichinger’s models. If the carbon exchange between DIC and solid carbonate occurs but without an exchange with soil CO2, the models reduces to Eichinger’s model. The carbon exchange between DIC and solid carbonate predominantly under closed system conditions occurs in the groundwater system studied; therefore, the Han and Plummer model will give a similar result to Eichinger’s model. Unlike the Han and Plummer’s model, assuming either open or closed systems, the IAEA model assumes mixed open and closed system evolutions [73]. That is, this model assumes that the isotopic composition of DIC evolved initially under completely open system condition and then evolved under closed system condition.
The measured and assumed parameters for estimating 14C0 using various adjustment models are summarized in Table S6. The measured parameters include temperature, pH, DIC, contents of dissolved CO2 (aq), HCO3 and CO3. The assumed parameters include 14C and δ13C values of soil gas CO2 and solid carbonate mineral. The 14C of soil CO2 (14Cg) and carbonate mineral (14Cs) are assumed to be 100 pMC and 0 pMC, respectively. The δ13C value for carbonate mineral (δ13Cs) is assumed to be 0‰. The vegetation in the Ordos Plateau is predominantly C3 species, and the soil profiles in Yulin, approximately 125 km away from study area, show that δ13C values range from −22‰ to −20‰ with an average of −21‰ [74,75]. Therefore, the δ13C value for soli gas CO213Cg) is assumed to be −21‰ in the study area. In addition, the calculated fraction factor for determination of 14C0 is also presented in Table S6. The εa/g and εb/g are 13C fractionate factors of dissolved CO2 and HCO3 with respect to gaseous CO2, respectively. The εa/s and εb/s are the 13C fractionate factors of dissolved CO2 and HCO3 with respect to carbonate mineral.
The measured 14CDIC values ranged from 1.4 to 89.6 pMC, with an average of 34.9 pMC (Table S7). Based on above these parameters, the calculated 14C0 and 14C age including the apparent and corrected ages are presented in Table S7. The models of Eichinger and of Han and Plummer yield similar ages, which are relatively lower than those determined by IAEA model. The results derived from these three models are overall comparable, suggesting that 14C method provides the relatively reliable and reasonable ages to represent the mean groundwater residence time. The average values of three models were considered as the final ages. For shallow groundwater in the recharge areas (type 1), the corrected 14C age is approximately 0 year BP (Before Present); for those in discharge areas (type 3 and 4), the corrected 14C age ranges from 3929 to 16,996 years BP with the mean value of 9450 years, which is consistent with the 3H content that is below the detection limit of 1.3 TU (Table S2). For deep groundwater (type 5), the corrected 14C age varies from 23,145 to 27,436 years BP with the average of 25,075 years, which is also in agreement with 3H content below the detection limit (Table S2). Summarizing, the age data further demonstrates that shallow groundwater in the recharge areas were recharged during the Holocene period, while those in the discharge areas and deep groundwater were recharged during the late Pleistocene period. That is, groundwater age agrees well with the δD and δ18O data discussed above. In addition, it is found that deep groundwater is about 2.7 times higher than shallow groundwater in the discharge areas. It is well-known that the groundwater age should increase along the groundwater flow direction. Therefore, such a vertical distribution of groundwater age further corroborates the hypothesis on the development of hierarchically nested groundwater flow systems. That is, groundwater age indicates that shallow and deep groundwater belong to local and regional flow system, respectively.
It is well-known that 14C ages are not identical to calendar age because atmospheric 14C concentration changes temporarily as a result of variation in the production rate, resulting from geomagnetic, solar modulation of the cosmic-ray flux and the carbon cycle [63,76]. In order to convert 14C ages into calendar ages relative to 1950 AD, many radiocarbon age calibration curves have been constructed using a variety of methods. In this study, the calibration curve IntCal13, representing the midlatitude Northern Hemisphere atmospheric reservoirs, was applied to estimate calendar age. Before it is applied, a conversion of conventional 14C ages calculated by using the modern 14C half-life (5730 years) to the Libby half-life (5568 years) is required because the calibration curves are obtained by first using the Libby half-life. The calibration curve IntCal13 and groundwater samples are shown in Figure 12, and the average calendar years of groundwater are also summarized in Table S7. The results of the calendar ages for groundwater show that shallow groundwater in the recharge areas is approximately 0 year BP, while those in discharge area ranges from 4.340 to 19,965 years BP with the mean value of 10,829 years; deep groundwater varies from 26,736 to 28,125 years BP with the average of 28,588 years.

5. Conclusions

Groundwater is the main water resource in the Ordos Plateau, which is located in arid-semiarid regions of China. In this study, environmental tracers including hydrochemistry and isotope, along with the inverse geochemical modeling technique were integrated successfully to obtain a comprehensive understanding on the origin and hydrochemical evolution of shallow and deep groundwater in the HTCH Lake Basin, NE Ordos Plateau. The major conclusions are summarized as follows:
1. The geochemical data support the hypothesis of hierarchically nested groundwater flow system and suggest the development of local and regional flow systems in the study area. Based on the isotopic (δ2H, δ18O, 3H and 14C) compositions, shallow and deep groundwater both originate primarily from local precipitation. Shallow local groundwater in the recharge areas and discharge areas was recharged during the Holocene and late Pleistocene periods, respectively; deep regional groundwater was recharged by precipitation in the late Pleistocene period. Therefore, from the viewpoint of sustainable use, shallow discharge zone groundwater and deep groundwater should not be exploited heavily.
2. The groundwater chemical pattern is closely related to the groundwater flow systems, and the hydrochemical evolution is primarily controlled by geochemical processes. The shallow local groundwater chemistry systematically changes along the flow path: The water type progressively evolves from type 1 (Ca·Mg–HCO3) in the recharge areas through type 2 (mixed cations-HCO3) to type 3 (Na–HCO3) and type 4 (Na–HCO3·SO4) in the discharge areas. This is primarily controlled by carbonates dissolution and cation exchange, which principally occur in the recharge areas and discharge areas of local flow system, respectively. In addition, human activities (e.g., agricultural fertilizer) also affect the chemical compositions of shallow recharge zone groundwater, particularly Cl and NO3 contents. Deep regional groundwater belongs to type 5 (Na–SO4), which is dominated by evaporite dissolution (e.g., gypsum and Glauber’s salt). Although the aquifer studied consists mainly of silicate-dominated Cretaceous sandstone, Sr isotope provides the compelling evidence that silicate weathering is not important in the study area.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/11/4/790/s1, Table S1: The well location, field parameters and major ion chemistry of the groundwater and lake water in the study area; Table S2: The minor ion chemistry, isotope chemistry, saturation indices and calculated 14C age of the groundwater in the study area; Table S3: The mean chemical and isotopic compositions along with the standard deviations (s.d.) in the 5 types of groundwater in the study area; Table S4: The mineralogical compositions of the aquifer core material in the study area; Table S5: A summary of the mass transfer for geochemical inverse modeling (The values are in mol/L): The phases and thermodynamic database are from PHREEQCI (version 2); Table S6: The measured and assumed values for the calculation of 14C0 for the groundwater in the study area; Table S7: The calculated 14C0 values and the corrected 14C ages using various models for the groundwater in the study area.

Author Contributions

M.L. performed the sampling and wrote the paper. Z.P. designed the research, and guided data interpretation and corrected the manuscript. L.Y., J.Z., S.Y., X.W. and T.G. helped interpret data. T.H. and Z.L. contributed to the mineralological analyses of the rock samples and to the related discussion.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC Grants 41727901 and 41672254), Hydrogeological Investigation at a 1:50,000 scale in the lake-concentrated areas of the Northern Ordos Basin (Grant DD20160293) and by the International Atomic Energy Agency through a regional cooperation project (IAEA/RCA RAS7030), all with Z.P. as CSI except 41672254 with T.H.

Acknowledgments

The authors thank Wuhui Jia, Tiangang Liu and Wenhao Xu for their assistance in the field work. We also thank the two anonymous reviewers for their constructive comments and suggestions that have helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders have no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology; Wiley: New York, NY, USA, 1998. [Google Scholar]
  2. Appelo, C.A.J.; Postma, D. Geochemistry, Groundwater and Pollution, 2nd ed.; Balkema: Leiden, The Netherlands, 2005. [Google Scholar]
  3. Edmunds, W.; Bath, A.; Miles, D. Hydrochemical evolution of the East Midlands Triassic sandstone aquifer, England. Geochim. Cosmochim. Acta 1982, 46, 2069–2081. [Google Scholar] [CrossRef]
  4. Herczeg, A.L.; Torgersen, T.; Chivas, A.; Habermehl, M. Geochemistry of ground waters from the Great Artesian Basin, Australia. J. Hydrol. 1991, 126, 225–245. [Google Scholar] [CrossRef]
  5. Fisher, R.S.; Mullican, I.; William, F. Hydrochemical Evolution of Sodium-Sulfate and Sodium-Chloride Groundwater Beneath the Northern Chihuahuan Desert, Trans-Pecos, Texas, USA. Hydrogeol. J. 1997, 5, 4–16. [Google Scholar] [CrossRef]
  6. Petrides, B.; Cartwright, I.; Weaver, T.R. The evolution of groundwater in the Tyrrell catchment, south-central Murray Basin, Victoria, Australia. Hydrogeol. J. 2006, 14, 1522–1543. [Google Scholar] [CrossRef]
  7. Zhang, G.; Deng, W.; Yang, Y.S.; Salama, R.B. Evolution study of a regional groundwater system using hydrochemistry and stable isotopes in Songnen Plain, northeast China. Hydrol. Process. 2007, 21, 1055–1065. [Google Scholar] [CrossRef]
  8. Zhu, G.F.; Li, Z.Z.; Su, Y.H.; Ma, J.Z.; Zhang, Y.Y. Hydrogeochemical and isotope evidence of groundwater evolution and recharge in Minqin Basin, Northwest China. J. Hydrol. 2007, 333, 239–251. [Google Scholar] [CrossRef]
  9. Currell, M.J.; Cartwright, I. Major-ion chemistry, δ13C and 87Sr/86Sr as indicators of hydrochemical evolution and sources of salinity in groundwater in the Yuncheng Basin, China. Hydrogeol. J. 2011, 19, 835. [Google Scholar] [CrossRef]
  10. Edmunds, W. Geochemistry’s vital contribution to solving water resource problems. Appl. Geochem. 2009, 24, 1058–1073. [Google Scholar] [CrossRef]
  11. Glynn, P.D.; Plummer, L.N. Geochemistry and the understanding of ground-water systems. Hydrogeol. J. 2005, 13, 263–287. [Google Scholar] [CrossRef]
  12. Herczeg, A.L.; Leaney, F. Environmental tracers in arid-zone hydrology. Hydrogeol. J. 2011, 19, 17–29. [Google Scholar] [CrossRef]
  13. Clark, I.D.; Fritz, P. Environmental Isotopes in Hydrogeology; CRC Press LLC: Boca Raton, FL, USA, 1997. [Google Scholar]
  14. Yin, L.; Hou, G.; Dou, Y.; Tao, Z.; Li, Y. Hydrogeochemical and isotopic study of groundwater in the Habor Lake Basin of the Ordos Plateau, NW China. Environ. Earth Sci. 2009, 64, 1575–1584. [Google Scholar] [CrossRef]
  15. Hou, G.C.; Zhang, M.S.; Liu, F. Groundwater Investigation and Research of Ordos Basin; Geological Publishing House: Beijing, China, 2008. [Google Scholar]
  16. Liu, F.; Song, X.; Yang, L.; Zhang, Y.; Han, D.; Ma, Y.; Bu, H. Identifying the origin and geochemical evolution of groundwater using hydrochemistry and stable isotopes in the Subei Lake basin, Ordos energy base, Northwestern China. Hydrol. Earth Syst. Sci. 2015, 19, 551–565. [Google Scholar] [CrossRef]
  17. Hou, G.C.; Yin, L.H.; Xu, D.D. Hydrogeology of the Ordos Basin, China. J. Groundw. Sci. Eng. 2017, 5, 104–115. [Google Scholar]
  18. Hou, G.C.; Liang, Y.P.; Su, X.S.; Zhao, Z.H.; Tao, Z.P.; Yin, L.H.; Yang, Y.C.; Wang, X.Y. Groundwater systems and resources in the Ordos Basin, China. Acta Geol. Sin. Engl. Ed. 2008, 82, 1061–1069. [Google Scholar]
  19. Yin, L.H.; Hou, G.C.; Tao, Z.P.; Li, Y. Origin and recharge estimates of groundwater in the ordos plateau, People’s Republic of China. Environ. Earth Sci. 2009, 60, 1731–1738. [Google Scholar] [CrossRef]
  20. Zhang, H.; Jiang, X.W.; Wan, L.; Ke, S.; Liu, S.-A.; Han, G.; Guo, H.; Dong, A. Fractionation of Mg isotopes by clay formation and calcite precipitation in groundwater with long residence times in a sandstone aquifer, Ordos Basin, China. Geochim. Cosmochim. Acta 2018, 237, 261–274. [Google Scholar] [CrossRef]
  21. Wang, H.; Jiang, X.W.; Wan, L.; Han, G.L.; Guo, H.M. Hydrogeochemical characterization of groundwater flow systems in the discharge area of a river basin. J. Hydrol. 2015, 527, 433–441. [Google Scholar] [CrossRef]
  22. Jiang, X.W.; Wan, L.; Wang, X.S.; Wang, D.; Wang, H.; Wang, J.Z.; Zhang, H.; Zhang, Z.Y.; Zhao, K.Y. A multi-method study of regional groundwater circulation in the Ordos Plateau, NW China. Hydrogeol. J. 2018, 26, 1657–1668. [Google Scholar] [CrossRef]
  23. Yin, L.H.; Hou, G.C.; Su, X.S.; Wang, D.; Dong, J.; Hao, Y.H.; Wang, X.Y. Isotopes (δD and δ18O) in precipitation, groundwater and surface water in the Ordos Plateau, China: Implications with respect to groundwater recharge and circulation. Hydrogeol. J. 2011, 19, 429–443. [Google Scholar] [CrossRef]
  24. Pan, G.; Li, X.; Zhang, J.; Liu, Y.; Liang, H. Groundwater-flow-system characterization with hydrogeochemistry: A case in the lakes discharge area of the Ordos Plateau, China. Hydrogeol. J. 2018, 27, 669–683. [Google Scholar] [CrossRef]
  25. Zhang, J.; Wang, W.; Wang, X.; Yin, L.; Zhu, L.; Sun, F.; Dong, J.; Xie, Y.; Robinson, N.I.; Love, A.J. Seasonal variation in the precipitation recharge coefficient for the Ordos Plateau, Northwest China. Hydrogeol. J. 2019, 27, 801–813. [Google Scholar] [CrossRef]
  26. Schwartz, F.W.; Domenico, P.A. Simulation of Hydrochemical Patterns in Regional Groundwater Flow. Water Resour. Res. 1973, 9, 707–720. [Google Scholar] [CrossRef]
  27. Plummer, L.N.; Parkhurst, D.L.; Thorstenson, D.C. Development of reaction models for ground-water systems. Geochim. Cosmochim. Acta 1983, 47, 665–685. [Google Scholar] [CrossRef]
  28. Plummer, L.N.; Busby, J.F.; Lee, R.W.; Hanshaw, B.B. Geochemical modeling of the Madison aquifer in parts of Montana, Wyoming, and South Dakota. Water Resour. Res. 1990, 26, 1981–2014. [Google Scholar] [CrossRef]
  29. Plummer, L.N.; Prestemon, E.C.; Parkhurst, D.L. An Interactive Code (NETPATH) for Modeling Net Geochemical Reactions along a Flow Path; US Geological Survey: Reston, VA, USA, 1991; Volume 94.
  30. Parkhurst, D.L.; Appelo, C. User’s guide to PHREEQC (Version 2): A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Water-Resour. Investig. Rep. 1999, 99, 312. [Google Scholar]
  31. Piper, A.M. A graphic procedure in the geochemical interpretation of water-analyses. Eos Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar] [CrossRef]
  32. Liu, F.; Song, X.; Yang, L.; Han, D.; Zhang, Y.; Ma, Y.; Bu, H. 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] [PubMed] [Green Version]
  33. Back, W. Origin of hydrochemical facies of ground water in the Atlantic Coastal Plain. In Proceedings of the 21st International Geological Congress, Copenhagen, Denmark, 1960; pp. 87–95. Available online: http://www.geology.wisc.edu/~andy/g929/Back_new.pdf (accessed on 12 April 2019).
  34. Chebotarev, I. Metamorphism of natural waters in the crust of weathering-1. Geochim. Cosmochim. Acta 1955, 8, 22–51. [Google Scholar] [CrossRef]
  35. Stuyfzand, P.J. Patterns in groundwater chemistry resulting from groundwater flow. Hydrogeol. J. 1999, 7, 15–27. [Google Scholar] [CrossRef]
  36. Toth, J. A theoretical analysis of groundwater flow in small drainage basins. J. Geophys. Res. 1963, 68, 4795–4812. [Google Scholar] [CrossRef]
  37. Freeze, R.A. Groundwater; Prentice Hall: Englewood Cliffs, NJ, USA, 1979. [Google Scholar]
  38. Tóth, J. Groundwater as a geologic agent: An overview of the causes, processes, and manifestations. Hydrogeol. J. 1999, 7, 1–14. [Google Scholar] [CrossRef]
  39. Craig, H. Isotopic variations in meteoric waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef] [PubMed]
  40. Yang, Y.C.; Shen, Z.L.; Weng, D.G.; Hou, G.C.; Zhao, Z.H.; Wang, D.; Pang, Z.H. Oxygen and hydrogen isotopes of waters in the Ordos Basin, China: Implications for recharge of groundwater in the north of Cretaceous Groundwater Basin. Acta Geol. Sin. Engl. Ed. 2009, 83, 103–113. [Google Scholar] [CrossRef]
  41. Gibbs, R.J. Mechanisms controlling world water chemistry. Science 1970, 170, 1088–1092. [Google Scholar] [CrossRef] [PubMed]
  42. Edmunds, W.M.; Carrillo-Rivera, J.J.; Cardona, A. Geochemical evolution of groundwater beneath Mexico City. J. Hydrol. 2002, 258, 1–24. [Google Scholar] [CrossRef]
  43. Edmunds, W.M.; Ma, J.; Aeschbach-Hertig, W.; Kipfer, R.; Darbyshire, D.P.F. Groundwater recharge history and hydrogeochemical evolution in the Minqin Basin, North West China. Appl. Geochem. 2006, 21, 2148–2170. [Google Scholar] [CrossRef]
  44. Cartwright, I.; Weaver, T.R. Hydrogeochemistry of the Goulburn Valley region of the Murray Basin, Australia: Implications for flow paths and resource vulnerability. Hydrogeol. J. 2005, 13, 752–770. [Google Scholar] [CrossRef]
  45. Foster, M.D. The origin of high sodium bicarbonate waters in the Atlantic and Gulf coastal plains. Geochim. Cosmochim.Acta 1950, 1, 33–48. [Google Scholar] [CrossRef]
  46. Varsányi, I.; Kovács, L. Chemical evolution of groundwater in the River Danube deposits in the southern part of the Pannonian Basin (Hungary). Appl. Geochem. 1997, 12, 625–636. [Google Scholar] [CrossRef]
  47. Somaratne, N.; Mustafa, S.; Lawson, J. Use of hydrochemistry, stable isotope, radiocarbon, 222Rn and terrigenic 4He to study the geochemical processes and the mode of vertical leakage to the Gambier Basin tertiary confined sand aquifer, South Australia. Water 2016, 8, 180. [Google Scholar] [CrossRef]
  48. Howard, K.; Beck, P. Hydrochemical interpretation of groundwater flow systems in Quaternary sediments of southern Ontario. Can. J. Earth Sci. 1986, 23, 938–947. [Google Scholar] [CrossRef]
  49. Schoeller, H. Qualitative evaluation of groundwater resources. In Methods and Techniques of Groundwater Investigations and Development; UNESCO: Paris, France, 1965; Volume 5483. [Google Scholar]
  50. Shand, P.; Darbyshire, D.P.F.; Love, A.J.; Edmunds, W.M. Sr isotopes in natural waters: Applications to source characterisation and water–rock interaction in contrasting landscapes. Appl. Geochem. 2009, 24, 574–586. [Google Scholar] [CrossRef]
  51. Capo, R.C.; Stewart, B.W.; Chadwick, O.A. Strontium isotopes as tracers of ecosystem processes: Theory and methods. Geoderma 1998, 82, 197–225. [Google Scholar] [CrossRef]
  52. Paces, J.B.; Wurster, F.C. Natural uranium and strontium isotope tracers of water sources and surface water–groundwater interactions in arid wetlands–Pahranagat Valley, Nevada, USA. J. Hydrol. 2014, 517, 213–225. [Google Scholar] [CrossRef]
  53. McNutt, R.H. Strontium isotopes. In Environmental Tracers in Subsurface Hydrology; Cook, P.G., Herczeg, A.L., Eds.; Kluwer Academic Publishers: Boston, MA, USA, 2000; pp. 233–260. [Google Scholar]
  54. Rao, W.B.; Han, G.L.; Tan, H.B.; Jiang, S. Chemical and Sr isotopic compositions of rainwater on the Ordos Desert Plateau, Northwest China. Environ. Earth Sci. 2015, 74, 5759–5771. [Google Scholar] [CrossRef]
  55. Dogramaci, S.; Skrzypek, G. Unravelling sources of solutes in groundwater of an ancient landscape in NW Australia using stable Sr, H and O isotopes. Chem. Geol. 2015, 393–394, 67–78. [Google Scholar] [CrossRef]
  56. Su, X.S.; Wu, C.Y.; Dong, W.H.; Hou, G.C. Strontium isotope evolution mechanism of the Cretaceous groundwater in Ordos Desert Plateau. J. Chengdu Univ. Technol. 2011, 38, 348–358. [Google Scholar]
  57. Rao, W.B.; Chen, J.; Tan, H.B.; Jiang, S.Y.; Su, J. Sr–Nd isotopic and REE geochemical constraints on the provenance of fine-grained sands in the Ordos deserts, north-central China. Geomorphology 2011, 132, 123–138. [Google Scholar] [CrossRef]
  58. Mayer, B. Potential and limitations of using sulphur isotope abundance ratios as an indicator for natural and anthropogenic environmental change. In Isotope Techniques in the Study of Environmental Change; International Atomic Energy Agency: Vienna, Austria, 1998; pp. 423–435. [Google Scholar]
  59. Marques, J.M.; Graça, H.; Eggenkamp, H.G.M.; Neves, O.; Carreira, P.M.; Matias, M.J.; Mayer, B.; Nunes, D.; Trancoso, V.N. Isotopic and hydrochemical data as indicators of recharge areas, flow paths and water–rock interaction in the Caldas da Rainha–Quinta das Janelas thermomineral carbonate rock aquifer (Central Portugal). J. Hydrol. 2013, 476, 302–313. [Google Scholar] [CrossRef]
  60. Dogramaci, S.S.; Herczeg, A.L.; Schiff, S.L.; Bone, Y. Controls on δ34S and δ18O of dissolved sulfate in aquifers of the Murray Basin, Australia and their use as indicators of flow processes. Appl. Geochem. 2001, 16, 475–488. [Google Scholar] [CrossRef]
  61. McIntosh, J.C.; Walter, L.M. Paleowaters in Silurian-Devonian carbonate aquifers: Geochemical evolution of groundwater in the Great Lakes region since the Late Pleistocene. Geochim. Cosmochim. Acta 2006, 70, 2454–2479. [Google Scholar] [CrossRef]
  62. Yang, Y.C.; Shen, Z.L.; Wen, D.G.; Hou, G.C.; She, H.Q.; Zhao, Z.H.; Wang, D.; Li, J.W. Distribution of δ34S and δ18O in SO42− in Groundwater from the Ordos Cretaceous Groundwater Basin and Geological Implications. Acta Geol. Sin. Engl. Ed. 2010, 84, 432–440. [Google Scholar] [CrossRef]
  63. Han, L.; Plummer, L. A review of single-sample-based models and other approaches for radiocarbon dating of dissolved inorganic carbon in groundwater. Earth-Sci. Rev. 2016, 152, 119–142. [Google Scholar] [CrossRef]
  64. Geyh, M.A. An overview of 14C analysis in the study of groundwater. Radiocarbon 2000, 42, 99–114. [Google Scholar] [CrossRef]
  65. Plummer, L.N.; Glynn, P.D. Radiocarbon dating in groundwater systems. In Isotope Methods for Dating Old Groundwater; Suckow, A., Aggarwal, P.K., Araguas-Araguas, L., Eds.; International Atomic Energy Agency: Vienna, Austria, 2013. [Google Scholar]
  66. Geyh, M.A. Dating of Old Groundwater-History, Potential, Limits and Future. In Isotopes in the Water Cycle: Past, Present and Future of a Developing Science; Aggarwal, P.K., Gat, J.R., Froehlich, K.F.O., Eds.; Springer: Dordrecht, The Netherlands, 2005; pp. 221–241. [Google Scholar]
  67. Vogel, J. Carbon-14 dating of groundwater. Isot. Hydrol. 1970, 1970, 225–240. [Google Scholar]
  68. Tamers, M. Validity of radiocarbon dates on ground water. Surv. Geophys. 1975, 2, 217–239. [Google Scholar] [CrossRef]
  69. Pearson, F.; Hanshaw, B. Sources of dissolved carbonate species in groundwater and their effects on carbon-14 dating. Isot. Hydrol. 1970, 1970, 271–285. [Google Scholar]
  70. Fontes, J.C.; Garnier, J.M. Determination of the initial 14C activity of the total dissolved carbon: A review of the existing models and a new approach. Water Resour. Res. 1979, 15, 399–413. [Google Scholar] [CrossRef]
  71. Eichinger, L. A contribution to the interpretation of 14C groundwater ages considering the example of a partially confined sandstone aquifer. Radiocarbon 1983, 25, 347–356. [Google Scholar] [CrossRef]
  72. Han, L.-F.; Plummer, L.N. Revision of Fontes & Garnier’s model for the initial 14C content of dissolved inorganic carbon used in groundwater dating. Chem. Geol. 2013, 351, 105–114. [Google Scholar]
  73. Salem, O.; Visser, J.; Dray, M.; Gonfiantini, R. Groundwater flow patterns in the western Libyan Arab Jamahiriya evaluated from isotopic data. In Groundwater Flow Patterns in the Western Libyan Arab Jamahiriya Evaluated from Isotopic Data; International Atomic Energy Agency: Vienna, Austria, 1980. [Google Scholar]
  74. Yang, S.; Ding, Z.; Li, Y.; Wang, X.; Jiang, W.; Huang, X. Warming-induced northwestward migration of the East Asian monsoon rain belt from the Last Glacial Maximum to the mid-Holocene. Proc. Natl. Acad. Sci. USA 2015, 112, 13178–13183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Huang, T.; Pang, Z.; Li, J.; Xiang, Y.; Zhao, Z. Mapping groundwater renewability using age data in the Baiyang alluvial fan, NW China. Hydrogeol. J. 2017, 25, 743–755. [Google Scholar] [CrossRef]
  76. Reimer, P.J.; Bard, E.; Bayliss, A.; Beck, J.W.; Blackwell, P.G.; Ramsey, C.B.; Buck, C.E.; Cheng, H.; Edwards, R.L.; Friedrich, M. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 2013, 55, 1869–1887. [Google Scholar] [CrossRef]
Figure 1. (a) The location map of the Ordos Basin in China; (b) a map showing three different groundwater systems in the Ordos Plateau and the location of study area; (c) a map showing the topography and groundwater sampling locations and the location of the hydrogeological cross section A–A″; and (d) the dotted lines represent the contours of water table elevation (m above sea level) in the study area, and the data in the lines is the specific value of the water table.
Figure 1. (a) The location map of the Ordos Basin in China; (b) a map showing three different groundwater systems in the Ordos Plateau and the location of study area; (c) a map showing the topography and groundwater sampling locations and the location of the hydrogeological cross section A–A″; and (d) the dotted lines represent the contours of water table elevation (m above sea level) in the study area, and the data in the lines is the specific value of the water table.
Water 11 00790 g001
Figure 2. (a) A schematic hydrogeological cross section of A–A″ in Figure 1c showing the hypothesis of hierarchically nested flow systems (local, intermediate and regional flow systems) controlled by wave topography: The blue, green and red solid lines with arrows indicate the local, intermediate and regional groundwater flow systems, respectively, and the pink dotted lines denote the boundaries of the different flow systems. (b) A schematic hydrogeological cross section of A′–A″ showing the sampling sites of all 38 groundwater samples and the boundaries of the different five types of groundwater: Types 1, 2, 3, 4 and 5 represent Ca(Mg)–HCO3, mixed cations-HCO3, Na–HCO3, Na–HCO3(SO4) and Na–SO4, respectively.
Figure 2. (a) A schematic hydrogeological cross section of A–A″ in Figure 1c showing the hypothesis of hierarchically nested flow systems (local, intermediate and regional flow systems) controlled by wave topography: The blue, green and red solid lines with arrows indicate the local, intermediate and regional groundwater flow systems, respectively, and the pink dotted lines denote the boundaries of the different flow systems. (b) A schematic hydrogeological cross section of A′–A″ showing the sampling sites of all 38 groundwater samples and the boundaries of the different five types of groundwater: Types 1, 2, 3, 4 and 5 represent Ca(Mg)–HCO3, mixed cations-HCO3, Na–HCO3, Na–HCO3(SO4) and Na–SO4, respectively.
Water 11 00790 g002
Figure 3. A Piper Plot illustrating the major ion chemistry of groundwater and lake water: The groundwater samples are divided into five geochemically distinct types: type 1 (Ca(Mg)–HCO3), type 2 (mixed cations-HCO3), type 3 (Na–HCO3), type 4 (Na–HCO3(SO4)) and type 5 (Na–SO4). The black solid lines in the diagrams are the groundwater chemical evolutionary trend. SG and DG in the legend represent shallow and deep groundwater, respectively.
Figure 3. A Piper Plot illustrating the major ion chemistry of groundwater and lake water: The groundwater samples are divided into five geochemically distinct types: type 1 (Ca(Mg)–HCO3), type 2 (mixed cations-HCO3), type 3 (Na–HCO3), type 4 (Na–HCO3(SO4)) and type 5 (Na–SO4). The black solid lines in the diagrams are the groundwater chemical evolutionary trend. SG and DG in the legend represent shallow and deep groundwater, respectively.
Water 11 00790 g003
Figure 4. The Box and Whisker plots of (a) total dissolved solids (TDS) content; (b) Na content; (c) Ca content; (d) Mg content; (e) Cl content; (f) SO4 content; (g) HCO3 content; and (h) NO3 content for the five types of groundwater.
Figure 4. The Box and Whisker plots of (a) total dissolved solids (TDS) content; (b) Na content; (c) Ca content; (d) Mg content; (e) Cl content; (f) SO4 content; (g) HCO3 content; and (h) NO3 content for the five types of groundwater.
Water 11 00790 g004
Figure 5. The plots of the chemical compositions of groundwater versus sampling depth: (a) TDS content; (b) Na+ content; (c) Ca2+ content; (d) Mg2+ content; (e) HCO3 content; (f) SO42 content; and (g) Cl content of groundwater.
Figure 5. The plots of the chemical compositions of groundwater versus sampling depth: (a) TDS content; (b) Na+ content; (c) Ca2+ content; (d) Mg2+ content; (e) HCO3 content; (f) SO42 content; and (g) Cl content of groundwater.
Water 11 00790 g005
Figure 6. The plot of the δ2H vs. δ18O values of groundwater investigated in this study as compared to the local meteoric water line (LMWL) and global meteoric water line (GMWL): The green pentagram represents the mean isotopic composition of modern rainfall in the Ordos Plateau (The data is from Reference [23]).
Figure 6. The plot of the δ2H vs. δ18O values of groundwater investigated in this study as compared to the local meteoric water line (LMWL) and global meteoric water line (GMWL): The green pentagram represents the mean isotopic composition of modern rainfall in the Ordos Plateau (The data is from Reference [23]).
Water 11 00790 g006
Figure 7. The plot of (a) the δ18O value and (b) the δD value in groundwater versus the depth of sampling: The dashed lines in Figure 7a,b denote the mean δ18O (−8.1‰) and δ2H (−58.5‰) values of modern rainfall in the Ordos Plateau, respectively.
Figure 7. The plot of (a) the δ18O value and (b) the δD value in groundwater versus the depth of sampling: The dashed lines in Figure 7a,b denote the mean δ18O (−8.1‰) and δ2H (−58.5‰) values of modern rainfall in the Ordos Plateau, respectively.
Water 11 00790 g007
Figure 8. The Gibbs diagram for groundwater in the study area: (a) logTDS vs. Na+/(Na+ + Ca2+) and (b) logTDS vs. Cl/(Cl + HCO3).
Figure 8. The Gibbs diagram for groundwater in the study area: (a) logTDS vs. Na+/(Na+ + Ca2+) and (b) logTDS vs. Cl/(Cl + HCO3).
Water 11 00790 g008
Figure 9. The plots of the correlations among different ions (and their ratios and indices) for the five water types: relationships between (a) Cl and Na; (b) (Ca + Mg) and (HCO3 + SO4); (c) Na and HCO3; (d) Na/Ca ratios (eq·eq−1) and TDS; (e) (Na − Cl) and (Ca + Mg); (f) CAI-1 and CAI-2 (CAI is the abbreviation for the Chloro-Alkaline indices); (g) Ca and SO4; and (h) Na and (SO4 + Cl).
Figure 9. The plots of the correlations among different ions (and their ratios and indices) for the five water types: relationships between (a) Cl and Na; (b) (Ca + Mg) and (HCO3 + SO4); (c) Na and HCO3; (d) Na/Ca ratios (eq·eq−1) and TDS; (e) (Na − Cl) and (Ca + Mg); (f) CAI-1 and CAI-2 (CAI is the abbreviation for the Chloro-Alkaline indices); (g) Ca and SO4; and (h) Na and (SO4 + Cl).
Water 11 00790 g009
Figure 10. (a) The plot of 87Sr/86Sr ratio versus 1/Sr of the groundwater: Three end-members are also shown. (b) The plot of 87Sr/86Sr ratio versus Sr concentration of groundwater.
Figure 10. (a) The plot of 87Sr/86Sr ratio versus 1/Sr of the groundwater: Three end-members are also shown. (b) The plot of 87Sr/86Sr ratio versus Sr concentration of groundwater.
Water 11 00790 g010
Figure 11. A plot of 14C versus δ13C (Han and Plummer diagram) for the groundwater samples in the study area: The δ13Ci are the calculated δ13CDIC according to Tamer model, which is about −11.1‰ at a temperature of 10 °C based on the assumed δ13C values listed in Table S5; δ13CE is the δ13C value of DIC in equilibrium with solid carbonate, which is about −0.15‰ at a temperature of 10 °C.
Figure 11. A plot of 14C versus δ13C (Han and Plummer diagram) for the groundwater samples in the study area: The δ13Ci are the calculated δ13CDIC according to Tamer model, which is about −11.1‰ at a temperature of 10 °C based on the assumed δ13C values listed in Table S5; δ13CE is the δ13C value of DIC in equilibrium with solid carbonate, which is about −0.15‰ at a temperature of 10 °C.
Water 11 00790 g011
Figure 12. The curve for the calibration of groundwater 14C ages: This curve is for the northern hemisphere IntCal13. The diagonal line shows where the curve would lie if the conventional 14C and calendar ages are the same.
Figure 12. The curve for the calibration of groundwater 14C ages: This curve is for the northern hemisphere IntCal13. The diagonal line shows where the curve would lie if the conventional 14C and calendar ages are the same.
Water 11 00790 g012

Share and Cite

MDPI and ACS Style

Lyu, M.; Pang, Z.; Yin, L.; Zhang, J.; Huang, T.; Yang, S.; Li, Z.; Wang, X.; Gulbostan, T. The Control of Groundwater Flow Systems and Geochemical Processes on Groundwater Chemistry: A Case Study in Wushenzhao Basin, NW China. Water 2019, 11, 790. https://doi.org/10.3390/w11040790

AMA Style

Lyu M, Pang Z, Yin L, Zhang J, Huang T, Yang S, Li Z, Wang X, Gulbostan T. The Control of Groundwater Flow Systems and Geochemical Processes on Groundwater Chemistry: A Case Study in Wushenzhao Basin, NW China. Water. 2019; 11(4):790. https://doi.org/10.3390/w11040790

Chicago/Turabian Style

Lyu, Min, Zhonghe Pang, Lihe Yin, Jun Zhang, Tianming Huang, Shuo Yang, Zhenbin Li, Xiaoyong Wang, and Tursun Gulbostan. 2019. "The Control of Groundwater Flow Systems and Geochemical Processes on Groundwater Chemistry: A Case Study in Wushenzhao Basin, NW China" Water 11, no. 4: 790. https://doi.org/10.3390/w11040790

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