Next Article in Journal
Photodegradation Behavior of Agricultural Antibiotic Oxytetracycline in Water
Next Article in Special Issue
Factors Controlling Natural Background Levels of Ammonium and Iodide in Shallow Groundwater of Coastal Aquifers, South China
Previous Article in Journal
Teaching Floods in the Context of Climate Change with the Use of Official Cartographic Viewers (Spain)
Previous Article in Special Issue
Key Factors Dominating the Groundwater Chemical Composition in a Grain Production Base: A Case Study of Muling–Xingkai Plain, Northeast China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anthropogenic Influences on the Hydrochemical Characteristics of the Groundwater in Xiamen City, China and Their Evolution

1
Fujian Provincial Key Laboratory of Water Cycling and Eco-Geological Processes, Xiamen 361000, China
2
Institute of Hydrogeology and Environmental Geology, CAGS, Shijiazhuang 050061, China
*
Authors to whom correspondence should be addressed.
Water 2022, 14(21), 3377; https://doi.org/10.3390/w14213377
Submission received: 14 September 2022 / Revised: 21 October 2022 / Accepted: 22 October 2022 / Published: 25 October 2022
(This article belongs to the Special Issue Groundwater Chemistry and Quality in Coastal Aquifers)

Abstract

:
This study analyzed the anthropogenic influences on the hydrochemical composition characteristics of the groundwater in Xiamen City, Fujian province, China, and their evolution. Based on the hydrochemical data of the groundwater of 1993 and 2019–2021, this study identified the indices of the anthropogenic influences using mathematical and statistical analysis methods, such as contrast coefficient, standard deviation, and Mahalanobis distance. The analytical results are summarized as follows: (1) the number of the indices affecting the groundwater quality in Xiamen increased from nine in 1993 to 15 in 2019, and the six increased indicators included NO3, Pb, NH4+, Al3+, NO2 and Cu (the contribution rates to poor-quality were 26.0%, 16.3%, 10.6%, 4.1%, 0.8% and 0.8%, respectively) which were related to the input of human activities. During this period, the number of hydrochemical types increased from 19 in 1993 to 28 in 2019, with a decrease in the water of the HCO3 type and an increase in the water of Cl and SO4 types; (2) In 2019, NO3 had higher content than SO42− in the groundwater and became a major anion, forming the water of NO3 type; (2) as indicated by the analytical results obtained using the Mahalanobis distance method, areas with strong anthropogenic influences include densely populated areas and areas with intensively distributed industrial enterprises, while anthropogenic influences are very weak in the northern forest land area.

1. Introduction

Safe drinking water, which is closely related to human health, food security, and the orderly development of the social economy, is essential for the sustainable development of a city [1,2]. With the acceleration of urbanization, the area of urban built-up areas and the number and density of the urban population have increased dramatically, leading to a sharp increase in the demand for water resources [3,4,5]. Groundwater is the main freshwater source for industrial and agricultural production, domestic drinking, and ecological water consumption and is the only source of water supply in arid and semi-arid regions [6,7]. The groundwater quality directly affects the quality of the human living environment and human health, and the groundwater environment is controlled by various factors such as meteorology, hydrogeology, and human activities [8,9,10,11,12]. Urbanization in China presents a sustainable and rapid growing trend, with the urbanization rate increasing from 18.57% in 1978 to 58.52% in 2017 [13]. Moreover, the ability of human beings to transform nature is continuously enhanced, and large-scale human activities have strong influences on climate [14,15,16], ecological environment [17,18], soil environment [19], and groundwater environment [20,21,22]. Many cities around the world, especially those in arid and semi-arid regions, are facing dual water stress exerted by groundwater shortage and the deterioration of water quality [23,24].
In 1984, the Commission on Groundwater Quality under the International Association of Hydrogeologists held a symposium entitled Impact of Agriculture on Groundwater in Ireland [25]. The international hydrological community has paid increasing attention to and studied the influence of human activities on the groundwater environment [24,26,27,28,29,30,31]. When the groundwater environment changes under anthropogenic influences, the hydrochemical characteristics of groundwater will also change accordingly, i.e., hydrochemical anomalies occur [32]. Since various geochemical indices frequently coexist in complex geological systems through mutual influences and constraints, the synergistic effects of various geochemical indices should be fully considered during the identification of geochemical anomalies, aiming to reduce the uncertainty of their identification performance. The multi-element geochemical anomalies can be presently identified using multiple methods, such as the hydrochemical diagrams (Durov diagram, ion ratio plot, et al.) method [33], coefficient of variation [34], affine equivariant [35], and multivariate statistical method [28]. The Mahalanobis distance is a preferred composite index used to identify multi-element geochemical anomalies. This index has been widely used in geochemical fields including geochemical prospecting and metallogenic prediction [36,37], the evaluation of geochemical anomalies of hydrocarbons [38], and the identification of geochemical anomalies of sediments [39]. However, it is rarely used in the identification of hydrochemical anomalies of groundwater.
Although Xiamen has abundant rainfall, it still suffers from a severe shortage of freshwater resources, which is a key problem restricting the sustainable economic development of the city. Groundwater accounts for only about 6% of Xiamen’s urban water supply [40], indicating great potential for groundwater exploitation. Moreover, groundwater is an important supplementary and backup emergency water source of Xiamen and its quality directly affects the emergency water supply ability and water supply safety of the city. This study identified the indices of anthropogenic influences using standard deviations, calculated the contrast coefficients of these indices and the Mahalanobis distances between these contrast coefficients, and determined the intensity of anthropogenic influences on the evolution of hydrochemical characteristics of the groundwater. This study qualitatively analyzed the anthropogenic influence on the changes in the groundwater chemistry of the study area since only 172 samples were collected from the large study area. The results of this study will provide some basis for the protection of groundwater and guide the management of groundwater resources to a certain extent.

2. Overview of the Study Area

2.1. Physical Geography and Social Economy

Xiamen is located in southeastern China, with geographical coordinates of 117°53′–118°26′ E and 24°23′–24°54′ N. This city lies in the central portion of the west coast of the Taiwan Strait and the center of the Golden Triangle of Southern Fujian Province. It covers a total area of 1700.61 km2, including approximately 390 km2 of seas, with a total coastline length of approximately 234 km. Its administrative region consists of land, seas, and islands. It generally dips to SE, with high and low hills, terraces, marine plains, and tidal flats distributed from NW to SE. Moreover, the Xiamen Island, the Gulangyu Islet, and seas lie in its south. Xiamen has a mild and humid subtropical marine monsoon climate, with an average multiyear temperature of 20.9 °C and average multiyear precipitation of 1100–2000 mm, which increases from SE to NW. Xiamen has developed water systems, a dense dendritic river network, and a runoff direction from NW to SE. The rivers in the city feature small runoff, short flows, narrow channels, and shallow riverbeds. Their water quantity varies greatly with seasons and has low sediment content. Xiamen enjoys abundant marine biological resources but lacks metal mineral resources and freshwater resources.

2.2. Geological and Hydrogeological Conditions

Xiamen is located in the coastal metamorphic zone of southeastern Fujian and is at the intersection of the large EW-trending Nanjing-Xiamen fault and the NES-trending Changdong-Nan’ao deep-seated fault. The exposed strata in the study area consist of the Upper Triassic Wenbinshan Formation, the Lower Jurassic Lishan and Fankeng formations, the Upper Jurassic Nanyuan Formation, and Quaternary strata including the Middle Pleistocene diluvium and alluvial-proluvial strata, the Upper Pleistocene alluvial-proluvial and marine strata, the Holocene alluvial-proluvial and marine strata, and the Quaternary eluvium.
The groundwater in Xiamen mainly occurs in aquifers of unconsolidated rock pores, weathered eluvial pores and fractures, and bedrock fissures (Figure 1). The pore water of unconsolidated rocks is hosted in alluvial-proluvial and marine deposits. The pore water in alluvial-proluvial deposits is distributed in the terraces on both sides of valleys and in piedmont proluvial fans. Its aquifers consist of sands, gravels, and pebbles and have a thickness of 3–9 m and a specific yield of <200 m3/d·m [41]. The depth of water level is 0.5–5.0 m. The pore water in marine deposits is distributed in the coastal and estuarine areas. Its aquifers consist of silty sands and medium-coarse-grained sands and have a small thickness. The pore-fissure water in weathered eluvial deposits is distributed in eluvial platforms and at piedmont slope toes. Its aquifers consist of detritus and breccias and generally have a thickness of 10–20 m and a specific yield of <10 m3/d·m. The bedrock fissure water is distributed in the fissures and faults of bedrock in mountainous areas, with the water yield controlled by structural faults. Its water content is extremely uneven, with a specific yield of <10 m3/d·m. The groundwater in Xiamen is primarily recharged by vertical infiltration of atmospheric precipitation and agricultural irrigation as well as the lateral flow of groundwater in bedrock mountainous areas during the rainy season. The annual average natural replenishment resource is 1.95 × 103 m3/a. The regional groundwater flows from the NW and NE to the southeast coast overall.

3. Materials and Methods

3.1. Data Sources

Two sets of hydrochemical data were used in this study. One set was the historical hydrochemical data of 1993 obtained from the study results of the study area. It included the data on 84 groundwater samples (19 samples of bedrock fissure water, 27 pore-fissure water samples, and 38 pore water samples), detailing the locations of sampling points, types of water-bearing media, hydrochemical types, pH, COD, the total dissolved solids (TDS) content, and the contents of HCO3, SO42−, Cl, K++Na+, Mg2+, NH4+, NO2, NO3, F, Fe, Al3+, Cu, Pb, Zn, As, and Cr6+. Unfortunately, the Ca2+ content was not recorded.
The other set of original data was 2019–2021 hydrochemical data of groundwater samples surveyed by the project team of this study. These samples were mostly collected from large-diameter wells or small pumping wells in villages and small towns, and 75% of the wells had a depth less than 50 m. The data on 172 groundwater samples were obtained in total, including 58 samples of bedrock fissure water, 72 pore-fissure water samples, and 42 pore water samples. Regarding sampling densities, seven and 8–18 samples were collected per 100 km2 in bedrock mountainous areas and plains, respectively. This set of data contained all the information that matches the historical data. Two 500 mL polyethylene bottles were used to store filtered (1.2-μm filter membrane) groundwater for the analysis of trace elements and major ions. The bottles were cleaned three times using the water to be sampled before sampling. The water sample in the bottle used to determine trace elements was acidized to pH < 2 using nitric acid, while the water sample in the bottle used to analyze major ions was not acidized. All samples were stored at 4 °C until laboratory procedures could be performed.

3.2. Tests and Analysis

All of the water samples were sent to the Groundwater, Mineral Water and Environmental Monitoring Center of the Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences to be tested in strict accordance with GB/T5750-2006 Standard Examination Methods for Drinking Water and GB 8538-2016 National Food Safety Standard–Methods for Examination of Drinking Natural Mineral Water. Metal ions and trace elements (K+, Na+, Mg2+, Fe, Cu, Pb, Zn, As, and Cr6+) were measured using inductively coupled plasma-mass spectrometry (ICP-MS) (Agilent 7500ce ICP-MS, Tokyo, Japan). SO42−, Cl, NH4+, and NO3 were measured using ion chromatography (IC) (Shimadzu LC-10ADvp, Japan). NO2 and HCO3 were determined using the UV spectrophotometry method, and acid-base titration method, respectively. Chemical oxygen demand (COD) and TDS content were measured using potassium dichromate titration and the gravimetric method, respectively. To ensure the accuracy and reliability of the test results, the anion-cation balance was analyzed according to the principle of ionic balance, with analysis errors less than 5%. All chemicals were purchased from Sinopharm Chemical Reagent Co., Ltd., Shanghai, China.

3.3. Multivariate Statistical Method

First, the anomaly thresholds or background values were determined according to the historical hydrochemical data. The indices of both anthropogenic influences and water-rock interactions were then identified using contrast coefficients and standard deviations. Finally, the Mahalanobis distances between the multivariate contrast coefficients of the indices of anthropogenic influences were calculated to analyze the changes in hydrochemical components and determine the intensity of anthropogenic influences on the hydrochemical characteristics of groundwater.

3.3.1. Identification of the Indicative Indices of Both Anthropogenic Influences and Water-Rock Interactions

Standard deviations were used to identify the indices of both anthropogenic influences and water-rock interactions.
A standard deviation (S) is the average of the distances between some data and their average, and a higher standard deviation means higher degrees of fluctuation and dispersion of a variable. The anthropogenic inputs in groundwater have a high degree of spatial dispersion. Therefore, the standard deviation was used to analyze the dispersion degree of groundwater components to determine the indices of the anthropogenic influences on the hydrochemical characteristics of groundwater. The standard deviation is expressed as:
S = 1 N i = 1 N ( X i X ¯ ) 2
where S is the standard deviation; N is the number of samples; Xi is the tested value; X ¯ is the arithmetic average.
The contrast coefficient (V) is a metric that measures the sharpness of anomalies and represents the fluctuation of the anomalies relative to the background values. To compare the anomaly intensity of a certain indicator element in different regions, the contrast coefficient can be used to eliminate the influences caused by the phenomenon that the indicator element has different background values in different regions. First, the ratio of each variable to the anomaly threshold or background value was calculated using Equation (2), and multiple variables with unequal averages were transformed into those with equal averages. The standard deviations of the contrast coefficients of various variables were then calculated and finally compared.
V = X i T
where V is the contrast coefficient; Xi is the tested value; T is the anomaly threshold or background value.
The anomaly threshold or background value (T) is a metric that reflects the geochemical background. The groundwater background values refer to the natural contents of physical and chemical elements in groundwater under uncontaminated conditions. However, a natural environment with no anthropogenic influence barely exists in the strict sense [42]. To analyze the anthropogenic influences on the hydrochemical characteristics of the groundwater in the study area in the past three decades, this study determined the range of background values using the quartile method based on the hydrochemical data of 1993 [43,44].

3.3.2. Intensity of Anthropogenic Influences

The Mahalanobis distance (MD) is a metric that measures multivariate distance. It is used for multivariate statistics of observed objects by characterizing the distance between the objects and their distribution center, during which covariance is considered [45]. This metric is frequently used to calculate the similarity of various indices, identify and quantify anomalies, and sometimes select representative data from a large quantity of measurement data [46]. This study analyzed the changes in hydrochemical components by calculating the Mahalanobis distances between the contrast coefficients of the indices of anthropogenic influences in order to further determine the intensity of anthropogenic influence on the hydrochemical characteristics of groundwater. The calculation formula for Mahalanobis distance is as follows:
D 2 = ( x i x ¯ ) S 1 ( x i x ¯ )
where D 2 is the Mahalanobis distance; x i is an index value of a sample; x ¯ is the arithmetic average of various indices; S is the standard deviation.
This study used the Kriging interpolation method of MapGis K9 to analyze the intensity grades of anthropogenic influences on the hydrochemical characteristics of the groundwater and then plotted the map of regional anthropogenic influence intensity.

4. Results and Discussion

4.1. Characteristics of Groundwater Quality Indicators

The descriptive statistics for the indicators analyzed in groundwater samples in Xiamen City are summarized in Table 1. Among anions in the groundwater of Xiamen, NO3 showed the highest median concentration of 42.23 mg/L, followed by Cl, SO42− and F. High concentrations of NH4+ occurred in groundwater in Xiamen, and up to 107.0 mg/L. Na+ ranged from 2.68 mg/L to 4043.8 mg/L, with the median of 32.55 mg/L. Among trace metal(loid)s in the groundwater of Xiamen, Fe showed the highest median concentration of 0.03 mg/L, followed by Zn, Pb, Cu, As, Cr6+ and Al3+. According to the standard for groundwater quality of China (GAQSIQPRC, 2017) [47], the percent of groundwater samples with the concentrations of indicators that exceeded the allowable values (PEV) for drinking purposes are also in Table 1. In Xiamen City, NO3 showed the highest PEV of 22.09%, followed by Fe, Pb, NH4+, F, Cl, Al3+, SO42−, NO2, Na+, Cu, and Zn, while the PEVs of the other two indicators were zero.

4.2. Hydrochemical Changes in Groundwater

4.2.1. Changes in Groundwater Chemical Composition

The evaluation indices used in this study were the water quality indices specified in GB/T14848-2017 Standard for groundwater quality. They were contained in the hydrochemical data of 1993 and consisted of 17 indices, namely pH, COD, TDS content, SO42−, Cl, Na+, NH4+, NO2, NO3, F, Fe, Al3+, Cu, Pb, Zn, As, and Cr6+. We next analyzed the changes in the chemical composition of the groundwater in Xiamen in the past three decades by comparing the number and contribution rates of the influencing factors of groundwater quality between 1993 and 2019–2021.
The following evaluation results of the groundwater quality of Xiamen in the past three decades were obtained based on the number of groundwater samples: (1) the proportion of Class I to III potable groundwater decreased from 40% in 1993 to 28.5% during 2019–2021; (2) the total proportion of poor-quality and undrinkable groundwater (Class IV and V) increased from 60% in 1993 to 71.5% in 2019–2021; (3) the proportion of Class IV water that can become portable water after proper treatment increased from 28.2% in 1993 to 51.2% in 2019–2021, and the proportion of non-portable Class V water decreased from 31.8% in 1993 to 20.3% in 2019–2021 (Table 2).
The number of the indices affecting the groundwater quality in Xiamen increased from nine in 1993 to 15 in 2019–2021 (Figure 2). The additional six indices during 2019–2021 consisted of NO3, Pb, NH4+, Al3+, NO2, and Cu (the contribution rates to poor-quality groundwater were 26.0%, 16.3%, 10.6%, 4.1%, 0.8% and 0.8%, respectively), among which NO3 and Pb had contribution rates second only to pH.

4.2.2. Changes in the Hydrochemical Types of Groundwater

Human activities have not only changed the groundwater quality but also affected and changed the natural hydrochemical field of groundwater due to a large number of anthropogenic inputs [48].
The groundwater of Xiamen had simple hydrochemical types in 1993. Specifically, 19 types were identified in 67 samples, and seven types were concentrated in 70% of the groundwater samples, namely Na-–HCO3, Ca·Na–HCO3, Na–HCO3·Cl, Na·Ca–HCO3, Na–Cl, Na·Ca–HCO3·Cl, and Na–Cl·HCO3. Regarding anions, the groundwater was dominated by the HCO3 type, which accounted for 46.3%; followed by the HCO3·Cl type, which accounted for 25.4%. Moreover, the groundwater of SO4·Cl and Cl·SO4 types accounted for a minimum proportion of 1.5% each. In terms of cations, the groundwater in Xiamen in 1993 was dominated by Na and Na·Ca types, which accounted for 41.8% and 25.4%, respectively. They were followed by the Ca·Na type, which accounted for 14.9%. Additionally, the water of Ca, Na·Mg, Ca·Mg, Ca·Na·Mg, and Na·Mg types accounted for 1.5%–6.05% each.
The hydrochemical types of groundwater in Xiamen increased during 2019–2021 compared to 1993. Specifically, 28 hydrochemical types were identified in 172 samples, and seven hydrochemical types were concentrated in 70% of the groundwater samples, namely Ca·Na–HCO3, Na·Ca–Cl, Na·Ca–HCO3, Ca·Na–HCO3·Cl, Na·Ca–HCO3·Cl, Na·Ca–Cl·HCO3, Ca–HCO3, and Na·Ca–Cl·SO4. Regarding anions, the groundwater was still dominated by the HCO3 type, whose proportion, however, decreased to 31.4%. By contrast, the Cl and SO42− contents in the groundwater increased, with the water of Cl type and the hydrochemical types dominated by Cl (Cl·HCO3, Cl·HCO3·SO4 type, and Cl·SO4 types) accounted for 16.9% and 17.4%, respectively. Moreover, the proportion of the water of SO4 type or types dominated by SO4 increased to 5.3%. Regarding cations, the Ca2+ content increased, the Mg2+ content decreased, and the groundwater was dominated by the Na·Ca type, whose proportion increased to 51.2%. This type was followed by the Ca·Na type, which accounted for 36.0%. Furthermore, the water of Ca and Na types had a low proportion of 6.4% each.
Groundwater in the pore-fissure aquifer showed higher Cl and SO42− concentrations than that in the bedrock fissure aquifers and pore aquifers (Figure 3). This may have occurred due to the greater number of industrial zones and population located in the area of the pore-fissure aquifer than in the areas of the other two types of aquifers, indicating more serious anthropogenic influences.
The increase in the absolute content of nitrates in groundwater affects the groundwater quality, while the increase in their relative content changes the hydrochemical characteristics of groundwater. As indicated by the hydrochemical test results of the groundwater in Xiamen in 2019, NO3 was a major anion in Xiamen’s groundwater, its relative content exceeded that of SO42− (Table 3), and thus more water of NO3 type exists in the groundwater. It is known that groundwater NO3 in Xiamen mainly originated from anthropogenic sources, including infiltration of domestic sewage and landfill leachate [40].

4.3. Anthropogenic Influences on the Hydrochemical Characteristics of Groundwater

4.3.1. Indicative Indices of the Anthropogenic Influences

Determination of the Background Values of Groundwater in Different Water-Bearing Media

This study statistically analyzed the data of 1993, obtaining the background values of the groundwater environment in aquifers of bedrock fissures, weathered pores and fissures, and pores in alluvial-proluvial deposits, respectively (Table 4). Furthermore, the background values of indices Cu, Pb, As, and Cr6+ were determined to be 0, since there were a few samples of these indices and they were not detected in the water of various water-bearing media.

Identification of the Indicative Indices of Anthropogenic Influences

To conduct a multivariate contrast, various variables were normalized first. In other words, the upper limits of the background values of various indices in Table 3 were used as the anomaly threshold or background value (T) to calculate the contrast coefficients (V) and the standard deviations of various indices of the groundwater of Xiamen in 2019.
Anthropogenic inputs in groundwater generally have high degrees of spatial dispersion and fluctuation [48]. Therefore, standard deviations can be used to distinguish between groundwater components originating from anthropogenic inputs and those from water-rock interactions. The contrast coefficients of pH and HCO3- had a standard deviation of 0.09 and 0.53, respectively, indicating very weak anthropogenic influences. By contrast, the contrast coefficients of Cl, NO3, and NH4+ had a standard deviation of 18.86, 35.83, and 42.76, respectively, indicating very strong anthropogenic influences. The standard deviations of the contrast coefficients of the 15 indices decreased in the order of pH, HCO3, COD, total hardness, Zn, Fe, NO2, SO42−, F, TDS, Mg2+, K++Na+, Cl, NO3, and NH4+. These 15 indices in the groundwater were divided into four categories according to their sources determined by the standard deviations of their contrast coefficients, as shown in Table 5.

4.3.2. Anthropogenic Influence Intensity

The Mahalanobis distance is a metric used for multivariate statistics of observed objects. This study combined five indices: Mg2+, K++Na+, Cl, NO3, and NH4+ subject to strong or very strong anthropogenic influences combined into a whole, then calculated the Mahalanobis distances of the contrast coefficients of various groundwater samples collected in 2019, and finally analyzed the anthropogenic influence intensity of various samples. Figure 4 shows the distribution of the Mahalanobis distances of 172 samples. The anthropogenic influence intensity of these samples was divided into four grades according to the Mahalanobis distances 0.25, 0.6, and 6 at inflexion points, which represented the boundaries of anthropogenic influence intensity grades 2, 3, and 4, respectively.
The intensity of anthropogenic influences on the hydrochemical composition of the groundwater in Xiamen was analyzed based on the Mahalanobis distance method. Figure 5 shows the types of land use and the intensity and distribution of anthropogenic influences on the hydrochemical composition of the groundwater in Xiamen. According to this figure, the groundwater subject to strong or relatively strong anthropogenic influences is mostly distributed on artificial surfaces and cultivated land, while the anthropogenic influences on the groundwater sites in forest land areas were mostly rated as weak.
Interpolation was conducted for the calculated Mahalanobis distances using the Kriging interpolation method of the MapGis software. The anthropogenic influence intensity was then divided into four grades according to the interpolated Mahalanobis distances, obtaining the map showing the zones of the anthropogenic influence intensity (Figure 6). According to this map, the zones subject to very strong anthropogenic influences include densely populated areas and areas with intensively distributed industrial enterprises, such as Xinxu Town in the Xiang’an District, Xike Street in the Tong’an District, Houxi Town and Xinglin Street in the Jimei District, Haicang Street in the Haicang District, and Dianqian Street in the Huli District. With low population density and high vegetation coverage, the northern bedrock mountainous area is a zone with very weak anthropogenic influences on groundwater.

5. Conclusions

1.
The number of indices affecting the groundwater quality in Xiamen increased from nine in 1993 to 15 in 2019, and the six increased indicators included NO3, Pb, NH4+, Al3+, NO2 and Cu. The number of hydrochemical types increased from 19 in 1993 to 28 in 2019. Moreover, the water of the HCO3 type has decreased, while the water of the Cl type (including types dominated by Cl) and the SO4 type (including types dominated by SO4) have increased, indicating that human activities are the powerful force driving the change in groundwater chemistry in Xiamen City.
2.
The increase in the relative contents of nitrates in the groundwater has changed the hydrochemical characteristics of the groundwater. Both the maximum and average relative content of NO3 exceed those of SO42−. Therefore, NO3 has become a major anion affecting the hydrochemical nomenclature of the groundwater in Xiamen, forming the water of NO3 type (according to Shukarev classification and Designation rules).
3.
As indicated by the analysis results of Mahalanobis distances, zones with very strong anthropogenic influences on groundwater include Houxi Town in Jimei District and Haicang Street in Haicang District, while zones with very weak anthropogenic influences on groundwater include the northern bedrock mountainous area. These results are consistent with the land use, industrial layout, and population distribution of Xiamen, indicating that the Mahalanobis distance method can be used to identify the hydrochemical anomalies of groundwater and determine the intensity of anthropogenic influences on groundwater and that the evaluation results obtained using this method is reasonable.

Author Contributions

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

Funding

This study was financially supported by China Geological Survey’s project (Grant No. DD20190303 and DD20221773).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Y.; Zheng, C.; Ma, R. Review: Safe and sustainable groundwater supply in China. Appl. Hydrogeol. 2018, 26, 1301–1324. [Google Scholar] [CrossRef]
  2. Turhan, A.; Ozmen, N. Influence of Chloride on Growth, Fruit Yield and Quality Parameters of Processing Pepper. J. Agric. Nat. 2021, 24, 1139–1144. [Google Scholar]
  3. Döll, P.; Müller Schmied, H.; Schuh, C.; Portmann, F.T.; Eicker, A. Global-scale assessment of groundwater depletion and related groundwater abstractions:Combining hydrological modeling with information from well observations and GRACE satellites. Water Resour. Res. 2014, 50, 5698–5720. [Google Scholar] [CrossRef]
  4. Mirdashtvan, M.; Najafinejad, A.; Malekian, A.; Sadoddin, A. Sustainable water supply and demand management in semi-arid regions: Optimizing water resources allocation based on RCPs scenarios. Water Resour. Manag. 2021, 35, 5307–5324. [Google Scholar] [CrossRef]
  5. Liu, C.L.; Zheng, J.H.; Li, Z.H.; Li, Y.S.; Hao, Q.C.; Li, J.F. Analysis on the situation and countermeasures of water resources supply and demand in the cities of small and medium-sized river basins along southeast coast of China-taking Xiamen City as an example. J. Groundw. Sci. Eng. 2021, 9, 350–358. [Google Scholar]
  6. Yang, H.F.; Meng, R.F.; Li, W.P.; Li, Z.Y.; Zhi, C.S.; Bao, X.L.; Li, C.Q.; Liu, F.T.; Wu, H.P.; Ren, Y. Groundwater resources of the Haihe River Basin and its development potential. Geol. China 2021, 48, 1032–1051. [Google Scholar]
  7. Nasiri, S.; Ansari, H.; Ziaei, A.N. Determination of water balance equation components in irrigated agricultural watersheds using SWAT and MODFLOW models: A case study of Samalqan plain in Iran. J. Groundw. Sci. Eng. 2022, 10, 44–56. [Google Scholar]
  8. Zhang, Y.L.; Ma, R.; Li, Z.H. Human health risk assessment of groundwater in Hetao Plain (Inner Mongolia Autonomous Region, China). Environ. Monit. Assess. 2014, 186, 4669–4684. [Google Scholar] [CrossRef]
  9. Tirkey, P.; Bhattacharya, T.; Chakraborty, S.; Baraik, S. Assessment of groundwater quality and associated health risks: A case study of Ranchi city, Jharkhand, India. Groundw. Sustain. Dev. 2017, 5, 85–100. [Google Scholar] [CrossRef]
  10. Adimalla, N.; Wu, J. Groundwater quality and associated health risks in a semi-arid region of south India: Implication to sustainable groundwater management. Hum. Ecol. Risk Assess. Int. J. 2019, 25, 191–216. [Google Scholar] [CrossRef]
  11. Adimalla, N. Controlling factors and mechanism of groundwater quality variation in semiarid region of South India: An ap-proach of water quality index (WQI) and health risk assessment (HRA). Env. Geochem. Health 2020, 42, 1725–1752. [Google Scholar] [CrossRef] [PubMed]
  12. Egbi, C.D.; Anornu, G.K.; 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]
  13. Li, X.Y.; Zhang, Y.L.; Li, Z.H.; Wang, R. Response of the groundwater environment to rapid urbanization in Hohhot, the provincial capital of western China. J. Hydrol. 2021, 603, 127033. [Google Scholar] [CrossRef]
  14. Mote, T.; Lacke, M.; Shepherd, J.M. Radar signatures of the urban effect on precipitation distribution: A case study for Atlanta, Georgia. Geophys. Res. Lett. 2007, 34, L20710. [Google Scholar] [CrossRef]
  15. Salerno, F.; Gaetano, V.; Gianni, T. Urbanization and climate change impacts on surface water quality: Enhancing the resilience by reducing impervious surfaces. Water Res. 2018, 144, 491–502. [Google Scholar] [CrossRef] [PubMed]
  16. Zhou, Q.Q.; Leng, G.Y.; Su, J.H.; Ren, Y. Comparison of urbanization and climate change impacts on urban flood volumes: Importance of urban planning and drainage adaptation. Sci. Total Environ. 2019, 658, 24–33. [Google Scholar] [CrossRef]
  17. Ariken, M.; Zhang, F.; Liu, K.; Fang, C.L.; Kung, H.T. Coupling coordination analysis of urbanization and eco-environment in Yanqi Basin based on multi-source remote sensing data. Ecol. Indic. 2020, 114, 106331. [Google Scholar] [CrossRef]
  18. Yu, B.B. Ecological effects of new-type urbanization in China. Renew. Sustain. Energy Rev. 2021, 135, 110239. [Google Scholar] [CrossRef]
  19. Shaharoona, B.; Al-Ismaily, S.; Al-Mayahi, A.; Al-Harrasi, N.; Al-Kindi, R.; Al-Sulaimi, A.; Al-Busaidi, H.; Al-Abri, M. The role of urbanization in soil and groundwater contamination by heavy metals and pathogenic bacteria: A case study from Oman. Heliyon 2019, 5, e01771. [Google Scholar] [CrossRef] [Green Version]
  20. Li, Z.H.; Zhang, Y.L.; Hu, B.; Wang, L.J.; Zhu, Y.C.; Li, J.F. Driving action of human activities on NO3—N pollution in confined groundwater of Togtoh County, Inner Mongolia. Acta Geosci. Sinica 2018, 39, 358–364. [Google Scholar]
  21. Kaur, L.; Rishi, M.S.; Siddiqui, A.U. Deterministic and probabilistic health risk assessment techniques to evaluate non-carcinogenic human health risk (NHHR) due to fluoride and nitrate in groundwater of Panipat, Haryana, India. Environ. Pollut. 2020, 259, 113711. [Google Scholar] [CrossRef] [PubMed]
  22. Zhu, L.; Yang, M.N.; Zhang, Y.X.; Liu, J.T. Research on distribution and influence factor of nitrate of shallow groundwater in Dongsheng and surrounding area. J. Water Resour. Water Eng. 2016, 27, 37–41,45. [Google Scholar]
  23. Kangabam, R.D.; Bhoominathan, S.D.; Kanagaraj, S.; Govindaraju, M. Development of a water quality index (WQI) for the Loktak Lake in India. Appl. Water Sci. 2017, 7, 2907–2918. [Google Scholar] [CrossRef] [Green Version]
  24. Cheng, X.; Chen, L.D.; Sun, R.H.; Jing, Y.C. Identification of regional water resource stress based on water quantity and quality: A case study in a rapid urbanization region of China. J. Clean. Prod. 2019, 209, 216–223. [Google Scholar] [CrossRef]
  25. Aldwell, C.R. Agriculture and groundwater must co-exist. Environ. Geol. Water Sci. 1986, 9, 1–2. [Google Scholar] [CrossRef]
  26. Wang, Y.X.; Ma, T.; Guo, Q.H.; Ma, R. Groundwater and environmental change. Earth Sci. Front. 2005, 12, 14–21. [Google Scholar]
  27. Wang, J.Z.; Zhang, G.H.; Nie, Z.L.; Yan, M.J.; Wang, Y. Disturbing degree of mankind activities to groundwater in the Hutuo River Valley area. Bull. Soil Water Conserv. 2010, 30, 65–69. [Google Scholar]
  28. Liu, L.; Zhou, X.; Ye, Y.H. Screening of characteristic indexes for shallow groundwater influenced by human activities using multivariate statistics. Resour. Surv. Environ. 2014, 35, 305–310. [Google Scholar]
  29. Ouyang, Q.C.; Lin, B.Q.; Chen, X.W. Effects of landuse changes on hydrological processes: An analysis based on SWAT model. J. Subtrop. Resour. Environ. 2016, 11, 65–70. [Google Scholar]
  30. Malki, M.; Bouchaou, L.; Hirich, A.; Brahim, Y.; Choukr-Allah, R. Impact of agricultural practices on groundwater quality in intensive irrigated area of Chtouka-Massa, Morocco. Sci. Total Environ. 2017, 574, 760–770. [Google Scholar] [CrossRef]
  31. Flem, B.; Reimann, C.; Fabian, K.; Birke, M.; Filzmoser, P.; Banks, D. Graphical statistics to explore the natural and anthropogenic processes influencing the inorganic quality of drinking water, ground water and surface water. Appl. Geochem. 2018, 88, 133–148. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, X.W.; He, J.T.; Peng, C.; Zhang, C.Y.; Ni, Z.H. Comparison of identification methods of main component hydrochemical anomalies in groundwater: A case study of Liujiang Basin. Environ. Sci. 2017, 38, 3225–3234. [Google Scholar]
  33. Peng, C.; He, J.T.; Liao, L.; Zhang, Z.G. Research on the influence degree of human activities on groundwater quality by the method of geochemistry: A case study from Liujiang Basin. Earth Sci. Front. 2017, 24, 321–331. [Google Scholar]
  34. Zhao, W.; Lin, J.; Wang, S.F.; Liu, J.L.; Chen, Z.R.; Kou, W.J. Influence of human activities on groundwater environment based on coefficient variation method. Environ. Sci. 2013, 34, 1277–1283. [Google Scholar]
  35. Rocke, D.M.; Woodruff, D.L. Identification of outliers in multivariate data. J. Am. Stat. Assoc. 1996, 91, 1047–1061. [Google Scholar] [CrossRef]
  36. Ji, H.J. Multivariate method for distinguishing geochemical background from geochemical anomalies. J. Chang. Univ. Earth Sci. 1988, 18, 311–320. [Google Scholar]
  37. Tian, L.M. Geochemical Data Processing and Targets Optimization in Eastern Kunlun Orogenic Belt, Qinghai Province; Wuhan China University of Geosciences: Wuhan, China, 2017. [Google Scholar]
  38. Tang, Y.P.; Liu, Y.L. The Application of Markov distance method to the anomaly Appraisal in oil-gas geochemical prospecting. Geophys. Geochem. Explor. 1998, 22, 231–233. [Google Scholar]
  39. Chen, X.N.; Wang, L.J.; Wang, H.P. Application of Mahalanobis distance anomaly in geochemical survey: Taking the 1:25000 geochemical survey of somewhere in Dulan County, Qinghai province as an example. J. Chongqing Univ. Technol. Nat. Sci. 2018, 11, 115–120. [Google Scholar]
  40. Li, Z.H.; Li, Y.S.; Hao, Q.C.; Li, J.F.; Zhu, Y.C. Characteristics, genesis and treatment measures of NO3-type groundwater in Xiamen, Fujian Province. Geol. China 2021, 48, 1441–1452. [Google Scholar]
  41. Li, Y.S.; Liu, C.L.; Hao, Q.C. Report on Comprehensive Geological Survey of Xiamen-Zhangzhou-Quanzhou City; Chinese Academy of Geological Sciences, The Institute of Hydrogeology and Environmental Geology: Shijiazhuang, China, 2021. (In Chinese) [Google Scholar]
  42. Qiu, H.X.; Huang, Q.Z. The concept of groundwater environment background and its determination. J. Ocean. Univ. Qingdao 1994, 16–19. [Google Scholar]
  43. Helsel, D.R. Less than obvious: Statistical treatment of data below the detection limit. Environ. Sci. Technol. 1990, 24, 1767–1774. [Google Scholar] [CrossRef]
  44. Lee, L.; Helsel, D. Baseline models of trace elements in major aquifers of the United States. Appl. Geochem. 2005, 20, 1560–1570. [Google Scholar] [CrossRef]
  45. Song, Y.H.; Li, Z.X.; Sun, L.H.; Jia, D.C.; Bu, X.J. Contrast between Mahalanob is distance and euclidean distance in geochemical exploration processing. Jinlin Geol. 2008, 27, 117–120. [Google Scholar]
  46. Jouan-Rimbaud, D.; Massart, D.L.; Saby, C.A.; Puel, C. Characterisation of the representativity of selected sets of samples in multivariate calibration and pattern recognition. Anal. Chim. Acta 1997, 350, 149–161. [Google Scholar] [CrossRef]
  47. General Administration of Quality Supervision Inspection and Quarantine of the People’s Republic of China (GAQSIQPRC); Standard for Groundwater Quality. Standards Press of China: Beijing, China, 2017.
  48. Cheng, D.H.; Chen, H.H.; He, J.T.; Lin, J.; Yang, D.G.; Yu, H. A study of indicators of anthropogenic influence and water-rock interaction in groundwater system in the urban region of Beijing. Hydrogeol. Eng. Geol. 2007, 34, 37–42. [Google Scholar]
Figure 1. Hydrogeological conditions and sampling points in Xiamen, (a) Hydrogeological conditions and sampling points; (b) Typical hydrogeological section.
Figure 1. Hydrogeological conditions and sampling points in Xiamen, (a) Hydrogeological conditions and sampling points; (b) Typical hydrogeological section.
Water 14 03377 g001
Figure 2. Main indices affecting groundwater quality and their contribution rates.
Figure 2. Main indices affecting groundwater quality and their contribution rates.
Water 14 03377 g002
Figure 3. Piper trilinear diagrams.
Figure 3. Piper trilinear diagrams.
Water 14 03377 g003
Figure 4. Mahalanobis distances of contrast coefficients of groundwater samples in 2019.
Figure 4. Mahalanobis distances of contrast coefficients of groundwater samples in 2019.
Water 14 03377 g004
Figure 5. Map showing the distribution of the intensity of anthropogenic influences on groundwater.
Figure 5. Map showing the distribution of the intensity of anthropogenic influences on groundwater.
Water 14 03377 g005
Figure 6. Map showing the zones of the intensity of anthropogenic influences on groundwater.
Figure 6. Map showing the zones of the intensity of anthropogenic influences on groundwater.
Water 14 03377 g006
Table 1. Descriptive statistics of concentrations of groundwater indicators.
Table 1. Descriptive statistics of concentrations of groundwater indicators.
ItemMinMedMaxSDDetection
Limits
Number of Non-DetectsAllowable ValuesPEV
(%)
SO42− (mg/L)0.3724.721171.0107.860.202501.74
Cl (mg/L)1.2935.197073.0548.160.102502.91
NO3 (mg/L)0.3942.232223.14177.530.2088.622.09
NO2 (mg/L)NDND5.200.550.0021103.291.16
F (mg/L)ND0.1911.901.040.1491.07.56
Na+ (mg/L)2.6832.554034.0305.490.01302001.16
NH4+ (mg/L)NDND107.011.450.041170.649.30
Al3+ (mg/L)NDND6.650.540.021460.22.33
Fe (mg/L)ND0.036.760.720.02490.313.95
Pb (mg/L)ND0.0010.8060.070.001660.0113.95
Cu (mg/L)NDND5.510.420.0101641.00.58
Zn (mg/L)ND0.017.740.600.002421.00.58
As (mg/L)NDND0.0090.0010.0011140.010
Cr6+ (mg/L)NDND0.0170.0020.0041610.050
SD = standard deviation; ND = below the limit of detection; Allowable Values = data from (GAQSIQPRC, 2017); PEV = percent of groundwater samples with the concentration of one index exceeded the allowable value.
Table 2. Comparison of the evaluation results of groundwater quality between 1993 and 2019–2021.
Table 2. Comparison of the evaluation results of groundwater quality between 1993 and 2019–2021.
YearNumber of SamplesProportion/%Water Inferior to Class III
Class IClass IIClass IIIClass IVClass VNumber of SamplesProportion/%
19938512.911.815.328.231.85160.0
20191721.24.722.751.220.312371.5
Table 3. Characteristic parameters of the relative contents of major anions in the groundwater of Xiamen.
Table 3. Characteristic parameters of the relative contents of major anions in the groundwater of Xiamen.
ItemYearSample NumberRangeMedAverageSD
HCO31993450.9–85.859.952.225.6
20191720.4–89.835.036.924.5
SO42−1993450.8–49.49.811.77.9
20191720.4–64.312.113.910.8
Cl1993456.5–97.521.931.823.2
20191722.8–88.225.727.315.9
NO31993450–22.52.54.35.1
20191720.1–84.720.221.816.8
Table 4. Background values of groundwater environment in different water-bearing media in Xiamen.
Table 4. Background values of groundwater environment in different water-bearing media in Xiamen.
Background Value
Scope
Bedrock
Fissure Water
Weathered
Pore Fissure Water
Pore Water in
Alluvial Proluvial Deposits
Item
pH6.57–7.046.42–6.936.05–7.09
COD0.41–1.320.90–1.220.57–1.33
Cl12.0–28.5427.81–76.0914.45–76.60
SO42−6.05–23.5418.97–40.715.64–42.81
HCO344.39–102.4938.04–218.1543.17–166.0
NO20.004–0.0480.021–0.1880–0.026
NO30.50–4.736.38–29.382.28–13.90
F0.04–0.200.04–0.200.04–0.20
K++Na+18.05–41.3230.36–114.5917.94–70.61
Mg2+1.32–5.434.18–17.632.19–6.77
NH4+0.02–0.180.03–0.280.01–0.10
Fe0.09–1.200.03–0.310.03–0.08
Al3+0–0.0260–0.0370
Zn0.048–0.4410.019–0.3900.022–0.350
TDS96.12–251.79180.4–379.55124.9–337.78
Table 5. Standard deviations of the contrast coefficients of indices and their indicative significance.
Table 5. Standard deviations of the contrast coefficients of indices and their indicative significance.
IndexStandard DeviationsSourcesAnthropogenic Influence Intensity
pH, HCO3≤1Water-rock interactionsVery weak
COD, total hardness, Zn, Fe, NO2, SO42−, F, and TDS content1–5Water-rock interactions, supplemented by anthropogenic inputsWeak
Mg2+, K++Na+5–10Anthropogenic inputs, supplemented by water-rock interactionsStrong
Cl, NO3, NH4+>10Anthropogenic inputsVery strong
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Li, J.; Zhu, Y.; Li, Y.; Hao, Q. Anthropogenic Influences on the Hydrochemical Characteristics of the Groundwater in Xiamen City, China and Their Evolution. Water 2022, 14, 3377. https://doi.org/10.3390/w14213377

AMA Style

Li Z, Li J, Zhu Y, Li Y, Hao Q. Anthropogenic Influences on the Hydrochemical Characteristics of the Groundwater in Xiamen City, China and Their Evolution. Water. 2022; 14(21):3377. https://doi.org/10.3390/w14213377

Chicago/Turabian Style

Li, Zhenghong, Jianfeng Li, Yuchen Zhu, Yasong Li, and Qichen Hao. 2022. "Anthropogenic Influences on the Hydrochemical Characteristics of the Groundwater in Xiamen City, China and Their Evolution" Water 14, no. 21: 3377. https://doi.org/10.3390/w14213377

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