Next Article in Journal
Is Urban Economic Output Decoupling from Water Use in Developing Countries?—Empirical Analysis of Beijing and Shanghai, China
Previous Article in Journal
Influence of Convex and Concave Curvatures in a Coastal Dike Line on Wave Run-up
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrochemical Evolution of Groundwater in a Typical Semi-Arid Groundwater Storage Basin Using a Zoning Model

1
Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130026, China
2
College of New Energy and Environment, Jilin University, Changchun 130021, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(7), 1334; https://doi.org/10.3390/w11071334
Submission received: 22 May 2019 / Revised: 17 June 2019 / Accepted: 24 June 2019 / Published: 28 June 2019
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Groundwater guarantees water resources and ecological environment security in semi-arid areas. Studying the chemical evolution of groundwater in semi-arid areas is of great significance to regional ecological environment protection and water resources management. The water storage basin is not only a space for groundwater storage and movement but also a place for water–rock–gas interaction and elemental migration, dispersion, and enrichment. Due to its unique climate and geological environment, the semi-arid water storage basins have different hydrochemical environments, forming a zonal hydrogeochemical character. In this study, a typical semi-arid water storage basin (west of Jilin Province) is taken as an example, through the cross section of the recharge–runoff–excretion zone. A three-level hydrogeochemical zoning model is constructed to reveal the hydrogeochemical evolution of the area. The model is divided into three layers from bottom to top. The first layer shows the geological and hydrogeological conditions, including the topography, lithology, geological time, and hydrodynamic characteristics of the study area. The second layer represents the hydrogeochemical processes, divided into the recharge zone, runoff zone, and discharge zone in the horizontal direction according to the hydrodynamic environment and hydrochemistry type. The hydrogeochemical action gradually changes from lixiviation to cation exchange, evaporation, and concentration, as the landform plays a key role in hydrochemistry formation in the discharge area. The third layer gives the characteristics of the groundwater chemical components, including chemistry type, total dissolved solids, main anion and cation, and characteristic element F. Qualitative and quantitative characterizations of hydrochemistry evolution by reverse simulation and hydrodynamic, hydrochemical and thermodynamic indicators are given.

1. Introduction

Groundwater is an essential part of the global water cycle, serving as the link between the lithosphere, the biosphere, and the atmosphere. The groundwater circulation involves the movement of solvents (water molecules) and solutes [1]. The formation of groundwater chemical characteristics is the result of the long-term interaction between groundwater and the surrounding environment, with a complex exchange of substances and energy, and the continuous migration, dispersion, and enrichment of elements. At present, the research on hydrogeochemical evolution mainly focuses on genesis and source analysis, the migration and enrichment of the feature elements (e.g., Fe, F, As), quantitative chemical simulation, groundwater pollution remediation, and the intrinsic relationship with groundwater circulation [2,3,4,5,6]. The study of regional hydrogeochemical evolution can help us reshape the regional hydrogeological history, clarify the origin and formation of groundwater, reveal the interaction mechanism between groundwater and the environment, and provide a scientific basis for environmental protection and water resource management [7,8].
It is well known that changes in and the distribution of climate, vegetation, and animals on the earth exhibit zonation to a certain extent, and the chemical composition of groundwater is no exception. The essence of hydrochemistry zoning is that the differentiation of the hydrogeochemical environment and the zoned scale is determined by the scale of the chemical environment. The scale can be a large continental structural scale. For example, in China, from southeast to northwest, groundwater of bicarbonate-type water, low in total dissolved solids (TDS), gradually turns into a complex SO42− type water, finally to become high-TDS Cl type water [9]. The scale can also be micro-geomorphic, such as groundwater in hills with high altitude and deep groundwater buried depth, which is generally HCO3-type water with low TDS, while groundwater in a low-lying fulje with a shallow groundwater buried depth is typically HCO3·Cl or Cl-type water with high TDS [10]. Usually, the large regional zoning depends on the climatic conditions, while the small-scale zoning depends on the topographical geological conditions.
The diversity of research subjects and zoning criteria determines the diversity of the zoning model. The authors of [11] divided the groundwater in Datong Basin into five hydrochemical zones: leaching zone, converging zone, enriching zone, oxidizing zone, and reducing zone according to the hydrogeochemical environment. The authors of [12] divided the piedmont plain of North China into the piedmont alluvial area, central alluvial-lake plain, and coastal plain, based on the topography. The authors of [13] divided the coastal area of Turkey into a freshwater holding area and a seawater intrusion area, according to the distribution of groundwater salinity. The authors of [4] divided the groundwater genesis of the Yinchuan Plain into a leaching zone, evaporation-concentration zone, and an evaporation-mixing zone, according to the degree of influence of natural and human irrigation activity on groundwater chemistry. The authors of [14] divided the Polish glacial-buried valley aquifer into a precipitation recharge area, lateral runoff recharge area, and leakage recharge area, based on different groundwater sources. However, hydrochemistry zonation can be inconspicuous because of local hydrogeochemical environmental anomalies [15,16].
These studies indicate that the chemical composition and formation of groundwater cannot be analyzed in isolation from the perspective of just chemistry. The evolutionary law must be revealed from the perspective of the interaction between water and the environment. However, most studies have only qualitatively analyzed groundwater chemical zoning. The authors of [1], taking the hydrogeological entity “hydrological geological storage unit” as the research object, used the combination of chemical thermodynamics and chemical kinetics to illustrate new hydrogeochemical differentiation theory in the water–gas–rock equilibrium system, and a new expression of Darcy’s law established by hydrochemistry indicators and a new method for calculating hydrogeological parameters and groundwater age based on hydrochemistry data was proposed. According to the theory, for typical hydrogeological water storage structural units, the difference in the scale, direction, and intensity of tectonic movement determines the difference between topography and stratum lithology, causing the groundwater dynamical characteristics, circulation conditions, and intensity of the water–rock–gas interaction to change from the edge to the center of the structural unit, which will affect the direction and progress of the chemical evolution of groundwater [1]. As a result, the chemical composition of the groundwater is spatially differentiated according to a certain law and exhibits horizontal zonation [1,4].
The western part of Jilin is a typical semi-arid water storage basin with less precipitation, strong evaporation, fragile ecological environment, and human life, depending heavily on groundwater [17]. Therefore, this study considers the western part of Jilin as the study area and establishes a zoning model of groundwater chemical evolution under the guidance of hydrogeochemical differentiation theory. The formation of major ions, and the migration and accumulation of trace elements were simulated by PHREEQC v3.4.0 (United States Geological Survey), combined with the interpretation of various parameters (e.g., hydrodynamic parameters, ion ratio coefficient, and thermodynamic index). The hydrogeochemical evolution of typical semi-arid water storage basins is further qualitatively and quantitatively described. For other semi-arid water storage basins, this study of groundwater chemistry evolution has a guiding significance.

2. Materials and Methods

2.1. Study Area

The study area is located in the southwestern part of the Songnen Plain in northeastern China, including Baicheng City, and western part of Songyuan City (Figure 1), with a total area of 41,140 km2 and a population of about 4,406,300. The study area is in the transition zone between the humid monsoon area and the arid inland (Figure 1a) and is a typical semi-arid inland basin. The annual average temperature is 3–6 °C, with average annual precipitation of 400–500 mm, mainly concentrated from June to September, with an average annual evaporation of 1500–2000 mm (Figure 2).
The surface water system in the basin is not well developed with only the Taoer River flowing through it, and the northern part of the basin is bordered by the Nen River and Songhua River. As shown in Figure 1b, the terrain is high in the east, south, and west, and low in the middle and north. The northwest is the hilly area of the Mt. Daxing’an fold belt, the southwest is the Songliao watershed, the central part is a broad alluvial and lacustrine plain, and the east is the tableland. The Quaternary strata cover the entire area, but the thickness varies, generally from 40 to 80 m [10]. The main landform types are the western Mountainous area, Taoer alluvial fan, Songnen low plain, tableland and valley plain (Figure 1c).
The study area is a typical semi-closed water storage basin structure, and the hydrogeological conditions are mainly controlled by the geomorphology and climate. The phreatic water and pore-fissure confined water occur in the Quaternary strata and the former is the object in this study. In the Taoer alluvial fan, the phreatic aquifer is composed of middle and upper Pleistocene gravel, with a thickness of 10–40 m and a groundwater buried depth of 10–20 m, which is recharged by meteoric water infiltration, river leakage, and lateral runoff from the western mountainous area [10]. In the low plain, the upper Pleistocene silty sand is dominant in the phreatic aquifer, followed by loess loam, with a thickness of 3–20 m and a groundwater buried depth less than 5 m [10,18]. The low and flat terrain of the low plain makes it conducive for the phreatic aquifer to receive rainfall infiltration recharge. Furthermore, the lateral runoff from the alluvial fan and tableland is also a source of recharge. Evaporation and artificial mining are the main methods of discharge. Controlled by topography, phreatic water flows from the east, south, and west to the middle and north, and flows into the Nen River and Songhua River at the northeastern boundary (Figure 3).

2.2. Sample Collection and Analysis

Groundwater sampling and water level survey were conducted from 3 to 24 November 2017. A total of 142 phreatic water samples (Figure 1c) were collected regularly from pumping wells with depths of 10–30 m (all water samples were from the Quaternary phreatic aquifer). The sampling was carried out 5–10 min after the start of the pumping, to drain the stagnant water in the suction pipe. The water samples were filtered through a 0.45 µm membrane and added to a 350 mL polyethylene plastic bottle, which was rinsed three times with deionized water. Three bottles of water were taken from the same well, and 10% HNO3 was added to one of the bottles to make the water sample pH less than 2 for cation analysis. Water temperature and pH were tested on-site using the calibrated HANNA (HI99131) portable pH/temperature analyzer, and the alkalinity was determined on-site by the Gran titration.
The water samples were tested by Pony Testing International Group in Changchun within one week. The water quality testing method is based on the Standard Examination Methods for Drinking Water (GB5750-2006). The main cations (K+, Na+, Ca2+, Mg2+, Fe3+) were tested by inductively coupled plasma atomic emission spectrometry (ICP-AES). F, Cl, NO3, and SO42− were measured by an ion chromatograph. NO2 and NH4+ were determined by the ultraviolet and visible spectrophotometers, and the TDS was measured by an electric blast drying oven and an electronic analytical balance (vapor-drying method). The reliability of the water sample analysis data was checked by the relative error of the anion and cation milliequivalent, and the error of all water samples was less than 5%. Statistics of groundwater chemical data in each zone are presented in Table 1.

2.3. Construction of the Hydrogeochemical Zoning Model

Simple and representative models are indispensable in disciplinary research. Therefore, the hydrogeological entity “groundwater storage basin” that exists in reality is selected as the research object [1]. The earth’s crust, under the geological action forced by internal and external power, forms folds, faults, and various landforms, which determines the distribution framework of the mountains and plains, and the dynamic environment of the water flow determines the spatial distribution of the particles inside the frame. The loosely deposited substance of different particle sizes is regularly differentiated and deposited, thereby controlling the spatial distribution of the aquifer and aquitard, and the structure of the groundwater system.
In a typical water storage structure, from the recharge area to the runoff area and discharge area, the permeability of the aquifer, the groundwater flow rate, the alternating characteristics of the groundwater circulation, the geological and hydrogeochemical environment, and the water–gas–rock interaction strength exhibit changes on a regular basis, leading to a clear zoning of hydrogeochemical characteristics, especially in semi-arid areas where precipitation and evaporation differ greatly.
Therefore, a hydrogeochemical zoning model of a semi-arid water storage basin based on typical profiles can be constructed (Figure 4). The real groundwater chemical formation process is approached by three levels: geological and geomorphological conditions, hydrochemical processes, and water chemical composition characteristics. AutoCAD2018 is used to describe and build the zoning model.

2.4. Hydrogeochemical Simulation

The quantitative analysis of groundwater chemical formation (e.g., rock dissolution and precipitation, evaporation, cation exchange) in the basin was carried out using hydrogeochemical reverse simulation. Considering that the water chemical evolution is mainly reflected in changes in the macroscopic components (Na/Mg/Ca and HCO3/Cl/SO4) and trace components (F), the representative mineral calcite/dolomite/gypsum and halite of the typical water storage structure system were selected as the simulated mineral phase. In addition, the following factors are considered: (1) clay minerals are widely distributed and cation exchange is prone to occur, (2) pore-phreatic aquifer is an open system, as CO2 participates in water–rock interaction, and (3) the strong evaporation in the basin, which can be indicated by the escape of H2O from the groundwater system. The possible mineral phases also include CaX2/NaX, CO2, and H2O. According to the groundwater flow direction and hydrogeological conditions, two simulation paths of I1–I5 and II1–II4 are selected (Figure 3). Reverse simulation and saturation index calculations were performed using PHREEQC 3.4.0 (United States Geological Survey) [19] developed by the USGS (U.S. Geological Survey).

3. Results

3.1. Hydrogeochemical Zoning Model

The variation in phreatic water chemical composition is consistent with the changes of the regional hydrogeological conditions. From the piedmont to the center of the basin, as the terrain becomes smooth, the groundwater chemical type gradually transitions from HCO3-Ca(Ca·Mg) to HCO3-Ca·Na(Na·Ca) to HCO3-Na or HCO3·Cl-Na type, showing good “geomorphology-hydrogeochemical process-groundwater chemical composition” zonal character. The range of recharge zones, runoff zones and discharge zones can be divided based on the geological features and the groundwater dynamic environment of the basin, combined with the distribution of groundwater chemistry types (Figure 3).
A hydrogeochemical zoning model (Figure 4) was established by the profile A-A’ (Figure 3) passing through two main recharge areas and discharge areas in the basin. The model was divided into three layers from bottom to top. The first layer represented geological and hydrogeological conditions, including the topography, lithology, geological time and hydrodynamic characteristics of the study area. The second layer represented hydrogeochemical processes, which was divided into the recharge zone/runoff zone and the discharge zone in the horizontal direction according to the hydrodynamic environment and hydrochemistry types, and the hydrogeochemical action gradually changes from lixiviation to cation exchange/evaporation and concentration. The third layer represented the characteristics of groundwater chemical components, including chemistry type, TDS, main anion and cation and the characteristic element F. Another analogous zoning model showing hydrogeochemical evolution controlled by micro-topography was built by profile P-P’ (Figure 5).

3.2. Hydrogeochemical Simulation

The results of the hydrogeochemical simulation are shown in Table 2. In the alluvial fan recharge zone (I1–I2), CO2 enters the groundwater, all minerals are in a dissolved state, and evaporation and cation exchange are relatively weak. When the groundwater flows through the runoff zone (I2–I3), calcite precipitates first (Equation (1)), with a precipitation amount of 0.94 mmol/L, while other minerals are dissolved. Evaporation is enhanced with the evaporation multiple being 1.67. Cation exchange is significant, and Ca2+ in groundwater is replaced by the Na+ adsorbed by the rock and soil (Equation (3)); the exchange amount is 1.31 mmol/L (calculated as Na+). When entering the discharge area (I3–I4), dolomite (Equation (2)) and calcite precipitate at the same time and decarbonation occurs; CO2 escapes from the groundwater; the evaporation and cation exchange is remarkable.
In the discharge area, calcite and dolomite precipitate while the fluorite/gypsum and halite are always in a dissolved state. Reverse cation exchange occurs (Equation (4)), and the exchange amount is 1.22 mmol/L (calculated as Na+). The mineral dissolution and precipitation of path II is similar to that of path I. The recharge area is dominated by mineral dissolution; calcite and dolomite precipitate simultaneously at the runoff zone while CO2 escapes from the groundwater. The evaporation is enhanced but the cation exchange is not significant; the reverse cation exchange occurs in the discharge zone.
Ca2+ + 2HCO3 → CaCO3↓ + CO2↑ + H2O
Ca2+ + Mg2+ + 4HCO3 → CaMg(CO3)2↓ + 2CO2↑ + 2H2O
Ca2+(l) + 2NaX(s)→2Na+(l) + CaX2(s)
2Na+(l) + CaX2(s) → Ca2+(l) + 2NaX(s)

3.3. Qualitative/Quantitative Indicator

The qualitative and quantitative characterization of zonation is carried out from two levels as hydrogeological and hydrochemical conditions; the results are shown in Table 3. Hydrogeological conditions determine the water flow environment, which is the key to control the hydrochemical evolution of groundwater. It is qualitatively described by lithology and recharge conditions, and quantitatively analyzed by elevation, groundwater buried depth, and hydrodynamic index (hydraulic gradient, coefficient of permeability, γCl/γCa2+).
The hydrochemical conditions are quantitatively characterized by the coefficients of the ion ratio, the saturation index, and the slope of γ(Na+ + K+ − Cl) and γ(Ca2+ + Mg2+) − (SO42 + HCO3). The coefficients of the ion ratio can be used to confirm the genesis of groundwater or the source of groundwater chemical composition [6,20]; the mineral saturation index is usually understood as a thermodynamic index describing the ability of groundwater to dissolve a mineral, which can be used to judge the reaction state between water and minerals [1,21]. The relationship of γ(Na+ + K+ − Cl) and γ(Ca2+ + Mg2+) − (SO42 + HCO3) can be used to indicate the extent of cation exchange; if the two are linear and the slope(kcation exchange) is close to −1, the cation exchange is significant [6]. The qualitative and quantitative description of each indicator in a different area (Figure 6 and Figure 7, and Table 3) explain the evolution and formation process of groundwater chemistry.

4. Discussion

4.1. Hydrogeochemical Zoning Model

4.1.1. Geological and Hydrogeological Conditions

The prototype of the Songnen Plain, the Delonggang Penelplane, was formed in the middle of the Middle Pleistocene, and the neotectonic movement formed it into the Songnen Plain known today [10,22]. The study area is located in the southern part of the Songnen Plain as a water-storage tectonic basin that opens to the north (Figure 1b). According to the difference in the form and intensity of the new tectonic movement, the study area is divided into five landforms: western mountainous area, the Taoer alluvial fan, the Songnen low plain, the tableland, and the valley plain (Figure 1c).
The Taoer alluvial fan is the transition zone from the uplift zone of the Mt. Daxing’an fold belt to the fault basin, and the upstream slope is steep (15‰), while the downstream is slow (7‰), with an elevation of 143–210 m. The aquifer consists of the middle Pleistocene alluvial gravel and sand gravel, which is overlying 0–2 m sandy loam, and the aquifer has good permeability, with a permeability coefficient of 50–200 m/day. The central low plain is a subsidence area, and the ground undulates in the form of microwaves, with an elevation of 130–150 m. It has long been accepted through sedimentary and accumulational evidence that the aeolian landform developed well. The upper part of the aquifer is loess loam while the lower part is fine sand or sandy loam, and the thickness is 5–25 m with low permeability (the permeability coefficient is only 0.2–0.3 m/day), and the underlying layer is thick silty clay. Due to the ground uplift, the eastern tableland has formed a ridge-like watershed with steep slopes and strong ground cutting, with an elevation of 173–290 m. The aquifer is mainly gravel in the north and Neogene sandstone in the south.
The buried depth and hydraulic gradient of groundwater in the basin are also controlled by the new tectonic movement. For the groundwater buried depth, generally the uplift area is larger than the fault area. The eastern tableland (10–40 m) and alluvial fan (7–15 m) have larger buried depths, while other areas are generally 1–5 m and the deeper depth can reach 5–10 m. The hydraulic gradient decreases from the uplift zone to the fault zone, the tableland (1.111 × 10−3) and the alluvial fan (0.667 × 10−3) are larger with great runoff conditions, and the low plain runoff is very sluggish (0.109 × 10−3). Atmospheric precipitation and leakage of the Taoer river are the main sources of recharge in the study area, and the discharge is mainly lateral runoff in the uplift area, while evaporation occurs in the fault area.
The neotectonic movement in the study area not only controls the formation of the landscape but also affects the dynamic environment of the groundwater, which changes the structure, lithology, burial, and circulation conditions of the aquifer system regularly from the edge of the basin to the center, and is a basis for regional groundwater hydrogeochemical evolution and differentiation.

4.1.2. Hydrogeochemical Processes

Under certain climatic and neotectonic conditions, the hydrogeochemical evolution of the groundwater system in a semi-arid basin presents typical zonal characteristics [1,5]. As the flow of water is the force of element migration, dispersion, and enrichment, starting from the groundwater dynamic environment, the recharge area (flow fast), runoff area (flow medium), and the discharge area (flow slow) are partitioned on the basis of the hydrochemistry to explain the process of hydrogeochemical evolution.

Recharge Area

The recharge area is located in the alluvial fan and eastern tableland. The lithology is mainly gravel or sand gravel, which makes the recharge area permeability coefficient 50–200 m/day in alluvial fan and 10–30 m/day in tableland, which is high. Due to the strong terrain cutting, the tableland is relatively steep, and its hydraulic gradient (1.111‰) is larger than that in alluvial fan (0.667‰), indicating that the hydrodynamic environment and runoff conditions are better in the recharge area, and the flow is dominated by horizontal movement.
The water chemistry in this environment is dominated by leaching [23]. The results of the hydrogeochemical simulation show that the water–rock reaction is mainly the dissolution of carbonate minerals. The good hydrodynamic environment and open system are conducive to the continuous entry of CO2, which strengthens its solvency. Long-term lixiviation causes the soluble minerals (chlorides/sulfates) in the aquifer to migrate into the water and be depleted, which makes the aquifer dominated by insoluble carbonates or silicate minerals [22,24]. Strong water alteration causes the chemical elements in the groundwater to remain in a migrated state, and the buried depth (7–40 m) is generally greater than the ultimate evaporation depth, causing the evaporation intensity of the phreatic water to be extremely weak (even if the pan evaporation is large) [25,26] and elements are difficult to accumulate. Therefore, the groundwater anion is dominated by HCO3, and the cation is dominated by Ca2+, forming HCO3-Ca type water with TDS less than 0.6 g/L, and the characteristic element F < 0.6 mg/L (Figure 4).

Runoff Area

As the water enters the fault basin and flows to the front edge of the alluvial fan and tableland, the terrain becomes slower and the hydraulic gradient becomes smaller (0.318‰ and 0.624‰ in the west and east, respectively). At the same time, the sediment changes from the alluvial phase to lacustrine phase, and the lithology of the aquifer becomes sandy loam, loess loam, or fine sand, with the permeability coefficient 3–5 m/day, which hinders the horizontal movement of groundwater, reduces the buried depth of groundwater (4–15 m), enhances evaporation and concentration, and accumulates chemical elements in the groundwater to some extent.
The results of simulation show that calcite and dolomite are saturated and precipitated successively in the runoff area, accompanied with the escape of CO2, but gypsum, halite, and fluorite are still in the dissolved state. The content of easily transportable elemental Na increased rapidly, while the increase in Cl content was slow (Figure 4), indicating that the increase in Na+ was caused not only by the dissolution of halite and evaporation [11,20]. As groundwater flows through the thick clay layer of the middle Pleistocene, it is prone to cation exchange, that is, Ca2+ in water is easily adsorbed by clay particles, and Na+ is desorbed into groundwater (Equation (3)). The simulation results show that the exchange amount is up to 1.82 mmol/L (calculated as Na+), while the eastern runoff area is smaller due to the thinner clay layer with a cation exchange amount of 0.28 mmol/L (calculated as Na+). Therefore, the water chemistry in the runoff area is relatively complicated, and there are some differences between east and west. Evaporation concentration and cation exchange cause a significant increase in readily soluble Na+, while the hydrochemical type changes from HCO3-Ca to HCO3-Ca·Na(Na·Ca) and HCO3-Na, with TDS reaching 0.3–1 g/L and F increasing to 0.5–2.5 mg/L.

Discharge Area

In the discharge area of the low plain, the terrain is low, and the groundwater runoff is sluggish, with a permeability coefficient of 0.2–0.3 m/day and a hydraulic gradient of 0.109‰. This causes the groundwater alternate speed to be very slow, and the movement is dominated by vertical movement while horizontal movement is extremely weak. Precipitation has become the main and the only source of recharge, while the amount of evaporation is 4.33 times that of precipitation (Figure 2), and the groundwater buried depth is generally shallow (usually less than 5 m), which makes evaporation and concentration the main discharge method and hydrochemical process.
As the concentrations of groundwater, calcite, and dolomite continue to precipitate, halite, gypsum, and fluorite are generally dissolved, and fluorites are over-saturated in five samples. The readily soluble elements, Na and Cl, peak in this zone. With the increase in Na+, reverse cation exchange occurs, that is, Na+ in water replaces Ca2+ adsorbed by rock (Equation (4)). The simulation results show that the exchange amount is 2.44 mmol/L and 0.86 mmol/L (calculated by Na+) in the west and east, respectively. Due to the uneven surface of the low plain, the aeolian landforms such as sandy dunes, sandy ridges, and sandy hills developed well, and they are distributed with the fulje. As shown in Figure 5, slow runoff and strong evaporation cause the chemical composition and salt transport of groundwater to be significantly affected by micro-geomorphology, and the changes in water chemistry are more complicated. Groundwater was buried deeper in the sandy hill with relatively high topography, and the effect of evaporation is relatively weak; the hydrochemical type is mostly HCO3-Na·Mg(Na·Ca) or HCO3-Na type with TDS of 0.8–1.5 g/L and F of 1–2.5 mg/L. Groundwater in the fulje is shallow and accompanied with strong evaporation, so the chemical type becomes HCO3·Cl-Na or HCO3·Cl-Na·Ca, with TDS of 1.7–4.3 g/L and F of 1.5–4 mg/L (the highest is 7.99 mg/L).

4.1.3. Characteristics of Groundwater Chemical Components

The neotectonic movement determines the spatial distribution of modern geomorphological landscapes of the basin [10]. The groundwater chemistry in different geomorphological units showing regular changes due to external forces (flowing water, sandstorm, climate). The dissolution or precipitation of the mineral and the accumulation or migration of the element caused spatial zoning of groundwater chemical components, which has the same evolutionary rules as many semi-arid basins [27].
The hydrochemistry type has typical horizontal zonation as the recharge area at the edge of the basin is mainly HCO3-Ca type, occasionally HCO3-Ca·Mg type, and changed from HCO3-Na·Ca(Ca·Na) to HCO3-Na through the runoff area. Because of the effects of evaporation, cation exchange, and micro-geomorphology, the hydrochemistry type is relatively complex, such that HCO3-Na·Mg(Na·Ca) and HCO3·Cl-Na (Na·Ca) type appear locally.
TDS also showed regular changes, as TDS in the recharge area was lower than 0.6 g/L, while in the discharge area it was generally higher than 0.8 g/L, and a few were higher than 1.5 g/L in the low-lying fulje under the micro-geomorphic effect. As Table 4 shows, 15.38%, 20.00%, and 39.53% of sample TDS in the recharge, runoff, and discharge areas, respectively, exceeded the standard of drinking water (the limit of TDS is 1 g/L in Sanitary Standard for Drinking Water of China).
Trace element F, which is closely related to human health, also exhibits change on a regular basis. In the recharge area, where the groundwater alteration is strong, F is easy to migrate and accumulate slowly, with a content less than 0.6 mg/L. When the groundwater reached the runoff zone, F begins to accumulate, as the content is greater than 0.5 mg/L. In the central discharge area of the basin, the strong evaporation and stagnation flow cause the elements to accumulate rapidly. The alkaline water environment (average of pH is 7.69) and the soda type of water chemistry are also conducive to the enrichment of F [28,29], finally forming high-fluorine water, which is not suitable for drinking. As Table 4 shows, 20.51%, 56.67%, 67.44% of sample F in the recharge, runoff, and discharge area, respectively, exceeded the standard of drinking water (the limit of F is 1 g/L in Sanitary Standard for Drinking Water of China).

4.2. Quantitative Analysis of Using Different Indicators

The differentiation of groundwater chemical characteristics in different zones can be explained and characterized by different indicators.

4.2.1. Hydrodynamic Indicator

As the carrier of dissolved minerals, water flow is the basic force for the migration, dispersion, and accumulation of elements in the groundwater. The change in the hydrodynamic conditions is caused by the friction of the surface of the particles when the water flows in the porous medium [1,30]. Therefore, the hydrodynamic index can reflect the time or extent of the interaction between the groundwater and porous medium.
The hydraulic gradient (0.667–1.111‰) and permeability coefficient (15–200 m/day) in the recharge area are much larger than that in the discharge area (hydraulic gradient: 0.109‰, permeability coefficient: 0.2–0.3 m/day) (Table 3), reflecting a relatively good hydrodynamic environment and strong water alteration. In addition, γCl/γCa2+ is one of the parameters that characterizes the hydrodynamics, as Cl is usually enriched in the stagnant region of groundwater and the cations in the rapid water flow are dominated by Ca2+. Therefore, high γCl/γCa2+ represents slow hydrodynamics condition. The increase in γCl/γCa2+ in the study area from the recharge zone to the runoff zone was slow, while it increased rapidly from the runoff zone to the discharge zone (Figure 6), reflecting the change in the hydrodynamic conditions from strong to weak, and the discharge zone was much weaker than that in the recharge and runoff area. Changes in regional hydrodynamic indicators and water chemistry indicate that the two are not independent but interrelated [1,31], especially in the semi-arid water storage basin, where a good hydrodynamic environment is characterized by elemental migration and low TDS water, while a poor hydrodynamic environment is characterized by elemental enrichment and higher TDS water.

4.2.2. Cation Exchange Index

Cation exchange is actually an adsorption–desorption process, and the occurrence possibility can be indicated by the relationship between γ(Na+ + K+ − Cl) and γ(Ca2+ + Mg2+) − (SO42 + HCO3) [20,32]. Except for the dissolution of halite, γ(Na+ + K+ − Cl) indicates the increase or decrease of Na+, while γ(Ca2+ + Mg2+) − (SO42 + HCO3) indicates the increase or decrease of Ca2+ and Mg2+, except for the dissolution of calcite, dolomite, and gypsum. When the cation exchange is significant, the relationship between the two is linear and the slope is close to −1. However, this parameter can only reflect the possibility of cation exchange but not the direction of exchange. As shown in Figure 6b, the western runoff zone (k = −0.8, R2 = 0.84) and the discharge zone (k = −1.05, R2 = 0.86) have a good linear relationship, indicating significant cation exchange, while the recharge zone and the eastern runoff zone do not.
According to [33], the clay minerals in the study area are mainly montmorillonite and illite, and the content of kaolinite is relatively small. As the montmorillonite (80–150 meq/100 g) and illite (10–40 meq/100 g) have the bigger CEC (Cation Exchange Capacity) and the CEC of kaolinite is only 3–15 meq/100 g [34], the clay minerals in the basin have strong cation exchange capacity. The alkaline pH (Table 1) in the basin can increase the variable negative charge in the clay and increase its CEC value [9,35]. Therefore, the strength of the cation exchange reaction in the basin is controlled by the distribution in the clay layer (Figure 4), which is more significant in the western runoff area and the middle discharge area with a thicker clay layer and consistent with the results of the cation exchange index (Figure 6b) and reverse simulation (Table 2).

4.2.3. Ion Ratio Coefficient

The ion ratio coefficient can be used to determine the origin of groundwater or the source or formation process of the groundwater chemical component [20,36]. The cation ratio γCa2+/γNa+, γMg2+/γNa+, and γCa2+/γMg2+ showed a regular downward trend from the recharge to the discharge zone, which indicated the gradual enrichment of Na+ and Mg2+ in groundwater. The migration order of the major cation is Na+ → Mg2+ → Ca2+, instead of the generally considered Na+ → Mg2+ → Ca2+ [34], which may be related to significant cation exchange. The authors of [37] also reached the parallel conclusion in the study of groundwater chemistry in the arid region.
Different from cations, anions γCl/γHCO3 and γSO42−/γHCO3 are all less than 1, and the trend from the recharge to the runoff zone shows a decreasing trend, indicating that the leaching effect is still significant, causing the increase in HCO3 to be greater than that in Cl or SO42−. When the groundwater enters the discharge area, the content of HCO3 reaches the peak without continuous increase (Figure 4), which is attributed to the saturation state of calcite and dolomite, while Cl and SO42− can continue to accumulate because of the unsaturated state of halite and gypsum, so that γCl/γHCO3 and γSO42−/γHCO3 are increased. The change in γCl/γSO42− indicates that the difference between the content of Cl and SO42− in the recharge and runoff regions is small, and the cumulative velocity of Cl is faster than that of SO42− in the discharge zone, but relatively speaking, HCO3 is still the dominant anion.
The value of γNa+/γCl in the recharge zone is greater than 1, that is, γNa+ is generally higher than γCl, which indicates that the dissolution of halite is not the only source of Na+, and the remaining Na+ may be derived from the weathering and dissolution of aluminosilicate, such as albite (Equation (5)) and sodium montmorillonite (Equation (6)). The abnormal increase in γNa+/γCl in the western runoff region may be the result of cation exchange (Equation (1)), while the decrease in the discharge zone may be affected by reverse cation exchange (Equation (4)).
4NaAlSi3O8 + 4CO2 + 22H2O→Al4(Si4O10)(OH)8 + 8H4SiO4 + 4Na+ + 4HCO3
3Na1/3Al7/3Si11/3O10 + 30H2O + 6OH→Na+ + 7Al(OH)4 + 11H4SiO4
The ratio γCa2+/γHCO3 is less than 1, indicating that the dissolution of dolomite is also an important source of Ca2+. The trend of decrease first then increase confirms the existence of the positive ion exchange in the runoff zone and the reverse cation exchange in the discharge zone. The value of γCa2+/γSO42− decreased significantly from the recharge zone to the discharge zone, indicating that the enrichment of SO42− is greater than that of Ca2+. As Figure 7 shows, gypsum is always in dissolved mode while calcite and dolomite undergo precipitation, indicating that the enrichment of Ca2+ may be limited by the cation exchange (Equation (3)) and the precipitation of calcite and dolomite (Equations (1) and (2)).

4.2.4. Thermodynamic Index

Weathering is the basis of groundwater chemical formation, which is reflected in the dissolution and precipitation of minerals in the water–rock interaction. The thermodynamic index of the ability of groundwater to dissolve minerals is usually expressed by the saturation index of minerals.
Figure 7 shows the tendency of the calcite and dolomite in the groundwater to gradually evolve from unsaturated to saturated from the recharge zone to the discharge zone. A considerable number of water samples have been found to be over-saturated in the recharge zone, indicating that the transition from unsaturated to saturated occurs in the recharge zone. In addition, there are still water samples in the runoff zone and the discharge zone that are unsaturated, which may be related to cation exchange and the microtopography effect.
For halite and gypsum, all water samples are in an unsaturated state, and the saturation index is between −3.24 and −0.18 and between −9.21 and −4.95, as it tends to increase from the recharge zone to the discharge zone, indicating that halite and gypsum are continuously dissolved. The unsaturation state of halite depends on its large solubility of chlorides, while the unsaturation state of gypsum is related to the alluvial and lacustrine aquifers, as there are relatively few sulfate minerals in the sedimentary facies [22,27].
For fluorite, only five water samples are saturated, and all of them are located at the discharge zone, where fluorite precipitation occurs. The elevation of these five sampling points is 125–136 m, which is significantly lower than the elevation of other sampling points (138–157 m) in the discharge area, that is, the sampling point is located at the low-lying area. Therefore, the dissolution and precipitation of fluorite, and the enrichment of F in the discharge area are mainly controlled by the microgeomorphology effect.
The calculation results of the saturation index of the study area are consistent with the results of the hydrogeochemical simulation. The state of minerals in the groundwater is also consistent with other semi-arid water storage basins, such as the Yinchuan Plain of China [4] and Weining Plain of China [37], Central Tunisia [6], and Hamersley Basin of northwest Australia [38]. Calcite and dolomite from the recharge zone to the discharge zone change from unsaturated to saturated, gypsum and rock salt are generally unsaturated, and fluorite is saturated only in the discharge area.

5. Conclusions

(1)
Taking the western part of Jilin Province as an example, a hydrogeochemical evolutionary zoning model of a typical semi-arid water storage basin was established. The model is divided into three layers from bottom to top. The first layer represents the geological and hydrogeological conditions, including the topography, lithology, geological time, and hydrodynamic characteristics. The second layer reveals the hydrogeochemical processes, divided into the recharge zone, runoff zone, and discharge zone in the horizontal direction. The third layer is the characteristics of the groundwater chemical components, including chemistry type, TDS, main anion and cation, and characteristic element F. This zoning model shows that the hydrogeochemical action gradually changes from lixiviation to cation exchange, evaporation controls the hydrogeochemical evolution, and landform plays a key role in hydrochemistry formation in the discharge area.
(2)
The quantitative analysis was carried out by hydrogeochemical reverse simulation. The results showed that from the recharge to the discharge zone, calcite and dolomite change from unsaturated to saturated, synchronized with the entry and escape of CO2 from the groundwater. Fluorite, gypsum, and halite are always in a dissolved state. The evaporation is gradually enhanced; the cation exchange changed from Na+ released into the water to Ca2+ released into the water. The differences in the lithology of the formation make the exchange amount quite different in the runoff area of the east and west.
(3)
The hydraulic gradient, permeability coefficient, and γCl/γCa2+ indicate that the change in the hydrodynamic environment worsens from the recharge zone to the discharge zone. The cation exchange index proves that the cation exchange is more significant in the western runoff zone and the discharge zone, while γCa2+/γNa+, γMg2+/γNa+, and γCa2+/γMg2+ show that the migration order of cations from the recharge zone to discharge zone is Na+→Mg2+→Ca2+. The γHCO3/γCl. γSO42−/γCl, γSO42−/γHCO3 results show that from the recharge to the runoff zone, the content of HCO3 shows a decreasing trend while reaching the peak without continuous increase in the discharge area, which is attributed to the saturation state of calcite and dolomite. Meanwhile, HCO3 is still the dominant anion in groundwater, and γNa+/γCl, γCa2+/γHCO3, and γCa2+/γSO42 indicate the source of Na+ and Ca2+, demonstrating the positive cation exchange in the western runoff zone and the reverse cation exchange in the discharge zone. The saturation indexes indicate that the calcite and dolomite from the recharge zone to the discharge zone change from unsaturated to saturated, gypsum and rock salt are generally unsaturated, and fluorite is saturated only in the discharge area, which is consistent with the results of the reverse simulation and consistent with the other semi-arid reservoirs, as well as with the results of hydrogeochemical simulation and other semi-arid water storage basins.

Author Contributions

M.L., S.H., and Y.C. processed the data and analyzed the results; M.L. wrote the manuscript; Y.C, X.L, and C.X. reviewed the manuscript and made helpful suggestions; M.L. revised the manuscript.

Funding

The study was financially supported by Natural Science Foundation of China (No. 41572216), the China Geological Survey Shenyang Geological Survey Center “Changji Economic Circle Geological Environment Survey” project (121201007000150012), the Provincial School Co-construction Project Special-Leading Technology Guide (SXGJQY2017-6), and the Jilin Province Key Geological Foundation Project (2014–13).

Acknowledgments

We would like to thank the anonymous reviewers and the editor.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cao, Y.Q.; Hu, K.R.; Li, Z.S. Groundwater Chemical Kinetics and EcoEnvironmental Zonation; Science Press: Beijing, China, 2009. [Google Scholar]
  2. Wang, W.; Song, X.; Ma, Y. Characterization of controlling hydrogeochemical processes using factor analysis in Puyang Yellow River irrigation district (China). Hydrol. Res. 2017, 48, 1438–1454. [Google Scholar] [CrossRef]
  3. Llewellyn, G.T.; Dorman, F.; Westland, J.L.; Yoxtheimer, D.; Grieve, P.; Sowers, T.; Humstonfulmer, E.; Brantley, S.L. Evaluating a groundwater supply contamination incident attributed to Marcellus Shale gas development. PNAS 2015, 112, 6325–6330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jing, X.-y.; Yang, H.; Cao, Y.; Wang, W. Identification of indicators of groundwater quality formation process using a zoning model. J. Hydrol. 2014, 514, 30–40. [Google Scholar] [CrossRef]
  5. Jia, Y.; Guo, H.; Jiang, Y.; Wu, Y.; Zhou, Y. Hydrogeochemical zonation and its implication for arsenic mobilization in deep groundwaters near alluvial fans in the Hetao Basin, Inner Mongolia. J. Hydrol. 2014, 518, 410–420. [Google Scholar] [CrossRef]
  6. Farid, I.; Trabelsi, R.; Zouari, K.; Abid, K.; Ayachi, M. Hydrogeochemical processes affecting groundwater in an irrigated land in; Central Tunisia. Environ. Earth Sci. 2013, 68, 1215–1231. [Google Scholar] [CrossRef]
  7. Yan, J.J.; Xu, F.; Kang, C.J.; Higano, Y. Effective stockbreeding biomass resource use and its impact on water environment from the viewpoint of sustainable development. J. Dev. Sustainable Agric. 2010, 5, 147–150. [Google Scholar]
  8. Wang, W.K.; Yang, Z.Y.; Cheng, D.H.; Wang, W.M. Methods of Ecology-Oriented Groundwater Resource Assessment in Arid and Semi-Arid Area. J. Jilin Univ. 2011, 41, 159–167. (In Chinese) [Google Scholar]
  9. Shen, Z. Fundamentals of Hydrogeochemistry; Geological Publishing House: Beijing, China, 1986. (In Chinese) [Google Scholar]
  10. Wang, Z.X.; Su, Q.S.; Cao, Y.Q. Groundwater and Quaternary Geology in Baicheng Area; Geological Publishing House: Beijing, China, 1984. (In Chinese) [Google Scholar]
  11. Guo, H.; Wang, Y. Geochemical characteristics of shallow groundwater in Datong basin, Northwestern China. J Geochem Explor 2005, 87, 109–120. [Google Scholar] [CrossRef]
  12. Xing, L.; Guo, H.; Zhan, Y. Groundwater hydrochemical characteristics and processes along flow paths in the North China Plain. J. Asian Earth Sci. 2013, 70, 250–264. [Google Scholar] [CrossRef]
  13. Somay, M.A.; Gemici, Ü. Assessment of the Salinization Process at the Coastal Area with Hydrogeochemical Tools and Geographical Information Systems (GIS): Selçuk Plain, Izmir, Turkey. Water Air Soil Pollut. 2009, 201, 55–74. [Google Scholar] [CrossRef]
  14. Dragon, K.; Górski, J. Identification of hydrogeochemical zones in postglacial buried valley aquifer (Wielkopolska Buried Valley aquifer, Poland). Environ. Geol. 2009, 58, 859–866. [Google Scholar] [CrossRef]
  15. Kløve, B.; Ala-Aho, P.; Bertrand, G.; Gurdak, J.J.; Kupfersberger, H.; Kværner, J.; Muotka, T.; Mykrä, H.; Preda, E.; Rossi, P. Climate change impacts on groundwater and dependent ecosystems. J. Hydrol. 2014, 518 (Pt B), 250–266. [Google Scholar] [CrossRef]
  16. Huang, G.; Sun, J.; Zhang, Y.; Chen, Z.; Liu, F. Impact of anthropogenic and natural processes on the evolution of groundwater chemistry in a rapidly urbanized coastal area, South China. Sci. Total Environ. 2013, 463, 209–221. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, B.; Song, X.; Zhang, Y.; Han, D.; Tang, C.; Yu, Y.; Ma, Y. Hydrochemical characteristics and water quality assessment of surface water and groundwater in Songnen plain, Northeast China. Water Res. 2012, 46, 2737–2748. [Google Scholar] [CrossRef] [PubMed]
  18. Jilin-University. The report of comprehensive hydrogeological atlas compilation in Jilin province; Jilin-University: Chuangchun, China, 2018. [Google Scholar]
  19. Parkhurst, D.L.; Appelo, C.A.J. User’s guide to PHREEQC (Version 2): A Computer Program for Speciation, Batch-Reaction, One-Dimensional Transport, and Inverse Geochemical Calculations; US Department of the Interior, US Geological Survey: Washington, DC, USA, 1999.
  20. García, G.; del, V.; Hidalgo, M.; Blesa, M. Geochemistry of groundwater in the alluvial plain of Tucumán province, Argentina. Hydrogeol. J. 2001, 9, 597–610. [Google Scholar]
  21. Razowska, L. Changes of groundwater chemistry caused by the flooding of iron mines (Czestochowa region, Southern Poland). J. Hydrol. 2001, 244, 17–32. [Google Scholar] [CrossRef]
  22. 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. Processes 2010, 21, 1055–1065. [Google Scholar] [CrossRef]
  23. Devic, G.; Djordjevic, D.; Sakan, S. Natural and anthropogenic factors affecting the groundwater quality in Serbia. Sci. Total Environ. 2014, 468, 933–942. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, Y.; Shvartsev, S.L.; Su, C. Genesis of arsenic/fluoride-enriched soda water: A case study at Datong, northern China. Appl. Geochem. 2009, 24, 641–64912. [Google Scholar] [CrossRef]
  25. Rose, D.A.; Konukcu, F.; Gowing, J.W. Effect of watertable depth on evaporation and salt accumulation from saline groundwater. Soil Res. 2005, 43, 565–573. [Google Scholar] [CrossRef]
  26. Yang, T.; Zhang, Q.; Wang, W.; Yu, Z.; Chen, Y.D.; Lu, G.; Hao, Z.; Baron, A.; Zhao, C.; Chen, X.; et al. Review of Advances in Hydrologic Science in China in the Last Decades: Impact Study of Climate Change and Human Activities. J. Hydrol. Eng. 2013, 18, 1380–1384. [Google Scholar] [CrossRef]
  27. Liao, Z.S.; Lin, X.Y. Chemical Characteristics and Variations of Groundwater Quality in Songnen Basin. Earth Sci. 2004, 29, 96–102. [Google Scholar]
  28. Brahman, K.D.; Kazi, T.G.; Afridi, H.I.; Naseem, S.; Arain, S.S.; Ullah, N. Evaluation of high levels of fluoride, arsenic species and other physicochemical parameters in underground water of two sub districts of Tharparkar, Pakistan: A multivariate study. Water Res. 2013, 47, 1005–1020. [Google Scholar] [CrossRef] [PubMed]
  29. Zabala, M.E.; Manzano, M.; Vives, L. Assessment of processes controlling the regional distribution of fluoride and arsenic in groundwater of the Pampeano Aquifer in the Del Azul Creek basin (Argentina). J. Hydrol. 2016, 541, 1067–1087. [Google Scholar] [CrossRef]
  30. Appelo, C.A.J.; Postma, D. Geochemistry, Groundwater and Pollution; Balkema: Rotterdam, The Netherland, 2005; pp. 256–270. [Google Scholar]
  31. Morin, K.A.; Cherry, J.A.; Dave, N.K.; Lim, T.P.; Vivyurka, A.J. Migration of acidic groundwater seepage from uranium-tailings impoundments, 1. Field study and conceptual hydrogeochemical model. J. Contam. Hydrol. 1988, 2, 271–303. [Google Scholar] [CrossRef]
  32. Hidalgo, M.C.; Cruz-Sanjulián, J. Groundwater composition, hydrochemical evolution and mass transfer in a regional detrital aquifer (Baza basin, southern Spain). Appl. Geochem. 2001, 16, 745–758. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, W.H. A Study on the Moisture Content Migration and Characteristics of Frost Heaving of Saline Soil in the Western of Jilin Province; Jilin University: ChangChun, China, 2011. (In Chinese) [Google Scholar]
  34. Caruccio, F.T. The Properties of Groundwater; John Wiley and Sons Inc.: New York, NY, USA, 1982. [Google Scholar]
  35. Matschonat, G.; Vogt, R. Effects of changes in pH, ionic strength, and sulphate concentration on the CEC of temperate acid forest soils. Eur. J. Soil Sci. 2010, 48, 163–171. [Google Scholar] [CrossRef]
  36. Starr, R.C.; Gillham, R.W. Controls on denitrification in shallow unconfined aquifers. In Contaminant Transport in Groundwater. Proc. of the International Symposium on Contaminant Transport in Groundwater; CRC Press: Boca Raton, FL, USA, 1989; pp. 51–56. [Google Scholar]
  37. Jiang, L.; Peicheng, L.I.; Guo, J. Hydrochemical Characteristics and Evolution Laws of Groundwater in Typical Oasis of Arid Areas on the West of Helan Mountain. J. Earth Sci. Environ. 2009, 31, 285–290. [Google Scholar]
  38. Dogramaci, S.; Grzegorz, D.; Wade, G.; Pauline, F. Stable isotope and hydrochemical evolution of groundwater in the semi-arid Hamersley Basin of subtropical northwest Australia. J. Hydrol. 2012, 475, 281–293. [Google Scholar] [CrossRef]
Figure 1. Location of the study area. (a) Songnen Plain in China and China climate zonation, (b) the terrain of the Songnen Plain and the study area in the Songnen Plain, (c) landform division and sampling point distribution map of the study area.
Figure 1. Location of the study area. (a) Songnen Plain in China and China climate zonation, (b) the terrain of the Songnen Plain and the study area in the Songnen Plain, (c) landform division and sampling point distribution map of the study area.
Water 11 01334 g001
Figure 2. Monthly mean value of precipitation and evaporation in the study area from 1951 to 2017.
Figure 2. Monthly mean value of precipitation and evaporation in the study area from 1951 to 2017.
Water 11 01334 g002
Figure 3. Groundwater chemical zonation/flow field and reverse chemical simulation path in the study area.
Figure 3. Groundwater chemical zonation/flow field and reverse chemical simulation path in the study area.
Water 11 01334 g003
Figure 4. Zoning model for the hydrogeochemical evolution of the study area. The color of stratum: Blue for aquifer, green for aquitard and yellow for aquiclude. Legend: (1) Holocene, (2) upper Pleistocene, (3) middle Pleistocene lower Cretaceous Nenjiang Formation, (4) lower Pleistocene, (5) Neogene system Taikang Formation, (6) lower Cretaceous Nenjiang Formation, (7) lower Cretaceous Yaojia Formation, (8) loess loam, (9) fine sand, (10) silty loam, (11) sandy gravel, (12) mudstone, (13) sandstone, (14) glutenite, (15) unconfined water level, (16) angular unconformity interface, (17) F, (18) TDS, (19) HCO3, (20) SO42−, (21) Cl, (22) Mg2+, (23) Na+, (24) Ca2+.
Figure 4. Zoning model for the hydrogeochemical evolution of the study area. The color of stratum: Blue for aquifer, green for aquitard and yellow for aquiclude. Legend: (1) Holocene, (2) upper Pleistocene, (3) middle Pleistocene lower Cretaceous Nenjiang Formation, (4) lower Pleistocene, (5) Neogene system Taikang Formation, (6) lower Cretaceous Nenjiang Formation, (7) lower Cretaceous Yaojia Formation, (8) loess loam, (9) fine sand, (10) silty loam, (11) sandy gravel, (12) mudstone, (13) sandstone, (14) glutenite, (15) unconfined water level, (16) angular unconformity interface, (17) F, (18) TDS, (19) HCO3, (20) SO42−, (21) Cl, (22) Mg2+, (23) Na+, (24) Ca2+.
Water 11 01334 g004
Figure 5. Hydrogeochemical evolution controlled by micro-topography. (1) aeolian deposit of Holocene, (2) upper part of Upper Pleistocene (lacustrine sediment), (3) lower part of Upper Pleistocene (lacustrine sediment), (4) fine sand, (5) loess loam, (6) sandy loam, (7) unconfined water level, (8) groundwater flow direction, (9) salinity transport direction, and (10) evaporation.
Figure 5. Hydrogeochemical evolution controlled by micro-topography. (1) aeolian deposit of Holocene, (2) upper part of Upper Pleistocene (lacustrine sediment), (3) lower part of Upper Pleistocene (lacustrine sediment), (4) fine sand, (5) loess loam, (6) sandy loam, (7) unconfined water level, (8) groundwater flow direction, (9) salinity transport direction, and (10) evaporation.
Water 11 01334 g005
Figure 6. Quantitative index in different partitions (a) γCl/γCa2+, (b) slope of γ(Na+ + K+ − Cl) and γ(Ca2+ + Mg2+) − (SO42 + HCO3), (c) cation ratio coefficient, (d) anion ratio coefficient; (e) provenance index
Figure 6. Quantitative index in different partitions (a) γCl/γCa2+, (b) slope of γ(Na+ + K+ − Cl) and γ(Ca2+ + Mg2+) − (SO42 + HCO3), (c) cation ratio coefficient, (d) anion ratio coefficient; (e) provenance index
Water 11 01334 g006
Figure 7. Saturation index in different partitions: (a) Calcite (b) Dolomite (c) Gypsum (d) Halite (e) Fluorite.
Figure 7. Saturation index in different partitions: (a) Calcite (b) Dolomite (c) Gypsum (d) Halite (e) Fluorite.
Water 11 01334 g007
Table 1. Statistics of groundwater chemical data in each hydrogeochemical zone (in mg/L, figures in parentheses represent mean values).
Table 1. Statistics of groundwater chemical data in each hydrogeochemical zone (in mg/L, figures in parentheses represent mean values).
ParametersWestern Recharge
Area (n = 19)
Western Runoff Area (n = 46)Central Discharge Area (n = 43)Eastern Runoff
Area (n = 13)
Eastern Recharge Area (n = 21)
pH7.2–7.8 (7.57)7.1–8.2 (7.66)7.1–8.3 (7.69)7.3–8 (7.68)7.5–8 (7.73)
TDS188–710 (413.29)362–1020 (577.81)529–3500 (1109.55)271–773 (495.42)183–527 (355.67)
K+0.29–1.94 (0.78)0.29–2.32 (1.04)0.52–29.3 (2.19)0.43–1.99 (0.83)0.17–0.89 (0.48)
Na+13.1–51.5 (26.65)23.2–221 (84.4)28.2–710 (199.66)13.3–156 (59.39)6.07–48.1 (18.18)
Ca2+31.6–115 (72.02)19.9–141 (73.03)19.5–341 (101.34)44.3–103 (69.08)30.6–166 (83.34)
Mg2+7.69–42.7 (18.07)12–79.2 (32.05)8.89–182 (57.61)13.1–82.3 (32.68)4.57–21.3 (12)
Fe3+0.01–1.39 (0.17)0.01–43.2 (3.71)<0.01–21.5 (1.88)<0.01–0.65 (0.26)<0.01–1.35 (0.18)
Cl9–84.9 (25.4)2.99–313 (40.09)7.98–638 (148.4)7.08–121 (40.77)3.79–68 (21.28)
SO42−12–56.7 (32.93)1.78–178 (41.3)0.88–336 (77.85)1.73–52.2 (17.67)1.17–61.9 (17.78)
HCO3126–530 (257.35)210–1170 (512.19)289–1480 (598.07)276–848 (457.83)146–374 (271.87)
F0.3–0.68(0.58)0.53–1.74 (1.15)1.03–3.07 (1.62)0.6–2.5 (1.29)0.21–0.86 (0.63)
NO3-N0.59–38.9 (15.63)<0.01–57 (3.03)<0.01–210 (28.9)<0.01–5.17 (1.16)<0.01–30.4 (10.52)
NO2-N<0.001–0.03 (0.01)<0.001–0.06 (0.01)<0.001–0.16 (0.02)<0.001–0.19 (0.03)<0.001–0.26 (0.03)
NH4-N0.03–0.34 (0.11)0.03–2.95 (0.61)0.02–3.38 (0.4)0.03–0.54 (0.22)0.03–0.37 (0.11)
Table 2. The result of hydrogeochemical simulation.
Table 2. The result of hydrogeochemical simulation.
ItemsPathway I
(Alluvial Fan → Low Plain)
Pathway II
(Tableland → Low Plain)
I1–I2I2–I3I3–I4I4–I5II1–II2II2–II3II3–II4
Evaporation Multiple1.001.671.611.201.032.141.17
Mineral dissolution and precipitation 1Calcite0.64−0.94−0.61−0.670.49−1.06−1.96
Fluorite0.010.0010.010.010.030.0020.003
Dolomite0.330.41−0.74−1.090.22−0.46−1.51
Gypsum0.130.020.310.210.060.510.01
Halite0.150.250.473.260.081.320.53
CO2(g)1.010.85−4.05−3.570.65−0.14−1.03
H2O(g) 20.00−37.27−33.88−11.28−1.91−63.3−9.25
Cation
Exchange 3
CaX2−0.23−0.65−0.911.22−0.25−0.140.43
NaX0.461.311.82−2.440.500.28−0.86
1 The positive value in the column of mineral dissolution and precipitation represents dissolution, and the negative value indicates the precipitation; 2 The unit of H2O is mol/L, and the rest is mmol/L; 3 A positive value of NaX represents the reaction of Equation (3), and a negative value indicates the reaction of Equation (4).
Table 3. Qualitative and quantitative indexes of hydrogeochemical zonation.
Table 3. Qualitative and quantitative indexes of hydrogeochemical zonation.
IndexesZones
Western Recharge AreaWestern Runoff AreaCentral Discharge AreaEastern Runoff AreaEastern Recharge Area
GeomorphologyLandformTaoer alluvial fanLow plainTableland
Elevation/m143–210135–150140–150 (hills)
130–140 (fulje)
155–230173–290
Qualitative description
of aquifer
Lithology
of aquifer
Gravel/sand gravel, with 0–2 m sandy loam covering on the surfaceThe upper part is sandy loam and the lower part is fine sandThe upper part is loess loam while lower part is fine sand or sandy loamThe upper part is loess loam or fine sand and the lower part is gravelThe northern part is a gravel and the southern part is Neogene sandstone
Quantitative
description
of aquifer
Hydraulic
gradient/10−3
0.6670.3180.1090.6241.111
Buried depth/m7–154–12<54–1510–40
Permeate
coefficient/m·day−1
50–2003–50.2–0.33–415–20
γCl/γCa2+0.320.380.860.360.21
Cation exchangekcation exchange−1.341
(R2 = 0.610)
−0.796
(R2 = 0.845)
−1.048
(R2 = 0.863)
−0.660
(R2 = 0.648)
−1.939
(R2 = 0.738)
Ion ratioγCa/γNa3.2731.2850.9312.1436.089
γMg/γNa1.4060.9150.7741.2711.442
γCa/γMg2.4191.4441.2651.5974.439
γCl/γHOC30.2100.2180.3800.1590.230
γSO4/γHCO30.2100.2300.2300.0500.130
γNa/γCl2.1874.3202.6202.6202.154
γCa/γSO45.0003.7602.9905.4406.020
γCa/γHCO30.9190.4950.6450.4921.130
Saturation
index
SIcalcite0.1200.5390.6310.3930.425
SIdolomite−0.1240.9651.2310.6380.245
SIgypsum−1.692−1.875−1.686−2.209−1.851
SIhalite−7.744−7.169−6.454−7.516−7.775
SIfluorite−1.459−1.093−0.910−0.728−1.173
Hydrochemical characteristicsGeochemical reactionMainly lixiviationMainly lixiviation, cation exchange and evaporation enhancedMainly evaporation and cation exchangeMainly lixiviation and evaporation enhancedMainly lixiviation
TDS0.2–0.60.4–10.5–3.50.3–0.80.2–0.5
Characteristic element: F0.3–0.70.5–21–30.6–2.50.2–0.9
Hydrochemical typeHCO3-CaHCO3-Ca·Na (Na·Ca)HCO3-Na (Na·Mg)HCO3·Cl-NaHCO3-Ca·Na (Na·Ca)HCO3-Ca
Table 4. The over-standard of TDS and F in each zone.
Table 4. The over-standard of TDS and F in each zone.
ItemsRecharge Area
(n = 39)
Runoff Area
(n = 60)
Discharge Area
(n = 43)
Total
(n = 142)
Average of pH7.597.677.697.65
TDSNumber of over standard6121735
Percentage15.38%20.00%39.53%26.65%
FNumber of over standard8342971
Percentage20.51%56.67%67.44%40.80%

Share and Cite

MDPI and ACS Style

Li, M.; Liang, X.; Xiao, C.; Cao, Y.; Hu, S. Hydrochemical Evolution of Groundwater in a Typical Semi-Arid Groundwater Storage Basin Using a Zoning Model. Water 2019, 11, 1334. https://doi.org/10.3390/w11071334

AMA Style

Li M, Liang X, Xiao C, Cao Y, Hu S. Hydrochemical Evolution of Groundwater in a Typical Semi-Arid Groundwater Storage Basin Using a Zoning Model. Water. 2019; 11(7):1334. https://doi.org/10.3390/w11071334

Chicago/Turabian Style

Li, Mingqian, Xiujuan Liang, Changlai Xiao, Yuqing Cao, and Shuya Hu. 2019. "Hydrochemical Evolution of Groundwater in a Typical Semi-Arid Groundwater Storage Basin Using a Zoning Model" Water 11, no. 7: 1334. https://doi.org/10.3390/w11071334

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