Next Article in Journal
National Surveillance of Injury in Children and Adolescents in the Republic of Korea: 2011–2017
Previous Article in Journal
Influence of Officially Ordered Restrictions during the First Wave of COVID-19 Pandemic on Physical Activity and Quality of Life in Patients after Kidney Transplantation in a Telemedicine Based Aftercare Program—A KTx360° Sub Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Groundwater Chemical Characteristics and Controlling Factors in a Region of Northern China with Intensive Human Activity

1
School of Chemical Engineering, Zhengzhou University, Zhengzhou 450001, China
2
School of Civil Engineering, Nanyang Institute of Technology, Nanyang 473004, China
3
Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang 050061, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(23), 9126; https://doi.org/10.3390/ijerph17239126
Submission received: 14 October 2020 / Revised: 3 December 2020 / Accepted: 4 December 2020 / Published: 7 December 2020

Abstract

:
The identification of groundwater chemical characteristics and the controlling factors is of major concern in water environment studies. In this study, we identified the groundwater chemical characteristics, evolution laws and main control factors in a region severely affected by human activities using hydrochemical and multivariate statistical techniques. The results showed that the concentrations of NO3 and TH were the primary pollution factors in the region with intensive human activity because of high concentration and over the standard rates. The major types of groundwater chemistry were HCO3·SO4–Ca and HCO3·SO4–Ca·Mg. The sulfate-type water was as high as 75.0%, 69.2% and 41.2% in the three hydrogeological units. In addition, there were Cl-type and Na-type waters, indicating that the groundwater in this area has been significantly affected by human activities. A principal component analysis (PCA) indicated that the three factors affecting groundwater hydrochemistry in the study area are domestic sewage and fertilizer, water–rock interactions and industrial wastewater. Therefore, we suggest that the government and water environment management departments should prevent the discharge of domestic and industrial wastewater without standardized treatment first in order to effectively prevent the further deterioration of groundwater quality in this area.

1. Introduction

Groundwater is one of the most vital resources for human living and production, especially in arid and semi-arid regions [1]. In recent years, with the rapid growth of China’s economy and population, processing capacity of sewage and garbage has not kept pace, and has led to unregulated discharge of sewage (especially in rural areas), the excessive application of pesticide and fertilizer, overexploitation of groundwater, and so on [2,3,4]. This has caused the deterioration of groundwater quality and changed groundwater chemical characteristics. It is well known that changes in the groundwater chemistry not only disrupt many ecological processes [5], but also affect groundwater use. Therefore, fully understanding the evolution law of groundwater chemistry is vital for the sustainable use of groundwater [6,7].
Previous studies found that the evolution of groundwater chemistry is mainly controlled by natural factors (including hydrogeological conditions, aquifer lithology, interactions of groundwater with rock, etc.) and human activities (such as domestic sewage, agricultural fertilizer and overexploitation of groundwater) [8,9]. In recent years, researchers have devoted to study the effect factors on groundwater chemistry [10,11,12]. Huang et al. [8] found that the main factors that control groundwater quality were seawater intrusion, water–rock interactions, lateral flow of river water, sewage discharge and agricultural pollution. Gu et al. [13] identified that the factors influencing groundwater quality were natural mineral dissolution and mine water mixing, domestic sewage and septic tank contamination and agricultural fertilizer contamination.
Selecting the correct analysis method is important for identifying the main factors that control groundwater hydrochemistry. Principal component analysis (PCA) is a mathematical procedure that uses an orthogonal transformation to convert a set of observations of possibly correlated variables into a set of linearly uncorrelated variables called principal components [3]. Therefore, we could identify the unobserved latent controlling factors of groundwater chemistry. Nowadays, PCA has been certified as a powerful tool for determining the controlling factors of groundwater hydrochemistry [14,15] and identifying pollution sources of groundwater quality [16,17].
The Hutuo River alluvial-pluvial fan is located in the western portion of the North China Plain. It is an important source of water in Shijiazhuang, and is one of the important receiving areas of the South–North Water Transfer Project. Thus, the groundwater quality directly affects the drinking water security of the residents in the area. However, with rapid urbanization, industrialization and population growth, the groundwater chemistry has been seriously affected by human activities (such as unregulated discharge of domestic sewage and industrial wastewater and overexploitation of groundwater) [18]. Therefore, it is significant for deepening the understanding of groundwater hydrochemical evolution to study the controlling factors of groundwater hydrochemical evolution in this region. The main aims of this study were to: (1) evaluate the characteristics of groundwater hydrochemistry in a region with intensive human activity; (2) analyze the major ion sources using hydrochemistry and multivariate statistical techniques; (3) identify main factors of controlling groundwater chemistry using PCA. The results could be helpful in developing effective water quality protection strategies and could deepen the understanding of groundwater hydrochemical evolution law in a region with intensive human activity.

2. Materials and Methods

2.1. Study Area

The study area is located in the Piedmont inclined plain at the eastern foot of Taihang Mountain, starting from Gangnan Reservoir in the west (the upstream recharge area of Hutuo River alluvial-proluvial fan) to Gaocheng in the east, and extending 5 km from Hutuohe River to the north and south sides. The geographical coordinates are 113°51′05″–114°56′00″ and 37°51′00″–38°24′40″. The administrative area includes Shijiazhuang City, Gaocheng District, Luquan District, Lingshou County, Zhengding County, Wuji County and Pingshan County, with a total area of 2442 km2 (Figure 1) and a total population of about 5.5 million. It has a temperate semi-humid and semi-arid monsoon climate, with average temperature and rainfall of 13.3–15.0 °C and 400–750 mm, respectively [1]. The topography is high in the west and low in the east. The landform is an alluvial-pluvial fan group formed by alluvial-pluvial action, and is a sloping plain formed by river and lake facies deposition. The western part is a hilly area with low and medium mountains, and the eastern part is flat North China Plain. The terrain slope is steeper in the west (5–6‰) and other areas are gentler (less than 2‰). The land use types in the study area include farmland (40%), city (11.60%), surface water (2.35%), town (13.09%), and wasteland (mainly including mountains and wasteland) (34.34%) (Figure 1).
The aquifer in the study area belongs to the Quaternary thick aquifer system of Hebei Plain. The stratigraphic lithology is mainly composed of mild clay, sandy loam, sandy mild clay and pebbles, gravel and sand with different particle sizes [19]. The western part of the aquifer is made of a thick layer of sandy gravel; the sandy gravel in the eastern part becomes finer and eventually transitions to mainly coarse sand (Figure 2). The aquifer in the study area is rich in water, and the flow direction of groundwater is the same as that of Hutuo River, which flows from the southwest to northeast. The groundwater mainly consists of bedrock fissure water and loose rock pore water, and the main recharge sources are rainfall vertical infiltration, surface water lateral seepage and flood lateral seepage in rainy season [20]. The study area can be divided into three hydrogeological units, namely the fissure pore water unit in the plain of Hutuo River valley between Gangnan and Huangbizhuang reservoir (HRVP), the hydrogeological unit at the top of Hutuo River alluvial-pluvial fan (THRVPF) and the hydrogeological unit in the central of Hutuo River alluvial-pluvial fan (CHRVPF). The HRVP is located in the transition zone from the mountain area to the plain area. The groundwater depth is shallow (between 2 and 21 m), and it has good conductivity and water richness. In the THRVPF area, the groundwater depth is about 12–40 m and the water conductivity is good. Currently, the main exploitation aquifer is the second water group due to the first aquifer group has been drained. In the CHRVPF area, the groundwater depth is about 40–50 m and the lithology of the aquifer is mainly coarse sand. The third and fourth water groups are the current exploitation aquifer. In recent years, due to the needs of local social and economic development, the groundwater has been intensively exploited, and two large reservoirs have been built in the upstream of the study area, resulting in the recharge source of groundwater being blocked. The depth of groundwater dropped rapidly. At present, one of the largest descending funnels in China has been formed [20].

2.2. Groundwater Sampling and Laboratory Analyses

In October 2015, 46 groups of groundwater samples were collected in the Hutuo River alluvial-pluvial fan area (Figure 1). The groundwater was taken from civil wells and agricultural irrigation wells, with a buried depth of 4.0–50.0 m. One 1.5 L and one 500 mL polyethylene plastic bottle were used for the analysis of anions and cations, respectively. The water samples for cation analysis were added to HNO3 to adjust the pH < 2. Before sampling, rinse the sampling bottles with purified water and original water sample for three times, respectively, and put the samples into the refrigerator for test.
The following 16 water chemical parameters were analyzed: pH, total dissolved solids (TDSs), potassium (K+), sodium (Na+), calcium (Ca2+), magnesium (Mg2+), nitrate (NO3), nitrite (NO2), ammonia (NH3+), bicarbonate (HCO3), carbonate (CO32), chloride (Cl), sulfate (SO42−), iron (Fe), manganese (Mn) and total hardness (TH).
The values of pH were measured in the field using a HQ40D multiparameter instrument (United States of America). Concentrations of the cations and trace elements (Fe and Mn) in the water samples were measured using an inductively coupled plasma atomic emission spectrometer (Agilent 7500ce ICP-MS, Tokyo, Japan), while anion analyses were carried out using spectrophotometry (Perkin-Elmer Lambda 35, Waltham, MA, USA). TDSs were measured using gravimetric methods, HCO3− was measured using acid–base titration and TH was determined using the ethylene diamine tetraacetic acid titration method. All analyses were carried out at the Groundwater Mineral Water and Environmental Monitoring Center at the Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences.

2.3. Data Analysis

Since the data obtained in this study did not satisfy the normal distribution, we used the non-parametric Kruskall–Wallis and Mann–Whitney U test methods to analyze the significant differences of the concentration of groundwater hydrochemical parameters in different hydrogeological units and land use types. The Spearman correlation analysis was used to identify the correlation between hydrochemical indicators. The principal component analysis was used to identify the main factors of controlling groundwater hydrochemistry. R (3.6.2) (Core Team, Vienna, Austria) and origin (9.0)(OriginLab Corporation, Northampton, MA, USA) were used for data analysis.

2.4. Multivariate Data Analysis

Contrast Coefficient Variance

Variance is used to describe the dispersion degree of a random variable to its mathematical expectation (mean). Since the absolute value ranges of the chemical components of groundwater may be very different, it is not possible to compare the variances. Therefore, the ratio of each variable (Xi) of each sample to the mean value of the variable (Xm) is calculated first, which is defined as the contrast coefficient v, Vi = Xi/(Xm) Since the average value of the contrast coefficient is 1, each variable is converted into another variable with the same mean (its contrast coefficient), and then the variance of the contrast coefficient of each variable is calculated to compare the variance.
The calculation formula of contrast coefficient variance (Vσi2) is as follows:
V σ i 2   =   1 n i = 1 n ( Vi Vmi ) 2
where: Vi is the contrast coefficient of the i-th chemical component of groundwater; Vmi is the average value of the contrast coefficient of the i-th chemical component of groundwater.

3. Results

3.1. Main Ionic Characteristics

The groundwater hydrochemical characteristics of the Hutuo River alluvial-pluvial fan are shown in Table 1. The pH of the groundwater in the study area is 7.38–8.40, and the average values in three hydrogeological units (valley plain area, top of alluvial-pluvial fan and central of alluvial-pluvial fan) are 7.73, 7.80 and 7.63, respectively. There is no excessive water sample, and the groundwater is weak alkaline water. The TDSs value is 152.2–1913 μs/cm. The mean values and over-standard rates of the three hydrogeological units are 950.4 μs/cm and 25%, 720.4 μs/cm and 11.5%, 412.8 μs/cm and 0%, respectively. The main cation in groundwater is Ca2+, with a concentration between 20.6 and 345.4 mg/L. The average values of Ca2+ in three hydrogeological units are 211.2, 143.8 and 69.8 mg/L, respectively. The order of cation concentration is Ca2+ > Na+ > Mg2+ >K+. The concentration of NH4+ is below the detection limit (0.04 mg/L) in three hydrogeological units in the study area. The main anion of groundwater is HCO3. Its concentration range is between 141.0 and 429.3 mg/L. The mean concentration of HCO3 in three hydrogeological units is 315.7, 311.7 and 249.3 mg/L, respectively. The order of anion concentration is: HCO3 > SO42− > NO3 > Cl > NO2 in the valley plain area, and HCO3 > SO42− > Cl > NO3 > NO2 in the top and middle of alluvial-pluvial fan. CO32− was not detected in the three hydrogeological units in the study area. It is worth noting that the concentration of NO3 is 1.76–465.0 mg/L, and the average concentrations in the three hydrogeological units are 176.2, 75.8 and 9.10 mg/L, with over-the-standard (the Grade III standard for groundwater quality in China [21]) rates of 62.5%, 34.6% and 0%, respectively. It is obvious that groundwater in the valley plain area and the top of the alluvial-pluvial fan has been seriously polluted by NO3. The concentration range of TH is 76.1–1058 mg/L. The average concentrations of TH in the three hydrogeological units are 650.6, 507.3 and 278.0 mg/L, respectively. The over-the-standard (the Grade III standard for groundwater quality in China [21]) rates are 87.5%, 73.1% and 0%, which is the most important pollution factor in the groundwater in this area. In addition, the concentration range of Fe is 0.01–1.313 mg/L. The average concentrations of Fe in the three hydrogeological units are 0.148, 0.067 and 0.214 mg/L, respectively, and the over-the-standard rates are 12.5%, 3.8% and 8.3%, respectively.
The spatial variation of hydrochemical indicators is relatively obvious: Hutuo River valley plain area > top of alluvial-pluvial fan > central of alluvial-pluvial fan, which may be related to the groundwater level depth, the lithology of aeration zone and the intensity of human activities [1]. From the valley plain area to the middle of the alluvial-pluvial fan, the buried depth of the groundwater level increases gradually, the particles in the aeration zone become thinner and the permeability becomes poor, which makes it difficult for pollutants to enter the groundwater and the concentrations of the hydrochemical indexes are relatively low.

3.2. Hydrochemical Type

It can be seen from Figure 3 that the main hydrochemical types of groundwater in the study area are HCO3·SO4-Ca and HCO3·SO4-Ca·Mg, accounting for 21.7% and 28.3%, respectively. The proportions of SO4-type water in the three hydrogeological units (valley plain area, top and middle of alluvial-pluvial fan) are 75.0%, 69.2% and 41.2%, respectively. In addition, Cl-type water also accounts for certain proportions in the three hydrogeological units, which are 50.0%, 26.9% and 8.3%, respectively. In this study, the proportions of HCO3-Ca(Mg)-type in the three hydrogeological units are 0%, 11.5% and 41.7%, respectively.

3.3. Analysis of Hydrochemistry

A Gibbs diagram is often used to reflect the main control effects of groundwater hydrochemical types (rock weathering, evaporation concentration and precipitation) [22]. The water samples that are greatly affected by atmospheric precipitation are located in the lower right corner of Gibbs figure; the samples that are mainly controlled by rock weathering are located in the middle left corner and the samples; those that are controlled by evaporation crystallization are located in the upper right corner [23]. It can be seen from Figure 4 that most of the water samples from different hydrogeological units fall in the rock weathering area, indicating that rock weathering is the main control factor of groundwater hydrochemistry in this area, and it is less affected by evaporative concentration and atmospheric precipitation.

4. Discussion

4.1. Groundwater Chemical Characteristics in a Region with Intensive Human Activity

Groundwater chemical characteristics are influenced by long-term geological history, seepage migration process and human activities [24]. Nowadays, China has experienced rapid urbanization and industrialization, and the influence of human activities on groundwater chemistry is increasing year by year. In the Hutuo River alluvial-pluvial fan region, the groundwater has been intensively exploited and has formed one of the largest descending funnels in China [20]. In addition, the main land use type is farmland in this area—thus, agricultural fertilization and pesticide spraying had serious impact on groundwater chemistry [25]. According to our research, TH and NO3 of the groundwater had significantly exceeded the Grade III standard for groundwater quality in China. Thus, TH and NO3 of groundwater in the study area have been significantly affected by human activities, such as domestic sewage and agricultural fertilizer [25,26]. This result is consistent with previous research. Zhang et al. (2019) found that NO3 was the main impact indicators for poor-quality groundwater in urbanized areas [27]. González Pérez et al. (2020) also found that intensive human activities have serious affected the level of NO3 in groundwater [28]. Egbi et al. (2020) also confirmed that the sources of NO3 in groundwater primarily originate from human contributions in the Lower Volta River basin of Ghana [29]. Therefore, TH and NO3 were the most important parameters, indicating that groundwater is affected by human activities.
SO42− is also a vital parameter that reflected the effect of human activities on groundwater chemistry [30]. Indeed, in this study region, especially in the HRVP area, the mean concentration and over-the-standard (the Grade III standard for groundwater quality in China) rate of groundwater SO42− were 176.7% and 12.5%, respectively. Previous research has confirmed that SO42− was an important indicator in highly anthropogenic regions, which can evaluate groundwater vulnerability and pollution risk [31]. Zhang et al. (2020) also found that the contribution proportion of sewage to SO42− in groundwater was 35.5%–42.7% in the Hutuo River basin [18]. Therefore, SO42− is also a vital indicator representing the effects of human activity.
In addition, the main hydrochemical type of groundwater in this study area was HCO3·SO4-Ca(Mg). However, previous studies have found that, in 1950, the hydrochemical type of groundwater in the Shijiazhuang area was mainly HCO3-Ca(Mg) [24]. In recent years, the hydrochemical type of groundwater has changed to HCO3·SO4-Ca(Mg), and there is Cl-type water and Na-type water. It can be seen that the groundwater in this area has been significantly affected by human activities, which could influence the safety of the drinking water.

4.2. Analysis of Major Ion Sources

4.2.1. Correlation Analysis

If there is a strong correlation between the chemical components in the water environment, it means that they have a common material source. Based on the correlation analysis of groundwater quality parameters in the study area (Table 2), we found that pH had a significant correlation with other indicators. However, TH and TDS showed significant positive correlations with Na+, Ca2+, Mg2+, NO3, Cl, HCO3 and SO42−. It indicated that these components may come from common material sources, such as domestic sewage, agricultural fertilizer, industrial wastewater [25,29]. Because it was in the rainy season, during the process of rainwater recharge, a large number of pollutants from the surface were carried and infiltrated into the aquifer, which promoted water–rock interactions, and increased the concentration of groundwater ions [1]. K+ was only positively correlated with Na+. Na+, Ca2+, Mg2+, Cl, NO3, HCO3 and SO42− all showed significant positive correlations, indicating that they came from the same source, such as agriculture runoff, urban runoff and water–rock interactions [8,25]. The correlation between Fe and Mn was very significant, but they did not show correlation with other indexes, which proved that they have a common source.

4.2.2. Source Judgment by Using Contrast Coefficient Variance

The factors that control the chemical composition of groundwater include natural conditions (such as water–rock interactions) and human activities. Under natural conditions, the chemical composition of groundwater will not change much in a complete groundwater system. However, the hydrochemical composition of a groundwater system affected by human activities often shows a great degree of change [32]. Due to the different impact characteristics on the chemical components of groundwater, the contrast coefficient variance is used in this study to distinguish whether the formation of the chemical components of groundwater is mainly dominated by the evolution of natural conditions or by human activities, and to compare the degree of influence of each component by human activities.
On the whole, the variance values of the contrast coefficients of pH, TDS, TH, Ca2+, Mg2+ and HCO3 in the three hydrogeological units are small (Table 3). They are mainly controlled by the evolution of natural hydrochemistry. Because of the shallow groundwater depth, the Hutuo River valley plain unit is greatly affected by human activities and the variance value of the contrast coefficient of each component in groundwater is greater than that of the top and middle hydrogeological units of the Hutuo River alluvial-pluvial fan. It showed that: in the valley plain area, σ2 (K+), σ2(Na+), σ2(Cl), σ2(NO3), σ2(Fe) and σ2(Mn) have large values (σ2 > 0.5), which indicated that these components are greatly affected by human activities. In the top of the alluvial-pluvial fan, the values of σ2(Na+), σ2(Cl), σ2(Mn) and σ2(Fe) are higher (σ2 > 0.5); in the middle of the alluvial-pluvial fan, the values of σ2(Na+), σ2(Cl), σ2(NO3) and σ2(Fe) exceeded 0.5, demonstrating that these components are significantly affected by human activities. In the three hydrogeological units, the variance value of the SO42− contrast coefficient is moderately low—0.199, 0.245 and 0.310, respectively—illustrating that it is less affected by human activities and mainly controlled by the evolution of natural hydrochemistry.

4.3. Controlling Factors of Groundwater Chemical

In order to prevent groundwater pollution, the main influencing factors must be identified accurately. In this study, we selected 13 water quality parameters (pH, TDS, TH, K+, Na+, Ca2+, Mg2+, NO3, Cl, SO42−, HCO3, Fe and Mn), and identified the main control factors of groundwater pollution in the Hutuo River alluvial-pluvial fan area by PCA. Before principal component analysis, a Kaiser–Meyer–Olkin (KMO) test and Barlett spherical test were carried out on all the data. The results showed that the KMO value was 0.713 and Barlett spherical test value was 963.07 (p < 0.001), indicating that the original data were suitable for PCA. Based on the eigenvalue being greater than 1, three main control factors (Figure 5) that caused groundwater pollution in the Hutuo River alluvial-proluvial fan area were identified. These factors can explain 74.01% of all variables, indicating that they basically contain all the information of the above 13 indicators.
The contribution rate of the first main factor (FC1) was TH, Ca2+, NO3, TDS, Na+, K+ and Cl, showing a strong positive correlation with FC1, and SO42− showed a moderate positive correlation with PC1. It can be seen that PC1 reflects that groundwater is mainly affected by human activities. The possible sources of Cl in nature were animal manure, domestic sewage, rising saltwater and road snow melting salt [33]. Since there was no rising saltwater that intruded into the aquifer in the study area, and it was in autumn without snowfall, it can be surmised that neither rising saltwater nor snowmelt salt were the main sources of Cl in the groundwater in the area. Therefore, the higher concentration of Cl is most likely to come from animal manure and domestic sewage. The main sources of NO3 included domestic sewage, agricultural fertilizers, manure and atmospheric deposition [34]. Previous studies have found that the nitrate pollution of groundwater is mainly related to the use of agricultural fertilizers, discharge of domestic sewage and sewage irrigation [35,36]. During the investigation and sampling, it was found that in the villages at the top of the Hutuo River alluvial-pluvial fan, a large amount of domestic sewage was directly discharged into the Hutuo River channel. In addition, the lithology of the stratum at the top of the fan was relatively coarse, and the sewage was easily infiltrated into the aquifer. As a result, domestic sewage was the main influencing factor of the groundwater in the study area. Furthermore, the main land use type was agricultural land (38.62%), so agricultural fertilizers were also an important source. In conclusion, we can conclude that the first factor (PC1) controlling groundwater chemistry in this area is human activity (domestic sewage and agricultural fertilizer).
The contribution rate of the second main factor (FC2) was 13.78%. HCO3 and Mg2+ showed a strong positive correlation with PC2, pH showed a medium negative correlation with PC2 and SO42− showed a medium positive correlation with PC2. In general, HCO3, Mg2+ and SO42− in groundwater were mainly from natural sources (dissolution of dolomite and limestone) [33]. Based on the previous analysis of the contrast coefficient, HCO3, Ph, Mg2+ and SO42− in the groundwater of the study area were mainly controlled by the chemical evolution of natural water. In recent years, with the development of industry and agriculture, the groundwater in the Hutuo River alluvial-pluvial fan area has been intensively exploited. In addition, two large reservoirs have been built in the upstream of the study area, which cut off the recharge source of groundwater in the area, and continuously reduced the groundwater level. A depression funnel was formed in some parts, breaking the original balance of the groundwater dynamic field, thus triggering a series of physicochemical reactions. The equilibrium is conducive to the dissolution of carbonate minerals (calcite and dolomite), and then the concentration of HCO3, Ca2+ and Mg2+ increased [1]. Therefore, the second factor controlling groundwater in this area is a natural factor (water–rock interaction).
The contribution rate of the third main factor (FC3) is 8.77%. Fe and Mn have strong positive correlations with PC3. In the Hutuo River alluvial-pluvial fan area, the average concentrations of Fe in the groundwater of three hydrogeological units are 0.148, 0.067 and 0.214 mg/L, respectively, and there were four samples that exceeded the national groundwater quality standard of 0.3 mg/L. Previous studies have found that high concentrations of Fe in the water environment were mainly from industrial pollution [37,38]. There was a steel factory upstream of the study area, and the sampling point near it (5) has the highest Fe concentration (0.552 mg/L). The sampling point near Gaocheng industrial park (35) downstream of the study area had an Fe concentration of 1.313 mg/L, which further confirms our hypothesis. Therefore, PC3 shows that the groundwater in the study area is polluted by industrial sewage. Thus, the third factor controlling the groundwater in the study area is the human activity (industrial sewage).

5. Managerial Implications and Future Research

The results indicate that high concentrations of NO3 and TH in groundwater are present in the region with intensive human activity. Local governments and environmental managers should pay due attention to NO3, which can cause not only eutrophication and toxic algal blooms [39], but also increase the risk of diseases, such as methemoglobinemia, blue baby syndrome and stomach cancer [40]. Therefore, groundwater in the region must be treated before drinking. In addition, the main factors controlling the groundwater hydrochemistry in this area are domestic sewage, industrial sewage and agricultural fertilizer. Thus, it is suggested that the government and water environment management department should first stop the substandard discharge of domestic and industrial sewage, and promote optimization of fertilization strategies (such as formula fertilization based on soil testing) in order to effectively prevent the further deterioration of groundwater quality in this area.
In this study, we have identified the main pollutants and factors controlling groundwater hydrochemistry in the region with intensive human activity. However, we have not quantitatively identified the source of the various pollutants. Nowadays, a quantitative traceability technique combined with isotope and source apportionment models (SIAR model) has been widely used in water environments. For example, previous researchers combined the isotope NO315N and δ18O) and SIAR model to identify the NO3 sources [41], and used the isotope SO42−34S and δ18O) and SIAR model to identify the SO42− sources [42]. Therefore, in order to formulate effective controlling measures for controlling groundwater pollution, it is necessary to quantitatively identify the sources of pollutants in groundwater in this region.

6. Conclusions

In this research, we applied hydrochemistry (Piper and Gibbs diagrams) and multivariate statistical techniques (variance of contrast coefficient, correlation analysis and PCA) to study the hydrochemistry characteristics and the main control factors of groundwater in the region with intensive human activity—the Hutuo River alluvial-pluvial fan of northern China. The results showed that the groundwater quality was relatively poor, among which the concentration range of NO3 and TH are 1.76–465.0 mg/L and 76.1–1058 mg/L, respectively. The over-the-standard rates in three hydrogeological units are 62.5%, 34.6% and 0%, and 87.5%, 73.1% and 0%, respectively, which were the most important parameters indicating that groundwater is affected by human activities. Rock weathering is the main control function of groundwater hydrochemistry in this area. The main types of groundwater hydrochemistry are HCO3·SO4–Ca and HCO3·SO4–Ca·Mg, accounting for 21.7% and 28.3%, respectively. SO4-type water in three hydrogeological units is as high as 75.0%, 69.2% and 41.2%, respectively. In addition, there are Cl- and Na-type waters, indicating that the groundwater in this area has been significantly affected by human activities. The variance analysis of contrast coefficient shows that K+, Na+, NO3, Cl, Fe and Mn are greatly influenced by human activities, and the other indexes are mainly influenced by natural causes.
The results of PCA show that the three factors controlling groundwater hydrochemistry in the study area are domestic sewage and chemical fertilizer (human activity), water–rock interactions (natural factor) and industrial sewage (human activity). In summary, the main factors controlling the groundwater hydrochemistry in this area are domestic sewage, industrial sewage and agricultural fertilizer. Therefore, in order to effectively prevent the further deterioration of groundwater quality in this area, it is suggested that the government and the water environment management department should first stop the substandard discharge of domestic and industrial sewage and promote optimization of fertilization strategies. In addition, it is necessary to quantitatively identify the sources of pollutants in the future in order to formulate effective controlling measures for controlling groundwater pollution.

Author Contributions

Conceptualization, Q.Z. and C.R.; methodology, Q.Z.; software, C.R.; validation, Q.Z.; formal analysis, C.R.; investigation, C.R., and Q.Z.; resources, C.R.; data curation, C.R.; writing—original draft preparation, C.R., and Q.Z.; writing—review and editing, Q.Z.; visualization, C.R.; supervision, Q.Z.; project administration, Q.Z.; funding acquisition, Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Projects of China Geological Survey, China, grant number DD20190331, the National Key R&D Program of China, grant number 2018YFC1803001, the Belt and Road Fund on Water and Sustainability, China, grant number U2019nkms01, and Open Foundation of the State Key Laboratory of Urban and Regional Ecology of China, grant number SKLURE2019-2-3.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Q.; Wang, H.; Wang, Y.; Yang, M.; Zhu, L. Groundwater quality assessment and pollution source apportionment in an intensely exploited region of northern China. Environ. Sci. Pollut. Res. 2017, 24, 16639. [Google Scholar] [CrossRef]
  2. Jiang, Y.; Wu, Y.; Groves, C.; Yuan, D.; Kambesis, P. Natural and anthropogenic factors affecting the groundwater quality in the Nandong karst underground river system in Yunan, China. J. Contam. Hydrol. 2009, 109, 49. [Google Scholar] [CrossRef]
  3. Qin, R.; Wu, Y.; Xu, Z.; Xie, D.; Zhang, C. Assessing the impact of natural and anthropogenic activities on groundwater quality in coastal alluvial aquifers of the lower Liaohe River Plain, NE China. Appl. Geochem. 2013, 31, 142. [Google Scholar] [CrossRef]
  4. Xue, D.; Botte, J.; Baets, B.D.; Accoe, F.; Nestler, A.; Taylor, P.; Cleemput, O.V.; Berglund, M.; Boeckx, P. Present limitations and future prospects of stable isotope methods for nitrate source identification in surface- and groundwater. Water Res. 2009, 43, 1159. [Google Scholar] [CrossRef]
  5. Moore, W.S. The subterranean estuary: A reaction zone of ground water and sea water. Mar. Chem. 1999, 65, 111–125. [Google Scholar] [CrossRef]
  6. Tziritis, E.; Skordas, K.; Kelepertsis, A. The use of hydrogeochemical analyses and multivariate statistics for the characterization of groundwater resources in a complex aquifer system. A case study in Amyros River basin, Thessaly, central Greece. Environ. Earth Sci. 2016, 75, 339. [Google Scholar] [CrossRef]
  7. Pazand, K.; Khosravi, D.; Ghaderi, M.R.; Rezvanianzadeh, M.R. Identification of the hydrogeochemical processes and assessment of groundwater in a semi-arid region using major ion chemistry: A case study of Ardestan basin in Central Iran. Groundw. Sustain. Dev. 2018, 6, 245–254. [Google Scholar] [CrossRef]
  8. Huang, G.; Sun, J.; Zhang, Y.; Chen, Z.; Liu, F. Impact of anthropogenic and natural processes on the evolution of groundwater chemistry in a rapidly urbanized coastal area, South China. Sci. Total Environ. 2013, 463, 209. [Google Scholar] [CrossRef] [PubMed]
  9. Lin, C.Y.; Abdullah, M.H.; Praveena, S.M.; Yahaya, A.H.B.; Musta, B. Delineation of temporal variability and governing factors influencing the spatial variability of shallow groundwater chemistry in a tropical sedimentary island. J. Hydrol. 2012, 432, 26. [Google Scholar] [CrossRef]
  10. Li, X.; Tang, C.; Cao, Y.; Li, D. A multiple isotope (H, O, N, C and S) approach to elucidate the hydrochemical evolution of shallow groundwater in a rapidly urbanized area of the Pearl River Delta, China. Sci. Total Environ. 2020, 724, 137930. [Google Scholar] [CrossRef]
  11. Güler, C.; Kurt, M.A.; Alpaslan, M.; Akbulut, C. Assessment of the impact of anthropogenic activities on the groundwater hydrology and chemistry in Tarsus coastal plain (Mersin, SE Turkey) using fuzzy clustering, multivariate statistics and GIS techniques. J. Hydrol. 2012, 414–415, 435–451. [Google Scholar]
  12. Li, Z.; Yang, Q.; Yang, Y.; Ma, H.; Wang, H.; Luo, J.; Bian, J.; Martin, J.D. Isotopic and geochemical interpretation of groundwater under the influences of anthropogenic activities. J. Hydrol. 2019, 576, 685–697. [Google Scholar] [CrossRef]
  13. Gu, H.; Chi, B.; Li, H.; Jiang, J.; Qin, W.; He, W. Assessment of groundwater quality and identification of contaminant sources of Liujiang basin in Qinhuangdao, North China. Environ. Earth Sci. 2015, 73, 6477. [Google Scholar] [CrossRef]
  14. Boateng, T.K.; Opoku, F.; Acquaah, S.O.; Akoto, O. Groundwater quality assessment using statistical approach and water quality index in Ejisu-Juaben Municipality, Ghana. Environ. Earth Sci. 2016, 75, 489. [Google Scholar] [CrossRef]
  15. Herojeet, R.; Rishi, M.S.; Lata, R.; Sharma, R. Application of environmetrics statistical models and water quality index for groundwater quality characterization of alluvial aquifer of Nalagarh Valley, Himachal Pradesh, India. Sustain. Water Resour. Manag. 2016, 2, 39. [Google Scholar] [CrossRef] [Green Version]
  16. Gholizadeh, M.H.; Melesse, A.M.; Reddi, L. Water quality assessment and apportionment of pollution sources using APCS-MLR and PMF receptor modeling techniques in three major rivers of South Florida. Sci. Total Environ. 2016, 566–567, 1552. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, Q.; Wang, X.; Wan, W.; Hou, P.; Li, R.; Ouyang, Z. The Spatial-Temporal Pattern and Source Apportionment of Water Pollution in a Trans-Urban River. Pol. J. Environ. Stud. 2015, 24, 841. [Google Scholar]
  18. Zhang, Q.Q.; Wang, H.W.; Lu, C. Tracing sulfate origin and transformation in an area with multiple sources of pollution in northern China by using environmental isotopes and Bayesian isotope mixing model. Environ. Pollut. 2020, 265, 115105. [Google Scholar] [CrossRef]
  19. Foster, S.; Garduno, H.; Evans, R.; Olson, D.; Tian, Y.; Zhang, W.; Han, Z. Quaternary aquifer of the north china plain—Assessing and achieving groundwater resource sustainability. Hydrogeol. J. 2004, 12, 81–93. [Google Scholar] [CrossRef]
  20. Li, Y.; Zhang, Z.; Fei, Y.; Chen, H.; Qian, Y.; Yu, D. Investigation of quality and pollution characteristics of groundwater in the Hutuo River Alluvial Plain, North China Plain. Environ. Earth Sci. 2016, 75, 581. [Google Scholar] [CrossRef]
  21. Ministry of Natural Resources of the People’s Republic of China. Standard for Groundwater Quality, (GB/T14848-2017); General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China: Beijing, China, 2017. (In Chinese)
  22. Li, P.; Wu, J.; Qian, H. Assessment of groundwater quality for irrigation purposes and identification of hydrogeochemical evolution mechanisms in Pengyang County, China. Environ. Earth Sci. 2013, 69, 2211. [Google Scholar] [CrossRef]
  23. Tang, J.P.; Zhang, Q.; Hu, Y. Groundwater chemical characteristics and analysis of their controlling factors in an alluvial fan of Jianjiang River. Environ. Sci. 2019, 40, 3089. (In Chinese) [Google Scholar]
  24. Zhang, Q.; Miao, L.; Wang, H.; Hou, J.; Li, Y. How Rapid Urbanization Drives Deteriorating Groundwater Quality in a Provincial Capital of China. Pol. J. Environ. Stud. 2020, 29, 441. [Google Scholar] [CrossRef]
  25. Zhang, Q.; Wang, H. Assessment of sources and transformation of nitrate in the alluvial-pluvial fan region of north China using a multi-isotope approach. J. Environ. Sci. 2020, 89, 9. [Google Scholar] [CrossRef] [PubMed]
  26. Meghdadi, A.; Javar, N. Quantification of spatial and seasonal variations in the proportional contribution of nitrate sources using a multi-isotope approach and Bayesian isotope mixing model. Environ. Pollut. 2018, 235, 207. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, F.E.; Huang, G.X.; Hou, Q.X.; Liu, C.Y.; Zhang, Y.; Zhang, Q. Groundwater quality in the Pearl River Delta after the rapid expansion of industrialization and urbanization: Distributions, main impact indicators, and driving forces. J. Hydrol. 2019, 577, 124004. [Google Scholar] [CrossRef]
  28. González Pérez, D.M.; Martín Martín, J.M.; Guaita Martínez, J.M.; Sáez-Fernández, F.J. An Analysis of the Cost of Water Supply Linked to the Tourism Industry. An Application to the Case of the Island of Ibiza in Spain. Water 2020, 12, 2006. [Google Scholar] [CrossRef]
  29. Egbi, C.D.; Anornu, G.A.; Ganyaglo, S.Y.; Appiah-Adjei, E.K.; Li, S.L.; Dampare, S.B. Nitrate contamination of groundwater in the Lower Volta River Basin of Ghana: Sources and related human health risks. Ecotoxicol. Environ. Saf. 2020, 191, 110227. [Google Scholar] [CrossRef]
  30. Ma, M.D.C.; Esteller, M.V.; Exposito, J.L.; Hirata, R. Impacts of urbanization on groundwater hydrodynamics and hydrochemistry of the Toluca Valley aquifer (Mexico). Environ. Monit. Assess. 2014, 186, 2979. [Google Scholar]
  31. Busico, G.; Kazakis, N.; Colombani, N.; Mastrocicco, M.; Voudouris, K.; Tedesco, D. A modified SINTACS method for groundwater vulnerability and pollution risk assessment in highly anthropized regions based on NO3 and SO42− concentrations. Sci. Total Environ. 2017, 609, 1512–1523. [Google Scholar] [CrossRef] [PubMed]
  32. Meng, F.A.; Liang, X.J.; Hao, Y.; Wang, Y.H.; Lou, Y.; Li, H.W. Dynamic Characteristics of Groundwater in Taoer River Fan. Water Sav. Irrig. 2016, 4, 65. (In Chinese) [Google Scholar]
  33. Liu, C.Q.; Li, S.L.; Lang, Y.C.; Xiao, H.Y. Using delta15N- and delta18O-values to identify nitrate sources in karst ground water, Guiyang, southwest China. Environ. Sci. Technol. 2006, 40, 6928. [Google Scholar] [CrossRef] [PubMed]
  34. Jin, Z.; Xue, Q.; Chen, L.; Jin, M.; Li, F. Using dual isotopes to evaluate sources and transformations of nitrate in the West Lake watershed, eastern China. J. Contam. Hydrol. 2015, 177, 64. [Google Scholar] [CrossRef] [PubMed]
  35. Bohlke, J.K. Groundwater recharge and agricultural contamination. Hydrogeol. J. 2002, 10, 153. [Google Scholar] [CrossRef]
  36. Chen, J.; Tang, C.; Sakura, Y.; Yu, J.; Fukushima, Y. Nitrate pollution from agriculture in different hydrogeological zones of the regional groundwater flow system in the North China Plain. Hydrogeol. J. 2005, 13, 481. [Google Scholar] [CrossRef]
  37. Bu, H.; Tan, X.; Li, S.; Zhang, Q. Temporal and spatial variations of water quality in the Jinshui River of the South Qinling Mts., China. Ecotoxicol. Environ. Saf. 2010, 73, 907. [Google Scholar] [CrossRef]
  38. Juahir, H.; Zain, S.M.; Yusoff, M.K.; Hanidza, T.T.; Armi, A.M.; Toriman, M.E.; Mokhtar, M. Spatial water quality assessment of Langat River Basin (Malaysia) using environmetric techniques. Environ. Monit. Assess. 2011, 173, 625. [Google Scholar] [CrossRef] [Green Version]
  39. Nestler, A.; Berglund, M.; Accoe, F.; Duta, S.; Xue, D.; Boeckx, P.; Taylor, P. Isotopes for improved management of nitrate pollution in aqueous resources: Review of surface water field studies. Environ. Sci. Pollut. Res. 2011, 18, 519–533. [Google Scholar] [CrossRef]
  40. Past en-Zapata, E.; Ledesma-Ruiz, R.; Harter, T.; Ramı´rez, A.I.; Mahlknecht, J. Assessment of sources and fate of nitrate in shallow groundwater of an agricultural area by using a multi-tracer approach. Sci. Total Environ. 2014, 470–471, 855–864. [Google Scholar] [CrossRef] [Green Version]
  41. Guo, Z.F.; Yan, C.Z.; Wang, Z.S.; Xu, F.F.; Yang, F. Quantitative identification of nitrate sources in a coastal peri-urban watershed using hydrogeochemical indicators and dual isotopes together with the statistical approaches. Chemosphere 2020, 243, 125364. [Google Scholar] [CrossRef]
  42. Torres-Martínez, J.A.; Mora, A.; Knappett, P.S.K.; Ornelas-Soto, N.; Mahlknecht, J. Tracking nitrate and sulfate sources in groundwater of an urbanized valley using a multi-tracer approach combined with a Bayesian isotope mixing model. Water Res. 2020, 182, 115962. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Monitoring sites of groundwater in Hutuo River alluvial-pluvial fan.
Figure 1. Monitoring sites of groundwater in Hutuo River alluvial-pluvial fan.
Ijerph 17 09126 g001
Figure 2. Hydrogeological cross-section of Hutuo River pluvial fan.
Figure 2. Hydrogeological cross-section of Hutuo River pluvial fan.
Ijerph 17 09126 g002
Figure 3. Piper diagram of groundwater in Hutuo River alluvial-pluvial fan.
Figure 3. Piper diagram of groundwater in Hutuo River alluvial-pluvial fan.
Ijerph 17 09126 g003
Figure 4. Gibbs diagram for the groundwater in Hutuo River alluvial-pluvial fan.
Figure 4. Gibbs diagram for the groundwater in Hutuo River alluvial-pluvial fan.
Ijerph 17 09126 g004
Figure 5. Loadings of 13 selected variables on varimax rotated factors in the Hutuo River alluvial-pluvial fan.
Figure 5. Loadings of 13 selected variables on varimax rotated factors in the Hutuo River alluvial-pluvial fan.
Ijerph 17 09126 g005
Table 1. Descriptive statistics of groundwater chemical parameters in three hydrogeological units in Hutuo River alluvial-pluvial fan.
Table 1. Descriptive statistics of groundwater chemical parameters in three hydrogeological units in Hutuo River alluvial-pluvial fan.
ParameterspHTDSK+Na+Ca2+Mg2+HCO3ClSO42−NO3NO2FeMnTH
(I) Hutuo River
valley plain
unit (n = 8)
Min7.64500.70.610.1135.210.23187.051.0962.468.80.0020.0170.001379.8
Max7.88191322.0180.8345.453.1381.4322.4322.3465.00.0070.5520.0121058
Mean7.73950.45.8050.9211.729.7315.7121.5176.7176.20.0030.1480.004650.6
SD0.08455.67.6656.168.315.870.099.378.7139.70.0020.1800.004219.7
Exceed standard rate (%)025012.512.562.5012.5087.5
(II) Top of the
Hutuo River
pluvial fan
(n = 26)
Min7.422410.279.0853.511.7171.715.921.85.040.0020.0100.001182.2
Max8.4015087.05307.2224.179.1429.3343.5318.2216.00.0410.3240.013814.2
Mean7.80720.42.0047.3143.836.0311.789.9151.475.80.0050.0670.003507.3
SD0.22280.71.2658.041.614.668.872.674.945.00.0090.0840.004157.7
Exceed standard rate (%)011.53.83.83.834.603.8073.1
(III) Center of
the Hutuo River
pluvial fan
(n = 12)
Min7.38152.21.2617.720.65.95141.01.065.271.760.0020.0100.00176.1
Max7.97626.52.56139110.141.5334.2135.3168.621.80.0041.3130.040445.9
Mean7.63412.82.1336.569.825.2249.342.681.49.100.0020.2140.005278.0
SD0.18124.90.5037.122.310.856.1335.045.37.350.0010.3990.01196.0
Exceed standard rate (%)00000008.300
Standard6.5–8.5100020025025088.63.290.30.1450
Note: Mean: average value; SD: standard deviation; Min: minimum value; Max: maximum value; Standard is grade III standard for groundwater quality in China (GB/T14848–2017).
Table 2. Spearman’s correlation matrix of water quality of groundwater in the Hutuo River alluvial-pluvial fan.
Table 2. Spearman’s correlation matrix of water quality of groundwater in the Hutuo River alluvial-pluvial fan.
ParameterspHTDSK+Na+Ca2+Mg2+ClSO42−HCO3NO3FeMnTH
pH1.000
TDS−0.0791.000
K+−0.0660.2151.000
Na+−0.0880.687 **0.320 *1.000
Ca2+−0.0670.917 **0.0600.439 **1.000
Mg2+−0.1770.752 **0.1960.540 **0.648 **1.000
Cl−0.1080.863 **0.1700.706 **0.741 **0.598 **1.000
SO42−0.0150.826 **0.2310.709 **0.746 **0.634 **0.681 **1.000
HCO3−0.1650.823 **0.2570.656 **0.742 **0.808 **0.631 **0.670 **1.000
NO30.0430.868 **0.2220.370 *0.909 **0.585 **0.683 **0.649 **0.645 **1.000
Fe0.1480.0800.1110.0510.0760.0670.063−0.0060.1270.0301.000
Mn0.0910.0320.205−0.0180.048−0.012−0.018−0.0070.1570.0290.737 **1.000
TH−0.0980.956 **0.1300.519 **0.964 **0.803 **0.770 **0.756 **0.820 **0.889 **0.0830.0401.000
Note: * p < 0.05; ** p < 0.01.
Table 3. Variance analysis of contrast coefficient of groundwater hydrochemical compositions in the Hutuo River alluvial-pluvial fan.
Table 3. Variance analysis of contrast coefficient of groundwater hydrochemical compositions in the Hutuo River alluvial-pluvial fan.
ParametersHRVPTHRAPFCHRAPF
pH0.0000.0010.001
TH0.1140.0970.119
TDS0.2300.1520.092
K+1.7440.3970.056
Na+1.2161.5031.033
Ca2+0.1040.0840.102
Mg2+0.2850.1650.183
Cl0.6680.6530.675
SO42−0.1990.2450.310
HCO30.0490.0490.051
NO30.6290.3530.652
Fe1.4741.5623.480
Mn1.0111.5410.292
Note: HRVP: Hutuo Rvier Valley Plain; THRAPF: Top of Hutuo River alluvial-pluvial fan; CHRAPF: Center of Hutuo River alluvial-pluvial fan.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ren, C.; Zhang, Q. Groundwater Chemical Characteristics and Controlling Factors in a Region of Northern China with Intensive Human Activity. Int. J. Environ. Res. Public Health 2020, 17, 9126. https://doi.org/10.3390/ijerph17239126

AMA Style

Ren C, Zhang Q. Groundwater Chemical Characteristics and Controlling Factors in a Region of Northern China with Intensive Human Activity. International Journal of Environmental Research and Public Health. 2020; 17(23):9126. https://doi.org/10.3390/ijerph17239126

Chicago/Turabian Style

Ren, Chaobin, and Qianqian Zhang. 2020. "Groundwater Chemical Characteristics and Controlling Factors in a Region of Northern China with Intensive Human Activity" International Journal of Environmental Research and Public Health 17, no. 23: 9126. https://doi.org/10.3390/ijerph17239126

APA Style

Ren, C., & Zhang, Q. (2020). Groundwater Chemical Characteristics and Controlling Factors in a Region of Northern China with Intensive Human Activity. International Journal of Environmental Research and Public Health, 17(23), 9126. https://doi.org/10.3390/ijerph17239126

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