Fluoride and nitrogen contamination is a global concern and has been a serious problem in agricultural areas. This study aims to identify the source of fluoride and nitrogen in the groundwater and assess groundwater quality and human health risks in the Guanzhong Plain, northwest China. The results showed that the concentrations were 0.15–4.74 mg/L for F, 0.02–89.89 mg/L for NO3-N, and BDL-2.40 mg/L for NH4+-N in groundwater. Distinct area-dependent distributions of fluoride and nitrogen were observed in the study region. Higher F and NO3-N concentrations in groundwater were detected in the northern part, and higher NH4+-N levels were observed in the southern part. Water–rock interaction and agricultural activities were the controlling factors for fluoride and nitrogen distribution in groundwater. About 80% of samples are considered to have good water quality with WQI < 100. Exposure to fluoride and nitrogen through drinking should require more attention. The total non-carcinogenic risks through oral ingestion of groundwater were 0.22–3.19 for adults and 0.51–7.44 for children, respectively. The order of pollutants in the groundwater in terms of their hazard to residents was F > NH4+-N > NO3-N > NO2-N. The findings of this study could provide more insights into groundwater management.

  • Distinct area-dependent distributions of fluoride and nitrogen were observed in the study region.

  • The primary factor affecting groundwater chemistry was the rock–water interaction.

  • This study provides an explanation for fluoride and nitrogen contamination by hydrogeochemistry.

  • The health risk of oral ingestion of groundwater is greater for children than for adults.

Over 2 billion individuals reside in regions experiencing significant water stress, a figure anticipated to rise. Furthermore, more than 1 billion people are deprived of access to clean and potable water (Tzanakakis et al. 2020). Groundwater serves as a vital resource, constituting the main supply for drinking, agricultural, industrial, and ecological uses (Jha et al. 2020; Chen et al. 2020), especially in arid and semi-arid regions, which comprise about one-third of the Earth land surface (Qian et al. 2020; Scanlon et al. 2023; Zou et al. 2024).

Understanding how hydrochemical alteration occurs in groundwater is essential for maintaining its quality (Ha et al. 2022). In hydrochemistry, anions and cations are widely regarded as crucial indicators for assessing water quality and evaluating the impact of human activities on groundwater (Gao et al. 2023). Groundwater hydrochemistry is influenced by hydrogeochemical processes, including precipitation, evaporation, weathering, adsorption and exchange reactions, and human activities. (Khan et al. 2021). The quality of groundwater has a significant impact on the health of both humans and plants (Kazi et al. 2009; Sun et al. 2024).

Groundwater with high fluoride and nitrogen levels is a global threat to drinking water safety (Andreah et al. 2023). Research indicates that more than 200 million people worldwide drink water with fluoride levels exceeding the acceptable threshold (Shaji et al. 2024). The primary source of fluoride is naturally occurring, originating from rocks and minerals that contain fluoride in their composition (Vithanage & Bhattacharya 2015; Maurya et al. 2020). Consuming groundwater with fluoride levels above the World Health Organization limit of 1.5 mg/L for drinking water can lead to widespread dental and skeletal fluorosis, affecting children's tooth development and causing leg pain, stress fractures, as well as potentially severe stomach irritation and ulcers (Chen et al. 2021; Dobrinas et al. 2022). The rise in nitrogen levels of groundwater can be linked to factors such as higher fertilizer usage, extensive irrigation, shallow groundwater tables, inadequate sanitation infrastructure, and lack of vegetative cover (Chen et al. 2016). Continuous consumption of nitrogen-rich water can lower oxygen capacity, causing methemoglobinemia in infants, and is linked to gastric cancer, multiple sclerosis, non-Hodgkin lymphoma, and thyroid enlargement (Zhang et al. 2018). In this regard, it is crucial to assess the risks associated with groundwater. The Human Health Risk Assessment (HHRA) by the United States Environmental Protection Agency (USEPA) (USEPA 2004) has been widely employed to evaluate the health risks associated with groundwater quality.

Groundwater is a vital resource for drinking, agriculture, industry, and ecosystems of the Guanzhong Plain in northwest China (Ren et al. 2021). However, high-fluoride groundwater caused the prevalence of dental fluorosis and skeletal fluorosis in the study region (Chen et al. 2021). Zhang et al. (2019a) identified that 24.61% of groundwater samples were above the nitrate limit for drinking in the plain due to agricultural activities. Due to the complex interplay of natural and anthropogenic factors, the mechanism of synchronous controls on groundwater quality and health risks is not yet fully understood.

Therefore, the main objective of this study were to (1) identify the hydrochemical properties of groundwater, (2) evaluate the groundwater quality for drinking water and irrigation purposes, and (3) identify the potential human health risks due to nitrogen and fluoride exposure in groundwater. The findings in this study would aid in gaining a comprehensive understanding of groundwater hydrochemistry and provide scientific support for decision-makers.

Study area

The study area is located in the middle of the Guanzhong Plain, which is located in Shaanxi Province, China (33°58′12″N–34°28′56″N, 107°58′59″E–108°39′45″E, Figure 1). The study area covers 2,793 km2, including Yangling, Wugong, Xingping, Zhouzhi, and Huyi Counties. The altitude ranges from a high of 1,918 m in the Qinling Mountains on the south side, gradually decreasing to 272 m toward the north. The Qinling Mountains have coarse granular particles and provide favorable conditions for capturing atmospheric precipitation and recharging streams. In the northern part, there is a 100-m-thick loess layer consisting of finer particles. The average annual temperature is 13.1 °C (Xu et al. 2019). Additionally, the annual average precipitation and evaporation are 530–700 mm and 1,000–1,200 mm, respectively (Chen et al. 2021). The Wei River flows through the study area from west to east.
Figure 1

Location map of the study area and groundwater sampling location.

Figure 1

Location map of the study area and groundwater sampling location.

Close modal

The Guanzhong Plain is responsible for 50% of the cultivated land in Shaanxi Province. It is known as the ‘Shaanxi Granary’ and the main food crops are wheat and maize (Zhang et al. 2023). Irrigation accounts for 60% of water use in the plain, with groundwater being the primary source, especially during dry seasons (Tang et al. 2015). However, the excessive use of chemical fertilizers in agriculture poses a great threat to groundwater quality. The fertilizer application rate in Shaanxi Province is 487.1 kg/ha (Meng-Yao et al. 2017), which is significantly higher than the mean national level 225 kg/ha (Ischinger 1979; Cai et al. 2018; Guo & Wang 2021). This indicates that excessive use of fertilizers, especially nitrogen, can lead to contamination of surface and groundwater through runoff or leaching.

Groundwater sampling and analysis

A total of 60 groundwater samples were collected from the study area in April 2019. To ensure the groundwater's representativeness, wells were flushed for 15 min prior to sampling. The groundwater was preserved in 500-mL polyethylene plastic bottles. Before sampling, the bottles were rinsed at least three times with the groundwater to be tested, ensuring no excess air remained inside, and the bottles were securely sealed. The pH and total dissolved solids (TDSs) were measured in the field using portable instruments (DDBJ-350F, INESA Scientific Instrument Co., Ltd., China). Na+ and K+ were assessed using flame atomic absorption spectrophotometry (SP-3520AA, Shanghai Spectrum Instrument Co., Ltd, China), while Mg2+, Ca2+ and TH concentrations were determined through EDTA titration. , Cl and F levels were evaluated using ion chromatography (Dionex 120, Dionex, Sunnyvale, USA), was measured by potentiometric titrimetric methods. Spectrophotometry was used to measure the contents of and , and (722N, Shanghai Precision & Scientific Instrument Co., LTD, China). The samples were used for quality assurance/quality control (QA/QC) analysis. All groundwater sampling, preservation, and testing procedures are adhered to the Chinese Groundwater Quality Standards (General Administration of Quality Supervision & Quarantine of China 2017), ensuring the quality and consistency of the groundwater data and applying the ion balance method to verify the validity of the groundwater quality analysis results. The charge balance checking shows ±10%, also receivable in the study.

All statistical analyses were performed using SPSS software (version 26). A Mann–Whitney U test, an efficient nonparametric method, was conducted to assess significant differences in fluoride and nitrogen concentrations between samples from the north and south banks of the Wei River. The null hypothesis assumed no variance between these two groups. This study employed ordinary kriging interpolation in ArcGIS (version 10.8.2) to examine the spatial variability of ions in groundwater.

Groundwater quality assessment

The Water Quality Index (WQI) is a numerical system used to evaluate water quality. This concept was initially introduced by Horton in 1965 (Horton 1965). Typically, the development of a WQI encompasses the following four stages: (i) choosing a group of parameters that accurately represent the water quality based on its intended use; (ii) converting each parameter into a set of rating scores; (iii) allocating weights to these parameters; and (iv) computing the index score by averaging these weighted scores (Lumb et al. 2011):
(1)
where Wi is the relative weight; ωi is the weight of each parameter; n is the number of parameters. The relative weight (Wi) values for each parameter have been computed and are listed in Table 2:
(2)
where Qi is the quality rating; Ci is the concentration of each chemical parameter in each groundwater sample; Si refers to the standard concentration for each chemical parameter in groundwater, as per the World Health Organization guidelines (WHO 2022); Cip is the optimal distilled value for a parameter, which is a pH of 7 and a concentration of zero for all other parameters:
(3)
where SIi is the sub-index of th parameter:
(4)
Table 1

Calculation of human health risk assessment parameters

ChildrenAdult
IR (L/day) 0.85 1.5 
ED (year) 30 
EF (days/year) 365 365 
BW (kg) 15 61.75 
ET (days) 2,190 10,950 
ChildrenAdult
IR (L/day) 0.85 1.5 
ED (year) 30 
EF (days/year) 365 365 
BW (kg) 15 61.75 
ET (days) 2,190 10,950 
Table 2

Results of hydrochemical analyses of groundwater samples in the study region

IndexUnitNBSBMean value for all samplesSDWHO guidelineChinese guidelinewiWi
MinMaxMeanMinMaxMean
pH – 7.7 8.8 8.3 7.1 8.5 8.1 8.1 0.3 6.5–8.5 6.5–8.5 0.06 
TDS mg/L 333 3,930 1,507 223 2,299 760 984 613.1 1,000 1,000 0.10 
K+ mg/L 0.59 50.20 6.52 0.91 290.00 10.79 9.51 37.6 10 – 0.04 
Na+ mg/L 59 770 279 434 90 147 136.2 200 – 0.04 
Ca2+ mg/L 5.4 206 52.5 17.8 159 63.9 60.5 26.7 75 – 0.06 
Mg2+ mg/L 10.6 175 76.1 5.2 167 36.0 48 40.1 30 – 0.04 
Cl mg/L 13.7 642 172 3.5 514 66 97.7 111.7 250 250 0.08 
 mg/L 20.3 1,071 263 4.8 665 101 149.9 171.3 250 250 0.04 
 mg/L 207 882 555 120 768 360 418.9 167.9 600 – 0.04 
 mg/L 0.04 0.55 0.13 BDL 2.40 0.27 0.23 0.3 – 0.5a 0.08 
 mg/L 0.05 89.89 18.88 0.02 26.78 5.61 89.9 15.3 10 20 0.08 
 mg/L BDL 0.01 0.0016 BDL 0.05 0.0034 0.0029 0.01 0.9 1a 0.10 
F mg/L 0.76 4.74 1.71 0.15 1.78 0.78 1.06 0.8 1.5 0.10 
TH mg/L 93.2 1,058 442 69.5 825 312 351 195.5 200 450 0.10 
IndexUnitNBSBMean value for all samplesSDWHO guidelineChinese guidelinewiWi
MinMaxMeanMinMaxMean
pH – 7.7 8.8 8.3 7.1 8.5 8.1 8.1 0.3 6.5–8.5 6.5–8.5 0.06 
TDS mg/L 333 3,930 1,507 223 2,299 760 984 613.1 1,000 1,000 0.10 
K+ mg/L 0.59 50.20 6.52 0.91 290.00 10.79 9.51 37.6 10 – 0.04 
Na+ mg/L 59 770 279 434 90 147 136.2 200 – 0.04 
Ca2+ mg/L 5.4 206 52.5 17.8 159 63.9 60.5 26.7 75 – 0.06 
Mg2+ mg/L 10.6 175 76.1 5.2 167 36.0 48 40.1 30 – 0.04 
Cl mg/L 13.7 642 172 3.5 514 66 97.7 111.7 250 250 0.08 
 mg/L 20.3 1,071 263 4.8 665 101 149.9 171.3 250 250 0.04 
 mg/L 207 882 555 120 768 360 418.9 167.9 600 – 0.04 
 mg/L 0.04 0.55 0.13 BDL 2.40 0.27 0.23 0.3 – 0.5a 0.08 
 mg/L 0.05 89.89 18.88 0.02 26.78 5.61 89.9 15.3 10 20 0.08 
 mg/L BDL 0.01 0.0016 BDL 0.05 0.0034 0.0029 0.01 0.9 1a 0.10 
F mg/L 0.76 4.74 1.71 0.15 1.78 0.78 1.06 0.8 1.5 0.10 
TH mg/L 93.2 1,058 442 69.5 825 312 351 195.5 200 450 0.10 

aIn the absence of WHO guidelines, Chinese guidelines are adopted.

BDL, below detection limit.

Table 3

Groundwater quality classification based on the WQI

WQIWater qualityNo. of sample (%)
NBSBTotal
<50 Excellent 1 (5.6) 13 (31) 14 (23.3) 
50–100 Good 9 (50) 25 (59.5) 34 (56.7) 
100–200 Average 6 (33.3) 4 (9.5) 10 (16.7) 
200–300 Fair 2 (11.1) – 2 (3.3) 
>300 Poor – – – 
WQIWater qualityNo. of sample (%)
NBSBTotal
<50 Excellent 1 (5.6) 13 (31) 14 (23.3) 
50–100 Good 9 (50) 25 (59.5) 34 (56.7) 
100–200 Average 6 (33.3) 4 (9.5) 10 (16.7) 
200–300 Fair 2 (11.1) – 2 (3.3) 
>300 Poor – – – 

Groundwater quality within the study area can be categorized into five distinct classes based on the WQI, WQI < 50, 50 ≤ WQI < 100, 100 ≤ WQI < 200, 200 ≤ WQI < 300, WQI ≥ 300 were regarded as (1) excellent; (2) good; (3) average; (4) fair; and (5) poor.

Human health risk assessment

A health risk evaluation is essential for determining the probability of adverse health effects in individuals who are exposed to chemical contaminants in groundwater. In this study, we utilized the empirical models established by the USEPA to assess the non-carcinogenic impacts of contaminants. The average daily intake of each element was determined using the subsequent equation (USEPA 1989):
(5)
where EDI is the exposed dose through intake of groundwater (mg/(kg day)); C is the concentration of a particular pollutant in groundwater (mg/L); IR is the groundwater ingestion rate (L/day); ED is the exposure duration (year); EF is the exposure frequency (days/year); BW is the average body weight (kg); ET is the average exposure time (days). Zhang et al. (2019b) conducted a HHRA of groundwater quality in the Jinghui Canal located in the Guanzhong Plain, utilizing the parameters he developed. The parameters is presented in Table 1.
The non-carcinogenic risk associated with pollutants in groundwater exposure was evaluated using the Health Quotient (HQ) for pollutants (, , , and F), as detailed in the corresponding equation (USEPA 1989):
(6)
where RfD is the reference dose for non-carcinogenic risk characterization, (mg/(kg day)). The RfD for , , , and F is 1.6, 0.16, 0.02, 0.06 mg/(kg day), respectively (USEPA 2024).
The HQ represents the potential risk of non-carcinogenic effects due to pollutants in groundwater. An HQ < 1 suggests that the groundwater is safe for consumption, while an HQ > 1 signifies a potential risk of adverse health impacts from consuming the groundwater. Assuming that the various pollutants in groundwater have no synergistic or antagonistic effects on human health, and only exhibit additive relationships, the total non-carcinogenic risk index should be the sum of the risk indices caused by multiple pollutants of concern. The calculation expression for this is presented in Equation (9). The Hazard Index (HI) serves as an indicator of the cumulative risk posed by , , , and F in groundwater, integrating the potential impacts of these pollutants. The higher the HI value, the higher the non-carcinogenic risk for human to the groundwater:
(7)

Distribution of hydrochemistry variables

Table 2 summarizes hydrochemistry from sample analyses. The groundwater in the study area exhibited predominantly mild alkalinity, with pH ranging from 7.1 to 8.8, averaging at 8.1. The TDS of groundwater was from 223 to 3,930 mg/L, with an average value of 984 mg/L. This study analyzed of hydrochemistry variables for groundwater samples in the south bank of the Wei River (NB) and the south bank of the Wei River (SB). The average value of TDS in the NB was 1,507 mg/L, which was greater than the average value of 760 mg/L in the SB. In the NB, 72.2% of the samples exceeded the limit value of TDS for drinking water (Ministry of Health of the PR China 2022; WHO 2022), and only 16.7% of the samples exceeded the limit value in the SB. The higher TDS in the NB may be due to the stronger water–rock interaction and evaporation concentration (Zhang et al. 2016).

The spatial variation of cations and anions values are shown in Figure 2(c–i). In the NB, the mean value of cation sequence was Na+ (279 mg/L) > Mg2+ (76.1 mg/L) > Ca2+ (52.5 mg/L) > K+ (6.52 mg/L). In the SB, the mean value of cation sequence was Na+ (90 mg/L) > Ca2+ (63.9 mg/L) > Mg2+ (36 mg/L) > K+ (10.79 mg/L). In both NB and SB, the abundance of anions followed by . Moreover, in the NB and SB samples ranged from 207 to 882 mg/L and 120–768 mg/L, with a mean of 555 mg/L (NB) and 360 mg/L (SB), respectively. The concentration varied from 20.3 to 1,071 mg/L (NB), and from 4.8 to 665 mg/L (SB), with a mean of 263 mg/L (NB) and 101 mg/L (SB). The Cl concentration varied from 13.7 to 642 mg/L (NB), and from 3.5 to 514 mg/L (SB), with mean of 172 mg/L (NB) and 66 mg/L (SB). A large amount of soft water was distributed in the SB, especially the TH value of groundwater near the Huyi was less than 300 mg/L.
Figure 2

Spatial variation of (a) TDS; (b) pH; (c) K+; (d) Na+; (e) Ca2+; (f) Mg2+; (g) Cl; (h) ; (i) ; (j) ; (k) ; (l) ; (m) F; and (n) and TH (the red line indicates the limit values of the WHO).

Figure 2

Spatial variation of (a) TDS; (b) pH; (c) K+; (d) Na+; (e) Ca2+; (f) Mg2+; (g) Cl; (h) ; (i) ; (j) ; (k) ; (l) ; (m) F; and (n) and TH (the red line indicates the limit values of the WHO).

Close modal

Fluoride and nitrogen contamination in groundwater

The values of fluoride in the SB ranged from 0.15 to 1.78 mg/L, averaged at 0.78 mg/L. The higher fluoride concentrations were found in the NB. The concentrations ranged from 0.76 to 4.74 mg/L, averaged at 1.71 mg/L. The results of the Mann–Whitney U test showed significant differences in fluoride concentrations in groundwater between the north and south banks of the Wei River (p = 0.000 < 0.01). Approximately 20.97% of all samples exceeded the WHO (2022) recommended reference value of 1.5 mg/L, while about 35.48% surpassed the national threshold of 1.0 mg/L (General Administration of Quality Supervision & Quarantine of China 2017). The groundwater samples in the northern part of Wugong and the northeastern part of Xingping had high fluoride levels above the safe limits (Figure 2(m)).

The average concentration of in the NB was 18.88 mg/L, which was significantly higher than that in the SB, which was 5.61 mg/L. The Mann–Whitney U test indicated significant differences in concentrations in groundwater between the north and south banks of the Wei River (p = 0.007 < 0.01). 50% of the samples in the NB, and 21.43% in the SB had high nitrate values exceeding the WHO reference value of 10 mg/L WHO (202). These groundwater samples were mainly distributed in western part of the study region (Figure 2(k)). The value of was between BDL and 0.05 mg/L in the study area, with an average of 0.003 mg/L. All the samples were below the limitation of 1.0 mg/L for WHO and Chinese standards. The concentration of in NB was 0.13 mg/L, which was lower than the concentration of 0.27 mg/L in SB. The Mann–Whitney U test indicated significant differences in concentrations in groundwater between the north and south banks of the Wei River (p = 0.022 < 0.05). A small portion of the points exceeding the standard levels were primarily located at the confluence of the Hei River and the Wei River.

Groundwater quality assessment

The WQI serves as an effective method for evaluating groundwater quality, simplifying extensive water quality data into a singular index to reflect the overall condition of groundwater (Gazzaz et al. 2012). By giving different weights (1–5) to the groundwater quality evaluation factors, the importance of each evaluation factor to the overall groundwater quality can be reflected. TDS, TH, nitrite, and fluoride are given higher weights and are considered to be an important part of the overall water quality (Xiao et al. 2023).

As shown in Table 3, the calculation results showed that 14 samples were excellent, most of which were located in Huyi. The 34 samples belonged to good, most of which were located in the piedmont of the northern foot of the Qinling Mountains in the SB. On the whole, 80% of the samples in the study area had WQI values less than 100, and 50 and 59.5% of the water samples in the NB and SB belonged to the good category, respectively. However, the water samples in the SB had a more prominent percentage of 31% in excellent category, representing the best water quality, compared with 5.6% of NB. Water samples with fair quality were generally found in the north bank of the Wei River. For instance, water sample 30 had the highest WQI at 272.5, with concentrations of TDS, Ca2+, Na+, and TH at 3,930, 206, 470, and 1,059 mg/L, respectively. In contrast, the water quality on the south bank of the Wei River in the study area was generally good, and this groundwater was suitable for drinking (Figure 3).
Figure 3

Water Quality Index (WQI) calculation zoning map of the study area.

Figure 3

Water Quality Index (WQI) calculation zoning map of the study area.

Close modal

Factors controlling groundwater hydrochemistry

The major ionic composition of groundwater is represented by the Piper diagram (Piper 1944). The groundwater of study area predominantly featured the HCO3–Ca·Mg type, as shown in Figure 4, shifting to HCO3–SO4–Ca·Na·Mg near Wugong in the NB and to HCO3–Na–Ca near Huyi in the SB.
Figure 4

Piper diagram of groundwater samples.

Figure 4

Piper diagram of groundwater samples.

Close modal
Based on the Gibbs diagram, the main governing processes of groundwater hydrochemistry include evaporation, precipitation, and water–rock interaction (Marandi & Shand 2018). The main driving forces in the hydrogeochemical evolution in the study area were rock weathering and evaporation dominance, as shown in Figure 5. Most of the groundwater samples located in the central parts of the sub-diagrams suggested that the chemical evolution of groundwater in the study area was mainly influenced by rock weathering and dissolution. According to Gao (2021), the groundwater depth in the north bank of the Wei River was less than 10 m. This indicates that evaporation had a stronger effect in the north bank, which is consistent with the observation.
Figure 5

Gibbs diagram of groundwater samples.

Figure 5

Gibbs diagram of groundwater samples.

Close modal
Figure 6(a) shows that almost all samples in the study area fall above the 1:1 line. It suggests that the study area is enriched in sodium ions, and that cation exchange interactions and weathering of the silicate minerals are the main factors in the enrichment of sodium ions (Adimalla et al. 2020). When gypsum dissolution is the only source of Ca2+ an d in groundwater, the ratio of should be 1:1. According to Figure 6(b), groundwater samples from SB in the study area are over-enriched in Ca2+, which may be due to the dissolution of carbonate minerals, such as calcite and dolomite (Morse & Arvidson 2002). Most of the groundwater samples from NB have a significant preponderance of in the bivariate diagram.
Figure 6

Plots of (a) Na+ vs. Cl; (b) vs. Ca2+; (c) ; (d) ; (e) ; (f) CAI-I vs. CAI-II; (g) ; (h) F vs. (Na+/ Ca2+); and (i) F vs. pH.

Figure 6

Plots of (a) Na+ vs. Cl; (b) vs. Ca2+; (c) ; (d) ; (e) ; (f) CAI-I vs. CAI-II; (g) ; (h) F vs. (Na+/ Ca2+); and (i) F vs. pH.

Close modal

When carbonates (e.g., calcite and dolomite) and sulfates (e.g., gypsum) are the dominant hydrologically acting base mineral types in the groundwater system, converge to a 1:1 trend line (Capaccioni et al. 2001). Figure 6(c) shows that almost all of the groundwater samples from NB fall above the 1:1 line, and the excess and may be due to weathering of silicate rock and cation exchange. Among (Figure 6(d)), the groundwater in the NB has more significant bicarbonate ions, indicating that the effect of calcite dissolution is secondary.

In Figure 6(e), most of the groundwater samples in the study area fall between the 1:1 and 2:1 line, indicating that dolomite dissolution is the most important source of calcium and magnesium ions.

The Chlor-Alkali Index is helpful to identify cation exchange direction. The negative values indicate Na+ and K+ increase, and Mg2+ and Ca2+ decrease, while positive values suggest the reverse process (Adimalla & Li 2019). Figure 6(f) illustrates that CAI-I values range from −37.83 to 0.17, while CAI-II values vary from −1.10 to 0.05. Predominantly, the indices for the groundwater samples are negative. This negative tendency significantly contributes to explaining the observed high Na+ and low Ca2+ concentrations in the samples, as delineated in Reaction 1.

The fine particles in the study region can provide a beneficial condition for cation exchange. In the NB, extensive distributions of loess and sandy deposits are prevalent. In contrast, in the SB, due to the proximity to the Qinling Mountains, the sediment particles are larger than the loess in the NB. Additionally, stronger surface runoff from the mountainous areas has eroded away the fine particles from the surface sediments:
(8)
(9)
(R1)

Fluoride and nitrogen contamination in groundwater

Fluoride naturally occurs in groundwater as a result of the dissolution of fluoride-bearing minerals, such as fluorite (CaF2) (Li et al. 2019). Excessive fluoride content in drinking water can produce dental fluorosis, and long-term intake of high-fluoride drinking water can also lead to fluorosis (Chen et al. 2021). Fluoride enrichment in shallow groundwater is influenced by a variety of factors. Under arid conditions, reduced groundwater flow extends the residence time with geological formations and enhances fluoride enrichment through evaporation (Brunt et al. 2004). Additionally, an environment with higher pH levels contributes to the enhancement of fluoride concentrations; the pH values for most of the samples were found to be between 7.8 and 8.8 (Figure 6(g)). The exchange of ions between Ca2+ and Na+ can further facilitate the dissolution and movement of CaF2 as Reaction 3 (Chen et al. 2021). This process leads to an increase in fluoride levels in groundwater. The Figure 6(h) shows the observed positive relationship between F and the Na+/ Ca2+ ratio suggests that cation exchange plays a role in the enrichment of F in groundwater. Elevated levels of can accelerate the dissolution of CaF2 (Figure 6(i)), leading to high F in groundwater (Reaction 4):
(R2)
(R3)
(R4)
Human activities are considered the primary factors influencing the concentration of nitrogenous pollutants in groundwater, such as the use of chemical fertilizers, aquaculture, and agricultural grazing (Laftouhi et al. 2003). The study area, with its long history as an agricultural irrigation district, identifies the use of chemical fertilizers as a significant source of nitrogenous pollutants in groundwater. The primary nitrogen compound in chemical fertilizers is urea (CO (NH2)2), which is hydrolyzed by the enzyme urease and subsequently converted into nitrate in the unsaturated zone of groundwater as Reactions (5–7) (Stadler et al. 2008). The hydrolysis of urea results in a temporary increase in pH, which facilitates the formation of ammonium (Widory et al. 2004), it can be verified by Figure 2(b) and 2(k). In the NB, Wugong and Xingping mainly produce wheat, while Huyi and Zhouzhi in the SB focus on fruit cultivation, especially kiwi fruits. Some scholars indicated that the application rate of nitrogen fertilizer was 750 kg/hm in the kiwifruit land, which was much higher than that of wheat 199.2 kg/hm (Lu et al. 2016; Zhao et al. 2016). Given intensive irrigation activities, nitrogen can infiltrate into groundwater system:
(R5)
(R6)
(R7)

Assessment of the hazard to human health

The HHRA developed by the USEPA (USEPA 1989) is extensively applied to evaluate the health risks associated with groundwater quality. The non-carcinogenic risk coefficient (HQ) ranges of , , and F in children were 0–6.79, 0–0.02, 0–3.18, 0.14–4.48, respectively. The non-carcinogenic risk coefficient (HQ) ranges of , , and F in adults were 0–2.91, 0–0.01, 0–1.36, 0.06–1.92, respectively. The non-carcinogenic risks of and have the same spatial distribution characteristics. The higher risk areas are distributed in the east of Zhouzhi and the west of Huyi. In the study region, the cumulative non-carcinogenic risk coefficient (HI) of adult through oral intake ranged from 0.22 to 3.19. The mean value was 0.86, which was less than the recommended reference value of 1.0 recommended by USEPA. The cumulative non-carcinogenic risk (HI) of oral intake in children was between 0.51 and 7.44 with the mean value of 2.0. The HI for 88.33% of groundwater samples in the study area exceeded the USEPA-recommended reference value of 1.0. The non-carcinogenic risk level for children is about 2.33 times that of adults. The HHR results of Zhang et al. (2019a) for groundwater in the Guanzhong Plain also found that children are exposed to greater potential risks. Therefore, the high non-carcinogenic risk associated with oral intake in the study area warrants close attention, particularly for children (Figure 7).
Figure 7

Non-carcinogenic human health risk assessment (the red line signifies the critical HQ and HI values).

Figure 7

Non-carcinogenic human health risk assessment (the red line signifies the critical HQ and HI values).

Close modal

The contribution rates of , , and F to the cumulative non-carcinogenic risk (HI) of oral intake in the study area were 17.01%, 0.05%, 32.82% and 50.12%, respectively. The order of non-carcinogenic risk of each index for human body was , and the cumulative contribution rate of F and to the cumulative non-carcinogenic risk (HI) of oral intake was as high as 82.94%. It can be considered that the concentration of F and plays a controlling role in the non-carcinogenic health risk in the study area. In particular, the spatial distribution of the cumulative non-carcinogenic risk coefficient HI is similar to the spatial distribution of fluoride. This confirms the findings of Chen et al. (2021) of higher health risks associated with fluoride on the northern bank of the Wei River in the Guanzhong Plain. Higher HI and higher HQ-F are distributed near Wugong and Xingping on the north bank of the Wei River. Excessive intake of fluorine will interfere with the activity of various enzymes in the body, destroy the metabolic balance of calcium and phosphorus, and cause fluorosis of teeth, bone and joint deformation.

The study detailed analysis of groundwater in the region has unveiled critical insights, underscoring the urgency for targeted environmental and health interventions. Groundwater in the study area shows mild alkalinity. TDS averaged 1,507 mg/L in the NB, higher than the value of SB (760 mg/L). In the NB, 72.2% of samples surpassed the drinking water TDS limit, versus 16.7% in the SB. Overall, the majority of groundwater samples in the study area were characterized by the HCO3–Ca ·Mg water type. Near Wugong County, located in the NB, the predominant water type was HCO3·SO4–Ca ·Na·Mg. Meanwhile, in the northwest of Huyi County in the SB, the dominant water type was HCO3–Na–Ca.

Fluoride in groundwater primarily originates from the natural weathering of fluoride-bearing minerals such as fluorite (CaF₂). The concentration of fluoride is influenced by several factors including the interaction between water and rock, higher pH levels that enhance fluoride solubility, ion exchange processes involving calcium and sodium, and evaporation concentration in arid conditions. In areas where fluoride levels exceed standards within the research region, policies can be implemented to adopt water treatment measures and increase regular monitoring to reduce the impact of fluoride on human health. Nitrogen contamination, mainly in the form of nitrates, originates from agricultural activities, especially the use of nitrogen-based fertilizers, as well as livestock farming and wastewater discharge. The concentration of nitrogen in groundwater is controlled by factors such as excessive fertilizer application, irrigation practices that promote leaching, shallow groundwater tables, inadequate sanitation infrastructure, and the biochemical processes of nitrification and denitrification in the soil and water system. Policy improvements to address these issues could include stricter regulations on fertilizer use, enhanced wastewater treatment standards, and initiatives to restore and protect soil and water ecosystems.

The study reveals 80% of samples having a WQI below 100, indicating suitability for drinking. A significant portion of good-quality samples is found near the Qinling foothills, suggesting lateral recharge into the area. Groundwater in the SB outperforms groundwater in the NB in excellent water quality percentage (31 vs. 5.6%), likely due to the Qinling Mountains continuous supply, indicating that most groundwater in the area is above average quality, with SB having notably better quality.

This study highlights that children in a specific area face higher non-carcinogenic health risks from groundwater contaminants such as , , , and F, compared to adults. The cumulative health risk (HI) for children often exceeds the USEPA safety threshold, with 88.33% of samples above the reference value, indicating significant health risks. Fluoride and ammonia nitrogen are the primary risk contributors, making up 82.94% of the cumulative HI. High-risk areas include Wugong and Xingping.

The findings of the present study can provide insights into the effect of natural factors and human activities on fluoride and nitrogen contamination in groundwater and provide a basis for further water management at a typical agricultural region.

We thank the associate editor and the reviewers for their useful feedback that improved the quality of this paper.

This research is supported by grants from the National Natural Science Foundation of China (42272289 and 42341102), the Natural Science Foundation of Shaanxi Province (2021JQ-862), Youth Talents Promotion Plan of Shaanxi Association for Science and Technology (20220711), and the Fundamental Research Funds for the Central Universities in Chang'an University (300102293205, 300102294905, 300102294906, and 300102294902).

J.Y. contributed to manuscript writing, sample collection, data analysis, result interpretation, and data generation. J.W. contributed to content review and data analysis. H.X contributed to data collection and for preparing cartography. Z.X. contributed to data collection and sample collection. Y.Z. contributed to language editing. J.C. contributed to funding acquisition and content review

Data cannot be made publicly available; readers should contact the corresponding author for details.

The authors declare there is no conflict.

Adimalla
N.
,
Dhakate
R.
,
Kasarla
A.
&
Taloor
A. K.
(
2020
)
Appraisal of groundwater quality for drinking and irrigation purposes in Central Telangana, India
,
Groundwater for Sustainable Development
,
10
,
100334
.
Brunt
R.
,
Vasak
L.
&
Griffioen
J.
(
2004
)
Fluoride in Groundwater: Probability of Occurrence of Excessive Concentration on Global Scale
.
Igrac
.
Chen
J.
,
Qian
H.
,
Gao
Y.
,
Wang
H.
&
Zhang
M.
(
2020
)
Insights into hydrological and hydrochemical processes in response to water replenishment for lakes in arid regions
,
Journal of Hydrology
,
581
,
124386
.
Dobrinas
S.
,
Soceanu
A.
,
Manea
N.
,
Sirbu
A.
,
Dumitrescu
C. I.
,
Popescu
V.
,
Birghila
S.
,
Matei
N.
&
Popovici
I. C.
(
2022
)
Health risk assessment of fluoride exposure due to groundwater consumption in Romania
,
Journal of Water and Health
,
20
,
1380
1392
.
Gao
Y.
(
2021
)
Spatio-Temporal Evolution of Hydrochemical Components and Hunman Health Risk Assessment of Groundwater in Guanzhong Plain
.
PhD
,
Chang'an University
.
Gao
H.
,
Yang
L.
,
Song
X.
,
Guo
M.
,
Li
B.
&
Cui
X.
(
2023
)
Sources and hydrogeochemical processes of groundwater under multiple water source recharge condition
,
Science of The Total Environment
,
903
,
166660
.
Gazzaz
N. M.
,
Yusoff
M. K.
,
Aris
A. Z.
,
Juahir
H.
&
Ramli
M. F.
(
2012
)
Artificial neural network modeling of the water quality index for Kinta River (Malaysia) using water quality variables as predictors
,
Marine Pollution Bulletin
,
64
,
2409
2420
.
General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China
(
2017
)
Standards for Groundwater Quality (GB/T 14848-2017)
.
Standards Press of China Beijing
.
Guo
Y.
&
Wang
J.
(
2021
)
Spatiotemporal changes of chemical fertilizer application and its environmental risks in China from 2000 to 2019
,
International Journal of Environmental Research and Public Health
,
18
,
11911
.
Ha
Q. K.
,
Ngoc
T. D. T.
,
Vo
P. L.
,
Nguyen
H. Q.
&
Dang
D. H.
(
2022
)
Groundwater in Southern Vietnam: Understanding geochemical processes to better preserve the critical water resource
,
Science of the Total Environment
,
807
,
151345
.
Horton
R. K.
(
1965
)
An index number system for rating water quality
.
Journal of the Water Pollution Control Federation
,
3
,
300
306
.
Ischinger
L. S.
(
1979
)
Nonpoint source pollution
,
Fisheries
,
4
,
50
52
.
Kazi
T. G.
,
Arain
M. B.
,
Jamali
M. K.
,
Jalbani
N.
,
Afridi
H. I.
,
Sarfraz
R. A.
,
BAIG
J. A.
&
Shah
A. Q.
(
2009
)
Assessment of water quality of polluted lake using multivariate statistical techniques: A case study
,
Ecotoxicology and Environmental Safety
,
72
,
301
309
.
Khan
F.
,
Krishnaraj
S.
,
Raja
P.
,
Selvaraj
G.
&
Cheelil
R.
(
2021
)
Impact of hydrogeochemical processes and its evolution in controlling groundwater chemistry along the east coast of Tamil Nadu and Puducherry, India
,
Environmental Science and Pollution Research
,
28
,
18567
18588
.
Laftouhi
N.-E.
,
Vanclooster
M.
,
Jalal
M.
,
Witam
O.
,
Aboufirassi
M.
,
Bahir
M.
&
Persoons
É
. (
2003
)
Groundwater nitrate pollution in the Essaouira Basin (Morocco)
,
Comptes Rendus. Géoscience
,
335
,
307
317
.
Lu
Y.-L.
,
Zhang
T.-T.
,
Zhang
X.-J.
,
Gao
J.-B.
,
Chen
Z.-J.
&
Zhou
J.-B.
(
2016
)
Evaluation of current fertilization status in kiwifruit orchards on the northern slope of Qinling Mountains: A case study of Yujiahe catchment, in Zhouzhi county
,
Journal of Plant Nutrition and Fertilizer
,
22
,
380
387
.
Lumb
A.
,
Sharma
T. C.
&
Bibeault
J.-F.
(
2011
)
A review of genesis and evolution of Water Quality Index (WQI) and some future directions
,
Water Quality, Exposure and Health
,
3
,
11
24
.
Marandi
A.
&
Shand
P.
(
2018
)
Groundwater chemistry and the Gibbs Diagram
,
Applied Geochemistry
,
97
,
209
212
.
Maurya
J.
,
Pradhan
S. N.
,
Seema
&
Ghosh
A.
(
2020
)
Evaluation of ground water quality and health risk assessment due to nitrate and fluoride in the Middle Indo-Gangetic plains of India
,
Human and Ecological Risk Assessment: An International Journal
,
27
,
1349
1365
.
Meng-Yao
H.
,
Zhang
L.
,
Zhi-Wen
W.
,
Dian-Lin
Y.
,
Li-Li
W.
,
Wei-Ming
X.
&
Jian-Ning
Z.
(
2017
)
Estimation of fertilizer usage from main crops in China
,
Journal of Agriculture Resources and Environment
,
34
,
360
.
Ministry of Health of the People's Republic of China, S.O.P.
(
2022
)
Standards for Drinking Water Quality (GB 5749–2022)
.
Beijing
:
China Standard Press
.
Morse
J. W.
&
Arvidson
R. S.
(
2002
)
The dissolution kinetics of major sedimentary carbonate minerals
,
Earth-Science Reviews
,
58
,
51
84
.
Piper
A. M.
(
1944
)
A graphic procedure in the geochemical interpretation of water-analyses
,
Eos, Transactions American Geophysical Union
,
25
,
914
928
.
Ren
X.
,
Li
P.
,
He
X.
,
Su
F.
&
Elumalai
V.
(
2021
)
Hydrogeochemical processes affecting groundwater chemistry in the central part of the Guanzhong Basin, China
,
Archives of Environmental Contamination and Toxicology
,
80
,
74
91
.
Scanlon
B. R.
,
Fakhreddine
S.
,
Rateb
A.
,
De Graaf
I.
,
Famiglietti
J.
,
Gleeson
T.
,
Grafton
R. Q.
,
Jobbagy
E.
,
Kebede
S.
&
Kolusu
S. R.
(
2023
)
Global water resources and the role of groundwater in a resilient water future
,
Nature Reviews Earth Environment
,
4
,
87
101
.
Shaji
E.
,
Sarath
K.
,
Santosh
M.
,
Krishnaprasad
P.
,
Arya
B.
&
Babu
M. S.
(
2024
)
Fluoride contamination in groundwater: A global review of the status, processes, challenges, and remedial measures
,
Geoscience Frontiers
,
15
,
101734
.
Stadler
S.
,
Osenbrück
K.
,
Knöller
K.
,
Suckow
A.
,
Sültenfuß
J.
,
Oster
H.
,
Himmelsbach
T.
&
Hötzl
H.
(
2008
)
Understanding the origin and fate of nitrate in groundwater of semi-arid environments
,
Journal of Arid Environments
,
72
,
1830
1842
.
Tzanakakis
V. A.
,
Paranychianakis
N. V.
&
Angelakis
A. N.
(
2020
)
Water supply and water scarcity
,
Water
,
12
,
2347
.
USEPA
(
1989
)
Risk Assessment Guidance for Superfund, volume1: Human Health Evaluation Manual (Part A)
.
Washington DC
:
Office of emergency and remedial response, USEPA
.
USEPA
(
2004
)
Risk Assessment Guidance for Superfund, Volume I: Human Health Evaluation Manual (Part E, Supplemental Guidance for Dermal Risk Assessment)
.
EPA/540/R/99
.
Washington DC
.
USEPA
(
2024
)
Integrated Risk Information System (IRIS)
.
Vithanage
M.
&
Bhattacharya
P.
(
2015
)
Fluoride in the environment: Sources, distribution and defluoridation
,
Environmental Chemistry Letters
,
13
,
131
147
.
WHO
(
2022
)
Guidelines for Drinking-Water Quality: Incorporating the First and Second Addenda
, 4th edn.
Geneva
:
World Health Organization
.
Widory
D.
,
Kloppmann
W.
,
Chery
L.
,
Bonnin
J.
,
Rochdi
H.
&
Guinamant
J.-L.
(
2004
)
Nitrate in groundwater: An isotopic multi-tracer approach
,
Journal of Contaminant Hydrology
,
72
,
165
188
.
Xiao
S.
,
Fang
Y.
,
Chen
J.
,
Zou
Z.
,
Gao
Y.
,
Xu
P.
,
Jiao
X.
&
Ren
M.
(
2023
)
Assessing the hydrochemistry, groundwater drinking quality, and possible hazard to human health in Shizuishan Area, Northwest China
,
Water
,
15
,
1082
.
Xu
P.
,
Zhang
Q.
,
Qian
H.
,
Li
M.
&
Hou
K.
(
2019
)
Characterization of geothermal water in the piedmont region of Qinling Mountains and Lantian-Bahe Group in Guanzhong Basin, China
,
Environmental Earth Sciences
,
78
,
1
17
.
Zhang
X.
,
Qian
H.
,
Wu
H.
,
Chen
J.
&
Qiao
L.
(
2016
)
Multivariate analysis of confined groundwater hydrochemistry of a long-Exploited sedimentary basin in northwest China
,
Journal of Chemistry
,
2016
,
1
15
.
Zhang
Q.
,
Xu
P.
&
Qian
H.
(
2019a
)
Assessment of groundwater quality and human health risk (HHR) evaluation of nitrate in the central-Western Guanzhong Basin, China
,
International Journal of Environmental Research and Public Health
,
16
,
4246
.
Zhang
Y.
,
Xu
B.
,
Guo
Z.
,
Han
J.
,
Li
H.
,
Jin
L.
,
Chen
F.
&
Xiong
Y.
(
2019b
)
Human health risk assessment of groundwater arsenic contamination in Jinghui irrigation district, China
,
Journal of Environmental Management
,
237
,
163
169
.
Zhao
H.-B.
,
Wang
Z.-H.
,
Gao
Y.-J.
&
Zhang
W.-F.
(
2016
)
Investigation and evaluation of household wheat fertilizer application in Shaanxi Province
,
Journal of Plant Nutrition and Fertilizer
,
22
,
245
253
.
Zou
Z.-H.
,
Zhao
B.
,
Chen
J.
&
Qian
H.
(
2024
)
Assessment of water resources carrying capacity and temporal-spatial variation characteristics of nine provinces along the Yellow River
,
Journal of Earth Sciences and Environment
,
46
(
3
),
351
363
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).