Next Article in Journal
Speciation Distribution Characteristic and Ecological Risk of Heavy Metals in Surface Sediments of Cascading Hydropower Dams in Lancang River
Next Article in Special Issue
Chemical and Isotopic Features of a High pCO2 Natural Mineral Water from Essentuki Field (Caucasian Mineral Water Region, Russia)
Previous Article in Journal
Effect of Attached Growth on Treatment Performance in Waste Stabilization Ponds
Previous Article in Special Issue
Comparison of Microbial Profiling and Tracer Testing for the Characterization of Injector-Producer Interwell Connectivities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genesis Mechanism and Health Risk Assessment of High Boron Water in the Zhaxikang Geothermal Area, South Tibet

1
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu 610059, China
2
Institute of Sedimentary Geology, Chengdu University of Technology, Chengdu 610059, China
3
College of Energy Resources, Chengdu University of Technology, Chengdu 610059, China
4
Chengdu Center, China Geological Survey, Chengdu 610081, China
5
Jiangsu Provincial Key Laboratory of Environmental Engineering, Jiangsu Provincial Academy of Environmental Science, Nanjing 210036, China
6
China Institute of Geo-Environment Monitoring, Beijing 100081, China
7
School of Engineering and Technology, China University of Geosciences, Xueyuan Road 29, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(20), 3243; https://doi.org/10.3390/w14203243
Submission received: 7 September 2022 / Revised: 5 October 2022 / Accepted: 6 October 2022 / Published: 14 October 2022
(This article belongs to the Special Issue Hydrochemistry and Isotopes in Groundwater Investigations)

Abstract

:
The natural discharge of geothermal water containing harmful components affects the water quality of the surrounding environment and brings security risks to drinking water safety. The geothermal water in Tibet is characterized by high boron content, but the water pollution caused by the discharge of this high boron geothermal water is not clear. In this study, we collected geothermal water and surface water from the Zhaxikang geothermal system in southern Tibet to investigate the causes of high boron geothermal water and the water pollution of water quality by its discharge. The results indicate that the hydrochemical type of geothermal water was HCO3-Cl-Na, while that of cold spring water, mine water, river water, and lake water was SO4-HCO3-Ca-Mg. Hydrogen and oxygen isotopes show that the recharge source of cold groundwater was mainly snow-melting water and meteoric water, while in addition to that, there is magmatic water for hot springs. The boron content of geothermal water in the study area is as high as 42.36 mg/L, far exceeding the World Health Organization limit for drinking water (0.5 mg/L). The analysis of ion components and PHREEQC modeling indicated that the dissolution of silicate minerals and cation exchange controlled the composition of groundwater, and the boron in groundwater mainly came from the volatilization of magmatic components and the leaching of shallow sediments. The entropy weight water quality index was used to evaluate the water quality of the study area; about 42.9% of the groundwater samples are of good quality and can be used for drinking, mainly cold water that has not been mixed with geothermal water in the upstream. With the discharge of geothermal water into the river (with a mix ratio of ~20%), the downstream water quality gradually deteriorated. The health risk assessment of drinking water in the study area showed that the hazard index (HI) of drinking water in the mixed area was higher than 1 (with an average of 1.594 for children and 1.366 for adults), indicating that children are at a higher health risk than adults. Geothermal water with high boron content has been found all over the world, and the adverse effects of its natural drainage cannot be ignored.

1. Introduction

The Qinghai–Tibet Plateau is well known as the “Asian Water Tower” (AWT) and is the source of the Yarlung Zangbo (Brahmaputra), Indus, Ganges, Mekong, Yangtze, and Yellow rivers. Its water environment is essential to the security of drinking water for nearly 1.4 billion people in more than a dozen countries in Asia [1]. In addition, the Qinghai–Tibet Plateau is an important part of the Mediterranean–Himalayan Geothermal Belt, with lots of geothermal systems and numerous surface thermal manifestations [2,3,4,5,6,7]. Additionally, under the trend of energy transformation, geothermal energy as clean energy will undoubtedly usher in prosperity, but some potential water security risks cannot be ignored [8]. Compared with other types of natural water bodies (such as shallow groundwater and surface water), geothermal water with higher temperatures can usually leach more chemical components from its host medium. The typical harmful components include CO2 and H2S, as well as mercury, arsenic, fluorine, lead, chromium, boron, and radon [9]. Sources of these toxic elements are related to hydrochemical processes such as water–rock interaction, mixing, and degassing [10,11]. Under the influence of various geological or human factors, such as hot spring water discharge and geothermal water recharge, geothermal water is discharged into shallow aquifers, surface water, and the surrounding environment, bringing hidden dangers to the safety of surrounding drinking water [12,13]. Additionally, the rich geothermal resources in the Qinghai–Tibet Plateau will be accelerated exploitation driven by the long-term goal of carbon neutrality [14]. Therefore, the security of the water environment in the AWT and the assessment of potential health risks have become more significant as the anthropogenic accelerate the discharge of harmful components from geothermal water.
Geothermal fluids are usually a two-phase mixture with complex compositions, and the characteristic toxic and harmful substances carried by different geothermal fields are different [15]. Those toxic and harmful components will have a certain impact on the water quality around the drainage area. At present, the research on water quality in the Qinghai–Tibet Plateau mainly focuses on the effects and characteristics of geothermal water on surface water. Tian et al. [16] investigated the effects of arsenic, selenium, and fluoride on the quality of natural river water. Li et al. [17] found heavy metals such as arsenic, chromium, and lead in several rivers on the Qinghai–Tibet Plateau. Wang et al. [18] used the pollution evaluation model to study the heavy metal pollution of soil and groundwater in Tibet. However, the impact of high boron geothermal water in Tibet on the water environment has received little attention [19,20].
High boron groundwater has been reported worldwide, such as in Greece [21], Bangladesh and the US [22], Italy [23], Ghana [24], and China [25]. Boron concentrations in groundwater are generally low; however, geological processes such as geothermal activity [25], seawater intrusion into the aquifer [23], and human activities [26,27] would increase its concentration. The standard limit for boron in drinking water is 0.5 mg/L due to potential health risks [28]. High concentrations of B in drinking water can lead to low fertility and infertility, low fetal weight, and acute neurological and metabolic disorders [24]. Therefore, water safety risks associated with the discharge of high boron geothermal water cannot be ignored.
In China, high boron geothermal water is mainly distributed in the Yunnan–Tibet tropics [29]. In remote areas, residents often take shallow groundwater and river water as drinking water, while hot spring water is often discharged into river water without treatment. Therefore, we believe that the discharge of geothermal water with a high concentration of elements such as B and As will affect the surrounding water quality and bring hidden dangers to the drinking water safety of local residents.
The Zhaxikang geothermal system was chosen for this research. The main sources of drinking water in this area are river water and shallow groundwater, but the health risks of these water resources have not been reported. Based on the hydrochemical and isotopic characteristics of hot springs, cold springs, river water, cold groundwater, and lake water, the main aims of this study are (1) to identify and characterize the main components in hot springs, cold groundwater, and surface water; (2) to identify the source of major contaminants; (3) to use entropy weight water quality index (EWQI) to assess the water quality of the river and groundwater, and (4) to evaluate the potential health risks of boron in geothermal water to local drinking water. The results are of great significance for the security and management of water resources in this region.

2. Materials and Methods

2.1. Study Area

The study area is located in the southern Tibetan and the northeastern part of the Himalayan mountains (E 91°55′–92°05′, N 28°15′–28°25′, Figure 1b). The Zhaxikang catchment area spans Longzi and Cuona counties in the Shannan district. Zhaxikang’s terrain is high in the west and low in the east, with elevations ranging from 4300 m to 5600 m. This region is characterized by a plateau temperate semiarid monsoon climate.
The study area is located at the intersection of the nearly east–west strike Zhegucuo–Longzi thrust belt and the nearly north–south extension Cuona–Woka rift, which is sandwiched by the Yarlung Zangbo Suture Zone and the southern Tibet Detachment System (Figure 1a) [30]. A series of near north–south trending extensional and torsional faults are arranged in an echelon in the study area (Figure 1c), and they are inferred to be the main channels for geothermal water to rise to the surface. Zhaxikang is an important metallogenic area in southern Tibet [31] and a large Sb-Pb-Zn-Ag deposit in the North Himalayas of southern Tibet [32]. Under the influence of artificial mining and geological processes, As and other metal elements will be released into the surrounding environment [33].
There are clusters of springs in an area of about 400 m2 [34]. The temperature of the hot spring is from 43 °C to 70 °C, with a total flow rate of about 2 L/s. There is a slight sulfur odor, and the hot water eventually flows into the Langmai River. Zhaxikang thermal area’s reservoir temperature ranges between 185 °C and 205 °C [34].

2.2. Sampling and Analysis

In September 2018, water samples were collected in the study area. To analyze the spatial evolution of water quality in the Zhaxikang catchment, sampling points were placed at various potential groundwater and surface water interaction areas, as well as downstream of sewage discharge from mining areas that may cause pollution. For the mixing area between geothermal water and river water, sampling occurred upstream, in the mixing area, and downstream along the river. A total of 21 water samples were collected, including 11 surface water samples (4 from lakes in the southwest of the study area and the rest from various sections of the Langmai River), 4 hot spring water samples, 3 cold spring water samples, and 3 mine water samples (mine water at a depth of 200–400 m, representing cold groundwater). The Zhaxikang hot spring flows through the valley and then forms a thermal water pool that drains into the Langmai River. The distribution of sampling points for different types of water is shown in Figure 1c.
Each water sample was collected in a polyethylene bottle, which was rinsed three times with deionized water before sampling. The water sample for anion analysis was filtered through a 0.45 μm microporous membrane and stored in bottles. The sample used for cation analysis was also filtered and acidified with 6 N HNO3 until pH < 2 and sealed for storage [35,36]. The water samples for hydrogen and oxygen isotope analysis were unfiltered and stored directly in polyethylene sampling bottles. In the field, a portable multiparameter (Hach HD40Q, Loveland, USA) was used to determine the temperature (°C), pH, total dissolved solids (TDS), oxidation-reduction potential (ORP), and electrical conductivity [37].
The water samples were tested and analyzed at the Institute of Geographic Sciences and Natural Resources, CAS (Chinese Academy of Sciences) within one month after the field works. Ion chromatography (ICS-1100, ThermoFisher Scientific, Waltham, USA) was used to measure anion concentrations, with an A SUPP5 150 anion separation column, conductance detector Suppressed CD (J002), and the eluents were 3.2 mmol/L Na2CO3 and 1.0mmol/L NaHCO3, with a detection limit of 0.01 mg·L−1. Inductively Coupled Plasma-Atomic Emission Spectrometry was used to determine the cations (ICP-AES, ThermoFisher Scientific, Waltham, USA). The cationic separation column was Metrosep C4-150/4.0, and the suppressed conductivity detector CD (J002), the leachate is 1.7 mmol/L HNO3 and 0.7 mmol/L C7H5NO4, and the detection limit was 0.01 mg·L−1. Atomic fluorescence spectroscopy was used to measure Fe3+, Cu2+, Pb2+, Cr3+, Ni2+, Cd2+, and As5+, with a detection limit of 0.06 mg/kg. Hydrogen and oxygen stable isotopes were determined using a laser water isotope analyzer (L1102-I, Picarro, Santa Clara, USA) and standard samples corresponding to the Vienna Standard Mean Ocean Water (VSMOW). The detection accuracy for δ18O and δD is ±0.1‰ and ±1‰, respectively. The ion charge balance error (%CBE) of each sample was calculated, and the results show that the %CBE of the cold groundwater and surface water samples is less than 5%, and that of hot spring water samples %CBE is less than 10%. This confirms the reliability of the ion analysis [38].

2.3. Data Processing and Analysis

This study uses the Entropy-weighted water quality index (EWQI) to describe the water quality characteristics of the study area. The calculation steps and results of EWQI are shown in Supplementary Tables S1 and S2 [39,40]. As a convenient and effective water quality assessment method, EWQI has been widely used by scholars all over the world (Supplementary Table S3) [39].
To assess the health risk of ground and surface water mixed with geothermal water, this study quantified the health risks from drinking water intake and dermal contact using the empirical model proposed by the United States Environmental Protection Agency [41], and the adverse health effects were assessed at various ages (children and adults). The detailed calculation process and result can be found in Supplementary Table S4.
In this study, the B element was used as a tracer to calculate the mixing ratio of the river water near the upstream hot spring cluster. See Supplementary Table S4 for the calculation method. The activity of dissolved ions and the saturation index (SI) for various mineral phases were computed to analyze the equilibrium status of groundwater by using aqueous geochemical modeling software (PHREEQC, Version 3.6, Reston, USA), based on the WATEQ4F database [42]).

3. Results and Discussion

3.1. Hydrochemical Characteristics

Table 1 and Table 2 show the statistical data of the physical and chemical parameters of 21 water samples in the Zhaxikang area and the standard acceptable limits for the drinking water standards of the WHO. It can be seen from Table 2 that the concentration of various ions in water exceeds the drinking water limits in the study area (Figure 2a,b).
The spring water is characterized by high concentrations of boron and bicarbonate (Figure 2a). This is a common feature of Tibetan hot springs [29]. The ORP value of hot spring water is negative, which is lower than that of other types of water, indicating reductive conditions for the deep hot water. Compared with cold spring water, hot spring water can leach out more surrounding rock components and has a higher ion concentration.
The positive ORP value of the lake water indicates an oxidizing environment. Additionally, the average temperature of the lake is 15 °C. Compared with other types of surface water, lake water has weaker discharge conditions and stronger evaporation.
The average pH value of the mine well water is 8.09, slightly alkaline. Sample M02, with an ORP value of −41.0 mV, is located downstream of the mining area, indicating a decrease in groundwater flow rate. The higher Fe concentration of M02 may be due to the less precipitation process under reduction conditions.
The variation in the hydrochemical parameters of the river water samples is large but regular, indicating the mixing of river water from other sources. Along the flowing direction from R04 to R05 (Figure 1c), the water temperature gradually increases, and the concentration of B also increases significantly. This is consistent with the results of the influence of groundwater containing arsenic, fluorine, and boron on drinking water in Argentina studied by Hudson-Edwards et al. [45].

3.2. Evolution Mechanism of Hydrochemical Components

3.2.1. Source of Geothermal Water

The relationship between hydrogen and oxygen isotopes in groundwater can be used to infer its origin [34,46]. The relationship between δ18O and δD in the study area is shown in Figure 3.
Surface water and groundwater are commonly distributed near the Global Meteoric Water Line (GMWL), with meteoric water as the recharge source [47]. The isotope data representing ice–snow meltwater collected from the Nyenqing–Tanggula Mountain at an altitude of approximately 5500 m are also shown in Figure 3. [43], and magmatic water with δD of 20‰ ± 10‰ and δ18O of 10‰ ± 2‰ from Giggenbach [44]. All of the hot spring samples were distributed within the mixing area of ice–snow melting water and magma water, indicating that they were recharge sources of geothermal water. The δ18O and δD values of S05 and S02 are lower and closer to GMWL, and the other spring water samples have higher values, indicating that S05 and S02 are supplied by a greater proportion of meteoric water, whereas other springs have experienced greater deep source mixing during migration [48]. According to Figure 3, ice melting and local meteoric waters are the recharge source of river water and mine water. The high δ18O and δD values of water collected from the lake, especially L02 and L03, indicate that the lake water occurred with strong evaporation and formed an evaporation line with R2 = 0.9913.

3.2.2. Hydrochemical Types and Water–Rock Interaction

According to the piper ternary diagram (Figure 3), the hot spring water in the study area is located in zone III of the diagram and has the hydrochemical type HCO3-Cl-Na. Hydrochemical types of mine water and river water are SO4-HCO3-Ca-Mg (Figure 4).
The hot spring water samples fall in region D of the cationic triangle on the left, indicating that the cations are dominated by Na+. The lake and river water samples fall in zones A and B, suggesting that the cations in the surface water are mainly nondominant and calcium types. It can be seen that HCO3 is the main anion, and the lake water was dominated by sulfuric acid type, which could be affected by evaporation.
The ion concentration in groundwater is determined by the hydrogeochemical processes that occur in the aquifer system [8,49]. Gibbs [50] classified the natural factors that control the chemical characteristics of groundwater as precipitation, evaporation, and water–rock interactions. The Gibbs diagram of the Zhaxikang area is shown in Figure 4.
In the study area, except for a few points that fall outside the anion ratio, almost all samples are controlled by water–rock action within the Gibbs map (Figure 5a). In the Gibbs cation ratio (Figure 5b), most samples are controlled by water–rock interaction, and some hot spring samples fall near the region controlled by evaporation. As a result, the chemical evolution of groundwater is primarily controlled by water–rock interaction and partly by evaporation. Furthermore, the increase in the Na+/(Na+ + Ca2+) ratio with increasing TDS indicates that the groundwater may be altered by cation exchange between Na+ and Ca2+ [24].
Through the above analysis, the types of water in the study area are HCO3-Cl-Na and SO4-HCO3-Ca-Mg. Their water chemical composition is controlled by water–rock interaction and evaporation, and it is necessary to further analyze the water–rock reaction process they experience.

3.2.3. Mineral Saturation States and Partial Pressure of CO2 (g)

Figure 6 shows the mineral saturation index (SI) of different types of groundwater and surface water in the study area along the flow direction. In general, the evaporite minerals (gypsum, anhydrite, and halite) are all undersaturated (SI < 0), indicating that the water in Tibet can still continue dissolving evaporite minerals. The climate in Tibet is dry, with strong weathering capacity and evaporation. Therefore, it is possible to provide more sources of evaporite minerals to river water [51].
For carbonate minerals, the analysis results of all water samples show a supersaturated state (SI > 0) and show a certain regularity. Along the flow direction of the river water, different types of groundwater discharge into the river water. The SI value of river water relative to carbonate minerals gradually increased, indicating that its saturation state gradually increased, which also shows that the impact of groundwater discharge cannot be ignored.
Figure 7 shows the relationship between HCO3 + CO32− and the logarithm of partial pressure of CO2 (g) (log (pCO2 (g))) for the groundwater samples. Most surface water and groundwater fall within the soil (pCO2 (g) = 10−1.5 atm) and atmospheric (pCO2 (g) = 10−3.5 atm) CO2 (g) range. Surface water is very close to equilibrium with atmospheric CO2 (g), indicating the fastest recharging waters [52]. The hot spring samples are in equilibrium with values above the upper limit of CO2 (g) in the soil region (pCO2 (g) = 10−1.5 atm), indicating that they are influenced by volcanic emissions of CO2 (g). Their pCO2 (g) values ranged from 10−1.13 to 10−0.66 atm.
The CO2 in the open system will be continuously replenished, and the pCO2 of the groundwater in the study area is greater than the atmospheric value, indicating that the groundwater system is an open system. Additionally, the study area is located in the recharge–recharge runoff area, with strong water alternation, which is conducive to the migration of ions and the continuous replenishment of atmospheric CO2. pCO2 also controls the concentration of Ca2+ and HCO3 in water by controlling the solubility of carbonate minerals.

3.2.4. Source of Chemical Components

Because of water–rock interactions, silicate minerals such as plagioclase and albite are common sources of Na and K in groundwater [53]. The process of Na+ replacing Ca2+ and Mg2+ in the cation exchange process may also be a major factor in cation change [54]. The interrelationships between major elements can be used to identify the water–rock processes that groundwater encounters [55].
The Na+-Cl diagram that cold spring, mine water, and surface water are all located under the curve of y = x (Figure 8a), and halite dissolution can be excluded. The Cl in groundwater generally comes from weathering, the leaching of rocks or sediments, or volatilization from deep magma. The hot springs have offset the curve of y = x (Figure 8a), indicating volatiles from deep magma [55]. The extra Na+ in cold water is most likely the result of silicate weathering and/or cation exchange processes.
If the relationship between Ca2+ and SO42− in water is on the y = x curve (Figure 8b), it indicates that they are primarily derived from gypsum dissolution. However, the groundwater samples in the study area are basically below the curve of y = x, which is inconsistent with the dissolution mechanism of gypsum, and cation exchange seems to be a better explanation for this phenomenon. In addition, the data of groundwater samples fall on the three regions divided by the curves of y = x and y = 2x from the relationship of Ca2+ and HCO3 (Figure 8c). Hot springs and some cold springs are above y = 2x, indicating that HCO3 comes from somewhere other than dissolution, such as deep metamorphic degassing or magma volatilization. Those between the two curves are thought to be the result of calcite dissolution.
Furthermore, the main rock types involved in the water–rock reaction can be identified using the Mg2+/Na+ vs. Ca2+/Na+ and HCO3/Na vs. Ca2+/Na+ (Figure 8d,e) [56]. According to Figure 8d,e, mine water is located in the transition zone of silicate and carbonate dissolution, while hot springs and some cold springs are located in the transition zone of silicate and evaporite. In Figure 8f, most of the samples are located above the isocline, which indicates that the Mg2+ and Ca2+ in groundwater are mainly from aluminosilicate dissolution, and carbonate dissolution may also exist. Excess Na+ in groundwater could be caused by silicate weathering. According to Figure 8g, the geothermal water samples mostly fall in the upper part of the region, which is affected by carbonate decomposition, while some samples of surface water and cold spring water fall in the lower part, indicating that they are also affected by evaporative decomposition. In addition, diagram Na+ + K+ − Cl vs. Ca2+ + Mg2+ − (HCO3 + SO42−) can also indicate the degree of cation exchange Figure 8g.
The chemical composition of Zhaxikang water was analyzed by PHREEQC, and Figure 9A–C shows the projections of all water samples in the stability diagrams for Na-, K-, and Ca-feldspars. It can be used to determine the feldspar weathering stage and the corresponding secondary minerals in the study area.
It can be seen from Figure 9 that the water samples in the study area are in a state of dissolution with respect to feldspar. The weathering process of albite in groundwater and surface water has passed from the gibbsite stage to the kaolinite stage, and the hot spring water samples have reached the stage of Na-montmorillonite (Figure 9A). The K-feldspar has not yet reached saturation and is in a state of dissolution, and the weathering has reached the kaolinite stage (Figure 9B). Weathering of anorthite is also at the kaolinite stage (Figure 9C).
The groundwater in the study area is supersaturated relative to quartz and undersaturated relative to amorphous silica (Figure 9). Under the undersaturated state, the solubility of amorphous silica controls the upper limit of H4SiO4 activity in the water sample. In shallow water areas, CO2 gas participates in the hydrolysis of silicate minerals to form gibbsite, which in turn forms kaolinite, then up to the montmorillonite stage and produces amorphous silica. The weathering of feldspar and its secondary minerals requires the consumption of H2CO3. While H2CO3 is related to pCO2, the content of Ca2+ and HCO3 in water is proportional to Pco2, indicating that Ca2+ and HCO3 are formed by the dissolution of carbonate and feldspar minerals by CO2, which makes the two ions accumulate continuously.
In conclusion, the ions in the water samples are mainly controlled by cation exchange and silicate weathering, and the related ions in hot spring water samples are mainly from sodium montmorillonite. The non-aluminosilicate minerals that control the dissolution and balance of chemical composition in water are mainly calcite, and the aluminosilicate minerals that participate in the water–rock reaction are mainly feldspar, followed by kaolinite and montmorillonite.

3.2.5. Source of B in The Groundwater

Boron is the main component of magmatic rock or fine-grained soil and often forms borates in water environments [57], and Its mobility is closely related to the pH of the fluid [58]. B concentration in geothermal water in the study area is relatively high, ranging from 32 to 66 mg/L, with an average value of 43 mg/L. When they are discharged into a river, the concentration of B downstream can be as high as 9 mg/L, exceeding the 0.5 mg/L drinking water standard. According to the analysis results of PHREEQC, B mainly exists in the geothermal water in the form of H3BO3 and H2BO3. The geothermal water environment in the Zhaxikang area is characterized by high temperature, high TDS, and strong reducibility, which is conducive to boron enrichment (Supplementary Figure S1) [59].
Pearson correlation analysis was performed on 13 hydrochemical parameters in Zhaxikang to explain the influencing factors of groundwater quality, and the value of the phase relationship is shown in Table 3 [60]. B is well correlated with most ions, and Cl (r = 0.99) and HCO3 (r = 0.63) are well correlated, which may be due to the dissolution of carbonate minerals [59]. pH is negatively correlated with bicarbonate and B, indicating that an acidic environment enhances bicarbonate and B enrichment.
B can be parsed into groundwater by flushing clay minerals with freshwater [22]. However, cation exchange is also accompanied by the process of freshwater flushing, resulting in Ca depletion and Na enrichment in groundwater, which is consistent with the results of the study area. Furthermore, B may be derived from deep magma volatilization in the geothermal area [25]. A portion of B enters the river directly from hot spring discharge, while another portion is adsorbed on the surface of clay or muddy minerals in near-surface sediments. It is clear from the good relationship between B, Na, and Cl that B is controlled by adsorption/dissolution in groundwater.

3.3. Relationship between the Surface Water and Groundwater

The study area is dominated by inorganic pollutants from geothermal water. The mixing of cold groundwater and deep hot water is the main reason for the differences in the properties and composition of hot springs. Surface water and groundwater from a complete water cycle are controlled by the geological structure and other factors.
River water can reduce the concentration of harmful elements through physical precipitation and ion exchange. The changes of Cd, Pb, B, V, and Fe along the river direction are shown in Figure 10a. B and V will increase significantly in the downstream river as the flowing water passes through the discharge area of the hot spring, resulting in a change in the river water quality of the river. According to the calculation of the mixing ratio of surface water B, the mixing ratio of geothermal water to river water is about 20.3%. Although the river water has a certain self-purification ability, the ratio downstream is still 14.3%. The concentration of B (B = 9.46 mg/L) still exceeded the standard drinking water limit (Figure 10b).
The elements such as B and As carried by geothermal water are significantly higher than those carried by surface water. In most geothermal areas of Tibet, such as Yangbajing, Yangyi, Tasha, Quzhuomu, Semi, and Mansarovar, water from geothermal wells or springs is discharged into rivers nearby directly. As a result, geothermal water will certainly affect the quality of the surrounding surface water. Furthermore, as hot springs rise from the subsurface to the surface, many toxic and harmful elements accumulate in the shallow sediments, and they are released into nearby rivers and lakes through processes such as rainfall, weathering, and mineral alteration, contaminating surface water quality on a long-term scale.

3.4. Health Risk Assessment

Boron is an essential trace element in plants, but it can cause damage to the human liver and kidney in high concentrations [61]. Human immune function may also be compromised if excessive amounts of boron are consumed in drinking water over a long period [62]. Most of the rural areas in Tibet are built along rivers, and the polluted water will inevitably affect the drinking water of residents in Zhaxikang [1]. The quality of drinking water in southern Tibet is not properly disinfected, which poses a substantial risk to drinking water safety [63].
The water quality in the Zhaxikang catchment ranges from excellent to extremely poor, with the EWQI value being closely related to high TDS, sulfate, nitrate, and B values (see Supplementary Tables S2 and S3). The average regional EWQI was 315.9, with values ranging from 14.8 to 1373.6. (Figure 11). The spatial distribution of water quality in Zhaxikang is a large variation. It results in geothermal water with very poor water quality being discharged into the river. The majority of the samples (42.9%) were classified as excellent, primarily upstream river water and shallow groundwater, while extremely poor samples accounted for 38.1% (eight samples), primarily hot spring water, polluted river, and part of cold spring. Two samples (9.5%), one sample (4.8%), and one sample (4.8%) were classified as poor, medium, and good quality, respectively, mostly mixed water and shallow groundwater.
According to the results of EWQI (Figure 11), there are two end elements of water quality in the Zhaxikang area, drinkable water and polluted water, and the mixing of them is the primary cause of the deterioration of water quality in the Zhaxikang area.
Dilution is the main self-purification process of rivers. In addition, the concentration of B in the river downstream from the geothermal water discharge increased slightly from R01 to R02 and then decreased to R03 (Figure 10b). The increase in B concentration at R02 could be explained by the fact that the point is located near a fault, and there could be hidden hot spring seepage at the river’s bed. According to Figure 10b, B concentration is significantly reduced due to the action of self-purification of river water; however, its value is still far above the limits of drinking water.
After the geothermal water is discharged into the river, all three mixed samples (R01, R02, and R03) have a HI (Hazard Index) value higher than 1, indicating that the noncarcinogenic risk is unacceptable [24,64]. The noncarcinogenic risk factor for children ranged from 0.002 to 1.844, with an average of 0.375, whereas the HQ (Hazard Quotient) for adults ranged from 0.002 to 1.580 (see Supplementary Table S5). Generally, children are at a higher health risk than adults [18,65]. Based on the perspective of spatial analysis of geothermal water pollution, the HQIngestion and HQDermal for the downstream water samples (R01 and later) after geothermal water mixing are both higher than 1, indicating that geothermal water mixing is the primary source of drinking water safety hazards in the Zhaxikang area.

3.5. Contamination Model of Drinking Water in Zhaxikang Area

Geothermal water is a potentially serious source of contamination, and water quality will be significantly affected as geothermal water is discharged into the river directly. Water quality is affected by a variety of factors, including topography, bedrock properties, soil characteristics, and meteoric precipitation [66].
This study identified the main water pollution sources, mixing patterns, and health risks in the Zhaxikang area using hydrochemistry, isotope analysis, and risk assessment. Figure 12 shows the mixed contamination model of geothermal water to surface water and shallow groundwater in the Zhaxikang area. Rainwater and snow meltwater are the main recharge sources of cold groundwater in this area, while thermal groundwater also has a magmatic water source. The deep hot water will mix with the shallow groundwater along the faults during the upward migration process, changing the physical and chemical properties of the shallow groundwater [37]. Additionally, the existence of deposit veins will also affect the surrounding groundwater properties [33].
Given the drinking water situation in Tibet, the presence of pollution sources in river water may pose some health risks to nearby residents’ drinking water sources, including drinking water safety hazards and skin contact hazards. Geothermal water is drained into the river with excessive ion concentration (B, HCO3, As) in the form of hot springs and seepage on the riverbed. Relatively speaking, the seepage in the river is a relatively hidden pollution source, which is not easy to detect. On the other hand, mining would also release arsenic, which can also contaminate shallow groundwater and surface water [33]. Arsenic concentrations in lake water and mine water, for example, exceed limits in the study area. The residents of the area are at risk from the high boron water produced by the mixing of geothermal water when these shallow groundwater and river waters are used for drinking water.

3.6. Prospect of Water Security in Geothermal Development

The Qinghai–Tibet Plateau region has the most abundant and concentrated geothermal resources in China [67,68,69,70]. Geothermal exploitation and utilization in the study area will enter a new era under the carbon-neutral framework. The pollution of geothermal water to surface water will become a major issue of security for this region’s water resources.
The Zhaxikang thermal field has not been developed and utilized on a large scale, and the harmful elements mainly come from the discharge of hot springs and mining wastewater. However, as the development and utilization of geothermal energy accelerate, the discharge and treatment of geothermal wastewater must be addressed in the future. For example, there were no endemic diseases in the Yangbajing area during geothermal exploration, but the residents experienced an increase in dental fluorosis, baldness, and other diseases after geothermal wastewater was discharged into the Zangbuqu River without proper treatment [12]. The toxic and harmful elements carried by geothermal water cause the water quality of the river to exceed the drinking water standard, even after the self-purification effect of the river.
As a result, to protect Zhaxikang’s water environment, the problem of geothermal wastewater discharge and purification must be clarified. Reinjecting geothermal wastewater into the reservoir is an environmentally friendly and effective way to avoid surface pollution and maintain the geothermal fluid pressure that is necessary for geothermal recovery. On the other hand, it is necessary to strengthen monitoring of the geothermal mining industries and other industries’ pollutant output to the environment, as well as to impose appropriate restrictions and regulations. In addition, local governments and relevant organizations also need to educate the residents about public health, strengthen the protection of children, and enhance awareness and knowledge of drinking water treatment.

4. Conclusions

Based on the hydrochemical and isotopic characteristics of hot springs, cold springs, river water, cold groundwater, and lake water in the Zhaxikang catchment, this paper reveals the water–rock interaction mechanism controlling the water chemical characteristics and B pollution source, as well as evaluates the potential water environmental health risk of geothermal water discharge. The main conclusion could be as follows:
  • The geothermal water is characterized by high boron content in the study area, and the chemical type of hot spring water is HCO3-Cl-Na, which is related to the metamorphic CO2 degassing and mixing process in the geothermal area.
  • Hydrogen and oxygen isotope, Gibbs diagram, Pearson correlation, and ion correlation were used to explore the evolution mechanism of groundwater quality. Ice melting water, meteoric water, and magmatic water are all sources of geothermal water. The groundwater is an open system, and the major components are controlled by water–rock processes such as silicate mineral dissolution and cation exchange.
  • B mainly comes from deep volatiles, cation exchange, and leaching in shallow sediments.
  • According to the EWQI results, the water quality of the river gradually deteriorated after flowing through the geothermal area and mixed with geothermal water (with a mix ratio of ~20%). Additionally, there are significant non-carcinogenic health risks if drinking contaminated water, and children are at a higher health risk than adults.
In the context of energy conservation and emission reduction, the exploitation of the rich geothermal resources in the Qinghai–Tibet Plateau will be accelerated. Many toxic and harmful elements of geothermal water discharged into the river will be a threat to the regional water environment and health risks. Therefore, attention should be paid to environmental monitoring during geothermal exploitation and its impact on the water environment and the drinking water safety of residents.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14203243/s1, Figure S1: Relationship of the B concentration with other measured hydrochemical parameters; Table S1: Information entropy and entropy weight of parameters; Table S2: Assessment results according to the calculated EWQI; Table S3: Water quality class based on EWQI; Table S4: Model Parameters for the calculation of exposure dose; Table S5: Assessment results of health risks based on drinking water intake and dermal contact. References [71,72,73,74,75,76,77,78,79,80,81,82,83,84] are cited in the Supplementary Materials.

Author Contributions

Y.W. provided supervision in the presented work and involved in the writing and editing of the manuscript. L.L. (Liang Li) conducted the investigation and data collection, developed the suitable methodology, and was involved with the writing of the manuscript. H.G., L.L. (Lianghua Lu) and L.L. (Luping Li) helped in methodology selection and visualization, and participated in editing of the earlier versions of the manuscript. J.P. and F.C. helped in figure preparation and edited and wrote the early version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (42002300, 41907174), the China Geological Survey’s project: DD20211381, the National Key Research and Development Program of China (2018YFC0406506), and the Open Fund (PLC2020032) of the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Chengdu University of Technology).

Institutional Review Board Statement

No applicable.

Informed Consent Statement

No applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to the editors and the three anonymous reviewers. Their comments are quite constructive and helpful for us to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qu, B.; Zhang, Y.; Kang, S.; Sillanpää, M. Water quality in the Tibetan Plateau: Major ions and trace elements in rivers of the “Water Tower of Asia”. Sci. Total Environ. 2019, 649, 571–581. [Google Scholar] [CrossRef] [PubMed]
  2. Liao, Z. Thermal Springs and Geothermal Energy in the Qinghai-Tibetan Plateau and the Surroundings; Springer: Singapore, 2018; p. 311. [Google Scholar] [CrossRef]
  3. Tian, J.; Pang, Z.; Guo, Q.; Wang, Y.; Li, J.; Huang, T.; Kong, Y. Geochemistry of geothermal fluids with implications on the sources of water and heat recharge to the Rekeng high-temperature geothermal system in the Eastern Himalayan Syntax. Geothermics 2018, 74, 92–105. [Google Scholar] [CrossRef]
  4. Tian, J.; Pang, Z.; Wang, Y.; Guo, Q. Fluid geochemistry of the Cuopu high temperature geothermal system in the eastern Himalayan syntaxis with implication on its genesis. Appl. Geochem. 2019, 110, 104422. [Google Scholar] [CrossRef]
  5. Wang, Y.; Pang, Z.; Hao, Y.; Fan, Y.; Tian, J.; Li, J. A revised method for heat flux measurement with applications to the fracture-controlled Kangding geothermal system in the Eastern Himalayan Syntaxis. Geothermics 2019, 77, 188–203. [Google Scholar] [CrossRef]
  6. Wang, Y.; Pang, Z. Heat flux measurements and thermal potential of the Garze geothermal area in the eastern Himalayan Syntaxis. J. Volcanol. Geotherm. Res. 2021, 420, 107407. [Google Scholar] [CrossRef]
  7. Wang, Y.; Pang, Z. Heat flux in volcanic and geothermal areas: Methods, principles, applications and future directions. Gondwana Res. 2022. [Google Scholar] [CrossRef]
  8. Zhao, Y.; Han, J.; Zhang, B.; Gong, J. Impact of transferred water on the hydrochemistry and water quality of surface water and groundwater in Baiyangdian Lake, North China. Geosci. Front. 2021, 12, 101086. [Google Scholar] [CrossRef]
  9. Sharifi, R.; Moore, F.; Keshavarzi, B. Mobility and chemical fate of arsenic and antimony in water and sediments of Sarouq River catchment, Takab geothermal field, northwest Iran. J. Environ. Manag. 2016, 170, 136–144. [Google Scholar] [CrossRef] [PubMed]
  10. Guo, Q.; Planer-Friedrich, B.; Liu, M.; Yan, K.; Wu, G. Magmatic fluid input explaining the geochemical anomaly of very high arsenic in some southern Tibetan geothermal waters. Chem. Geol. 2019, 513, 32–43. [Google Scholar] [CrossRef]
  11. Villalba, E.; Tanjal, C.; Borzi, G.; Páez, G.; Carol, E. Geogenic arsenic contamination of wet-meadows associated with a geothermal system in an arid region and its relevance for drinking water. Sci. Total Environ. 2020, 720, 137571. [Google Scholar] [CrossRef]
  12. Guo, Q.; Wang, Y.; Liu, W. B, As, and F contamination of river water due to wastewater discharge of the Yangbajing geothermal power plant, Tibet, China. Environ. Earth Sci. 2008, 56, 197–205. [Google Scholar] [CrossRef]
  13. MacDonald, A.; Bonsor, H.C.; Ahmed, K.M.; Burgess, W.G.; Basharat, M.; Calow, R.C.; Dixit, A.; Foster, S.S.D.; Gopal, K.; Lapworth, D.J.; et al. Groundwater quality and depletion in the Indo-Gangetic Basin mapped from in situ observations. Nat. Geosci. 2016, 9, 762–766. [Google Scholar] [CrossRef] [Green Version]
  14. Wang, Y.; Li, L.; Wen, H.; Hao, Y. Geochemical evidence for the nonexistence of supercritical geothermal fluids at the Yangbajing geothermal field, southern Tibet. J. Hydrol. 2022, 604, 127243. [Google Scholar] [CrossRef]
  15. Hao, Y.; Pang, Z.; Tian, J.; Wang, Y.; Li, Z.; Li, L.; Xing, L. Origin and evolution of hydrogen-rich gas discharges from a hot spring in the eastern coastal area of China. Chem. Geol. 2020, 538, 119477. [Google Scholar] [CrossRef]
  16. Tian, Y.; Yu, C.; Zha, X.; Wu, J.; Gao, X.; Feng, C.; Luo, K. Distribution and Potential Health Risks of Arsenic, Selenium, and Fluorine in Natural Waters in Tibet, China. Water 2016, 8, 568. [Google Scholar] [CrossRef] [Green Version]
  17. Li, C.; Kang, S.; Chen, P.; Zhang, Q.; Mi, J.; Gao, S.; Sillanpää, M. Geothermal spring causes arsenic contamination in river waters of the southern Tibetan Plateau, China. Environ. Earth Sci. 2014, 71, 4143–4148. [Google Scholar] [CrossRef]
  18. Wang, X.; Dan, Z.; Cui, X.; Zhang, R.; Zhou, S.; Wenga, T.; Yan, B.; Chen, G.; Zhang, Q.; Zhong, L. Contamination, ecological and health risks of trace elements in soil of landfill and geothermal sites in Tibet. Sci. Total Environ. 2020, 715, 136639. [Google Scholar] [CrossRef] [PubMed]
  19. Liao, Z.; Zhao, P. Geotropics of Yunnan and Tibet: Geothermal Resources and Typical Geothermal Systems; Chinese Academy of Sciences: Beijing, China, 1999. (In Chinese) [Google Scholar]
  20. Lü, Y.; Zheng, M.; Zhao, P.; Xu, R. Geochemical processes and origin of boron isotopes in geothermal water in the Yunnan-Tibet geothermal zone. Sci. China Earth Sci. 2014, 57, 2934–2944. (In Chinese) [Google Scholar] [CrossRef]
  21. Voutsa, D.; Dotsika, E.; Kouras, A.; Poutoukis, D.; Kouimtzis, T. Study on distribution and origin of boron in groundwater in the area of Chalkidiki, Northern Greece by employing chemical and isotopic tracers. J. Hazard. Mater. 2009, 172, 1264–1272. [Google Scholar] [CrossRef] [PubMed]
  22. Ravenscroft, P.; McArthur, J. Mechanism of regional enrichment of groundwater by boron: The examples of Bangladesh and Michigan, USA. Appl. Geochem. 2004, 19, 1413–1430. [Google Scholar] [CrossRef]
  23. Pennisi, M.; Bianchini, G.; Muti, A.; Kloppmann, W.; Gonfiantini, R. Behaviour of boron and strontium isotopes in groundwater–aquifer interactions in the Cornia Plain (Tuscany, Italy). Appl. Geochem. 2006, 21, 1169–1183. [Google Scholar] [CrossRef]
  24. Zango, M.S.; Sunkari, E.D.; Abu, M.; Lermi, A. Hydrogeochemical controls and human health risk assessment of groundwater fluoride and boron in the semi-arid North East region of Ghana. J. Geochem. Explor. 2019, 207, 106363. [Google Scholar] [CrossRef]
  25. Liu, M.; Guo, Q.; Luo, L.; He, T. Environmental impacts of geothermal waters with extremely high boron concentrations: Insight from a case study in Tibet, China. J. Volcanol. Geotherm. Res. 2020, 397, 106887. [Google Scholar] [CrossRef]
  26. Nigro, A.; Sappa, G.; Barbieri, M. Boron isotopes in groundwater: Evidence from contamination and interaction with terrigenous–evaporitic sequence, east-central Italy. Geochem. Explor. Environ. Anal. 2018, 18, 343–350. [Google Scholar] [CrossRef]
  27. Vengosh, A. The isotopic composition of anthropogenic boron and its potential impact on the environment. Biol. Trace Element Res. 1998, 66, 145–151. [Google Scholar] [CrossRef] [PubMed]
  28. WHO. Guidelines for Drinking-Water Quality: Incorporating First Addendum, World Health Organization. 2017. Available online: https://www.who.int/publications-detail-redirect/9789241549950 (accessed on 8 July 2022).
  29. Tong, W.; Liao, Z.; Liu, S. Thermal Springs in Tibet; Science Press: Beijing, China, 2000. [Google Scholar]
  30. Duan, J.; Tang, J.; Lin, B. Zinc and lead isotope signatures of the Zhaxikang PbZn deposit, South Tibet: Implications for the source of the ore-forming metals. Ore Geol. Rev. 2016, 78, 58–68. [Google Scholar] [CrossRef]
  31. Qing, C.; Ding, J.; Li, G.; Zhou, Q.; Lai, Y. Element Distribution of the Surrounding Rock of Zhaxikang Deposit in Longzi Country, Tibet. Metal Mine 2014, 43, 97–100. (In Chinese) [Google Scholar]
  32. Wang, D.; Sun, X.; Zheng, Y.; Wu, S.; Xia, S.; Chang, H.; Yu, M. Two pulses of mineralization and genesis of the Zhaxikang Sb–Pb–Zn–Ag deposit in southern Tibet: Constraints from Fe–Zn isotopes. Ore Geol. Rev. 2017, 84, 347–363. [Google Scholar] [CrossRef]
  33. Gu, H.; Ni, H.; Wang, Y.; Liu, Y.; Li, D.; Zhang, Z. Hydrogeochemical Characteristics and Impact of Arsenic Released from a Gold Deposit in Tibet. Mine Water Environ. 2020, 39, 746–757. [Google Scholar] [CrossRef]
  34. Wang, Y.; Gu, H.; Li, D.; Lyu, M.; Lu, L.; Zuo, Y.; Song, R. Hydrochemical characteristics and genesis analysis of geothermal fluid in the Zhaxikang geothermal field in Cuona County, southern Tibet. Environ. Earth Sci. 2021, 80, 415. [Google Scholar] [CrossRef]
  35. Bhutiani, R.; Kulkarni, D.B.; Khanna, D.R.; Gautam, A. Water Quality, Pollution Source Apportionment and Health Risk Assessment of Heavy Metals in Groundwater of an Industrial Area in North India. Expo. Health 2016, 8, 3–18. [Google Scholar] [CrossRef]
  36. Wongsasuluk, P.; Chotpantarat, S.; Siriwong, W.; Robson, M. Heavy metal contamination and human health risk assessment in drinking water from shallow groundwater wells in an agricultural area in Ubon Ratchathani province, Thailand. Environ. Geochem. Health 2014, 36, 169–182. [Google Scholar] [CrossRef] [PubMed]
  37. Guo, Q.; Liu, M.; Li, J.; Zhang, X.; Guo, W.; Wang, Y. Fluid geochemical constraints on the heat source and reservoir temperature of the Banglazhang hydrothermal system, Yunnan-Tibet Geothermal Province, China. J. Geochem. Explor. 2017, 172, 109–119. [Google Scholar] [CrossRef]
  38. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology; Wiley: New York, NY, USA, 1998. [Google Scholar]
  39. Adimalla, N. Groundwater quality delineation based on fuzzy comprehensive assessment method (FCAM): A case study. Arab. J. Geosci. 2020, 13, 1256. [Google Scholar] [CrossRef]
  40. Chen, F.; Yao, L.; Mei, G.; Shang, Y.; Xiong, F.; Ding, Z. Groundwater Quality and Potential Human Health Risk Assessment for Drinking and Irrigation Purposes: A Case Study in the Semiarid Region of North China. Water 2021, 13, 783. [Google Scholar] [CrossRef]
  41. USEPA. User’s Guide: Human Health Risk Assessment; United States Environmental Protection Agency: Washington, DC, USA, 2008.
  42. Parkhurst, D.L.; Appelo, C.A.J. Description of Input and Examples for PHREEQC Version 3—A Computer Program for Speciation, Batch-reaction, One-dimensional Transport, and Inverse Geochemical Calculations. In Techniques and Methods; U.S. Geological Survey: Reston, VA, USA, 2013. [Google Scholar]
  43. Guo, Q.; Wang, Y.; Liu, W. O, H, and Sr isotope evidences of mixing processes in two geothermal fluid reservoirs at Yangbajing, Tibet, China. Environ. Earth Sci. 2010, 59, 1589–1597. [Google Scholar] [CrossRef]
  44. Giggenbach, W. Isotopic shifts in waters from geothermal and volcanic systems along convergent plate boundaries and their origin. Earth Planet. Sci. Lett. 1992, 113, 495–510. [Google Scholar] [CrossRef]
  45. Hudson-Edwards, K.A.; Archer, J. Geochemistry of As-, F- and B-bearing waters in and around San Antonio de los Cobres, Argentina, and implications for drinking and irrigation water quality. J. Geochem. Explor. 2012, 112, 276–284. [Google Scholar] [CrossRef] [Green Version]
  46. Lu, L.; Pang, Z.; Kong, Y.; Guo, Q.; Wang, Y.; Xu, C.; Gu, W.; Zhou, L.; Yu, D. Geochemical and isotopic evidence on the recharge and circulation of geothermal water in the Tangshan Geothermal System near Nanjing, China: Implications for sustainable development. Appl. Hydrogeol. 2018, 26, 1705–1719. [Google Scholar] [CrossRef]
  47. Craig, H. Isotopic Variations in Meteoric Waters. Science 1961, 133, 1833–1834. [Google Scholar] [CrossRef]
  48. Wen, R.; Tian, L.; Weng, Y.; Liu, Z.; Zhao, Z. The altitude effect of δ 18O in precipitation and river water in the Southern Himalayas. Chin. Sci. Bull. 2012, 57, 1693–1698. [Google Scholar] [CrossRef] [Green Version]
  49. Shah, B.; Kansara, B.; Shankar, J.; Soni, M.; Bhimjiyani, P.; Bhanushali, T.; Shah, M.; Sircar, A. Reckoning of water quality for irrigation and drinking purposes in the konkan geothermal provinces, Maharashtra, India. Groundw. Sustain. Dev. 2019, 9, 100247. [Google Scholar] [CrossRef]
  50. Gibbs, R.J. Mechanisms Controlling World Water Chemistry. Science 1970, 170, 1088–1090. [Google Scholar] [CrossRef]
  51. Gou, Q.Q.; Han, Z.W.; Wang, G.H. Research Progress on Soil Salinization in Arid Irrigated Area of Northwestern China. Chin. Agric. Sci. Bull. 2011, 29, 246–250. (In Chinese) [Google Scholar]
  52. Bennett, G.; Van Reybrouck, J.; Shemsanga, C.; Kisaka, M.; Tomašek, I.; Fontijn, K.; Kervyn, M.; Walraevens, K. Identification of low fluoride areas using conceptual groundwater flow model and hydrogeochemical system analysis in the aquifer system on the flanks of an active volcano: Mount Meru, Northern Tanzania. Sci. Total Environ. 2021, 814, 152682. [Google Scholar] [CrossRef] [PubMed]
  53. Bretzler, A.; Osenbrück, K.; Gloaguen, R.; Ruprecht, J.S.; Kebede, S.; Stadler, S. Groundwater origin and flow dynamics in active rift systems—A multi-isotope approach in the Main Ethiopian Rift. J. Hydrol. 2011, 402, 274–289. [Google Scholar] [CrossRef]
  54. Farooqi, A.; Masuda, H.; Firdous, N. Toxic fluoride and arsenic contaminated groundwater in the Lahore and Kasur districts, Punjab, Pakistan and possible contaminant sources. Environ. Pollut. 2007, 145, 839–849. [Google Scholar] [CrossRef] [PubMed]
  55. Zhang, Y.; Dai, Y.; Wang, Y.; Huang, X.; Xiao, Y.; Pei, Q. Hydrochemistry, quality and potential health risk appraisal of nitrate enriched groundwater in the Nanchong area, southwestern China. Sci. Total Environ. 2021, 784, 147186. [Google Scholar] [CrossRef] [PubMed]
  56. Gaillardet, J.; Dupré, B.; Louvat, P.; Allègre, C. Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers. Chem. Geol. 1999, 159, 3–30. [Google Scholar] [CrossRef]
  57. Palmucci, W.; Rusi, S. Boron-rich groundwater in Central Eastern Italy: A hydrogeochemical and statistical approach to define origin and distribution. Environ. Earth Sci. 2014, 72, 5139–5157. [Google Scholar] [CrossRef]
  58. Gemici, U.; Tarcan, G. Distribution of boron in thermal waters of western Anatolia, Turkey, and examples of their environmental impacts. Environ. Earth Sci. 2002, 43, 87–98. [Google Scholar] [CrossRef]
  59. Liu, M.; Guo, Q.; Wu, G.; Guo, W.; She, W.; Yan, W. Boron geochemistry of the geothermal waters from two typical hydrothermal systems in Southern Tibet (China): Daggyai and Quzhuomu. Geothermics 2019, 82, 190–202. [Google Scholar] [CrossRef]
  60. Das, S.; Nag, S.K. Application of multivariate statistical analysis concepts for assessment of hydrogeochemistry of groundwater—A study in Suri I and II blocks of Birbhum District, West Bengal, India. Appl. Water Sci. 2015, 7, 873–888. [Google Scholar] [CrossRef] [Green Version]
  61. Karro, E.; Uppin, M. The occurrence and hydrochemistry of fluoride and boron in carbonate aquifer system, central and western Estonia. Environ. Monit. Assess. 2013, 185, 3735–3748. [Google Scholar] [CrossRef]
  62. Murray, F. A Human Health Risk Assessment of Boron (Boric Acid and Borax) in Drinking Water. Regul. Toxicol. Pharmacol. 1995, 22, 221–230. [Google Scholar] [CrossRef] [PubMed]
  63. Zhang, X.; Li, X.; Budu; Liu, Y. Analysis on sanitary monitoring of drinking water quality in Linzhi region of Tibet in 2011. Tibet. J. Med. 2013, 1, 62–64. (In Chinese) [Google Scholar]
  64. Adimalla, N.; Qian, H. Groundwater quality evaluation using water quality index (WQI) for drinking purposes and human health risk (HHR) assessment in an agricultural region of Nanganur, south India. Ecotoxicol. Environ. Saf. 2019, 176, 153–161. [Google Scholar] [CrossRef] [PubMed]
  65. Durowoju, O.S.; Ekosse, G.-I.E.; Odiyo, J.O. Occurrence and Health-Risk Assessment of Trace Metals in Geothermal Springs within Soutpansberg, Limpopo Province, South Africa. Int. J. Environ. Res. Public Health 2020, 17, 4438. [Google Scholar] [CrossRef] [PubMed]
  66. Ravikumar, P.; Mehmood, M.A.; Somashekar, R.K. Water quality index to determine the surface water quality of Sankey tank and Mallathahalli lake, Bangalore urban district, Karnataka, India. Appl. Water Sci. 2013, 3, 247–261. [Google Scholar] [CrossRef] [Green Version]
  67. Guo, Q.; Pang, Z.; Wang, Y.; Tian, J. Fluid geochemistry and geothermometry applications of the Kangding high-temperature geothermal system in eastern Himalayas. Appl. Geochem. 2017, 81, 63–75. [Google Scholar] [CrossRef]
  68. Luo, J.; Pang, Z.; Kong, Y.; Wang, Y. Geothermal potential evaluation and development prioritization based on geochemistry of geothermal waters from Kangding area, western Sichuan, China. Environ. Earth Sci. 2017, 76, 343. [Google Scholar] [CrossRef]
  69. Dong, B.; Dong, S.; Wang, Y.; Wen, F.; Yu, C.; Zhou, J.; Song, R. Detecting Geothermal Resources in a Plateau Area: Constraints From Land Surface Temperature Characteristics Using Landsat 8 Data. Front. Earth Sci. 2022, 10, 785900. [Google Scholar] [CrossRef]
  70. Hu, J.; Tian, Y.; Li, Z.; Jiang, G.; Zuo, Y.; Zhang, C.; Wang, Y.; Wang, Y.; Hu, S. Geothermo Tool: An open-source software for basic geothermal calculations and plots. Geothermics 2022, 106, 102551. [Google Scholar] [CrossRef]
  71. Li, P.; Qian, H.; Wu, J. Groundwater Quality Assessment Based on Improved Water Quality Index in Pengyang County, Ningxia, Northwest China. E-J. Chem. 2010, 7, 451304. [Google Scholar] [CrossRef] [Green Version]
  72. Ugran, V. Groundwater fluoride contamination and its possible health implications in Indi taluk of Vijayapura District (Karnataka State), India. Environ. Geochem. Health 2017, 39, 1017–1029. [Google Scholar] [CrossRef] [PubMed]
  73. Sathe, S.; Mahanta, C.; Subbiah, S. Hydrogeochemical Evaluation of Intermittent Alluvial Aquifers Controlling Arsenic and Fluoride Contamination and Corresponding Health Risk Assessment. Expo. Health 2021, 13, 661–680. [Google Scholar] [CrossRef]
  74. Zhang, Y.; Wu, J.; Xu, B. Human health risk assessment of groundwater nitrogen pollution in Jinghui canal irrigation area of the loess region, northwest China. Environ. Earth Sci. 2018, 77, 273. [Google Scholar] [CrossRef]
  75. USEPA. Human Health Evaluation Manual, Supplemental Guidance: Update of Standard Default Exposure Factors; United States Environmental Protection Agency: Washington, DC, USA, 2014.
  76. USEPA. Integrated Risk Information System; United States Environmental Protection Agency: Washington, DC, USA, 2012.
  77. Li, P.; Tian, R.; Liu, R. Solute Geochemistry and Multivariate Analysis of Water Quality in the Guohua Phosphorite Mine, Guizhou Province, China. Expo. Health 2019, 11, 81–94. [Google Scholar] [CrossRef]
  78. Wu, J.; Sun, Z. Evaluation of Shallow Groundwater Contamination and Associated Human Health Risk in an Alluvial Plain Impacted by Agricultural and Industrial Activities, Mid-west China. Expo. Health 2016, 8, 311–329. [Google Scholar] [CrossRef]
  79. Edet, A.; Worden, R. Monitoring of the physical parameters and evaluation of the chemical composition of river and groundwater in Calabar (Southeastern Nigeria). Environ. Monit. Assess. 2009, 157, 243–258. [Google Scholar] [CrossRef] [PubMed]
  80. Davis, S.; Whittemore, D.; Fabryka-Martin, J. Uses of Chloride/Bromide Ratios in Studies of Potable Water. Ground Water 1998, 36, 338–350. [Google Scholar] [CrossRef]
  81. Rautio, A.; Korkka-Niemi, K. Chemical and isotopic tracers indicating groundwater/surface-water interaction within a boreal lake catchment in Finland. Hydrogeol. J. 2015, 23, 687–705. [Google Scholar] [CrossRef]
  82. Dun, Y.; Tang, C.; Shen, Y. Identifying interactions between river water and groundwater in the North China Plain using multiple tracers. Environ. Earth Sci. 2013, 72, 99–110. [Google Scholar] [CrossRef]
  83. Lamontagne, S. River infiltration to a subtropical alluvial aquifer inferred using multiple environmental tracers. Water Resour. Res. 2015, 51, 4532–4549. [Google Scholar] [CrossRef] [Green Version]
  84. Crandall, C.; Katz, B.; Hirten, J. Hydrochemical evidence for mixing of river water and groundwater during high-flow conditions, lower Suwannee River basin, Florida, USA. Hydrogeol. J. 1999, 7, 454–467. [Google Scholar] [CrossRef]
Figure 1. (a) Geographical distribution map of the study area; (b) tectonic and topography and (c) geology and sampling points distribution map of the study area. TB: Tarim Basin, SGT: Songpan–Ganzi Terrain, QT: Qiangtang Terrain, LT: Lasha Terrain, HM: Himalayan Mountains.
Figure 1. (a) Geographical distribution map of the study area; (b) tectonic and topography and (c) geology and sampling points distribution map of the study area. TB: Tarim Basin, SGT: Songpan–Ganzi Terrain, QT: Qiangtang Terrain, LT: Lasha Terrain, HM: Himalayan Mountains.
Water 14 03243 g001
Figure 2. Box plots of chemical parameters for different water types in the study area.
Figure 2. Box plots of chemical parameters for different water types in the study area.
Water 14 03243 g002
Figure 3. Relationship of δD vs. δ18O in surface water, mine water, cold spring, and hot spring [43,44].
Figure 3. Relationship of δD vs. δ18O in surface water, mine water, cold spring, and hot spring [43,44].
Water 14 03243 g003
Figure 4. Piper diagram of river water, lake water, spring water and mine water.
Figure 4. Piper diagram of river water, lake water, spring water and mine water.
Water 14 03243 g004
Figure 5. Gibbs diagrams of groundwater samples. (a) ρNa+/ρ(Na+ + Ca2+); (b) ρCl/ρ(Cl + HCO3).
Figure 5. Gibbs diagrams of groundwater samples. (a) ρNa+/ρ(Na+ + Ca2+); (b) ρCl/ρ(Cl + HCO3).
Water 14 03243 g005
Figure 6. Saturation index (SI) of potentially relevant minerals for the water samples in the study area.
Figure 6. Saturation index (SI) of potentially relevant minerals for the water samples in the study area.
Water 14 03243 g006
Figure 7. Relationship between HCO3 + CO32− and log (pCO2 (g)) for the water samples in the study area.
Figure 7. Relationship between HCO3 + CO32− and log (pCO2 (g)) for the water samples in the study area.
Water 14 03243 g007
Figure 8. Relationship diagrams of (a) Cl vs. Na+, (b) SO42− vs. Ca2+, (c) HCO3 vs. Ca2+, (d) (Mg2+/Na+) vs. (Ca2+/Na+), (e) (HCO3/Na+) vs. (Ca2+/Na+), (f) HCO3 + SO42− vs. Ca2+ + Mg2+, (g) Na+ + K+ − Cl vs. Ca2+ + Mg2+ − (HCO3 + SO42−).
Figure 8. Relationship diagrams of (a) Cl vs. Na+, (b) SO42− vs. Ca2+, (c) HCO3 vs. Ca2+, (d) (Mg2+/Na+) vs. (Ca2+/Na+), (e) (HCO3/Na+) vs. (Ca2+/Na+), (f) HCO3 + SO42− vs. Ca2+ + Mg2+, (g) Na+ + K+ − Cl vs. Ca2+ + Mg2+ − (HCO3 + SO42−).
Water 14 03243 g008
Figure 9. (AC): Stability diagrams of Na-, K-, and Ca-feldspars and their weathering products for water samples from the study area. The dashed lines represent saturation with respect to quartz and amorphous silica.
Figure 9. (AC): Stability diagrams of Na-, K-, and Ca-feldspars and their weathering products for water samples from the study area. The dashed lines represent saturation with respect to quartz and amorphous silica.
Water 14 03243 g009
Figure 10. (a) The covariant relationship of Cd, Pb, B, V and Fe along the river direction; (b) Variations of B concentrations of Langmai River waters along the flow direction.
Figure 10. (a) The covariant relationship of Cd, Pb, B, V and Fe along the river direction; (b) Variations of B concentrations of Langmai River waters along the flow direction.
Water 14 03243 g010
Figure 11. Correlation diagram between EWQI and TDS.
Figure 11. Correlation diagram between EWQI and TDS.
Water 14 03243 g011
Figure 12. Schematic contamination model of geothermal water to surface water and shallow groundwater in Zhaxikang Area.
Figure 12. Schematic contamination model of geothermal water to surface water and shallow groundwater in Zhaxikang Area.
Water 14 03243 g012
Table 1. Physical parameters and hydrogen and oxygen isotopes of water samples from the Zhaxikang groundwater system.
Table 1. Physical parameters and hydrogen and oxygen isotopes of water samples from the Zhaxikang groundwater system.
Sample IDWater TypepHTemp [°C]EC [μS/cm]ORP [mV]TDS [mg/L]δ2H [‰]δ18O [‰]
R01rw8.2219.6555921.5404.5−131−17.2
R02rw8.7917.84604.520455.5−133−17.1
R03rw8.7316.9461331.8470.9−141−18.2
R04rw8.25.5534327.8354.8−129−17.8
R05rw8.196.33302.935305.9−128−17.6
R06rw8.6815.08232.425.6186.3−128−17.4
R07rw8.817.17284.612.6217.6−129−17.2
L01lw8.9418.34281.123.2209.3−128−17.4
L02lw9.2414.898944.723.27206−76.7−8.26
L03lw8.7410.662630.725.42356−67−6.59
L04lw9.6516.1422.629.8331.2−99.3−11.4
M01mw8.4819.2793.240.3580.1−157−20.1
M02mw8.0515.65935.6−41740.4−154−21.1
M03mw7.7421.4718.438.7554−156−20.3
S01hsw6.8250.23523.1−10.61547−144−17.3
S02csw8.326.53240.523.3241.5−147−19
S03csw6.758.172035.415.61950−144−17
S04csw8.618.35115214.5857.7−109−11.5
S05hsw7.0349.972372.1−19.61044−147−19.5
S06hsw7.0143.334138.7−17.71992−142−16.3
S07hsw7.3569.514698−262.51653−142−16.6
Snow-melt water a -----−175.0−24.5
Magmatic water b -----−20 ± 1010 ± 2
Note: rw: river water, lw: lake water, mw: mine well water, csw: cold spring water, hsw: hot spring water; “-” indicates no value; a, [43]; b, [44].
Table 2. Statistical results of chemical parameters for different water types in the study area (unit: mg/L).
Table 2. Statistical results of chemical parameters for different water types in the study area (unit: mg/L).
ParametersWHO StandardsHot Spring Water (n = 4)Cold Spring Water (n = 3)Lake Water (n = 4)Mine Water (n = 3)River Water (n = 7)
MinMaxMeanMinMaxMeanMinMaxMeanMinMaxMeanMinMaxMean
K+≤20011.6835.4026.310.5426.1012.750.5454.2317.650.352.251.450.427.803.02
Na+≤200193.00695.00486.258.43584.00288.145.731028.00318.158.7528.3020.254.5869.5027.87
Ca2+≤20044.20151.90107.3058.61150.9090.667.5176.3135.7683.04149.30106.3943.8992.5861.60
Mg2+≤1507.4328.0717.7718.8147.8532.9113.082123.00671.7159.7184.4670.1912.7843.6319.54
Cl≤250128.00621.00433.557.60530.80258.830.11103.5031.351.6410.055.070.0069.8322.86
SO42−≤250175.30235.90201.0066.74163.20104.0969.618524.002654.55161.40427.20255.1372.15201.40128.77
HCO3≤120622.621401.001004.00186.791556.54840.5362.261867.85786.05311.31466.96394.3293.39280.18191.23
NO3≤200.005.821.550.1923.989.930.000.700.230.001.040.450.001.720.49
B<0.58.5566.3642.360.3557.0230.020.020.240.080.070.990.660.019.463.52
SiO2-54.54105.5072.237.9624.8213.961.515.883.1213.0925.6419.874.4720.4412.22
V (μg/L)-1.554.043.230.624.453.350.952.521.710.280.430.350.381.460.93
As (μg/L)≤500.036.962.860.035.921.990.03126.6834.560.034.502.980.032.220.34
Cd (μg/L)<50.360.460.410.450.390.420.420.540.460.480.570.530.350.580.42
Pb (μg/L)<100.640.700.660.620.650.630.590.690.710.650.680.660.640.660.65
Fe (μg/L)<30028.60108.8948.3725.0465.8847.523.85163.0957.0643.56103.8266.0331.6672.3151.15
pH6.5~8.56.827.357.05 *6.758.607.89 *8.749.659.14 *7.748.488.09 *8.198.808.52 *
TDS≤10001044.001992.001559.00241.50 1950.001016.40209.307206.002525.63554.00740.40624.83186.30470.90342.21
Note: * The pH is just a statistical mean value.
Table 3. Correlation coefficients of components for the collected sample.
Table 3. Correlation coefficients of components for the collected sample.
pHTDSK+Na+Ca2+Mg2+ClSO42−HCO3NO3AsCdB
pH1.00
TDS0.011.00
K+−0.330.881.00
Na+−0.350.860.991.00
Ca2+−0.71−0.180.040.101.00
Mg2+0.320.930.650.63−0.391.00
Cl−0.730.290.700.730.49−0.071.00
SO42−0.310.930.660.64−0.370.99−0.051.00
HCO3−0.420.840.930.950.160.590.700.591.00
NO3−0.460.090.220.280.44−0.090.47−0.100.401.00
As0.300.920.670.65−0.380.99−0.030.990.58−0.111.00
Cd0.160.270.060.030.040.35−0.250.340.140.040.341.00
B−0.720.180.610.640.51−0.180.99−0.170.630.48 *−0.14−0.271.00
Note: * represents the statistical mean value.
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, L.; Wang, Y.; Gu, H.; Lu, L.; Li, L.; Pang, J.; Chen, F. The Genesis Mechanism and Health Risk Assessment of High Boron Water in the Zhaxikang Geothermal Area, South Tibet. Water 2022, 14, 3243. https://doi.org/10.3390/w14203243

AMA Style

Li L, Wang Y, Gu H, Lu L, Li L, Pang J, Chen F. The Genesis Mechanism and Health Risk Assessment of High Boron Water in the Zhaxikang Geothermal Area, South Tibet. Water. 2022; 14(20):3243. https://doi.org/10.3390/w14203243

Chicago/Turabian Style

Li, Liang, Yingchun Wang, Hongyu Gu, Lianghua Lu, Luping Li, Jumei Pang, and Feifei Chen. 2022. "The Genesis Mechanism and Health Risk Assessment of High Boron Water in the Zhaxikang Geothermal Area, South Tibet" Water 14, no. 20: 3243. https://doi.org/10.3390/w14203243

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