Next Article in Journal
Application of Terrestrial Laser Scanning to Tree Trunk Bark Structure Characteristics Evaluation and Analysis of Their Effect on the Flow Resistance Coefficient
Previous Article in Journal
Water Resources Sustainability of Ulaanbaatar City, Mongolia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Groundwater Hydrochemical Zoning in Inland Plains and its Genetic Mechanisms

1
College of Hydraulic and Environmental Engineering, University of Jinan, Jinan 250022, China
2
Engineering Technology Institute for Groundwater Numerical Simulation and Contamination Control, Jinan 250022, China
3
Shandong Insititute of Geological Survey, Jinan 250014, China
4
Shandong Provincial Bureal of Geology and Mineral Resources, Jinan 250013, China
*
Author to whom correspondence should be addressed.
Water 2018, 10(6), 752; https://doi.org/10.3390/w10060752
Submission received: 3 May 2018 / Revised: 30 May 2018 / Accepted: 6 June 2018 / Published: 8 June 2018
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Pore water in inland plain areas, generally having poor water quality, contain complex hydrochemical properties. In order to examine groundwater chemical composition formation characteristics, groundwater in the Jiyang area of Lubei Plain was studied using stratified monitoring of drilling, analysis of water level and water quality, isotope analysis, ion ratio coefficient and isothermal adsorption experiments, hydrochemical characteristics, and analysis of variations in different shallow depths. Results show that: (1) Numerous hydrochemistry types are present in the diving. Along with the direction of groundwater flow, total dissolved solids (TDS) of diving in the study area generally increases and the hydrochemical type changes from the HCO3 type to the HCO3·SO4 type, Cl·HCO3 type and the Cl·SO4 type. (2) Shallow brackish water and freshwater in the horizontal direction are alternately distributed, and shallow brackish water is distributed in the area between old channels, showing sporadic spots or bands, whose hydrochemistry type is predominantly Cl·SO4-Na·Mg·Ca. (3) Affected by the sedimentary environment, hydrodynamic conditions and other factors; diving, middle brackish water and deep freshwater are vertically deposited in the study area. The dynamics of middle brackish water quality are stable due to the sedimentary environment and clay deposits. The hydrochemistry types of middle brackish water are mainly Cl·SO4-Mg·Na and SO4·Cl-Na·Mg, while the deep confined water is dominated by HCO3. (4) The optimal adsorption isotherms of Na+, Ca2+ and Mg2+ in groundwater from clay, with a thickness raging from 6–112 m, conformed to the Henry equation and the Langmuir equation. The retardation of Na+, Ca2+ and Mg2+ in groundwater differed with differing depths of the clay deposit. The trend of change in retardation strength correlates strongly with the TDS of groundwater. Groundwater in the inland plain area is affected by complicated hydrogeochemical processes; in addition, leaching, cation exchange, evaporative concentration and clay retardation have an important influence on the horizontal and vertical zonation of groundwater chemical components.

1. Introduction

The chemical composition of groundwater, especially brackish water in inland plains, is controlled by numerous factors such as meteoric, climatic, hydrological, geological background, hydrodynamic conditions and the hydrochemical field. The genetic type and evolution of groundwater are extremely complicated. Quaternary sediments in North China have been subjected to multiple transgressions and salinization effects, and they are widely distributed in shallow freshwater, mid-deep saline water and deep freshwater. The shallow pore water in the western part of North China is mainly affected by the salinization of the continent [1,2]. The relative lack of freshwater resources, insufficient utilization of inland saltwater bodies, and the deep over-exploitation of freshwater have resulted in land subsidence [3,4]. Since inland saltwater has been regarded as an unconventional and inferior water source, the current understanding of the alternating role of fresh and brackish water in inland plains is not well developed. Study on the characteristics, chemical composition and the causes of inland groundwater, which not only has important academic significance, but also important practical value of developing and utilizing groundwater resources rationally [5,6].
The chemical composition of groundwater is the product of the long-term interaction between groundwater and the surrounding environment [7,8]. During the evolution of a groundwater system, different hydrogeochemical processes result in groundwater having different chemical compositions [9], manifested as certain hydrochemical zonality [10,11], which is widely used in groundwater hydrochemistry and its genetic mechanisms [12,13]. In order to compare the formation processes of groundwater’s chemical composition, the Piper trilinear nomograph, the Shukalev classification and kurllov’s formula have been traditionally used [14,15]. Currently, statistical analysis, isotope analysis and ion proportion coefficient analysis are often used to characterize the variability in chemical compositions of groundwater [16,17]. For example, monitoring the variation of isotopic composition can distinguish their different origins [18,19]. The ion ratio coefficient can be used to analyze the hydrogeochemical effects experienced by groundwater [20,21], and to invert the formation of groundwater hydrochemistry [22,23].
Previous studies have shown that water-rock interaction, water evaporation, solute dissolution, evaporation concentration and mixing action have important effects on the chemical type of groundwater [24,25]. The aquifer medium can, to a certain extent, block the migration of ion components in groundwater, thus affecting the chemical composition [26,27]. The strength of the retardation is not only related to the nature of the media [28,29], but also to the properties of the components in the aqueous solutions [30,31]. The aquifer media can also adsorb groundwater contaminants [32,33], indicating that adsorption is one of the factors that affects its retarding properties [34,35]. However, studies on the retardation of the aquifer have predominantly focused on blocking pollutants in groundwater [36,37]. The media is less concerned by studies on conventional ion retardation in groundwater. Clay retardation can also be considered as an important hydrogeochemical process in groundwater systems. It is important to carry out relevant research for the rational development and application of groundwater [38].
Groundwater in inland plain areas mainly occur in newly accumulated loose sediments. The aquifer medium is mainly composed of clay and powdery sticky sand containing thin silt sand. The clay aquitard is an important part of the inland plains groundwater system [39,40], and it cannot ignore that clay retarding action has a significant impact on the evolution of hydrochemical composition [41]. Thus, taking the Jiyang area of Lubei Plain as a case, the zonal characteristics and genetic mechanisms of groundwater hydrochemical composition in an inland plain area are investigated using core drilling, long-term monitoring of the groundwater dynamics, stratified sampling, and routine groundwater ion tests. Ion ratio coefficient and isotope analysis are also used to provide a scientific basis for the evolution of brackish water.

2. Research Area Overview

The research area was located in Jiyang County and parts of Shanghe County, Linyi County, Qihe County (Figure 1). The study area, being hot and rainy in the summer and cold and dry in the winter, is characterized by a warm, temperate, humid continental monsoon climate. Affected by frequent diversions of the Yellow River, paleoclimate and other factors, this area has a complex sedimentary environment and contains alluvial and lacustrine [7]. Groundwater depth less than 500 m can be divided into three aquifer types. From top to bottom they are: (i) Shallow diving type (up to a depth of 40 m), having a silt and fine silt lithology; (ii) middle confined water type (60–200 m), having a sand, medium sand and fine sand lithology, shallow diving type and middle confined water type, separated by clay and silty clay aquitard; and (iii) deep confined water (roof depth of 200–300 m and a floor depth of 450–500 m), having a lithology comprising a predominantly silty clay, clay and a thin layer of fine sand (Figure 2).
The source of diving recharge in the study area is mainly atmospheric precipitation and surface water. The terrain in the study area descends from the southwest to the northeast, having a slope of 1/7000–1/8000 (Figure 1). Ground elevation is between 18–20 m. Because of the slow movement of groundwater and the small gradient, evaporation becomes the main means of discharge route [6]. The middle confined water and deep confined water were formed during the ancient geological period, and their recharge time and supply path were longer. The water quality of the middle confined water is poor, and the exploitation and utilization are less; the deep confined water is mainly excreted manually, and may be supplied by diving and over-flow of middle confined water.

3. Materials and Methods

3.1. Collection and Testing of Groundwater Samples

The hydrological survey was undertaken in the study area, resulting in the collection of 70 diving samples whose depth is less than 60m (from an area of 1076 km2 in Jiyang County), including 30 samples in the dry period, 40 samples in the wet period and 101 deep groundwater samples whose depth is greater than 200 m in and around Jiyang County. At the same time, 20 shallow brackish water samples and 40 middle brackish water samples whose depth is between 60–200 m were collected in 3.2 km2 typical saltwater areas (Figure 1). Before taking the water sample, we rinsed the washed 50 mL sample bottle with freshwater three times and then sealed the water. 14C, 3H, CFC-11, CFC-12 samples used a 60 mL glass bottle. Then, before sampling, the glass bottle was rinsed and stored in sealed conditions to avoid exposure to air. 14C, 3H, CFC-11, CFC-12 samples were sent to the Beta Analytic Radiocarbon Dating Laboratory to complete the test [7]. K+ and Na+ were measured using a flame photometer and Ca2+ and Mg2+ were measured using a Shimadzu-AA7000 atomic absorption spectrophotometer. SO42− was determined using the EDTA-Mg titration method. CO32−, HCO3 was determined using the acid-base titration method with methyl orange as the indicator. Titration of Cl used a potassium chromate indicator. Control samples and the relative balance of anion balance was within ±5%.

3.2. Dynamic Monitoring of Groundwater

In order to collect soil samples from different depths for adsorption analysis, coring techniques were used. After coring, six stratified observation holes were created (WX01-06; Table 1). Groundwater dynamic monitoring was undertaken over a two-year period to obtain the aquifer level and water quality data at different depths.
(1)
Stratified monitoring of the groundwater level was undertaken at different depths in the observation holes on a weekly basis.
(2)
Stratified monitoring of groundwater quality characteristics was undertaken at different depths using an Aqualtroll 600-type multi-parameter water quality detector. Conductivity (EC) changes from the top to the bottom of the profiles were recorded using a test interval of 1 m. Test accuracy was ±0.5% + 1 μs/cm. Weekly monitoring commenced in November, 2016.

3.3. Adsorption Isothermal Experiment

3.3.1. Experimental Procedure

To examine the adsorption characteristics of clays at different depths to Na+, Ca2+ and Mg2+ in groundwater, the indoor adsorption isothermal batch method [39,40] was used. For analysis, soil samples (3000 g) were taken from different depths. All debris was removed and the samples were air dried. Samples were then sieved through a 20 mesh and oven dried. Samples were ground (using a mortar) and passed through a 100 mesh sieve and retained for the adsorption experiment.
A 30 g sub-sample of clay from different depth aquifers, obtained by drilling and which represents different aquifer media, was placed into a 250 mL Erlenmeyer flask. A 150 mL volume different concentration of water samples were added to the flasks, respectively, after which the flask was placed into water bath thermostatic shaker and shaken for 3 h until adsorption reached saturation. Samples were centrifuged at 4000 r·min−1 for 1 min and filtered through a 0.45 μm microporous membrane. 2 mL of the filtrate was then used to test for Na+, Ca2+ and Mg2+ concentrations (using atomic absorption spectrophotometer-AA7000). The experiment was repeated three times to obtain effective data of 307 groups.

3.3.2. Calculation of the Retardation Coefficient

The adsorption capacity of clay on groundwater Na+, Ca2+ and Mg2+ is indicated by adsorption capacity Qe, and the adsorption efficiency is expressed by ε, the formula are (1) and (2):
Q e = ( C 0 C e ) × V m
ε = ( C 0 C e ) C 0 × 100 %
The Temkin equation, Langmuir equation, Freundlich equation and Henry equation were used to fit the adsorption isotherms of Na+, Ca2+ and Mg2+ in the groundwater. The Temkin, Freundlich and Langmuir equations belong to the monolayer adsorption. The Henry equation is linear. Equation is (3)–(6) [42,43]:
Q e = k × lg C e + a
Q e = K L × Q m × C e 1 + K L × C e
Q e = K F × C e n
Q e = K d × C e
Clay effective porosity, ne, was calculated according to Equation (7) and the retardation coefficient, Rd, was calculated according to Equation (8) [44,45]:
n e = e 0 1 + e 0
R d = 1 + 1 n n K d
where C0 is the initial concentration of conventional ions in groundwater, mg·L−1; Ce is the equilibrium concentration of adsorption, mg·L−1; Qe is the adsorption capacity, mg·kg−1; V is the volume of the test water sample, L; m is the volume of soil samples, kg; K is the parameter of the Temkin equation and describes the conventional ionic strength of soil in the sample; Qm and KL are the maximum adsorption capacity and adsorption energy-related parameters for the Langmuir adsorption equation, mg·kg−1 and L·mg−1, respectively; KF is the Freundlich adsorption equation; Kd is the partition coefficient of the Henry equation and describes the ratio of the conventional ion concentration in the clay to the concentration in the groundwater at adsorption equilibrium, L·kg−1; e0 is the porosity ratio, as determined by the geotechnical test; ne is the effective porosity; and Rd is the retardation coefficient which characterizes the medium’s ability to block solutes.

4. Results and Discussion

4.1. Groundwater Chemical Composition and Zoning Characteristics

4.1.1. Chemical Characteristics of Shallow Groundwater

Analysis of 70 groundwater samples in Jiyang County indicated that TDS of the shallow groundwater was between 0.55 g·L−1 and 13.10 g·L−1 (Table 2). Apart from individual water samples whose depths less than 10 m deep, there are almost no changes in the type of hydrochemistry and the content of various ions; although the total dissolved solids of water samples whose depths less than 10 m deep are larger during the dry season and smaller during the rainy season, there is no change in the type of hydrochemistry. Freshwater samples with TDS < 2 g·L−1 accounted for 74.3% of total samples, mainly accounting for HCO3-type water, such as the HCO3-Na (Mg·Na) type and the HCO3 (HCO3·Cl)-Na·Mg type. Brackish water with a TDS of 2–3 g·L−1 accounted for 15.7% of total samples. The major hydrochemical types were Cl·HCO3 (HCO3·Cl·SO4, HCO3·SO4·Cl)-Na·Mg type, HCO3-Na (Na·Mg) type and HCO3·SO4-Na type. Brackish water with a TDS of 3–5 g·L−1 accounted for 8.6% of water samples, with Cl·HCO3-Mg·Na·Ca as the main type and Cl·SO4-Na (Na·Mg) as the secondary type. Saltwater with a TDS > 5 g·L−1 contained Cl·SO4-Na·Mg or Cl·SO4-Mg·Na types. The shallow groundwater HCO3 coefficient of variation was less and the Cl, SO42− and Na+ coefficient of variation was larger, indicting complex chemical types.

4.1.2. Hydrochemical Types and Horizontal Zoning Characteristics

Shallow groundwater hydrochemical types reflect a certain level of zonation.
(1)
Groundwater flow from the Southwest to the Northeast. The longer migration route, the greater the TDS of the groundwater. In Section I, Fengjiajing Village in Taiping Town (a located in the Southwest)—Hujiacun (b located in the Central)—Tanjia Village in Renfeng Town (c located in the Northeast) (Figure 1) of Taiping Town, groundwater was changed from freshwater with a TDS < 1 g·L−1 to brackish water with a TDS > 3 g·L−1. Na+, Ca2+, Mg2+, Cl, HCO3, SO42− and TDS all recorded a gradual increase. The hydrochemistry types were HCO3-Na·Mg ⟶ HCO3·Cl ·SO4-Na·Mg ⟶ Cl·HCO3-Na·Mg (Figure 3).
(2)
Brackish water and freshwater show an alternate distribution in local areas. Groundwater with a high salinity had a sporadic or irregular strip-shaped distributing interfluve between old channels (Figure 4). There are numerous seasonal streams and diversion channels in the Jiyang surface water system. Freshwater with a TDS < 2 g·L−1 is mainly distributed along the Yellow River, along the Tuhai River, near the Dasi River and around the Heishi River. The main hydrochemistry type is HCO3 type. The brackish water-saltwater area is mainly distributed between river, showing north-easterly. This area recorded a general TDS > 3 g·L−1, with results occasionally being >5 g·L−1. TDS samples in the Old Tang Temple were recorded to be >10 g·L−1. In the direction of Section II, the hydrochemistry types were HCO3·Cl-Mg·Na·Ca, Cl·SO4—Mg·Na·Ca, HCO3·Cl-Mg·Na and HCO3·Cl-Mg·Na, Cl·SO4—Na·Mg type, formed a freshwater-saltwater-freshwater-saltwater pattern.

4.2. Groundwater Chemical Composition of the Vertical Zone Characteristics

4.2.1. Changes in the Chemical Composition of Saltwater with Depth

Taking Sungeng town as a case-study, the distribution area of shallow saltwater was 1.41 km2, pH ranged from 6.96–7.11, TDS > 10 g·L−1, and the order of cations (by average concentration) was: Na+ (2029.10 mg·L−1) > Ca2+ (798.65 mg·L−1) > Mg2+ (791.81 mg·L−1). The order of anions (by average concentration) was: Cl (3791.59 mg·L−1) > SO42− (3251.69 mg·L−1) > HCO3 (886.96 mg·L−1) (Table 3). Due to differences in burial conditions, the pH value of the middle saltwater ranged from 7.01–8.16, this being weakly alkaline. Average results for cations and anions were: Na+ (2108.01 mg·L−1) > Mg2+ (714.79 mg·L−1) > Ca2+ (670.10 mg·L−1) and SO42− (3266.29 mg·L−1) > Cl (2865.51 mg·L−1) > HCO3 (244.71 mg·L−1), respectively. Due to ground subsidence caused by deep groundwater exploitation, the deep production wells in Jiyang have been closed. Sampling and analysis in the surrounding counties indicates that TDS of deep pressure water is less than 3 g·L−1, and the highest cation concentration (Na+) was 345.29 mg·L−1. Results for HCO3, Cl and SO42− were 429.08 mg·L−1, 247.64 mg·L−1 and 205.99 mg·L−1, respectively. Therefore, vertical distribution of groundwater have obvious zonation characteristics in this area.

4.2.2. Groundwater Conductivity with Depth

Long-term monitoring results (2016–2017) indicate that groundwater from similar depths have similar conductivity values, and conductivity depths greater than 10 m are relatively stable. The variation amplitude between samples at the same depth was around 1.5 ms/cm (Figure 5a). There is a change in vertical conductivity at different depths, (Figure 5b). The rainfall recharge [7] had a significant effect on conductivity depth less than 10 m. With the increase of depth, the effect of low conductivity rainwater mixing was weaker, and with vertical conductivity gradually increasing from 8088 μs/cm, the rate of increase was 254 μs/cm per meter. Conductivity between 11–28 m and 28–45 m were relatively stable. At 45 m, conductivity abruptly increased to 20,699 μs/cm before decreasing. Conductivity from 70–76 m was again relatively stable (18,621 μs/cm). From 82–89 m, conductivity increased slightly with depth (to 19,758 μs/cm). After 100 m a final sudden decline was recorded, with conductivity stabilizing around 12,019 μs/cm.
According to the core observation, gypsum crystals (maximum crystal length of 5 cm) are occasionally present in clays and silty clays at a depth of 45 to 90 m, and the maximum crystal length is 5 cm. The graded and abrupt change of groundwater conductance in the vertical direction reflects that the formation of clay in different sedimentary environments has a significant influence on vertical zoning of groundwater conductivity.

4.2.3. Groundwater Hydrochemical Vertical Zoning Characteristics

According to the Shukalev classification, shallow saltwater was dominated by Mg2+, Na+ and Ca2+ salts and Cl and SO42− anions. The coefficients of variation of Mg2+, Na+ and Cl are small. The hydrochemistry type of the saltwater was Cl·SO4-Mg·Na·Ca type and the individual is Cl·SO4-Na·Mg·Ca type. Cations in the middle brackish water were mainly Na+ and Mg2+, with the main anions being SO42− and Cl. The coefficient of variation of SO42− was larger and the remainder of the conventional ion coefficient of variation was less. The main hydrochemistry types were SO4·Cl-Na·Mg and Cl·SO4-Mg·Na. Hydrochemistry types of deep confined water were the HCO3 type, such as HCO3-Na, HCO3·Cl (HCO3·SO4·Cl, HCO3·Cl·SO4)-Na and Cl·HCO3·SO4-Na, HCO3-Ca·Mg. The coefficient of variation for HCO3 and Na+ was small. The vertical distribution of groundwater (depth less than 500 m) can be divided into: shallow freshwater—middle saltwater—deep freshwater; shallow brackish water—middle saltwater—deep freshwater; and shallow saltwater—middle saltwater—deep freshwater. Shallow saltwater in Sungeng town often has no clear boundary with the mid-level saltwater, and the groundwater hydrochemical type is Cl·SO4 or SO4·Cl. Hydrological and geochemical characteristics of confined water in Yinchuan area are similar to those in the study area [18].

4.3. Inland Plain Groundwater Mechanisms

4.3.1. Groundwater Characteristics

The Piper trilinear nomograph (Figure 6) shows that there is a difference in the hydrochemical characteristics between shallow and deep confined water. Shallow groundwater receiving recharge of atmospheric precipitation is mainly dominated by leaching and the filtration effect. Saltwater samples are located in the I and IV districts of the Piper trilinear nomograph, reflecting evaporation, concentration and cation exchange during geological history. Deep confined water samples are mainly located in the II district, reflecting cation exchange, while shallow water samples in the lower half of the II district are only 2.8% (Figure 6).

4.3.2. Effect of Clay Retarding Action on Chemical Zoning in Water

(1) Isothermal adsorption equation
Clay plays a controlling role in the hydrodynamic field of a groundwater system. Clay formed in different depositional environments also has an impact on the chemical composition of groundwater. Adsorption experiments can therefore be used to investigate the effect of clay retardation on the chemical composition of water. Na+, Ca2+ and Mg2+ results from adsorption experiments were fitted with the Temkin, Langmuir, Freundlich and Henry adsorption isotherm equations. The optimal isothermal adsorption curves for Na+, Ca2+ and Mg2+ in groundwater with different depths clay were well fitted with the Henry and Langmuir equations (Table 4). Fitting results of the Temkin and Freundlich equations were not as suitable. Clay adsorption of Na+, Ca2+ and Mg2+ in groundwater has a monolayer or linear adsorption process [36].
The fitting results of the Langmuir adsorption isotherm equation show that the maximum adsorption capacity of clay to groundwater is different. For example, the order of the conventional ion adsorption capacity of groundwater with a depth of 6–6.2 m was Ca2+ > Mg2+ > Na+. The order of routine ion adsorption of clay to groundwater at a depth of 51–51.2 m was: Na+ > Mg2+ > Ca2+, and at a depth of 71–71.2 m the order of conventional ion adsorption of clay to groundwater was: Mg2+ > Ca2+ > Na+. At depths of 84–84.2 m, 102–102.2 m and 112–112.2 m, ion adsorption was in the order of: Mg2+ > Na+ > Ca2+ (Table 4). In addition, Henry equation Kd is the linear distribution coefficient [37]. The larger the Kd value, the stronger the clay adsorption capacity. The order of clay Kd values of the same depth fitted by the Henry model were in accordance with the order of the maximum adsorption capacity from the Langmuir equation. The maximum adsorption capacity of Na+, Ca2+ and Mg2+ is different in the same depth clays. The maximum adsorption capacity of the same ions is also different in various depth clays (Table 4).
(2) Clay blocking effect
The blocking strength of clay in Na+, Ca2+ and Mg2+ in groundwater was expressed as the retardation coefficient Rd. The retardation coefficients of Na+, Ca2+ and Mg2+ in groundwater, calculated in Equations (7) and (8), show that the order of the conventional ion blockages of groundwater with depths of 6 m, 84 m, 97 m and 102 m was Mg2+ > Na+ > Ca2+. Groundwater with depths of 48 m and 112 m had a conventional ion block strength order of: Na+ > Mg2+ > Ca2+. The order of the conventional ion block strength of groundwater with a depth of 51.2 m was Na+ > Ca2+ > Mg2+, and the order of the conventional ionic blockages of clay to groundwater at the depth of 71.2 m was: Ca2+ > Na+ > Mg2+ (Table 5). Results from these equations indicated that clays at the same depth have different retardation strengths for different conventional ions in groundwater. The blocking strength of clay at different depths for Ca2+, Mg2+ and Na+ in groundwater varied alternately, reflecting that the sedimentary environment has a certain impact on the quality of the groundwater [40].
By comparing the clay retardation strength of different depths with the corresponding chemical composition of the groundwater, the correlation between the total retardation of Na+, Ca2+, Mg2+ and groundwater TDS from corresponding depths was calculated (Figure 7). The relationship between total retardation and groundwater TDS was calculated using the following equation:
y = 4818.3x − 9496 (R2 = 0.9167)
where, x was total retardation strength (a dimensionless value); and y was total dissolved solids, mg·L−1.
The clay block in the hinterland plain area has some control over the vertical ion composition of groundwater, and that clay retardation plays a certain role in promoting Na·Mg; Mg·Na saltwater (Figure 7).

4.3.3. Groundwater Hydrogeochemistry Reflected in the Ion Ratio Coefficient

(1) The role of leaching processes in shallow freshwater
γNa/γCl can be used to indicate the hydrogeochemical parameters of Na+ enrichment in groundwater. Near the southern part of the study area and regions near the river, γNa/γCl from 1 to 2 accounts for 55.7% of the total, indicating that the groundwater is affected by the leaching effect [7,26]. The results of the study are similar to those in Northwestern Shandong [7].
(2) Strong saltwater cation exchange
As well as determining the intensity of ion exchange reactions in water chemistry, the γCa/γNa and γMg/γNa coefficients can also characterize the evolution of water quality and mineralization strength [7]. γNa/γCl of the shallow brackish water in the study area ranged from 0.5–0.6; γNa/γCl in the middle confined aquifer ranged from 0.45–0.8; and the ration in the deep confined water was greater than 3. The higher γNa/γCl corresponds to the lower γMg and γCa (Figure 8). γ Mg/γNa > 1, γCa/γNa < 1, and Mg2+ in the saltwater zone were dominant, followed by Na+ and Ca2+. Results indicate that the formation mechanism of deep confined freshwater, middle confined saltwater and shallow saltwater differed. The cation exchange activity in the saltwater area was stronger than that in the freshwater; Ca2+ and Mg2+ in middle and deep groundwater may have undergone cation exchange with Na+, which may have certain reference significance for the evolution of groundwater in other regions.
(3) Slow saltwater groundwater runoff rates
The γCa/γCl ratio can be used to describe hydrodynamic conditions [11]. When TDS < 1.2 g·L−1, γCa/γCl > 2 indicates favorable shallow groundwater runoff conditions. Saltwater TDS values in Sungeng town were >3 g·L−1, and the γCa/γCl ratio ranged from 0.2–0.4. Cl enrichment indicated that its hydrodynamic conditions were very poor, this being the conditions consistent with a small hydraulic gradient and slow runoff in the study area [7]. The coefficients of γHCO3/γSO4 and γHCO3/γCl in the saltwater area were both less than 1, indicating that the saltwater area is affected by evaporite dissolution [26], which is helpful to ionic composition analysis in other regions. The average value of the γHCO3/γCl coefficient of deep groundwater is 2.24, indicating that there is a difference in the anionic components between the deep confined water, middle confined water and the diving.

4.3.4. The Mechanism of Formation of Hydrochemistry

(1) The formation of shallow brackish water
The study area located on the Yellow River plain undergoes atmospheric precipitation infiltration and lateral recharge from the Yellow River. Due to the high terrain and water level, the hydraulic gradient is about 1/2000–1/3000, a larger hydraulic gradient leads to faster runoff, thus leaching and filtration are intense. The hydrochemical type in this area is HCO3-Ca type, which is good water quality. The plain away from the Yellow River is characterized by a complex sedimentary environment which has an increase in clay content. Along the direction of the groundwater, flow velocity decreases and the hydraulic gradient is less than 1/5000 (Figure 1). The concentrations of Na+, Mg2+, Cl and SO42− in the groundwater increased and the hydrochemical types changed from HCO3-Ca to HCO3·Cl-Na·Mg, HCO ·Cl·SO4-Na·Mg and Cl·HCO3-Na·Mg.
Diving in the Jiyang area shows a pattern of horizontal distribution of fresh and brackish water (Figure 4). The old channel rich in water is better and has good water quality. The interfluve between two old channels, having sediments consisting of mainly fine particles, distributed shallow saltwater showing a sporadic shape. This area undergoes atmospheric precipitation recharge and shallow freshwater recharge. The sediments in the interfluve are fine. The groundwater table has a shallow depth, poor runoff conditions and intense evaporation and concentration to form shallow brackish water. Prior to the 1960s, the groundwater head in the middle and deep levels was higher than the ground elevation. In the low-energy region of the regional flow system and the local flow system, brackish water in the mid-level moves upwards due to evaporation, thus forming a modern, shallow saltwater. Although the mixing effect influences the isotope test, it does not affect the analysis result. The shallow brackish water of CFC-11 and CFC-12 is therefore only 30–50 years old.
(2) The formation of the middle saltwater
Drilling data shows that the river and lake phase sediments (40–110 m) in this area mainly comprise clay, silty clay, clip silt, thin silt sand, and clay containing gypsum greater than 5 cm in diameter. 14C data indicate that the isotopic age of clay at the depth of 43–43.2 m is 23.79–24.085 ka BP, and that of clay at the depth of 90 m is 31.56–31.265 ka BP. The age of brackish water depth of 50 m measured by 14C was 6.43–6.55 ka BP. The 14C datings in the aquifer of figuig indicate that the age of groundwater ranges between 2495 and 10,696 years [14]. The δ18O content of this water was −8.4‰. Middle confined saltwater from a depth of 99–110 m was dated to 39–36 ka BP, this water having an δ18O content of −9.0~−9.2‰. This time period coincided with the Dali ice age in the Quaternary, resulting in a dry climate [46]. Evaporation and concentration led to an increase in salinity, thus the paleoclimate and sedimentary environment had an important influence on the formation of the middle saltwater.
(3) The dynamic stability of saltwater quality
Due to climatic changes since the Late Pleistocene, alternating glacial and interglacial episodes occurred in the study area. With the exploration depth being 321 m, sedimentary cycles resulting in fine to coarse grained deposits appeared several times (Figure 2). The sand layer was thin and discontinuous, whilst clay and silty clay were abundant. Long-term observation data showed that the water level at different depths had different dynamic characteristics. The dynamics of the shallow groundwater level was mainly controlled by meteorology and hydrology, whose level decreased in the dry season and increased in the wet season. The elevation of the water level was between 20.5–21.5 m, indicating seasonal variation. Middle saltwater had little development and utilization of value. The water level was relatively stable, with a depth of 18–18.5 m. Freshwater at a depth of 300–320 m had a conductivity of 1161–1308 um/cm, a 14C age of 41.74–43.5 ka BP, and δ18O content of −9.8~−9.9‰. In the late 1970s, the deep confined water level was about 26 m. Due to poor supply conditions of deep confined water, the water level continued to slowly decline due to over-exploitation (Figure 9), and the rate of decrease was about 0.1–0.5 m/a [47]. At present, the water level is 12.9–13.3 m, this being lower than the middle saltwater level and the diving level. It illustrates that the stable distribution of cohesive soil is a good regional weak aquitard.
Before the 1970s, the groundwater level in the middle and deep layers was higher than shallow groundwater. Over the last 30 years, the groundwater level in the mid-deep layers was lower than shallow groundwater due to the diversion of irrigation water from the Yellow River; however, water quality of the middle saltwater was still stable. Combining groundwater level, dynamic characteristics of water quality, groundwater age and sedimentary environment, groundwater in the study area (500 m) can be divided into three systems: a shallow partial flow system; an intermediate saltwater flow system; and a deep freshwater flow system (Figure 10). The stable distribution of clay aquitards weakens the vertical hydraulic exchange of saltwater. The combination of the sedimentary environment, clay retardation and hydrogeochemical affect the saltwater quality, and the evolution of saltwater water quality is extremely slow. The alternating flow between the three mobile systems is weak, thus water quality and the water level of the middle brackish water is dynamic and stable, which hass implications for the study and utilization of other similar inland saline waters. These factors have resulted in local populations being unable to exploit this water source for thousands of years, and people have been taking water from 3 km away for domestic use.

5. Conclusions

(1)
The influence of hydrodynamic conditions, sedimentary environment and topography has resulted in a complex groundwater state in the inland plain. TDS increases along the groundwater flow and with the horizontal distribution of freshwater and brackish water (saltwater), showing alternate distribution. Shallow saltwater is distributed in the interfluve between old channels as sporadic or irregular bands.
(2)
The optimal isothermal adsorption curves of Na+, Ca2+ and Mg2+ in groundwater over different depths accord with Henry’s equation and Langmuir equation, which is the monolayer chemical adsorption or linear adsorption process. The 6–112 m clay exhibits a large-small-large trend of Na+ and Ca2+ retardation in groundwater. The trend of Mg2+ was small-large-small, and the block strength of Na+ and Mg2+ in groundwater was higher than that of Ca2+, facilitating a Na-Mg or Mg-Na saltwater formation. The correlation between the block strength of clay for Na+, Ca2+ and Mg2+, and groundwater TDS was strong.
(3)
Water quality of middle saltwater was relatively stable. The main hydrochemical types were Cl·SO4-Na·Mg and SO4·Cl-Mg·Na, which belong to the intermediate circulation system. The dynamic stability of saltwater quality was affected by the stability of clay distribution, its weak permeability and blocking effect, and the sedimentary environment.
(4)
Groundwater in the study area is not only affected by complex hydrogeochemistry, such as the leaching effect, cation exchange and evaporation, it is also closely related to the retardation effect of clay. The inland plains area has horizontal and vertical hydrogeochemical zonation characteristics.

Author Contributions

L.X. and X.H. conceived and designed the experiments; L.H. performed the experiments; L.Y. and H.Z. analyzed the data; G.C. and J.X. contributed reagents/materials/analysis tools; L.X. and X.H. wrote the paper.

Funding

This research is supported by the National Natural Science Foundation of China (41772257, 41472216, 41727225).

Conflicts of Interest

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

References

  1. Mehta, S.; Fryar, A.E.; Brady, R.M.; Morin, R.H. Modeling regional salinization of the ogallala aquifer, Southern High Plains, TX, USA. J. Hydrol. 2000, 238, 44–64. [Google Scholar] [CrossRef]
  2. Weisbrod, N.; Dragila, M.I. Potential impact of convective fracture venting on salt-crust buildup and ground-water salinization in arid environments. J. Arid Environ. 2006, 65, 386–399. [Google Scholar] [CrossRef]
  3. Hunter, R.; Day, J.; Shaffer, G.; Lane, R.; Englande, A.; Reimers, R.; Kandalepas, D.; Wood, W.; Day, J.; Hillmann, E. Restoration and management of a degraded baldcypress swamp and freshwater marsh in Coastal Louisiana. Water 2016, 8, 71. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Xu, M.; Li, X.; Qi, J.; Zhang, Q.; Guo, J.; Yu, L.; Zhao, R. Hydrochemical characteristics and multivariate statistical analysis of natural water system: A case study in Kangding County, Southwestern China. Water 2017, 10, 80. [Google Scholar] [CrossRef]
  5. Dedzo, M.G.; Tsozué, D.; Mimba, M.E.; Teddy, F.; Nembungwe, R.M.; Linida, S. Importance of rocks and their weathering products on groundwater quality in Central-East Cameroon. Hydrology 2017, 4, 23. [Google Scholar] [CrossRef]
  6. An, L.S.; Zhao, Q.S.; Ye, S.Y.; Liu, G.Q.; Ding, X.G. Hydrochemical characteristics and formation mechanism of shallow groundwater in the Yellow River Delta. Environ. Sci. 2012, 33, 370–378. [Google Scholar]
  7. Xing, L.; Zhang, F.; Changsuo, L.I.; Yang, L.; Wang, L. Hydro-chemical characteristics of shallow groundwater in the North Shandong Plain. J. Irrig. Drain. 2015, 17, 50–56. [Google Scholar]
  8. Wang, X.S.; Wan, L.; Jiang, X.W.; Li, H.; Zhou, Y.; Wang, J.; Ji, X. Identifying three-dimensional nested groundwater flow systems in a Tóthian Basin. Adv. Water Resour. 2017, 108, 139–156. [Google Scholar] [CrossRef]
  9. Li, Q.; Zhou, J.; Zhou, Y.; Bai, C.; Tao, H.; Jia, R.; Ji, Y.; Yang, G. Variation of groundwater hydrochemical characteristics in the plain area of the Tarim Basin, Xinjiang Region, China. Environ. Earth Sci. 2014, 72, 4249–4263. [Google Scholar] [CrossRef]
  10. Yuan, J.; Xu, F.; Deng, G.; Tang, Y.; Li, P. Hydrogeochemistry of shallow groundwater in a karst aquifer system of Bijie City, Guizhou Province. Water 2017, 9, 625. [Google Scholar] [CrossRef]
  11. Najib, S.; Fadili, A.; Mehdi, K.; Riss, J.; Makan, A. Contribution of hydrochemical and geoelectrical approaches to investigate salinization process and seawater intrusion in the coastal aquifers of Chaouia, Morocco. J. Contam. Hydrol. 2017, 198, 24–36. [Google Scholar] [CrossRef] [PubMed]
  12. Liao, L.; Jean, J.S.; Chakraborty, S.; Lee, M.K.; Kar, S.; Yang, H.J.; Li, Z. Hydrogeochemistry of groundwater and arsenic adsorption characteristics of subsurface sediments in an Alluvial Plain, SW Taiwan. Sustainability 2016, 8, 1305. [Google Scholar] [CrossRef]
  13. Zhang, W.; Cao, Y.; Jiang, J.; Wang, W.; Duan, L. Numerical simulation of the formation of groundwater chemical components in the Mayukou area, Datong, China. In Proceedings of the International Symposium on Water-Rock Interaction, Kunming, China, 31 July–5 August 2007. [Google Scholar]
  14. Jilali, A.; Fagel, N.; Amar, M.; Abbas, M.; Zarhloule, Y. Hydrogeochemical processes constrained by multivariate statistical methods and isotopic evidence of groundwater recharge in the aquifer of Figuig, Eastern High Atlas of Morocco. Arab. J. Geosci. 2016, 9, 1–11. [Google Scholar] [CrossRef]
  15. Cao, W.; Yang, H.; Liu, C.; Li, Y.; Bai, H. Hydrogeochemical characteristics and evolution of the aquifer systems of Gonghe Basin, Northern China. Geosci. Front. 2018, 9, 907–916. [Google Scholar] [CrossRef]
  16. Awaleh, M.O.; Baudron, P.; Soubaneh, Y.D.; Boschetti, T.; Hoch, F.B.; Egueh, N.M.; Mohamed, J.; Dabar, O.A.; Masse-Dufresne, J.; Gassani, J. Recharge, groundwater flow pattern and contamination processes in an arid volcanic area: Insights from isotopic and geochemical tracers (bara aquifer system, Republic of Djibouti). J. Geochem. Explor. 2017, 175, 82–98. [Google Scholar] [CrossRef]
  17. Swaid, F.A.; Issa, A.A. Deep groundwater quality and its suitability for drinking and agricultural purposes, Wadibaye, Middle Region, Libya. In Sustainable Irrigation and Drainage V: Management, Technologies and Policies; WIT Press: Southampton, UK, 2014; pp. 183–196. [Google Scholar]
  18. Zhihong, L.I.; Fusheng, H.U.; Zhou, W.; Zhao, Z.; Wang, G.; Shi, Z. Hydrogeochemical characteristics and controlling factors of confined water in Yinchuan. Hydrogeol. Eng. Geol. 2017, 44, 31–39. [Google Scholar]
  19. Mustafa, O.; Merkel, B.; Weise, S. Assessment of hydrogeochemistry and environmental isotopes in karst springs of makook anticline, Kurdistan Region, Iraq. Hydrology 2015, 2, 48–68. [Google Scholar] [CrossRef]
  20. Sun, Y.Q.; Qian, H.; Wu, X.-H. Hydrogeochemical modeling of minerals in groundwater in Yinchuan Region. J. Saf. Environ. 2006, 6, 32. [Google Scholar]
  21. Banaś, K. The hydrochemistry of peatland lakes as a result of the morphological characteristics of their basins. Oceanol. Hydrobiol. Stud. 2013, 42, 28–39. [Google Scholar] [CrossRef] [Green Version]
  22. Kolpakova, M.N.; Naymushina, O.S. Zonal features of bog and lake water chemistry along a transect from boreal to arid landscapes in the south of Western Siberia, Russia. IOP Conf. Ser. Earth Environ. Sci. 2016, 44, 022025. [Google Scholar] [CrossRef]
  23. Chen, Q.; Mao, Y.; Cheng, L.I. Hydrochemistry characteristics of shallow groundwater in redbed of Luzhou, Sichuan. China Meas. Test 2013, 39, 42–45. [Google Scholar]
  24. Li, X.; Wu, H.; Qian, H.; Gao, Y. Groundwater chemistry regulated by hydrochemical processes and geological structures: A case study in Tongchuan, China. Water 2018, 10, 338. [Google Scholar] [CrossRef]
  25. Jiang, X.W.; Wan, L.; Wang, X.S.; Li, H.-L. Distribution of groundwater age in drainage basins. Hydrogeol. Eng. Geol. 2012, 4, 003. [Google Scholar]
  26. Feng, X.; Zhang, Y.Z. An analysis of the ions ratio coefficients of groundwater in Shenzhou area. China Rural Water Hydropower 2014, 4, 18–21. [Google Scholar]
  27. Zhang, G.X.; Wei, D.; Yan, H.E.; Salama, R. Hydrochemical characteristics and evolution laws of groundwater in Songnen Plain, Northeast China. Adv. Water Sci. 2006, 17, 20–28. [Google Scholar]
  28. Xing, W.; Ma, R.; Sun, Z.; Hu, Y.; Ge, M.; Chang, Q. Hydrochemical characteristics and water quality assessment of groundwater in the Dunhuang Basin, Northwestern China. Geol. Sci. Technol. Inf. 2016, 35, 196–202. [Google Scholar]
  29. Loni, O.A.; Zaidi, F.K.; Alhumimidi, M.S.; Alharbi, O.A.; Hussein, M.T.; Dafalla, M.; Alyousef, K.A.; Kassem, O.M.K. Evaluation of groundwater quality in an evaporation dominant arid environment; A case study from al asyah area in Saudi Arabia. Arab. J. Geosci. 2015, 8, 6237–6247. [Google Scholar] [CrossRef]
  30. Nakagawa, K.; Amano, H.; Kawamura, A.; Berndtsson, R. Classification of groundwater chemistry in Shimabara, using self-organizing maps. Hydrol. Res. 2017, 48, 840–850. [Google Scholar] [CrossRef]
  31. Sherwood, J.D.; Craster, B. Transport of water and ions through a clay membrane. J. Colloid Interface Sci. 2000, 230, 349–358. [Google Scholar] [CrossRef] [PubMed]
  32. Yamei, L.I.; Cheng, J.; Minmin, L.I.; Cui, L. Experimental study on migration and transformation in cohesive soil of salt water. Water Resour. Prot. 2015, 31, 88–93. [Google Scholar]
  33. Ishiguro, M.; Matsuura, T.; Detellier, C. Reverse osmosis separation for a montmorillonite membrane. J. Membr. Sci. 1995, 107, 87–92. [Google Scholar] [CrossRef]
  34. Saber, M.; Abdelshafy, M.; Faragallah, E.A.A.; Abd-Alla, M.H. Hydrochemical and bacteriological analyses of groundwater and its suitability for drinking and agricultural uses at Manfalut District, Assuit, Egypt. Arab. J. Geosci. 2014, 7, 4593–4613. [Google Scholar] [CrossRef]
  35. Yin, J.; Deng, C.; Yu, Z.; Wang, X.; Xu, G. Effective removal of lead ions from aqueous solution using nano illite/smectite clay: Isotherm, kinetic, and thermodynamic modeling of adsorption. Water 2018, 10, 210. [Google Scholar] [CrossRef]
  36. Zhang, Y.; Fu, C.; Mao, L.; Gong, X.L.; Li, X.-Q. Hydrochemical characteristics and formation mechanism of the groundwater in Yancheng, Jiangsu Province. Resour. Environ. Yangtze Basin 2017, 26, 598–605. [Google Scholar]
  37. Xu, W.; Zhang, J.; Liu, Y. Degradation of organic matter and ammonia nitrogen in constructed rapid infiltration system. Fresenius Environ. Bull. 2011, 20, 1685–1690. [Google Scholar]
  38. Maleki, S.; Karimi-Jashni, A. Effect of ball milling process on the structure of local clay and its adsorption performance for Ni(II) removal. Appl. Clay Sci. 2017, 137, 213–224. [Google Scholar] [CrossRef]
  39. Yechieli, Y.; Wood, W.W. Hydrogeologic processes in saline systems: Playas, sabkhas, and saline lakes. Earth Sci. Rev. 2002, 58, 343–365. [Google Scholar] [CrossRef]
  40. Yang, Y.; Xing, L.T.; Zhang, W.J.; Xiang, H.; Song, G.Z.; Chi, G.Y.; Hou, X.Y. The retardation of clay in shallow groundwater system of inland plains. Sci. Technol. Eng. 2017, 17, 50–56. [Google Scholar]
  41. Huang, L.; Wang, L.; Zhang, Y.; Xing, L.; Hao, Q.; Xiao, Y.; Yang, L.; Zhu, H. Identification of groundwater pollution sources by a SCE-UA algorithm-based simulation/optimization model. Water 2018, 10, 193. [Google Scholar] [CrossRef]
  42. Hui, B.; Ping, H. Experimental study of water and salt redistribution of clay soil in an opening system with constant temperature. Environ. Geol. 2008, 55, 717–721. [Google Scholar]
  43. Zhang, Z.; Qin, Y.; Wang, G.X.; Fu, X.H. Numerical description of coalbed methane desorption stages based on isothermal adsorption experiment. Sci. China Earth Sci. 2013, 56, 1029–1036. [Google Scholar] [CrossRef]
  44. Vimonses, V.; Lei, S.; Bo, J.; Chow, C.W.K.; Saint, C. Kinetic study and equilibrium isotherm analysis of Congo red adsorption by clay materials. Chem. Eng. J. 2009, 148, 354–364. [Google Scholar] [CrossRef]
  45. Nelson, C.; Mildenhall, D.; Todd, A.; Pocknall, D. Subsurface stratigraphy, paleoenvironments, palynology, and depositional history of the late neogene tauranga group at ohinewai, lower waikato lowland, South Auckland, New Zealand. N. Z. J. Geol. Geophys. 1988, 31, 21–40. [Google Scholar] [CrossRef]
  46. Sugita, F.; Gillham, R.W. Pore scale variation in retardation factor as a cause of nonideal reactive breakthrough curves: 3. Column investigations. Water Resour. Res. 1995, 31, 121–128. [Google Scholar] [CrossRef]
  47. Zhang, Z.; Shi, D.; Ren, F.; Yin, Z.; Sun, J.; Zhang, C. Evolution of quaternary groundwater system in North China plain. Sci. China 1997, 40, 276–283. [Google Scholar] [CrossRef]
Figure 1. The sketch map of the groundwater flow field and the detection point in the study area.
Figure 1. The sketch map of the groundwater flow field and the detection point in the study area.
Water 10 00752 g001
Figure 2. The shallow stratum histogram of a 321 m profile and groundwater conductivity at different depths.
Figure 2. The shallow stratum histogram of a 321 m profile and groundwater conductivity at different depths.
Water 10 00752 g002
Figure 3. Variation trend of the shallow groundwater component.
Figure 3. Variation trend of the shallow groundwater component.
Water 10 00752 g003
Figure 4. Sketch map of horizontal hydrochemical zonation of phreatic water in the study area.
Figure 4. Sketch map of horizontal hydrochemical zonation of phreatic water in the study area.
Water 10 00752 g004
Figure 5. The conductivity and dynamic diagram of groundwater of different depths (a) Dynamic curve of electrical conductivity at different depths; (b) Vertical conductivity distribution diagram of vertical depth).
Figure 5. The conductivity and dynamic diagram of groundwater of different depths (a) Dynamic curve of electrical conductivity at different depths; (b) Vertical conductivity distribution diagram of vertical depth).
Water 10 00752 g005
Figure 6. The Piper diagram for shallow unconfined water, middle confined water and deep confined water.
Figure 6. The Piper diagram for shallow unconfined water, middle confined water and deep confined water.
Water 10 00752 g006
Figure 7. Variation trend of ion content and retardation coefficient of clay at different depths.
Figure 7. Variation trend of ion content and retardation coefficient of clay at different depths.
Water 10 00752 g007
Figure 8. γMg, γCa and γNa/γCl for different types of groundwater.
Figure 8. γMg, γCa and γNa/γCl for different types of groundwater.
Water 10 00752 g008
Figure 9. Dynamic characteristics of groundwater level for different types of groundwater.
Figure 9. Dynamic characteristics of groundwater level for different types of groundwater.
Water 10 00752 g009
Figure 10. The sketch map of the vertical and horizontal zoning of water chemistry in the study area.
Figure 10. The sketch map of the vertical and horizontal zoning of water chemistry in the study area.
Water 10 00752 g010
Table 1. Basic parameters of typical test wells.
Table 1. Basic parameters of typical test wells.
Well NumberWell Diameter/mmWell Depth/mFilter Tube Depth/mGroundwater Level/mAquifer Lithology
WX0114683–823.59Silty clay
WX02146128–1221.03Silt
WX032355013–5020.99Silty clay
WX041409178–9019.16Silty clay
WX0514011399–10818.26silty sand
WX06180336274–33324.00Fine sand, fine sand, medium sand
Table 2. Unconfined common water chemical composition characteristic value (n = 70).
Table 2. Unconfined common water chemical composition characteristic value (n = 70).
IndexNa+Ca2+Mg2+HCO3SO42−ClTDS
Freshwater (n = 50)Maximum384.00198.36127.811009.28460.45358.991985.00
Minimum61.4050.7446.4858.564.0023.98549.89
Average167.9595.3976.90580.06174.85159.071254.23
Coefficient of variation (%)46.1629.8933.0733.3557.3755.1546.16
Brackish water (n = 15)Maximum688.45342.62226.691170.511165.65907.133973.94
Minimum223.5070.14101.1637.65178.04167.922014.23
Average469.83151.90161.63898.33548.04503.832733.57
Coefficient of variation (%)35.9848.9928.2520.6751.3048.6423.54
Saltwater (n = 5)Maximum2223.69646.68310.74540.086783.092738.7713,482.32
Minimum1244.00155.32470.961292.751470.891600.006255.17
Average1582.56360.53603.573191.943497.991177.339065.82
Coefficient of variation (%)20.4032.2416.6437.4619.0437.8023.45
Note: the unit of measurement for each mg·L−1.
Table 3. Characteristic values of conventional hydrochemical composition of groundwater (n = 141).
Table 3. Characteristic values of conventional hydrochemical composition of groundwater (n = 141).
IndexNa+Ca2+Mg2+HCO3SO42−ClPHTDS
Shallow salt water (n = 20)Maximum4201.201062.511022.001183.86494.005321.247.1116,984.20
Minimum1094.24264.47283.21645.541877.401217.006.995242.40
Average2029.10798.65791.81886.963251.693791.597.0410,077.96
Coefficient of variation (%)41.8225.6737.7912.2827.3232.981.0926.05
Middle confined water (n = 20)Maximum2250.16879.25853.25318.143266.293616.508.1613,482.32
Minimum1088.18395.33336.75192.37849.09849.097.015572.77
Average2108.01545.05560.93244.713266.292865.517.598991.46
Coefficient of variation (%)29.3225.3529.8317.9258.119.094.8731.04
Deep freshwater (n = 101)Maximum970.00125.16110.99540.081274.011036.32\3339.83
Minimum11.005.181.36186.2016.9721.52\372.72
Average345.2927.5923.74429.08205.99247.64\1285.15
Coefficient of variation (%)34.5991.5676.6726.1574.8759.03\28.21
Note: In addition to pH, the unit of measurement for each mg·L−1. Deep freshwater was not pH real-time monitored.
Table 4. The adsorption isotherm model fitting parameters of different clay depth adsorptions for Na+, Ca2+ and Mg2+ in the groundwater.
Table 4. The adsorption isotherm model fitting parameters of different clay depth adsorptions for Na+, Ca2+ and Mg2+ in the groundwater.
Clay Sampling Depth/mIon in GroundwaterLangmuir Equation ParametersHenry Equation Parameters
Qm/(mg·kg−1)KL/(L·mg−1)R2Kd/(L·mg−1)R2
6–6.2Na+4.60.00110.91640.23960.6311
Ca2+55.870.01500.74380.1020.4001
Mg2+2.720.02690.03942.92920.3851
48–48.2Na+1009.090.04530.7150.860 60.078 3
Ca2+466.670.00360.94080.83390.8511
Mg2+6000.00130.98530.84510.9637
51–51.2Na+23.920.00130.75741.07390.8907
Ca2+70.770.00160.86860.5430.7074
Mg2+549.450.00030.66050.81320.7468
71–71.2Na+130.210.00080.74640.48550.8182
Ca2+132.450.02200.92290.67380.9463
Mg2+194.550.00850.91160.43380.7336
84–84.2Na+9.630.00250.24210.63650.1848
Ca2+3.570.00600.99760.02620.826
Mg2+183.820.00040.55740.65700.9237
97–97.2Na+423.730.00060.92820.53680.9011
Ca2+40.680.00260.45540.16360.588
Mg2+182.820.00080.88100.68760.8955
102–102.2Na+359.710.00220.82670.45180.7132
Ca2+22.520.00790.63430.2070.8226
Mg2+1373.330.00030.98360.46060.9069
112–112.2Na+4.160.00130.76981.50950.9251
Ca2+4.160.00140.90310.21110.5134
Mg2+288.180.01780.94030.37740.8569
Table 5. The retardation coefficient of Na+, Ca2+ and Mg2+ in groundwater in cohesive soils of different depths.
Table 5. The retardation coefficient of Na+, Ca2+ and Mg2+ in groundwater in cohesive soils of different depths.
Clay Sampling Depth/mRd-NaRd-CaRd-MgRd-sum
61.2471.1054.0266.378
482.5942.5442.5657.703
512.9942.0082.5107.513
711.9022.2511.8065.959
842.1821.1012.2205.504
971.9971.3042.2775.578
1021.8391.3841.8555.079
1123.8031.3921.7016.896

Share and Cite

MDPI and ACS Style

Xing, L.; Huang, L.; Hou, X.; Yang, L.; Chi, G.; Xu, J.; Zhu, H. Groundwater Hydrochemical Zoning in Inland Plains and its Genetic Mechanisms. Water 2018, 10, 752. https://doi.org/10.3390/w10060752

AMA Style

Xing L, Huang L, Hou X, Yang L, Chi G, Xu J, Zhu H. Groundwater Hydrochemical Zoning in Inland Plains and its Genetic Mechanisms. Water. 2018; 10(6):752. https://doi.org/10.3390/w10060752

Chicago/Turabian Style

Xing, Liting, Linxian Huang, Xinyu Hou, Lizhi Yang, Guangyao Chi, Junxiang Xu, and Henghua Zhu. 2018. "Groundwater Hydrochemical Zoning in Inland Plains and its Genetic Mechanisms" Water 10, no. 6: 752. https://doi.org/10.3390/w10060752

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