Next Article in Journal
Child- and Parent-Related Correlates of Total and Prolonged Sedentary Time in 5- to 6-Year-Old Children
Previous Article in Journal
Evaluating the Coordination of Industrial-Economic Development Based on Anthropogenic Carbon Emissions in Henan Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Blocking Effect of Clay in Groundwater Systems: A Case Study in an Inland Plain Area

1
School of Water Conservancy and Environment, University of Jinan, Jinan 250022, China
2
Engineering Technology Institute for Groundwater Numerical Simulation and Contamination Control, Jinan 250022, China
3
Jinan Water Conservancy Construction Survey Design Research Institute, Jinan 250014, China
4
Shandong Provincial Bureau of Geology & Mineral Resources, Jinan 250013, China
5
Shandong Yu Qiang Geological Resource Exploration and Development Co., Ltd., Taian 271000, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(9), 1816; https://doi.org/10.3390/ijerph15091816
Submission received: 22 May 2018 / Revised: 11 July 2018 / Accepted: 2 August 2018 / Published: 22 August 2018
(This article belongs to the Section Environmental Science and Engineering)

Abstract

:
In order to increase understanding of the hydrogeochemical effects that influence changes in the quality of salt water, we investigated the distribution of saline and fresh water in an inland plain area and, in particular, the scarcity of fresh water resources. Taking the inland plain in Jiyang County as a specific case study, samples of undisturbed clay and underground saline water from different depths were collected to examine hydrogeological changes. A wide variety of methods was used to analyze the blocking effect of clay on the chemical characteristics of the groundwater. These include real-time monitoring for field water quality, tests for isothermal adsorption, a factor analysis model, physiochemical analysis, and correlation analysis. Our results show that the optimal adsorption isotherm of clay for Na+, Ca2+ and Mg2+ in groundwater conform to the established Henry and Langmuir equations for adsorption isotherms. The influence of clay mineral types and content on the blockage of Na+, Ca2+ and Mg2+ in groundwater samples were evident at different depths, with the clay adsorption capacity increasing in line with increases in the clay mineral content. Clay at different depths was found to have the strongest blocking effect on Na+ in groundwater, being systematically greater than its effect on Ca2+ and Mg2+. It is believed that the blocking effect of clay has an important influence on the hydrochemical zoning of groundwater in inland plains and the formation of saline water in groundwater systems. This study therefore provides concrete evidence in support of this supposed effect.

1. Introduction

Clay has a negative charge under normal or salt water conditions. In order to maintain its neutrality, clay adsorbs cations in groundwater. As Na+, Ca2+ and Mg2+ are the most common cations in groundwater, when the groundwater concentration reaches a certain value, they can be blocked by clay [1].
The strength of the blocking effect is not only related to the properties of the clay itself [2] and the hydrodynamic characteristics of the groundwater flow system [3,4], but also to the degree of mineralization of the groundwater [5,6]. In general, the ability of clay to block Na+, Ca2+ and Mg2+ in groundwater depends on the nature and mass fraction of the clay minerals. The larger the specific surface area of the clay minerals, the greater their adsorption and ability to block cations [7,8,9]. Clay particles with a large specific surface area form a double electrical layer on the surface of the clay when Na+, Ca2+ and Mg2+ interact with groundwater. The pH value of the groundwater, however, causes a change in the properties of this double electrical layer, which affects the adsorption capacity of the clay particles for Na+, Ca2+ and Mg2+ [8,9,10,11,12]. So, to give an example, once the mass concentration of Na+ (Ca2+ and Mg2+) in the groundwater exceeds 200 mg/L, the clay will begin to absorb the Na+ [13,14]. A large number of studies have also demonstrated that the retarding strength of clay on underground salt water decreases when there is an increase in effective porosity or the hydraulic gradient. Once the hydraulic gradient passes a certain threshold, shallow salty water will flow into deep freshwater [15,16,17,18]. During this leakage process, the retardation coefficient related to the strength of the clay blocking effect becomes inversely proportional to its own permeability coefficient [19,20,21].
There is a large body of literature dedicated to the water chemistry mechanisms associated with saline water leakage. Various methods have been used to examine the system in more detail, including factor analysis [22], the ion ratio method [23], isotope tests [24], trace element tests [25], and so on. Factor analysis is the principal approach used for the comprehensive evaluation of multi-factor systems. Factor analysis is able to fully encompass overlaps in information between different evaluation indicators. It is also able to bring about a comprehensive reduction of high-dimensional variables by retaining original information to the maximum extent. The weight of each index can then be objectively determined, thus offsetting any subjective arbitrariness.
There are many ways to study the blocking effect of clay on groundwater chemical components, principal amongst them being high-pressure permeation tests [26] and adsorption isothermal tests [27]. Adsorption isotherm experiments are relatively straightforward, with a shorter duration than pressure permeation tests and the adsorption isotherm equation and test data have a higher degree of fit. Some scholars have found a particular value in using a Langmuir formula to describe the adsorption process. This is because it defines an adsorption limit at a given adsorption site that will meet certain specified criteria, then estimates the total adsorption capacity of the clay for different adsorbates. This provides a solid foundation for the study of the clay adsorption mechanism. The process for clay adsorption of groundwater chemical components is mainly linear adsorption (what is known as the Henry model). As an example, the adsorption isotherms of clay for niobium and ammonia nitrogen conforming to the Henry model will have retardation coefficients of 2523.36 and 60.87, respectively [28,29]. For monolayer chemical adsorption, the adsorption isotherms of clay for organophosphorus and Mn2+ will conform to the Langmuir model [30,31,32]. For physical adsorption, the clay adsorption of ammonium ions [33,34] will conform to what is known as the Freundlich model [18,21]. On the basis of these various models, a significant number of previous investigations have shown that clay has a blocking effect on pollutants and have sought to highlight some of the key factors involved. However, there is a notable lack of discussion regarding the mechanisms whereby clay can serve to block Na+, Ca2+ and Mg2+ in groundwater.
The presence of finer sediment particles in inland plains, the high proportion of clay minerals, and stagnant groundwater runoff all contribute to a vertical alternation of saline and fresh water, with hydrochemically characteristic horizontal and vertical zoning [35]. Over the past few decades, excessive exploitation of deep fresh water in inland plains has caused shallow saline water to flow into deep aquifers, which in turn has led to their salinization. In that case, controlling and preventing salt water from continuing to invade deep fresh water in inland plains has become important to the sustainable use of groundwater resources. As an important aspect of the saline water leakage system in inland plains, clay plays a key role in the distribution of brackish water and the downward invasion of underground salt water. In order to analyze the influence of adsorption and the blocking effect of clay on the chemical components of groundwater, we used geological drilling and pumping tests to collect undisturbed clay and water samples from different depths in the inland plains of Jiyang county in China. Real-time monitoring of field water quality, adsorption isotherm experiments, factor analysis and correlation analysis were all then used to investigate the mechanisms whereby clay can serve to block major common cations. This, in turn, will serve to increase our understanding of the hydrogeochemical effects and development processes that influence the distribution of saline water.

2. Materials and Methods

2.1. Overview of the Study Area

The research took place in Jiyang County, Shandong Province, an inland region of China (Figure 1) with a warm temperate continental monsoon climate. The average annual temperature is 14.3 °C, the average annual rainfall is 665.7 mm and the average annual evaporation is 1700.2 mm. So its four seasons are distinct. The study area is mainly covered with loose Quaternary deposits, the other deposits being Neogene, Permian, Carboniferous, and Ordovician (Figure 2). The 500 m of shallow water in the area can mostly be found in the pores of the Quaternary and Neogene loose deposits. Salt water has a wide distribution in the mainland. Here, the degree of salinity for the salt water in the vertical direction, first of all increases, then decreases with an increase in depth. The change range is 5–9 g/L. In the plane, the water chemistry is varied and its composition is complex. The Total Dissolved Solids (TDS) is generally 2~5 g/L. The salt water is located in the runoff stagnation area between the ancient rivers.

2.2. Sample Collection

After taking into account current field hydrogeological conditions and site constraints, Wangxing Village, Sunjing Town was selected as a typical location for investigating salt water in the inland plains. Four boreholes with a range of 65 m × 45 m were selected to serve as the water level and water quality measurement boreholes (see Figure 1). The boreholes were numbered WX01, WX02, WX03 and WX04, with No. 4 being an existing borehole. Water filtration pipes for hydrogeological observations were placed within the boreholes at depths of 8–12 m, 79–94 m, 98–107.9 m, and 27–49.2 m. The lithology for the hydrogeological observation borehole WX01 was silty soil and silty clay soil. For borehole WX02, the lithology was silty clay and cohesive soil. The lithography for borehole WX03 was silty clay and clay, whilst the lithology for borehole WX04 was dominated by silt, fine sand, and medium sand. The basic parameters for the hydrogeological drillings are shown in Table 1.
The physical properties of the collected undisturbed soil (grain specific gravity, water content, density, and void ratio) were tested using a geotechnical test. The particle size of the soil was determined using a particle size sedimentation analysis test, and the X-Ray Powder Diffusion (XRD) analysis software, MDI Jade (San Francisco, CA, USA), was used for phase analysis and quantitative analysis of the clay. Adsorption and the blocking effect of the clay on Na+, Ca2+, and Mg2+ in the groundwater at different depths were determined by means of isothermal adsorption experiments. In accordance with technical regulations, undisturbed soil samples from depths of 6.0–6.2 m, 48.0–48.2 m, 51.0–51.2 m, 71.0–71.2 m, 84.0–84.2 m, 97.0–97.2 m, 102.0–102.2 m and 112.0–112.2 m were collected. The void ratio and clay mineral content of the undisturbed soil samples were measured, with the void ratio being between 0.496 and 0.849 and the clay mineral content being within the range of 21–46%. The coefficient of variation was 24.49% (see Table 2). The soil samples (3000 g) used for sorption isotherm tests were collected from each depth sample, from which the debris was removed. The samples were air dried in a cool ventilated area, sieved, ground, fully mixed, bottled and retained for adsorption tests.
Water samples for analysis were taken from the four hydrological observation boreholes, sealed and stored [36]. Sample preparation was completed within 7 days of returning the samples to the laboratory. The milliequivalent percentage of water chemical composition is shown in Table 3. The hydrochemistry type for the groundwater in boreholes WX01, WX03 and WX04 was ClSO4–NaMg. The hydrochemistry type for the groundwater in borehole WX02 was ClSO4–MgNa.

2.3. Real-Time Field Groundwater Quality Monitoring

Groundwater from the water filtration sections of the four hydrological observation boreholes was analyzed for water temperature (T), pH and electrical conductivity (EC) using an Aqualtroll600 multi-parameter water quality detector (In-Situ, Los Angeles, CA, USA). The test accuracy was ±0.1 °C, ±0.1 and 1 ± 0.1 μS/cm for the three physiochemical indicators, respectively. As conductivity reflects the number of ions in the water, this indicator can be used to characterize changes in the degree of mineralization [37].

2.4. R-Type Factor Analysis

Numerous studies have suggested that R-factor analysis is the best approach to establishing water chemistry [17], so this was used to infer the hydrochemistry for the groundwater at different depths. The test indicators were Na+, Ca2+, Mg2+, HCO3, Cl, SO42−, total hardness, pH and EC. After correlation analysis and factor analysis using SPSS 10.0 (San Francisco, CA, USA), the correlation coefficient matrix for each index was obtained. A main component analysis method was used to extract the eigenvalues and the main factor was selected. In order to accurately interpret the meaning represented by the main factor, a maximum variance rotation method was used to rotate the index matrix.

2.5. Indoor Adsorption Isothermal Test and Determination of the Retardation Factor

2.5.1. Indoor Adsorption

Indoor adsorption isothermal tests were performed using a batch method [38]. Clay samples from depths of 6, 48, 51, 71, 84, 97, 102 and 112 m were selected, washed and dried with ultra-pure water. The test water samples were diluted to 10 solutions with a gradient of conductivity of 2000, 4000, 6000, 8000, 10,000, 12,000, 14,000, 16,000, 17,000, and 188,400 us/cm. First of all, ten samples of 30 g of clay were transferred to a 250 mL stoppered erlenmeyer flask. The above-mentioned solutions were then added separately. The temperature of the solution was maintained at 18 ± 1 °C throughout the test and the pH was 7.5 ± 0.5. The erlenmeyer flask was then placed in a constant-temperature shaker bath (SHZ-82A, Shimadzu company, Tokyo, Japan) to oscillate for 3 h. Once adsorption had reached a point of saturation, samples were extracted. These samples were centrifuged at 4000 r/min, then filtered through a 0.45 μm microporous membrane. The filtrate was analyzed using a PXSJ-216F ion meter (Shanghai Yidian Scientific Instrument Company, Shanghai, China) to determine the mass concentration of Na+ in the samples. Then a Shimadzu-AA7000 (Shimadzu company, Tokyo, Japan.) atomic absorption spectrophotometer was used to measure the concentration of Ca2+ and Mg2+ in the filtrate. The sample analysis was repeated three times for each solution. Three hundred and seven groups of valid data were obtained by the indoor adsorption isothermal test. Adsorption isotherm analysis was then performed on the valid data.

2.5.2. Determination of the Retardation Factor

The adsorption capacity of the adsorbent (clay) for the adsorbates (Na+, Ca2+ and Mg2+; qe) was calculated as follows [39]:
q e = ( C 0 C e )   V M
The adsorption isotherm indicates the relationship between the adsorbent surface and the adsorbate in the solution. The isotherm was fitted using the Langmuir isotherm Equations (2) and the Henry isotherm Equations (3). The effectiveness of the fittings were evaluated using the determination coefficient R2 [39].
q e = q m K L C e 1 + K L C e
q e = K d C e
For this analysis, the linear distribution coefficient and the effective porosity of the clay were calculated first (4), then the retardation coefficient (Rd) (5) [40]:
n = e 0 1 + e 0
R d = 1 + 1 n n K d
In the above Equations (1)–(5), C0 is the initial mass concentration of the adsorbate in the solution (mg L−1); Ce is the adsorption equilibrium mass concentration (mg L−1) of the adsorbate in the solution; qe is the equilibrium adsorption capacity of the adsorbent (mg kg−1); V is the volume of the solution (L); M is the mass of the adsorbent (kg); qm is the maximum adsorption capacity of the adsorbent in the Langmuir adsorption equation (mg kg−1); KL is the adsorption equilibrium constant (L mg−1) of the adsorbent in the Langmuir adsorption equation; Kd is the partition coefficient in the Henry equation (L kg−1); e0 is the void ratio; n is the effective porosity; and Rd is the retardation coefficient.

3. Results and Analysis

3.1. The Effects of Clay Adsorption on the Chemical Characteristics of Groundwater

Three main factors (F1–F3) were identified through the factor analysis model (Table 4 and Table 5). The first main factor (F1) influencing groundwater in the aquifer at a depth of 8.0–12.0 m included Mg2+, Cl, SO42− and total hardness, with a variance contribution rate of 64.164%. The main factor load values of the four hydrochemical indicators were all above 0.8. The second main factor (F2) was mainly determined by the pH, with a variance rate of 31.235%. This result indicates that F1 and F2 reflect the effects of leaching and cation exchange adsorption, together with the influence of an acid-based evolution of the chemical composition of the groundwater during runoff.
For the groundwater at a depth of 13.0–49.2 m, F1 included seven water chemistry indicators (Na+, Ca2+, Mg2+, Cl, SO42−, total hardness and EC), with a variance contribution rate of 79.446%. The main factor load value was above 0.9 and EC and F1 had a negative correlation. In this case, F2 was mainly determined by HCO3 and the pH and its variance contribution rate was 16.941%. Both the load values of the two water chemistry indicators in F2 were above 0.5 and the pH was negatively correlated with F2. This indicates that F1 and F2 reflect the effects of mixing, cation exchange adsorption, carbonic acid evolution and the acid-based evolution of the chemical composition of the groundwater.
F1 in the groundwater aquifer at a depth of 79.0–94.0 m included Na+, SO42−, pH and EC. The variance contribution rate was 34.943% and the values of the four water chemical indicators in F1 were above 0.6. F2 was mainly determined by Mg2+ and Cl, with a variance contribution rate of 22.397%. The load values of the two water chemical indicators were both greater than 0.7 and there was a negative correlation between Cl and F2. The third main factor (F3) in this aquifer was mainly determined by Ca2+, with a variance contribution rate of 13.456%. The water chemistry indices in the main factor load values were all greater than 0.6. In this case, F1, F2 and F3 reflect the effects on the chemical composition of the groundwater of concentration, acid-based evolution, leaching, ion exchange adsorption and dissolution/precipitation of calcite and gypsum.
Groundwater from an aquifer depth of 98.0−107.9 m recorded an F1 that was dominated by Ca2+, SO42− and total hardness, with a variance contribution rate of 54.647%. The three water chemical indicators in the main factor load value were all greater than 0.7. F2 was mainly determined by Ca2+, HCO3 and EC, with a variance contribution rate of 14.607%. The water chemical index was greater than 0.45 and there was a negative correlation between Cl and F2. F3 was mainly determined by EC, its variance contribution rate being 12.138%. The water chemistry index in the main factor load value was greater than 0.7. Here then, F1, F2 and F3 are reflecting the effects on the chemical composition of the groundwater of dissolution/precipitation, leaching, ion exchange adsorption and the mixing of calcite and gypsum.

3.2. The Adsorption Model of Na+, Ca2+ and Mg2+ in Groundwater by Clay

The adsorption characteristics and adsorption isotherms of clay for Na+, Ca2+ and Mg2+ in groundwater (Figure 3) show that the adsorption capacity increases in line with an increase in the concentration of Na+, Ca2+ and Mg2+. It can also be seen that, as the concentration of the ions increases, the adsorption capacity continues to increase. This indicates a fast-paced process of change. The maximum adsorption capacity of clay for Na+, Ca2+ and Mg2+ in groundwater corresponds to the maximum adsorption capacity obtained from the Langmuir equation. The deviation was only 5–15% and the adsorption isotherm could be fitted to the Henry equation. As the adsorption isotherms of Na+, Ca2+ and Mg2+ for clay have a relatively high correlation (see Table 6) and conform to the Henry and Langmuir equations, this indicates either monomolecular chemical adsorption or linear adsorption processes.

3.3. Clay Adsorption Affecting the Formation of Na+, Ca2+ and Mg2+ in Groundwater

The maximum buffering capacity (MBC) of clay is the product of the maximum adsorption quantity (qm) and the adsorption equilibrium constant (KL) in the Langmuir sorption isotherm equation. The linear distribution coefficient (Kd) in the clay adsorption process is the ratio of the conventional ion concentration in the clay to the conventional ion concentration in the groundwater. These factors all reflect the adsorption capacity of clay. Based on the calculation and analysis of MBC and Kd (see Figure 4), the order of adsorption of clay on conventional ions in groundwater at a depth of 6 m and 48 m was: Mg2+ > Na+ > Ca2+. At depths of 51 m and 112 m the order was: Na+ > Ca2+ > Mg2+; and at depths of 71 m, 84 m, 97 m and 102 m the order was: Na+ > Mg2+ > Ca2+. The correlation between the linear distribution coefficient for the different depths of clay and the total content of clay minerals was as high as 0.8754 and the linear distribution coefficient increased with an increase in the clay mineral content (Figure 5). The difference in adsorption capacity may therefore be due to differences in the amount of illite and montmorillonite, the main clay minerals that were found to differ at different depths, together with differences in the type of ion exchange and cation exchange adsorption capacity [41]. To sum up, the adsorption ability of clay for Na+ in groundwater is relatively strong with a certain controlling effect on the formation of chemical components in the groundwater, whilst the clay adsorption capacity is affected by the quantity of clay minerals.

3.4. The Blocking Effect of Clay on the Vertical Migration of Na+, Ca2+ and Mg2+ in Groundwater

According to the retardation coefficient of clay for Na+, Ca2+ and Mg2+ in groundwater (see Table 7), the clay at certain depths has a different retarding strength for the different conventional ions in the groundwater. The order of blocking strength of clay on conventional ions at depths of 6 m, 48 m and 112 m was: Na+ > Mg2+ > Ca2+. At depths of 51 m and 71 m, the order was: Na+ > Ca2+ > Mg2+; and at depths of 84 m, 97 m and 102 m the order of blocking was: Mg2+ > Na+ > Ca2+. Different depths of clay have different degrees of blocking effect on Na+ (Ca2+ and Mg2+) in the groundwater. The retarding strength of clay for Na+, Ca2+ and Mg2+ in the groundwater decreased from 6 m to 102 m to 112 m. From depths of 6 m to 97 m to 112 m, and from 6 m to 102 m to 112 m, a trend of increase-decrease was recorded.
According to the real-time monitoring of groundwater quality from 31 December 2016 to 4 March 2017, the groundwater water quality in each hydrological observation borehole was relatively stable. The average conductivity in the four boreholes, at groundwater depths of 8.0–12.0 m, 13.0–49.2 m, 79.0–90.0 m and 98.0–107.9 m, was 8801.82 μS/cm, 14,603.57 μS/cm, 18,734.42 μS/cm and 11,314.42 μS/cm, respectively. The correlation coefficient of the retardation coefficient of clay for Na+, Ca2+ and Mg2+ in groundwater and the milliequivalent percentage of Na+, Ca2+ and Mg2+ in equivalent depths of groundwater were as high as 0.89078 (Figure 6a). The content of Na+, Ca2+ and Mg2+ in groundwater increases as the clay blocking strength increases. The correlation between the sum of the retardation coefficients of clay for Na+, Ca2+ and Mg2+ in groundwater and the corresponding depth’s electrical conductivity observed in real-time was as high as 0.89025 (Figure 6b). This finding indicates that the larger the sum of the retardation coefficients, the higher the conductivity. Thus, we can observe that the blocking effect of clay has an important influence on the vertical migration of Na+, Ca2+ and Mg2+ in groundwater, the type of water chemistry and water quality zoning. This further confirms that clay has a controlling effect on the downward movement of salt water.

4. Discussion

(1) It is a consensus that illite and montmorillonite in clay can adsorb certain ionic components in groundwater, such as Cu2+ [42], Pb2+, Zn2+, Cd2+, Cr3+ [43], cadmium [44], Na+, Ca2+ and Mg2+ [45,46]. These minerals, affecting the chemical properties of groundwater, are contained in the clay of this area, which is consistent with the conclusions of the factor analysis methods in this paper.
(2) Former researchers concluded different adsorption models, to name a few, Henry model, namely linear adsorption; Langmuir model, i.e., mono-layer chemical adsorption; Freundlich model, namely physical adsorption. The adsorption isotherm curve of clay for anthracene and ammonia nitrogen conforms to the Henry model [28,29]. The process of clay adsorption of organic phosphorus and Mn2+ in groundwater is consistent with the Langmuir model [30,31,32]. The adsorption isotherm curve of clay for ammonium ions accords with the Freundlich model [33,34]. It can be seen that different ions have different models. Therefore, it is reasonable to say that the adsorption model of clay in this area is different for Na+, Ca2+ and Mg2+ in groundwater.
(3) Different depths of clay have different degrees of retardation on Na+ (Ca2+ and Mg2+) in groundwater, and with the increase of depth, the block strength of clay in groundwater decreases, which is related to the nature of the clay and its age of sedimentation. The newly deposited clay has strong retardation ability. The aged clay deposited in geological history can be traced back to a long time, and the retardation ability is reduced due to the secondary fracture caused by clay consolidation. Within the tens of thousands of square kilometers of the North China Plain, such as Jiyang, Dezhou, Cangzhou, Tianjin, and Dongying, different groundwater types exist at different depths [47,48,49]. In the vertical direction, there is shallow fresh water—middle salt water—deep fresh water within the depth of 600 m [33]. The clay of less than 200 m is relatively new so that the block strength is stronger, and the groundwater upward conductivity is significantly zonal. The clay deposits at a depth of 300–600 m are ages ago, but the difference between 300–600 m is not salient, and only deep fresh water is present. Therefore, the retardation of clay is one of the factors affecting the vertical zoning of groundwater hydrochemistry.

5. Conclusions

In recent years, the downward intrusion of shallow saline groundwater in inland plains (especially the North China Plain) has disabled a large number of deep aquifer extraction wells because of the salinization of deep groundwater. Control and prevention of this downward intrusion is therefore key to the sustainable use of groundwater resources in the inland plains. Clay can play a critical role in controlling the downward intrusion of shallow saline groundwater, so understanding the mechanism involved in this downward intrusion is a fundamental prerequisite for being able to exercise such control. There is large body of literature dedicated to the blocking effects of clay, but most of it focuses on the factors affecting the clay blocking coefficient, such as effective porosity, permeability, concentration of seepage fluid, particle release and hydraulic gradient. There is an additional literature that focuses on the blocking effect of clay on inorganic substances and pollutants. However, research regarding how clay affects the downward intrusion of shallow saline groundwater and conventional ion blockage is still largely absent.
In this study, we analyzed the blocking effect of clay on Na+, Ca2+ and Mg2+ in groundwater using real-time monitoring of field water quality, adsorption isothermal experiments, factor analysis and correlation analysis. The results have shown that the formation of groundwater chemical components in inland plains is affected by cation exchange adsorption, acid-based evolution, mixing, carbonic acid evolution, concentration, dissolution/precipitation of calcite and gypsum, and leaching. The effects of precipitation and leaching on the quality of shallow saline water, however, are not obvious. It was also found that the optimal adsorption isotherms of clay for Na+, Ca2+ and Mg2+ in groundwater at different depths are in accordance with the Henry and Langmuir equations. The results show that clay has a blocking effect on Na+, which is consistent with the results of previous tests. However, the research results regarding the blocking effect of clay on Ca2+ and Mg2+ in groundwater are limited and make no comparison. Here we have found that the adsorption capacity of clay increases with an increase in the concentration of Na+, Ca2+ and Mg2+. We have also found that, as the concentration of ions increases, the adsorption capacity continues to increase, suggesting a rapid rate of change. The intensity of the blocking effect of clay on Na+ (Ca2+, Mg2+) varied according to depth and was controlled by the type and content of minerals, with clay having a strong capacity to retard Na+ at certain depths. We can conclude, in that case, that the blocking effect of clay influences the formation and zoning of water chemical components, further confirming that clay has a controlling effect on the downward movement of salt water.
However, the formation of groundwater chemistry components is not only related to hydrogeochemistry, but also closely related to palaeogeographic climate, hydrodynamic conditions and stratum lithology. This study has only focused on the blocking effect that clay can have on Na+, Ca2+ and Mg2+ in groundwater. The controlling effect of clay on conventional anions and other factors affecting the formation of water chemistry still require further research.

Author Contributions

L.X and Y.Y. conceived and designed this research; L.H., J.X. and W.Z. performed the experiments; G.C. and X.H. analyzed the data; L.X. and Y.Y. wrote the paper. All authors have read, edited and approved the final manuscript.

Funding

This research is supported by the National Natural Science Foundation of China (41772257, 41472216) and the Project of Shandong Province Higher Educational Science and Technology Program (J17KA191).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Majhi, A.; Monash, P.; Pugazhenthi, G. Fabrication and characterization of γ-Al2O3–clay composite ultrafiltration membrane for the separation of electrolytes from its aqueous solution. J. Membr. Sci. 2009, 340, 181–191. [Google Scholar] [CrossRef]
  2. Frank Jones, J.R. Influence of chemical composition of water on clay blocking of permeability. J. Pet. Technol. 1964, 16, 441–446. [Google Scholar] [CrossRef]
  3. Yu, D.; Jia, Z.; Wei, X.; Huafeng, S.U.; Guan, Z.; Zhou, J. Salt migration phenomena during phase transition in cooling water and its impact on leakage current. Proc. CSEE 2012, 32, 184–191. [Google Scholar]
  4. Wang, D.; Hou, Q.; Luo, Y.; Zhu, Y.; Fan, H. Blocking ability and flow characteristics of nitrogen foam stabilized with clay particles in porous media. J. Dispersion Sci. Technol. 2015, 36, 170–176. [Google Scholar] [CrossRef]
  5. Li, Y.; Cheng, J.; Li, M.M.; Cui, L.H. Experimental study on migration and transformation in cohesive soil of salt water. Water Resour. Prot. 2015, 31, 3, 88–93. [Google Scholar]
  6. Hagedorn, F.; Spinnler, D.; Siegwolf, R. Increased n deposition retards mineralization of old soil organic matter. Soil Biol. Biochem. 2003, 35, 1683–1692. [Google Scholar] [CrossRef]
  7. Kahle, M.; Kleber, M.; Jahn, R. Retention of dissolved organic matter by phyllosilicate and soil clay fractions in relation to mineral properties. Org. Geochem. 2004, 35, 269–276. [Google Scholar] [CrossRef]
  8. Jackson, M.L. An Introduction to Clay Colloid Chemistry, 1st Ed.; van Olphen, H., Ed.; Interscience Publishers: New York, NY, USA, 1963; Volume 27, p. 4. [Google Scholar]
  9. Churchman, G.J. Is the geological concept of clay minerals appropriate for soil science? A literature-based and philosophical analysis. Phys. Chem. Earth Parts A/B/C 2010, 35, 927–940. [Google Scholar] [CrossRef]
  10. Zuo, W.; Yejing, X.U.; Ying, W.; Mingyan, L.I.; University, H.U.; Amp, H.I. Experimental studies on the membrane behavior of the cohesive soil. Ion Exch. Adsorpt. 2013, 29, 551–560. [Google Scholar]
  11. Chang, F.R.C.; Sposito, G. The electrical double layer of a disk-shaped clay mineral particle: Effects of electrolyte properties and surface charge density. J. Colloid Interface Sci. 1996, 178, 555–564. [Google Scholar] [CrossRef]
  12. Abdou, F.M.; Elnennah, M.E.; Foda, M.S. Selective adsorption of potassium and calcium on clay minerals [Egypt]. Crop Sci. 1987, 48, 1292–1304. [Google Scholar]
  13. Tertre, E.; Savoye, S.; Hubert, F.; Prêt, D.; Dabat, T.; Ferrage, E. Diffusion of water through the dual porosity swelling clay mineral vermiculite. Environ. Sci. Technol. 2018. [Google Scholar] [CrossRef] [PubMed]
  14. Cao, W.B.; Wan, L.; Gong, B.; Wang, D.C. Experimental on osmosis through clayey soil under the condition of variable water levels. Hydrogeol. Eng. Geol. 2006, 2, 118–122. [Google Scholar]
  15. Wang, P.; Hu, F.S.; Han, Z.T.; Kong, X.K. Experimental study on dispersion of saline water flowing through clay-bearing soil. South North Water Transf. Water Sci. Technol. 2014, 12, 101–104. [Google Scholar]
  16. Hua, K.; Yang, S.; Wang, W.; Feng, L.; Li, Y.; Wang, Y.; Mu, G. Research about Adsorption Characteristics of Pb2+ on Sediments at Weihe river. In Proceedings of theInternational Symposium on Water Resource and Environmental Protection, Xi’an, China, 20–22 May 2011; pp. 1526–1528. [Google Scholar]
  17. Alabi, O. Permeability Characterization of Clay under Different Hydraulic Gradient; Springer Publishing: New York, NY, USA, 2017; Volume 24, pp. 29–48. [Google Scholar]
  18. Xing, L.; Zhang, F.J.; Li, C.S.; Yang, L.Z.; Wang, L.Y. Hydro-chemical characteristics of shallow groundwater in the north Shandong plain. J. Irrig. Drain. 2015, 34, 6, 90–94. [Google Scholar]
  19. Li, X.Q.; Zhu, Y.N.; Guang, Y.U.; Xie, X.; Song, S. Study on permeation and migration law of ammonia nitrogen in aquitard. J. China Hydrol. 2011, 31, 3, 71–75. [Google Scholar]
  20. Achari, G.; Joshi, R.C.; Bentley, L.R.; Chatterji, S. Prediction of the hydraulic conductivity of clays using the electric double layer theory. Can. Geotech. J. 1999, 36, 783–792. [Google Scholar] [CrossRef]
  21. Carvalho, M.N.; Abreu, C.A.M.D.; Benachour, M.; Sales, D.C.S.; Baraúna, O.S.; Sobrinho, M.A.D.M. Applying combined Langmuir-Freundlich model to the multi-component adsorption of BTEX and phenol on smectite clay. Adsorpt. Sci. Technol. 2012, 30, 691–700. [Google Scholar] [CrossRef]
  22. Alariqi, W.S.T.; Abduljalil, A.D.S. Assessment of hydrochemical quality of ground water under some urban areas within Sana'a secreteriat. Eclética Quím. 2018, 35, 77–84. [Google Scholar] [CrossRef]
  23. Batayneh, A.; Zaman, H.; Zumlot, T.; Ghrefat, H.; Mogren, S.; Nazzal, Y.; Elawadi, E.; Qaisy, S.; Bahkaly, I.; Altaani, A. Hydrochemical facies and ionic ratios of the coastal groundwater aquifer of Saudi gulf of Aqaba: Implication for seawater intrusion. J. Coast. Res. 2014, 30, 75–87. [Google Scholar] [CrossRef]
  24. Seibert, S.L.; Holt, T.; Reckhardt, A.; Ahrens, J.; Pollmann, T.; Giani, L.; Waska, H.; Böttcher, M.E.; Greskowiak, J.; Massmann, G. Hydrochemical evolution of a freshwater lens below a barrier island (Spiekeroog, Germany): The role of carbonate mineral reactions, cation exchange and redox processes. Appl. Geochem. 2018, 92, 196–208. [Google Scholar] [CrossRef]
  25. Tesfamichael, S.G.; Ndlovu, A. Utility of aster and landsat for quantifying hydrochemical concentrations in abandoned gold mining. Sci. Total Environ. 2018, 618, 1560. [Google Scholar] [CrossRef] [PubMed]
  26. Sunitha, N.; Chandrashekar, J.S.; Nandini, N. Groundwater quality in handigundi state forest block of Ramanagaram Taluk, Karnataka, India. Res. J. Pharm. Technol. 2012, 1, 2012. [Google Scholar]
  27. Kelleher, B.P.; Oppenheimer, S.F.; Han, F.X.; Willeford, K.O.; Simpson, M.J.; Simpson, A.J.; Kingery, W.L. Dynamical systems and phase plane analysis of protease−clay interactions. Langmuir 2003, 19, 9411–9417. [Google Scholar] [CrossRef]
  28. Vico, L.I.; Acebal, S.G. A langmuir modelation for Mn2+ adsorption by fibrous clay minerals. An. Asoc. Quim. Argent. 2000, 88, 45–52. [Google Scholar]
  29. Song, D.; Hongyan, L.I.; Xuqian, L.I.; Zhou, J.; Luo, Y. Study on transport and transformation rules of organophosphor in aquitard. Environ. Pollut. Control 2011, 33, 1–205. [Google Scholar]
  30. Zhou, W.; Ren, L.; Zhu, L. Reducement of cadmium adsorption on clay minerals by the presence of dissolved organic matter from animal manure. Environ. Pollut. 2017, 223, 247–254. [Google Scholar] [CrossRef] [PubMed]
  31. Zayed, A.M.; Selim, A.Q.; Mohamed, E.A.; Wahed, M.S.M.A.; Seliem, M.K.; Sillanpӓӓ, M. Adsorption characteristics of Na-A zeolites synthesized from Egyptian kaolinite for manganese in aqueous solutions: Response surface modeling and optimization. Appl. Clay Sci. 2017, 140, 17–24. [Google Scholar] [CrossRef]
  32. Luukkonen, T.; Sarkkinen, M.; Kemppainen, K.; Rämö, J.; Lassi, U. Metakaolin geopolymer characterization and application for ammonium removal from model solutions and landfill leachate. Appl. Clay Sci. 2016, 119, 266–276. [Google Scholar] [CrossRef]
  33. Wang, L.Y.; Xing, L.T. An experimental study of water-rock action in inland saltwater region. China Rural Water Hydropower 2013, 10, 119–121. [Google Scholar]
  34. Lebron, I.; Suarez, D.L.; Amrhein, C.; Strong, J.E. Size of mica domains and distribution of the adsorbed Na-Ca ions. Clays Clay Miner. 1993, 41, 380–388. [Google Scholar] [CrossRef]
  35. Li, C.; Xing, L.; Wang, L.; Zhang, F.; Bo, L. Mixing experimental study on water-rock interaction (WRI) in the inland shallow salt groundwater area. Open Chem. Eng. J. 2015, 9, 25–28. [Google Scholar]
  36. Zheng, Z.; Yong, Q.; Wang, G.X.; Xuehai, F.U. Numerical description of coalbed methane desorption stages based on isothermal adsorption experiment. Sci. China Earth Sci. 2013, 56, 1029–1036. [Google Scholar] [CrossRef]
  37. Liu, Y.; Chang, X.Y.; Hai ni, L.I.; Lin, T.T.; Gao, J.; Qiang, F.U. Relationships between salinity, conductivity and water content of soils in a pristine tidal flat in Jiangsu province. Water Saving Irrig. 2015, 8, 4–7. [Google Scholar]
  38. Omer, M.F.; Renuga, V.; Sripriya, N.; Prakash, N.K.U. Adsorption isothermal studies on the removal of bod and cod of a leather tannery effluent using clinoptilolite. Int. J. Chem. Res. 2017, 10, 496–500. [Google Scholar]
  39. Barreca, S.; Orecchio, S.; Pace, A. The effect of montmorillonite clay in alginate gel beads for polychlorinated biphenyl adsorption: Isothermal and kinetic studies. Appl. Clay Sci. 2014, 99, 220–228. [Google Scholar] [CrossRef] [Green Version]
  40. 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]
  41. Qian, J.; Wang, L.; Ma, L.; Lu, Y.H.; Zhao, W.; Zhang, Y. Multivariate statistical analysis of water chemistry in evaluating groundwater geochemical evolution and aquifer connectivity near a large coal mine, Anhui, China. Environ. Earth Sci. 2016, 75, 747. [Google Scholar] [CrossRef]
  42. Zhu, C.S.; Xu, G.X.; Wei, T.; Chen, W.B. Isotherm Adsorption of Cu(Ⅱ) Ions by Marine Sediments from Lianyungang. J. Huaihai Inst. Technol. (Nat. Sci. Ed.) 2008, 17, 49–52. [Google Scholar]
  43. Cruz-Guzmán, M.; Celis, R.; Hermosín, M.C.; Koskinen, W.C.; Nater, E.A.; Cornejo, J. Heavy Metal Adsorption by Montmorillonites Modified with Natural Organic Cations. Soil Sci. Soc. Am. J. 2006, 70, 215–221. [Google Scholar] [CrossRef] [Green Version]
  44. Keren, R. Effect of Exchangeable Ions and Ionic Strength on Boron Adsorption by Montmorillonite and Illite1. Clays Clay Miner. 1982, 30, 341–346. [Google Scholar] [CrossRef]
  45. Tehrani-Bagha, A.R.; Bahrami, H.; Movassagh, B.; Arami, M.; Amirshahi, S.H.; Menger, F.M. Dynamic Adsorption of Gemini and Conventional Cationic Surfactants onto Polyacrylonitrile. Colloids Surf. A 2007, 307, 121–127. [Google Scholar] [CrossRef]
  46. Rytwo, G.; Banin, A.; Nir, S. Exchange Reactions in the Ca-Mg-Na-Montmorillonite System. Clays Clay Miner. 2013, 442, 76–85. [Google Scholar] [CrossRef]
  47. Zhang, W.; Cao, Y.; Jiang, J.; Wang, W.; Duan, L.; 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 Water Rock Interaction, Proceedings of the 12th International Symposium on water-rock interaction, Kunming, China, 31 July–5 August, 2007; Taylor and Francis: London, UK, 2007; Volume 12, pp. 407–414. [Google Scholar]
  48. Scheiber, L.; Ayora, C.; Vázquez-Suñé, E.; Cendón, D.I.; Soler, A.; Custodio, E.; Baquero, J.C. Recent and Old Groundwater in the Niebla-Posadas Regional Aquifer (Southern Spain): Implications for Its Management. J. Hydrol. 2015, 523, 624–635. [Google Scholar] [CrossRef]
  49. Hamzaoui-Azaza, F.; Ketata, M.; Bouhlila, R.; Gueddari, M.; Riberio, L. Hydrogeochemical Characteristics and Assessment of Drinking Water Quality in Zeuss–Koutine Aquifer, Southeastern Tunisia. Environ. Monit. Assess. 2011, 174, 283–298. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Ijerph 15 01816 g001
Figure 2. Stratigraphic histogram for depths of 0–113 m.
Figure 2. Stratigraphic histogram for depths of 0–113 m.
Ijerph 15 01816 g002
Figure 3. Adsorption isotherm of clay for (a) Na+, (b) Ca2+ and (c) Mg2+ in groundwater at different depths. Note: Ce is the adsorption equilibrium mass concentration (mg L−1) of the adsorbate in the solution; qe is the equilibrium adsorption capacity of the adsorbent (mg kg−1).
Figure 3. Adsorption isotherm of clay for (a) Na+, (b) Ca2+ and (c) Mg2+ in groundwater at different depths. Note: Ce is the adsorption equilibrium mass concentration (mg L−1) of the adsorbate in the solution; qe is the equilibrium adsorption capacity of the adsorbent (mg kg−1).
Ijerph 15 01816 g003
Figure 4. The maximum buffering capacity (MBC) trend and the linear distribution coefficient of clay for Na+, Ca2+ and Mg2+ in groundwater. Note: Kd is the partition coefficient in the Henry equation (L kg−1).
Figure 4. The maximum buffering capacity (MBC) trend and the linear distribution coefficient of clay for Na+, Ca2+ and Mg2+ in groundwater. Note: Kd is the partition coefficient in the Henry equation (L kg−1).
Ijerph 15 01816 g004
Figure 5. Correlation between the linear distribution coefficient at different depths of clay and the total content of clay minerals.
Figure 5. Correlation between the linear distribution coefficient at different depths of clay and the total content of clay minerals.
Ijerph 15 01816 g005
Figure 6. (a) Correlation between the coefficient of retardation and the milliequivalent percentage of the ionic components; (b) Correlation between the sum of retardation coefficients and groundwater conductivity.
Figure 6. (a) Correlation between the coefficient of retardation and the milliequivalent percentage of the ionic components; (b) Correlation between the sum of retardation coefficients and groundwater conductivity.
Ijerph 15 01816 g006aIjerph 15 01816 g006b
Table 1. Basic parameters for the hydrogeological drilling.
Table 1. Basic parameters for the hydrogeological drilling.
Hydrological Observation Borehole NumberWX01WX02WX03WX04
Borehole depth/m12.049.295.0113.6
Filtration pipe depth/m8.0–12.013.0–49.279.0–94.098.0–107.9
Lithology at the depth of the Filtration pipeSilty soil, silty Clay clip siltSilty clay, silty Soil, silty sandSilty clay, ClaySilt, fine sand, Medium sand
Table 2. Basic physical and chemical properties of the clay samples.
Table 2. Basic physical and chemical properties of the clay samples.
Clay Sampling Depth/m6.0–6.248.0–48.251.0–51.271.0–71.284.0–84.297.0–97.2102.0–102.2112.0–112.2
Porosity Ratio0.6120.5350.5100.5120.5640.6450.4960.849
Mineral Amount/%Quartz55.034.3-36.247.452.621.137.7
Calcite0.012.6-16.78.31.154.925.7
Feldspar15.07.5-11.81715.10.115.5
Illite8.035.8-31.224.429.514.513.5
Montmorillonite12.02.0-0.40.40.48.90.4
Kaolinite10.07.8-3.72.51.30.57.2
Table 3. Milliequivalent percentage statistics for the experimental water chemical composition (n = 392).
Table 3. Milliequivalent percentage statistics for the experimental water chemical composition (n = 392).
Test CompositionCation Ratio/%Anion Ratio/%
K+Na+Ca2+Mg2+ClSO42−HCO3
WX010.0451.7111.5736.6857.7929.2312.98
WX020.0241.0915.0343.8649.8448.381.78
WX030.0441.5918.6639.7165.8131.482.71
WX040.0341.4627.131.4144.5747.38.13
Table 4. Eigenvalues and cumulative contributions for the factor analysis.
Table 4. Eigenvalues and cumulative contributions for the factor analysis.
Hydrological Observation BoreholeMain IngredientInitial Feature ValueRotation Square and Load
EigenvaluesVariance/%Cumulative Variance/%EigenvaluesVariance/%Cumulative Variance/%
WX0115.77564.16464.1645.77564.16464.164
22.81131.23595.3982.81131.23595.398
WX0217.15079.44679.4467.15079.44679.446
21.52516.94196.3881.52516.94196.388
WX0313.14534.94334.9433.14534.94334.943
22.01622.39757.3402.01622.39757.340
31.21113.45670.7971.21113.45670.797
WX0414.91854.64754.6474.91854.64754.647
21.31514.60769.2541.31514.60769.254
31.09212.13881.3931.09212.13881.393
Table 5. Rotational factor load matrix for the groundwater at different depths.
Table 5. Rotational factor load matrix for the groundwater at different depths.
Water Chemistry IndexWX01WX02WX03WX04
F1F2F1F2F1F2F3F1F2F3
Na+−0.1350.9730.9840.1740.908−0.0080.1570.850−0.121−0.004
Ca2+−0.9920.1090.967−0.0910.145−0.5040.6200.7010.425−0.357
Mg2+0.999−0.0170.9710.174−0.3230.8450.1050.790−0.1720.456
HCO3−0.9780.690−0.7800.5630.5240.3080.2950.671−0.6460.100
Cl0.9650.2570.9740.1900.242−0.730−0.317−0.730−0.0910.235
SO42−0.8630.4690.9460.2730.6540.2950.4100.9010.121−0.001
pH−0.3640.9210.141−0.9880.6890.462−0.2620.6880.510−0.237
EC0.203−0.825−0.9190.2120.911−0.250−0.168−0.0520.6100.765
Total Hardness0.9840.1790.9920.078−0.361−0.2380.5830.900−0.1340.224
Note: F1, F2, F3 represents main factors. EC represents conductivity.
Table 6. The adsorption isotherm model parameters of clay for Na+, Ca2+, Mg2+ in groundwater.
Table 6. The adsorption isotherm model parameters of clay for Na+, Ca2+, Mg2+ in groundwater.
Clay Sampling Depth/mIon in GroundwaterLangmuir Equation ParametersHenry Equation Parameters
qm/(mg·kg−1)KL/(L·mg−1)R2Kd/(L·mg−1)R2
6Na+4.600.00110.91640.29120.6311
Ca2+1.870.01500.74380.10200.4001
Mg2+2.720.02690.03940.23960.3851
48Na+1009.090.04530.71500.86060.0783
Ca2+466.670.00360.94080.83390.8511
Mg2+600.000.00130.98530.84510.9637
51Na+22.370.00140.75281.07390.8907
Ca2+909.090.00060.77180.91910.7074
Mg2+833.330.00030.88890.81320.7468
71Na+1000.000.00040.47640.48550.8182
Ca2+15,873.020.00010.92290.67380.9463
Mg2+909.090.00050.91160.37010.7336
84Na+277.780.00230.24210.72540.1848
Ca2+1250.000.00020.99760.02620.8260
Mg2+1250.000.00030.55740.65700.9237
97Na+416.670.00060.93820.53680.9011
Ca2+222.220.00070.45540.16360.5880
Mg2+769.230.00080.59400.68760.8955
102Na+400.000.00190.89420.45180.7132
Ca2+625.000.00030.63430.20700.8226
Mg2+11,235.960.00010.98360.46060.9069
112Na+192.310.00120.63521.50950.9251
Ca2+5000.000.00010.90310.21110.5134
Mg2+3333.330.00010.94030.37740.8569
Note: qm is the maximum adsorption capacity of the adsorbent in the Langmuir adsorption equation (mg kg−1); KL is the adsorption equilibrium constant (L mg−1) of the adsorbent in the Langmuir adsorption equation; R2 is correlation coefficient.
Table 7. The retardation coefficient of clay for Na+, Ca2+ and Mg2+ in groundwater.
Table 7. The retardation coefficient of clay for Na+, Ca2+ and Mg2+ in groundwater.
Sampling Depth/mRd-NaRd-CaRd-MgSum of Retardation Coefficients
61.24741.10534.02556.3782
482.01681.98521.98525.9873
515.59073.32092.228111.1397
713.06753.8692.84706.9365
843.46041.10123.53978.1013
972.81441.55293.32407.6913
1022.98591.90993.02457.9203
1124.87621.54201.96918.3873
Note: Rd-Na is the retardation coefficient of Na+. Rd-Ca is the retardation coefficient of Ca2+. Rd-Mg is the retardation coefficient of Mg2+.

Share and Cite

MDPI and ACS Style

Xing, L.; Huang, L.; Yang, Y.; Xu, J.; Zhang, W.; Chi, G.; Hou, X. The Blocking Effect of Clay in Groundwater Systems: A Case Study in an Inland Plain Area. Int. J. Environ. Res. Public Health 2018, 15, 1816. https://doi.org/10.3390/ijerph15091816

AMA Style

Xing L, Huang L, Yang Y, Xu J, Zhang W, Chi G, Hou X. The Blocking Effect of Clay in Groundwater Systems: A Case Study in an Inland Plain Area. International Journal of Environmental Research and Public Health. 2018; 15(9):1816. https://doi.org/10.3390/ijerph15091816

Chicago/Turabian Style

Xing, Liting, Linxian Huang, Yi Yang, Junxiang Xu, Wenjing Zhang, Guangyao Chi, and Xinyu Hou. 2018. "The Blocking Effect of Clay in Groundwater Systems: A Case Study in an Inland Plain Area" International Journal of Environmental Research and Public Health 15, no. 9: 1816. https://doi.org/10.3390/ijerph15091816

APA Style

Xing, L., Huang, L., Yang, Y., Xu, J., Zhang, W., Chi, G., & Hou, X. (2018). The Blocking Effect of Clay in Groundwater Systems: A Case Study in an Inland Plain Area. International Journal of Environmental Research and Public Health, 15(9), 1816. https://doi.org/10.3390/ijerph15091816

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