Hydrochemical Characteristics and Quality Evaluation of Irrigation and Drinking Water in Bangong Co Lake Watershed in Northwest Tibetan Plateau

: Bangong Lake is a narrow and long lake in the arid region of the plateau in northern Tibet. The salinity of the east of the lake is different from that the west, resulting in differences in the natural environment and human living conditions on each side. Watershed hydrochemical analysis and spatial statistical analysis can help to understand regional hydrochemical evolution and water quality evaluation. In this study, the hydrochemical characteristics of surface water (glacier, river, and lake) and groundwater in the Bangong Co Lake Watershed were investigated to reveal the relationships between various water bodies. The drinking water quality index (DWQI) and USSL classiﬁcation were applied to assess groundwater quality suitability for agricultural and drinking purposes. The hydrochemical characteristics show the differences among water bodies and their spatial distribution. The analyzed groundwater and surface water samples, such as river water and glaciers, were mainly Ca-HCO 3 -type and the lake water was mainly categorized as Na-Cl-type with some Na-HCO 3 -Cl type. The lake water’s chemical components are mainly affected by evaporative karst decomposition. The main mineralization process of groundwater and river water was related to the dissolution of reservoir minerals, such as dolomite and calcite, as well as halite. The drinking water quality index (DWQI) indicates that 79% of the groundwater samples in the study area were of good enough quality for drinking. In terms of irrigation water quality, the electrical conductivity (EC), calculated sodium adsorption ratio (SAR), and magnesium hazardous ratio (MHR) showed that more than 13% of the total samples were not suitable for irrigation. However, the USSL classiﬁcation indicated that glacier and river water are relatively suitable for irrigation. Additionally, some groundwater and lake water has very high alkalinity or salinity, which is alarming when considering them for irrigation.


Introduction
The Tibetan Plateau, known as the roof of the world and the water tower of Asia, is the birthplace of the Yellow River, the Yangtze River, the Ganges River, the Mekong River, the Indus River, and other rivers, and has the largest, highest, and most densely distributed alpine lake cluster in the world [1,2]. In the Tibetan Plateau, lakes play an important role in regulating the ecosystem and hydrochemistry and are a vital indicator of water-rock interaction and environmental change in the basins [3,4]. In the last few years, many academics have conducted investigations on lakes in the Tibetan Plateau, which mainly focused on hydrochemical characteristics and its spatial distribution [5,6], typical water pollutants [7], water quality assessment [8,9], and hydrochemistry evolution and sources [10,11]. Due to the vulnerability of the water environment to climate change and the scarcity of water resources in the alpine desert, water quality assessment and hydrochemical evolution in these areas have attracted increasing attention [12,13].
Natural factors, such as the geological structure, hydrogeological conditions, vegetation, climate, rock weathering, terrain, and human activities, greatly impact the chemical composition of water [14]. Sun and Jin pointed out the differences in chemical characteristics and evolution between tectonic and glacial lakes on the plateau based on the analysis of physicochemical parameters and ion concentrations of lakes in the Tibetan Plateau [6]. Li et al. discussed the relationship between lake change and climate in recent decades and clarified that lakes are important supporters of the water cycle and environmental change [15]. Based on spatial analysis of the heavy metal contents in ores, air pollutants, and lake water samples, Bazova and Moiseenko discussed the sources of heavy metals in lakes in northwestern Russia, and the results showed that the concentration of Cu, Ni, Mn, Zn, Pb, and As in lake water was caused by anthropogenic pollution, such as flue gas emissions of metallurgical enterprises [16]. Baranov et al. studied the chemical composition, natural geochemistry, and human factors of dissolved effluents in the northern Valdai Mountains using statistical methods [17]. The results showed that the soil mineral composition, rainfall intensity, and biogeochemical processes had a great influence on the chemical composition of water bodies in the aeration zone. Moreover, many scholars have gained a significant understanding of regional soil salinization, regional geological background, element migration and enrichment law, petrophysical weathering process, and other aspects by studying ion characteristics of watershed water bodies [18]. It can be seen that the hydrochemical characteristics of natural water bodies are one of the important indicators when studying regional water cycles and evolution.
The reserves and quality of water resources are important for the stable and sustainable development of human society and natural ecology. Good water quality is especially valuable in water shortage areas. High evaporation and low rainfall lead to salinization or heavy metal pollution in shallow groundwater or surface water in lake basins in the arid region of the plateau [19]. Previously, water quality assessments were carried out on many lake basins in the Tibetan Plateau, such as Qinghai Lake [20], the Ebinur Lake Watershed [21], Dal Lake in Kashmir Valley [22], Mapam Yumco [23], and Nam Co [24]. The results showed that the water quality of most lakes was relatively uniform and fluctuated in a small range. The results also indicated that lake water with high salinity had a great influence on the shallow surface water of the lake shore, and the salinity of groundwater had a great influence on plant growth, human life, animal husbandry development, and industrial activities in the region.
Bangong Co Lake is a long and narrow tectonic lake with a nearly east-west distribution. The western part of the Tibetan Plateau is characterized by a dry and cold climate with strong evaporation [25]. As the main source of replenishment for the Lake Bangong Co Watershed, glacial meltwater is insufficient to balance evaporation in the area. Likewise, affected by evaporation across the whole watershed and river water supply in the east of Bangong Co Lake, the salinity of the lake gradually increases from east to west [26]. The water circulation system of the plateau arid watershed is very fragile and sensitive, which makes the water resources in the study area vulnerable to the impacts of human activities, animal husbandry activities, and climate change.
Some scholars investigating lake groups on the Tibetan Plateau have pointed out that the lake water type on the east side of Bangong Co Lake is Cl-Na or HCO 3 ·Cl-Mg·Na [27,28]. Lin et al. researched information on toxic organic pollutants and metals of Bangong Co Lake, and their results indicate that polycyclic aromatic hydrocarbons (ΣPAH) and phthalic acid esters (ΣPAE) concentrations in lake water have no relationship with hydrochemical parameters, and organic pollutants have been the main source of domestic waste related to increasing regional human activities in recent years [29]. The domestic water and irrigation water in the study area are mainly shallow aquifers in the lakeshore zone and river valley or supplemented by snowmelt water and river water. However, there are few relevant studies on the various water bodies of Bangong Co Lake Watershed. Therefore, the hydrochemical characteristics and controlling factors of the river water, lake water, groundwater, and other water bodies in the watershed of Bangong Co Lake are still unknown. Additionally, there is no available information on water quality assessment in this region.
The present investigation was performed to explore the glaciers, rivers, lakes, and groundwater along the lakeshore in the Bangong Co Lake Watershed. Statistical analysis, a Gibbs diagram, a Piper diagram, and an ion ratio analysis were used to analyze the hydrochemical characteristics and formation mechanisms of different water bodies. The drinking water quality index (DWQI) and irrigation water quality parameters were used to assess the water quality and its suitability in the study area. The specific objectives of this work are as follows: (i) To explain the hydrochemistry and evolution of different water bodies in Bangong Lake Basin, (ii) to identify the various factors controlling the water chemistry in this area, and (iii) to assess the suitability of the water quality for drinking and irrigation.

Study Area
Bangong Co Lake Watershed is located in the northwest of Tibet, China. Based on ArcGIS hydrological analysis module and DEM data, the watershed area of Bangong Co Lake is about 33,000 km 2 , with the location between 32 • 40 -34 • 30 N and 78 • 10 -81 • 15 E, with an elevation range of 3736-6771 m. The watershed reaches the Kanakoram Mountains in the north and the Gangdises Mountains in the south. The highest point is located in the eastern mountains of Zecuo Lake, and the lowest is located in Bangong Co Lake, in the middle of the area (Figures 1 and 2). The geomorphology of the watershed is that of plateau lake basin geomorphic type, which has been fully eroded by snow and water for a long time. The watershed covers the Bangong Lake replenishment river basin and several small intermountain lake basins, such as Spangour Co, Zecuo, Rebangcuo, Aiyongcuo, and Sharda Co [30]. This region has a cold and dry climate, with an average annual temperature of 0.1-2 • C. The annual rainfall is about 70-80 mm, with high variability, and is mostly concentrated in July or August [31]. The main source of water in the study area is glacial meltwater. Due to gravity, the water mainly flows along the slope or gully to feed rivers or veins of rock and eventually drains into the lake. Because the sediments in river valleys in the middle or lower reaches and the sediments in lake shores are thick and loose, the river water is likely to penetrate, forming relatively stable shallow groundwater in the alluvial fan or valley area.
Geologically, the study area is located in the nearly east-west Bangong Co-Nujiang suture zone, with Qiangtang block to the north and Lhasa block to the south ( Figure 1) [32,33]. Qiangtang block in the region is mainly composed of carbonate rocks, sandstone, and slate rocks from Yanshiping formation, Tunlonggongba formation, and Nuoco formation, and Lhasa Terrane in the study area is mainly composed of ophiolite mélange, Cretaceous gabbros and granites, and sandwiched limestone and flysch sediments from Mugagangri Group and Langshan formation ( Figure 1) [31,34]. Quaternary sediments are mainly distributed on the sides of lake and river valleys. NWW-trending and NE-trending faults are mainly developed in the Bangong Co Lake area. Most faults align with the boundary of the lake basin in geomorphology, and the fault surface is shovel-type dip into the lake. Meanwhile, the NWW-trending faults are the leading faults in the region, belonging to the Bangong Co-Nujiang suture system, which has undergone multiple periods of subduction, collision, and other tectonic processes, and are nearly consistent with the overall tectonic line and the ophiolitic hybrid zone. The NE-trending faults are mostly developed in some north-south gullies and valleys and have the characteristics of water conduction after differentiation and denudation [35].
lies east and southeast of the lake, mainly growing plateau meadows and reeds. Plateau desert meadow is the main vegetation type in the Bangong Co Lake catchment area. The dry period is a maximum of nine months in this area, so the grass yield of the regional meadow is relatively low [36]. Animal husbandry activities are mainly concentrated in the rainy season, and animal husbandry activities are basically practiced in areas with abundant groundwater, such as alluvial fan front, lakeshore, and valley land.  There are many aquatic plants and fishes in Bangong Co Lake. Fishes mainly include Schizothoracinae and Triplophysa stenura, and the aquatic plants encompass Pyrrophyta, Bacillariophyta, and Chlorophyta [28]. Moreover, a small field of wetlands lies east and southeast of the lake, mainly growing plateau meadows and reeds. Plateau desert meadow is the main vegetation type in the Bangong Co Lake catchment area. The dry period is a maximum of nine months in this area, so the grass yield of the regional meadow is relatively low [36]. Animal husbandry activities are mainly concentrated in the rainy season, and animal husbandry activities are basically practiced in areas with abundant groundwater, such as alluvial fan front, lakeshore, and valley land.

Sampling and Measurement
In early June 2022, field investigation was carried out in Bangong Lake Basin, and a total of 60 water samples were collected, including 14 lake water samples, 19 groundwater samples, 24 river water samples, and 3 glacier water samples ( Figure 2). Among them, river samples were mainly distributed in Makazangbu and Doma Rivers, and groundwater samples were mainly collected from wells on alluvial fan bodies in lakeshore zone. The water depth of the collected wells was less than 15 m, and most of the glacier samples were new snow samples. Water samples were mostly collected at a depth of 0.1-0.5 m, filtered through 0.45 μm cellulose acetate membrane, and stored in 1500 mL cleaned polyethylene bottles. Before collection, the collection containers were rinsed with sam-

Sampling and Measurement
In early June 2022, field investigation was carried out in Bangong Lake Basin, and a total of 60 water samples were collected, including 14 lake water samples, 19 groundwater samples, 24 river water samples, and 3 glacier water samples ( Figure 2). Among them, river samples were mainly distributed in Makazangbu and Doma Rivers, and groundwater samples were mainly collected from wells on alluvial fan bodies in lakeshore zone. The water depth of the collected wells was less than 15 m, and most of the glacier samples were new snow samples. Water samples were mostly collected at a depth of 0.1-0.5 m, filtered through 0.45 µm cellulose acetate membrane, and stored in 1500 mL cleaned polyethylene bottles. Before collection, the collection containers were rinsed with sampling water 2 to 3 times. The lake water samples were taken from a location about 1-2 m away from the lake. In addition, water electrical conductivity (EC) and pH were measured in the field using a multi-parameter portable water analyzer (Manta, Eureka Inc., Austin, TX, USA), with a pH accuracy of ±0.1 pH and conductivity accuracy of ±1 µS/cm.
All collected samples were stored in 0-4 • C incubators until they were analyzed in the laboratory. Parallel samples were used to control the test's accuracy during the experiment. For any water sample, a part of it was acidified with 1:1 nitric acid and stored in a plastic bottle for cation analysis. The remaining part was placed in glass bottles and sealed with parafilm sealing film to analyze anions such as SO 4 2− and Cl − . The content of HCO 3 − was measured by acid-base titration, with an accuracy of ±0.05%. Water salinity (TDS) was determined using the dry weight method at 105 • C, with an accuracy of ±0.2%. Total hardness (TH) was determined using the EDTA method with an accuracy of ±0.05%. The main cation (Na + , K + , Ca 2+ , and Mg 2+ ) concentrations and trace heavy metal (Fe, Mn, Zn, Cu, and Cr) levels were determined using inductively coupled plasma emission spectrometer (Agilent5110, Agilent Technologies Inc., Santa Clara, CA, USA), and the testing error was less than 2%. The total As content was detected via atomic fluorescence (AFS-933, Beijing Jitian Instrument Co., LTD, Beijing, China) with an accuracy of ±0.5%. Anion concentration (SO 4 2− , Cl − , F − , and NO 3 − ) was determined via chromatography (ECO IC, Vantone China Co., LTD, Beijing, China), and the testing error was less than 0.5%. The ion balance of the major soluble cationic equivalent (TZ + = Na + + K + + 2Ca 2+ + 2Mg 2+ ) and soluble anionic equivalent (TZ − = Cl − + NO 3 − + HCO 3 − + F − + 2SO 4 2− ) of all samples was compared for ion balance test, and it was within ±5%. All water samples were finished by Xizang Shengyuan Environmental Engineering Co., LTD, Tibet, China.

Data Processing
SPSS v22.0 (SPSS Inc., Chicago, IL, USA) was used to determine the average (Mean), maximum (Max), minimum (Min), and coefficient of variation (CV). The Piper diagrams [37], principal component analysis, Gibbs diagrams [38], and ion ratio analysis [39] were implemented in Origin v2021 (OriginLab, Northampton, MA, USA). and Microsoft Excel v2013 (Microsoft Corp., Redmond, WA, USA) to explore the classification and main controlling factors of the nature water in Lake Bangong Co Watershed. International standards [40], drinking water quality index [41,42], USSL classification [43], and other irrigation water indicators were used to evaluate groundwater suitability for drinking and irrigation. ArcGIS v10.4 (Esri Inc., Redlands, CA, USA) and CorelDRAW Graphics Suite v18 (Corel Corp., Ottawa, ON, Canada) were used for mapping and modification.

Statistical Characteristics of Parameters
The parameters of each water body in Bangong Co Watershed are shown in Table 1. The pH of glaciers, lake water, river water, and shallow groundwater along the lakeshore in the study area ranged from 6.6 to 8.5. The pH of glaciers was relatively low (pH = 6.6~7.6), with an average value of 7.2, which is neutral. Other water bodies were weakly alkaline as a whole, among which the pH of Bangong Co Lake water was relatively high (pH = 8.1~8.5), with an average value of 8.4. The coefficient of variation (CV) for the water pH of each water body was <10%, which showed weak variability. Salinity (TDS) and electrical conductivity (EC) are parameters that can be used to determine the salinity of water [44]. The salinity (TDS) of glaciers, river water, groundwater, and lake samples ranged from 13 to 25, 81 to 438, 124 to 972, and 206 to 4091 mg/L, respectively. The salinity (TDS) values of 10 lake samples, mainly distributed in the western or middle part of the lake, exceeded the WHO's standard limit (1 g/L) [40]. The average salinity (TDS) of river water was 282 mg/L, which was also higher than the global average and values for other large rivers in China [45]. Similarly, the electrical conductivity (EC) of glacier water, river water, groundwater, and lake samples ranged from 23 to 43, 145 to 779, 220 to 1706, and 385 to 7861 µS/cm, respectively. Water hardness is a very important parameter in the growth and reproduction of aquatic organisms [46]. The total hardness (TH) of the glaciers, groundwater, river, and lake were in the range of 23-43 mg/L, 27-750 mg/L, 58-467 mg/L, and 258-2854 mg/L, respectively. With the decrease in altitude, the values regarding the total hardness (TH), electrical conductivity (EC), and salinity (TDS) increased gradually, which were the highest around Bangong Co Lake. The coefficient of variation (CV%) for total hardness (TH), electrical conductivity (EC), and salinity (TDS) of different water bodies were between 10% and 66%, which demonstrates medium variability. Based on the total hardness (TH) and salinity (TDS) [47], lake water is hard-brackish water, glacier water is soft-fresh water, and river water changes from soft-fresh water to hard-fresh water along the stream ( Figure 3). river water, groundwater, and lake samples ranged from 23 to 43, 145 to 779, 220 to 1706, and 385 to 7861 μS/cm, respectively. Water hardness is a very important parameter in the growth and reproduction of aquatic organisms [46]. The total hardness (TH) of the glaciers, groundwater, river, and lake were in the range of 23-43 mg/L, 27-750 mg/L, 58-467 mg/L, and 258-2854 mg/L, respectively. With the decrease in altitude, the values regarding the total hardness (TH), electrical conductivity (EC), and salinity (TDS) increased gradually, which were the highest around Bangong Co Lake. The coefficient of variation (CV%) for total hardness (TH), electrical conductivity (EC), and salinity (TDS) of different water bodies were between 10% and 66%, which demonstrates medium variability. Based on the total hardness (TH) and salinity (TDS) [47], lake water is hard-brackish water, glacier water is soft-fresh water, and river water changes from soft-fresh water to hard-fresh water along the stream ( Figure 3).

Ion Spatial Distribution
The ion diagrams ( Figure 4) depict the milligram equivalent and spatial variation in anions and cations in the different water bodies across the basin. In the water samples of groundwater, river water, and glaciers, the cations were dominated by Ca 2+ , and the order of abundance was Ca 2+ > Na + > K + > Mg 2+ . However, Na + was the main cation of lake water, and the cation content followed the order Na + > Mg 2+ > K + > Ca 2+ . The major anion in groundwater, river water, and glacier samples was HCO 3 − , and the order of the abundance in these samples was HCO 3 − > SO 4 2− > Cl − > NO 3 − , whereas the primary anion in lake samples was Cl − , and the anion content followed the order Cl − > HCO 3 − > SO 4 2− > NO 3 − . A strong coefficient of variation was found in Na + , Cl − , and SO 4 2− of groundwater, and NO 3 − of river and lake water, which was between 1.15 and 2.85, indicating that the spatial distribution of these ionic components was unstable. The coefficient of variation of most chemical parameters of different water bodies in the area was less than 1 (Table 1), which is a weak or medium variation. Spatially, the Doma River and Maka Zangbo River showed a gradually increasing trend of ions from upstream to downstream. The main ions concentration in the lake increased gradually from east to west ( Figure 4c). However, for lake samples L7 to L5 (Figures 2 and 4c), the ions and salinity (TDS) increased slowly compared to the other sections, which may be due to a water exchange barrier caused by the narrow body of the lake [50]. According to Figure 4d, there is no obvious spatial rule of groundwater in the lakeshore zone. which is a weak or medium variation. Spatially, the Doma River and Maka Zangbo River showed a gradually increasing trend of ions from upstream to downstream. The main ions concentration in the lake increased gradually from east to west (Figure 4c). However, for lake samples L7 to L5 (Figures 2 and 4c), the ions and salinity (TDS) increased slowly compared to the other sections, which may be due to a water exchange barrier caused by the narrow body of the lake [50]. According to Figure 4d, there is no obvious spatial rule of groundwater in the lakeshore zone. Al, Cr, Pb, Cd, Hg, Cu, and Ni were below the limit of quantification. Moreover, Fe, Mn, As, and F − were found at very insignificant concentrations. Table 1 shows the summary of the datasets of trace metals. Concentrations of Fe in various water bodies varied from ND to 0.224 mg/L, which was within the WHO's tolerable limit (0.3 mg/L). Mean- Al, Cr, Pb, Cd, Hg, Cu, and Ni were below the limit of quantification. Moreover, Fe, Mn, As, and F − were found at very insignificant concentrations. Table 1 shows the summary of the datasets of trace metals. Concentrations of Fe in various water bodies varied from ND to 0.224 mg/L, which was within the WHO's tolerable limit (0.3 mg/L). Meanwhile, concentrations of As in lake samples varied from 0.002 to 0.035 mg/L, with a mean concentration of 0.021 ± 0.012 mg/L, most of which exceeded the WHO's tolerable limits (0.01 mg/L). The results were similar to those of previous studies, which suggested that the As in Bangong Co Lake was mainly derived from natural sources such as rock weathering [29]. Mn concentrations in the river and groundwater ranged from ND to 0.398 mg/L and ND to 0.441 mg/L. 4.2% (n = 1) in river water samples and 15.8% (n = 3) in shallow groundwater, exceeding the WHO's limit for Mn content (0.1 mg/L). However, few lake water samples had Mn content exceeding the WHO's standard. The Mn exceedance points of groundwater are located near the shoreline of the lake, with the pH ranging from 6.83 to 7.72. The river water exceedance point is located at a spring recharge in the lower Doma River, and the pH is 7.19. Low water pH provides favorable conditions for Mn enrichment in groundwater. On the contrary, Mn and Fe ions easily form hydroxide precipitation in lake water with relatively high pH (8.06-8.52). The flu-lacustrine purplish neogene conglomerates, rich in Fe, Mn, and other elements, are scattered in some river valleys in the study area, which may provide Mn sources for moderate and acidic groundwater. F − concentrations of various water bodies were all below the acceptable limit according to WHO standards.

Piper Diagram and Hydrochemical Classification
Water hydrochemical characteristics can be classified using the Piper diagram, which shows scatter plots of cations (Na + , K + , Ca 2+ , and Mg 2+ ) and anions (HCO 3 − , Cl − , and SO 4 2− ) [37,42]. According to Figure 5, water samples in the Bangong Co Lake Watershed were classified into three water types. River samples, glaciers, and most of the groundwater corresponded to zone 2 of the Piper diagram, characterized as Ca-HCO 3 -type samples. The lake samples were distributed in zone 5, classified as saline water with high Cl − or SO 4 2− and Na + . Due to Cl − being the main cation of lake samples, lake water was classified as a Na-Cl-type sample. Only a few groundwater samples corresponded to zone 3, representing contaminated water with high Cl − and Ca 2+ concentrations. Eastern water samples were in the transition between the river point and the lake point, indicating that these water samples exhibit ion exchange and simple dissolution or mixing.

Gibbs diagrams
Based on an analysis of many global hydrochemical compositions of surface waters, such as large lakes, rivers, and precipitation water samples, diagrams were proposed by Gibbs to determine the features of ionic distribution in different natural water bodies [38]. Rock weathering, atmospheric precipitation, and evaporation-crystallization process were identified as three major factors controlling the surface water chemistry in Gibbs plots [11,51]. In this paper, a TDS-Na + /(Na + + Ca 2+ ) diagram and TDS-Cl − /(Cl − + HCO3 − ) diagram were used to distinguish the ionic characteristics of water bodies in the Bangong Co Lake Watershed ( Figure 6).
With regard to the river and groundwater samples, the cation ratios Na + /(Na + + Ca 2+ ) varied from 0.07 to 0.66, and the anion ratios Cl − /(Cl − + HCO3 − ) ranged from 0.02 to 0.43. Additionally, all the samples for both diagrams fell in the rock weathering domain. For the lake water samples, the ranges of the cation ratios Na + /(Na + + Ca 2+ ) and the anion ratios Cl − /(Cl − + HCO3 − ) were 0.66-0.99 and 0.02-0.52, respectively. The western lake samples mainly fell into the evaporation domain (Figure 6c) or near the evaporation dominance zone (Figure 6b), and the eastern lake samples fell in the transition zone between evaporation and rock weathering (Figure 6b,c). This indicates that the dominant ions of various water bodies were significantly affected by rock weathering and evaporationcrystallization may affect the western lake water and some groundwater. Duoma and Makazangbu rivers which were the main sources of replenishment, joined Bangong Co Lake on the east side, resulting in a characteristic transition between evaporation and rock weathering in the eastern lake. Glacier water is considered to be a precipitation sample in the Bangong Co Lake area. The glacier water samples are part of the end ele-

Gibbs Diagrams
Based on an analysis of many global hydrochemical compositions of surface waters, such as large lakes, rivers, and precipitation water samples, diagrams were proposed by Gibbs to determine the features of ionic distribution in different natural water bodies [38]. Rock weathering, atmospheric precipitation, and evaporation-crystallization process were identified as three major factors controlling the surface water chemistry in Gibbs plots [11,51]. In this paper, a TDS-Na + /(Na + + Ca 2+ ) diagram and TDS-Cl − /(Cl − + HCO 3 − ) diagram were used to distinguish the ionic characteristics of water bodies in the Bangong Co Lake Watershed ( Figure 6).
With regard to the river and groundwater samples, the cation ratios Na + /(Na + + Ca 2+ ) varied from 0.07 to 0.66, and the anion ratios Cl − /(Cl − + HCO 3 − ) ranged from 0.02 to 0.43. Additionally, all the samples for both diagrams fell in the rock weathering domain. For the lake water samples, the ranges of the cation ratios Na + /(Na + + Ca 2+ ) and the anion ratios Cl − /(Cl − + HCO 3 − ) were 0.66-0.99 and 0.02-0.52, respectively. The western lake samples mainly fell into the evaporation domain (Figure 6c) or near the evaporation dominance zone (Figure 6b), and the eastern lake samples fell in the transition zone between evaporation and rock weathering (Figure 6b,c). This indicates that the dominant ions of various water bodies were significantly affected by rock weathering and evaporation-crystallization may affect the western lake water and some groundwater. Duoma and Makazangbu rivers which were the main sources of replenishment, joined Bangong Co Lake on the east side, resulting in a characteristic transition between evaporation and rock weathering in the eastern lake. Glacier water is considered to be a precipitation sample in the Bangong Co Lake area. The glacier water samples are part of the end element of rock weathering and are slightly closer to the rain zone in Gibbs diagrams. Therefore, the dissolved ions in the rainwater of the Bangong Co Lake Watershed are slightly affected by ocean evaporation and are mainly controlled by the dissolution of atmospheric CaCO 3 particles. Therefore, the dissolved ions in the rainwater of the Bangong Co Lake Watershed are slightly affected by ocean evaporation and are mainly controlled by the dissolution of atmospheric CaCO3 particles.

Ion ratio analysis
The cations and anions dissolved by the chemical weathering of different rocks contribute to the combination of water in nature [52]. Therefore, the ion ratio characteristics can reflect the water-rock reaction and evolution process between natural water and rock strata [39]. As can be seen in the regional geological map (Figure 1), the main runoff areas of rivers are mainly developed strata, such as the Mugagangri Group, Langshan formation, Yanshiping formation, and Nuoco formation. The lithology of these strata is principally carbonate rock, sandstone, and conglomerate. Na + and K + are mainly derived from evaporite or silicate, and Ca 2+ and Mg 2+ are mainly derived from carbonate and silicate weathering and evaporative dissolution [53]. Thus, when the ratio of γ(Na + + K + ) and γ(Cl − ) is close to or above the 1:1 line, dissolution of evaporate is the leading role for Na + and K + in water chemical evolution. Figure 7a shows that the sample points of γ(Na + + K + ) and γ(Cl − ) are distributed around a straight line, with a slope of 1.63 and R 2 of 0.99. Most of the water bodies (except glaciers) are almost on the upper side of the 1:1 line, indicating that halite (NaCl) and potassium salt (KCl) are the main sources of Na + , K + , and Cl − , and the Cl − content is not enough to balance the content of Na + and K + in water. The excess Na + and K + may be derived from the weathering and dissolution of silicate rock.

Ion Ratio Analysis
The cations and anions dissolved by the chemical weathering of different rocks contribute to the combination of water in nature [52]. Therefore, the ion ratio characteristics can reflect the water-rock reaction and evolution process between natural water and rock strata [39]. As can be seen in the regional geological map (Figure 1), the main runoff areas of rivers are mainly developed strata, such as the Mugagangri Group, Langshan formation, Yanshiping formation, and Nuoco formation. The lithology of these strata is principally carbonate rock, sandstone, and conglomerate. Na + and K + are mainly derived from evaporite or silicate, and Ca 2+ and Mg 2+ are mainly derived from carbonate and silicate weathering and evaporative dissolution [53]. Thus, when the ratio of γ(Na + + K + ) and γ(Cl − ) is close to or above the 1:1 line, dissolution of evaporate is the leading role for Na + and K + in water chemical evolution. Figure 7a shows that the sample points of γ(Na + + K + ) and γ(Cl − ) are distributed around a straight line, with a slope of 1.63 and R 2 of 0.99. Most of the water bodies (except glaciers) are almost on the upper side of the 1:1 line, indicating that halite (NaCl) and potassium salt (KCl) are the main sources of Na + , K + , and Cl − , and the Cl − content is not enough to balance the content of Na + and K + in water. The excess Na + and K + may be derived from the weathering and dissolution of silicate rock.  Carbonates, evaporates, or silicates containing calcium and magnesium are the main sources of Ca 2+ and Mg 2+ ions in natural water [24]. Figure 7b shows the ratio of (Ca 2+ + Mg 2+ ) to (HCO3 − + SO4 2− ) in the water sample. The slope of the scattered fitting line of the water samples is 0.43, and the R 2 is 0.85, indicating that many sample points are below the 1:1 line, especially the lake points. In these water bodies, Ca 2+ and Mg 2+ contents are not high enough to reach equilibrium with HCO3 − and SO4 2− , and it is generally believed that the dissolution of evaporites that contain mirabilite (NaSO4), polyhalite (K2Ca2Mg[SO4]4·2H2O), etc., may provide abundant SO4 2− or HCO3 − for the water bodies.
The ratio of γ(Ca 2+ + Mg 2+ ) and γ(HCO3 − ) can further be used to investigate the source of Ca 2+ and Mg 2+ [6]. As can be seen in Figure 7c, these sample points fit a line with a slope of 1.01, and the R 2 is 0.85, which is close to 1:1 line. This indicates that Ca 2+ , Mg 2+ , and HCO3 − in river water, groundwater, and other water bodies are mainly derived from the dissolution of dolomite (CaMg(CO3)2) and calcite (CaCO3) rather than gypsum (CaSO4).
Generally, silicates are more difficult to weather than carbonates. Thus, the ratio of (Ca 2+ + Mg 2+ )/(Na + + K + ) or γ(HCO3 − )/γ(Na + + K + ) in the water can be used as a symbol for judging the main types of weathered rocks in natural water [24]. As shown in Figure 7d, there are obvious differences between lake water and other water points. The high ratios of γ(Ca 2+ + Mg 2+ )/γ(Na + + K + ) in the river water and groundwater indicate flow through the carbonate weathering regions (Figure 7d). Additionally, the lake samples were lo- Carbonates, evaporates, or silicates containing calcium and magnesium are the main sources of Ca 2+ and Mg 2+ ions in natural water [24]. Figure 7b shows the ratio of (Ca 2+ + Mg 2+ ) to (HCO 3 − + SO 4 2− ) in the water sample. The slope of the scattered fitting line of the water samples is 0.43, and the R 2 is 0.85, indicating that many sample points are below the 1:1 line, especially the lake points. In these water bodies, Ca 2+ and Mg 2+ contents are not high enough to reach equilibrium with HCO 3 − and SO 4 2− , and it is generally believed that the dissolution of evaporites that contain mirabilite (NaSO 4 ), polyhalite (K 2 Ca 2 Mg[SO 4 ] 4 ·2H 2 O), etc., may provide abundant SO 4 2− or HCO 3 − for the water bodies. The ratio of γ(Ca 2+ + Mg 2+ ) and γ(HCO 3 − ) can further be used to investigate the source of Ca 2+ and Mg 2+ [6]. As can be seen in Figure 7c, these sample points fit a line with a slope of 1.01, and the R 2 is 0.85, which is close to 1:1 line. This indicates that Ca 2+ , Mg 2+ , and HCO 3 − in river water, groundwater, and other water bodies are mainly derived from the dissolution of dolomite (CaMg(CO 3 ) 2 ) and calcite (CaCO 3 ) rather than gypsum (CaSO 4 ).
Generally, silicates are more difficult to weather than carbonates. Thus, the ratio of (Ca 2+ + Mg 2+ )/(Na + + K + ) or γ(HCO 3 − )/γ(Na + + K + ) in the water can be used as a symbol for judging the main types of weathered rocks in natural water [24]. As shown in Figure 7d, there are obvious differences between lake water and other water points. The high ratios of γ(Ca 2+ + Mg 2+ )/γ(Na + + K + ) in the river water and groundwater indicate flow through the carbonate weathering regions (Figure 7d). Additionally, the lake samples were located on the lower side of the 1:1 line, which confirmed that silicate dissolution or evaporate dissolution contributes to the main ion characteristics (Figure 7d).

Hydrochemical Modeling of Mineral Saturation Index (SI)
The saturation index (SI) of the main minerals in natural water indicates the equilibrium state between minerals and water and distinguishes the dissolution and precipitation reaction of minerals to simulate the trend of water chemistry in the environment. Phreeqc software is a hydrochemical simulation software developed by the United States District Survey, which can calculate the saturation index of minerals in groundwater. The saturation index (SI) is calculated as where IAP is the ionic activity product, and K is the equilibrium constant. SI > 0 indicates that the mineral is saturated, and precipitate occurs. SI = 0 indicates that the mineral reached the equilibrium state, relatively stable. SI < 0 indicates that the mineral is not saturated, and a dissolution reaction occurs [51].
The ion ratio analysis of water in the study area showed that halite and carbonate minerals were the main minerals of the water-rock reaction. The solubility relationship between dolomite, calcite, and halite is shown in Figure 8 The saturation index (SI) of the main minerals in natural water indicates the equilibrium state between minerals and water and distinguishes the dissolution and precipitation reaction of minerals to simulate the trend of water chemistry in the environment. Phreeqc software is a hydrochemical simulation software developed by the United States District Survey, which can calculate the saturation index of minerals in groundwater. The saturation index (SI) is calculated as

= lg
where IAP is the ionic activity product, and K is the equilibrium constant. SI > 0 indicates that the mineral is saturated, and precipitate occurs. SI = 0 indicates that the mineral reached the equilibrium state, relatively stable. SI < 0 indicates that the mineral is not saturated, and a dissolution reaction occurs [51].
The ion ratio analysis of water in the study area showed that halite and carbonate minerals were the main minerals of the water-rock reaction. The solubility relationship between dolomite, calcite, and halite is shown in Figure 8. There are some differences in the saturation index (SI) of minerals in various water bodies. The salt mines have the largest halite SI, varying from −3.96 to 0.27. The ranges of halite SI in the lake, groundwater, river water, and glaciers are −8. The results showed that the saturation state of calcite, dolomite, and halite was not reached in glacier meltwater and upstream river water bodies. With the contact between the river and the carbonate rock strata in the flowing area, calcite and dolomite in the river bodies in the middle and lower and Bangong Co Lake water reached saturation. Overall, the halite saturation index of Bangong Co Lake showed an increasing trend from east to west, but it did not reach the saturation state like the river water. This indicates that the calcite and dolomite gradually become saturated from unsaturated to saturated, while halite dissolves during water flow from upstream to downstream and into the lake. The data collected from the salt lakes in this area show that the saturation index of calcite, dolomite, and halite in some salt lakes is greater than 0, indicating the precipitation of calcite, dolomite, and halite.  Table 1). Figure 8. SI relationship between halite/dolomite (a) and halite/calcite (b) (the salt lake samples are data collected from the literature, as shown in Table 1).
The results showed that the saturation state of calcite, dolomite, and halite was not reached in glacier meltwater and upstream river water bodies. With the contact between the river and the carbonate rock strata in the flowing area, calcite and dolomite in the river bodies in the middle and lower and Bangong Co Lake water reached saturation. Overall, the halite saturation index of Bangong Co Lake showed an increasing trend from east to west, but it did not reach the saturation state like the river water. This indicates that the calcite and dolomite gradually become saturated from unsaturated to saturated, while halite dissolves during water flow from upstream to downstream and into the lake. The data collected from the salt lakes in this area show that the saturation index of calcite, dolomite, and halite in some salt lakes is greater than 0, indicating the precipitation of calcite, dolomite, and halite.
Previous studies [35,48] pointed out that the Bangong Co Lake basin experienced a period of basin shrinkage during the plateau's uplift. During this time, some independent intermountain lakes entered the evaporite-deposition period due to abrupt regional drought climate, forming the characteristics of dispersed salt lake salt mines in the area, such as the Jieze salt mine [48] and the Taicuo salt mine [49] (Figure 1). Through the collection of chemical composition data on the salt lake water body and salt mine soil layer (Table 1), it can be seen that there are rich minerals such as halite (NaCl), dolomite (CaMg(CO 3 ) 2 ), calcite (CaCO 3 ), and mirabilite (NaSO 4 ) in these lacustrine chemical sedimentary profiles. Under rain erosion or lateral river water intrusion, the interlacustrine group continuously replenishes rivers and lakes, providing abundant Na + , K + , Cl − , and SO 4 2− ions. Therefore, the rich carbonate strata, the rock salt minerals produced in the historical environment, and the arid climate environment in the study area lead to the diverse hydrochemical characteristics in the basin.

Principal Component Analysis (PCA)
Principal component analysis (PCA) can reduce the dimension of the data set by converting the original variables into new and uncorrelated variables generated by retaining original information [54]. Many studies have used the principal component analysis (PCA) technique to identify important water quality parameters [55,56]. Before component analysis, all data passed the Kaiser-Meyer-Olkin (KMO) value (0.704) and Bartlett's test of sphericity statistics (p < 0.05) by using Origin2022 to evaluate the feasibility of principal component analysis (PCA) for source apportionment [57].
By filtering principal components (PCs) with an eigenvalue greater than 1 [29], only two PCs were extracted from the scree plot and explained 89.0% of the total variance (Table 2 and Figure 9). K + , Na + , Mg 2+ , HCO 3 − , Cl − , SO 4 2− , TDS, As, and EC have a relatively high loading on the principal component PC1, accounting for 75.4% of the total variance. PC2 explained 13.6% of the total variance and has a strong positive loading for Ca 2+ and NO 3 − and weak negative loading for pH. Given the water sample loading distribution (Figure 9a), the loading of lake water points in PC1 is higher than that of other water bodies. The high salinity (TDS) and other ions in lake water are mainly caused by the dry environment and the lithology of the regional strata. Therefore, PC1 mainly indicates the characteristics of natural sources. Lin, Dong, and Wang pointed out that the concentration of heavy metals, such as As, Mn, and Ni, in the water body of Bangong Co Lake was lower than the background value of the water body in Tibetan areas, indicating that those heavy metals in the water were mainly derived from natural rock weathering [29]. This also verified that PC1 containing high loadings of As and TDS in this study mainly represented natural sources. PC2 comprised Ca 2+ , NO 3 − , and pH. According to the loading distribution of the sample points (Figure 9a), the points with high PC2 loadings are mainly groundwater and a few river water points, which are mainly located in the human living, farming area, and the valley area where the meadow grows in the alluvial fan of the Doma River. This indicates that PC2 mainly illustrates animal husbandry activities and human activities. The middle and lower reaches of the Doma River and Makazangbu River are important animal husbandry areas, agricultural irrigation areas, and human living areas, indicating that agricultural fertilization and human activities have a certain impact on the water environment.

Assessment of Groundwater Suitability for Drinking Purposes
A drinking water quality index (DWQI) method was used to assess the suitability of groundwater for drinking [41,42]. The drinking water quality index (DWQI) is a water quality assessment method proposed by Horton RK and widely used [58]. This technique used the weighted index method to change many water quality parameters into a single index that can be compared, effectively providing a comprehensive groundwater quality evaluation model [59]. The calculation of the drinking water quality index (DWQI) can be divided into several steps: (a) Assignment of each parameter's weight according to its relative importance (values from 2 to 5), (b) calculation of the relative weight for each

Assessment of Groundwater Suitability for Drinking Purposes
A drinking water quality index (DWQI) method was used to assess the suitability of groundwater for drinking [41,42]. The drinking water quality index (DWQI) is a water quality assessment method proposed by Horton RK and widely used [58]. This technique used the weighted index method to change many water quality parameters into a single index that can be compared, effectively providing a comprehensive groundwater quality evaluation model [59]. The calculation of the drinking water quality index (DWQI) can be divided into several steps: (a) Assignment of each parameter's weight according to its relative importance (values from 2 to 5), (b) calculation of the relative weight for each index (W i ), (c) the rate calculation of the quality parameter (Q i ), and (d) the water quality index computation of each sample. These data formulas are calculated via the following equations: where W i is the relative weight, w i is the weight of each parameter, n is the total number of parameters considered for the drinking water quality index (DWQI) calculation, and C i is the concentration of each parameter (mg/L). C ip is the ideal value of the parameter in pure water, and S i is the standard value of each parameter.
In this study, an analysis of eleven water quality parameters, namely, pH, TH, TDS, Cl − , SO 4 2− , NO 3 − , F − , Na + , As, and Mn, was used to evaluate the suitability of groundwater for drinking. The relative weight of each parameter and its weight used in these calculations are presented in Table 3. Based on drinking water quality index (DWQI) values [41,51,60], the groundwater quality status can be categorized into five types: Excellent water (<20), good water (≥20 and <40), poor water (≥40 and <80), very poor water (≥80 and <120) and water unsuitable for drinking (≥120). In the Bangong Co Lake Watershed, the computed DWQI values vary from 12 to 86, with an average value of 32 ± 20. According to the DWQI classification, 26% (n = 5) of the total groundwater samples fall under the excellent category, and 53% (n = 10) are classified as good water, respectively. About 21% (n = 4) of groundwater samples fell under moderate water quality (poor and very poor), and no sample was labeled as unsuitable for drinking ( Figure 10). Figure 10 shows the spatial distribution of the drinking water quality index (DWQI), which illustrates that the chemical characteristics of groundwater in the lakeshore zone do not have a regular distribution in the east-west direction. Meanwhile, wells located near the lakeshore have higher drinking water quality index (DWQI) values than those away from the lakeshore, indicating that groundwater near the lakeshore may be affected by lake water (Figure 10).  Figure 10. Assessment map of groundwater suitability in Bangong Lake Watershed.

Water Evaluation for Irrigation
In the plateau and semi-arid region of Bangong Co Lake Watershed, agriculture and livestock husbandry are the main occupations, and the main crops in the region are barley, rape, wheat, and peas [61,62]. Agriculture is the basic department of regional economic and social development. In addition, afforestation is an initiative advocated by the local government to improve the ecological environment [61]. Due to the fact that water quality plays an important role in crops, the survival of tree planting, and soil characteristics in the study area, groundwater irrigation suitability evaluation is necessary. To better evaluate the suitability of regional surface water and groundwater, this study uses the USSL classification proposed by Richards [43] and several other methods, such as electrical conductivity (EC), adsorption ratio (SAR), permeability index (PI), and magnesium hazardous ratio (MHR) [43,63,64]. Table 4 shows the calculation methods and sources of these parameters. High salinity and high sodium concentration in irrigation water are the main causes of soil salinization, which affects the growth of plants and crops [66]. The sodium adsorption ratio (SAR), recommended by Richards [43], is one of the important indices used to calculate the harm of sodium in irrigation water. In the Bangong Co Lake Watershed, the sodium adsorption ratio (SAR) value ranged from 0.03 to 16.72, with a mean value of 3.05 ± 4.92. Based on the classification of sodium adsorption ratio (SAR) values [43,67], water samples in the Bangong Co Lake Watershed are classified into good (n = 51, 85%), poor (n = 1, 2%), and unsuitable (n = 8, 13%). The water points with sodium adsorption ratio (SAR) values exceeding the limits are located in the west and middle parts of Bangong Co Lake. Electrical conductivity (EC), the salinity parameter of irrigation water, ranges from 23 to 7861 μS/cm with an average value of 1493 ± 2179 μS/cm. According to electrical conductivity (EC) value classification, water samples in the Bangong Co Lake

Water Evaluation for Irrigation
In the plateau and semi-arid region of Bangong Co Lake Watershed, agriculture and livestock husbandry are the main occupations, and the main crops in the region are barley, rape, wheat, and peas [61,62]. Agriculture is the basic department of regional economic and social development. In addition, afforestation is an initiative advocated by the local government to improve the ecological environment [61]. Due to the fact that water quality plays an important role in crops, the survival of tree planting, and soil characteristics in the study area, groundwater irrigation suitability evaluation is necessary. To better evaluate the suitability of regional surface water and groundwater, this study uses the USSL classification proposed by Richards [43] and several other methods, such as electrical conductivity (EC), adsorption ratio (SAR), permeability index (PI), and magnesium hazardous ratio (MHR) [43,63,64]. Table 4 shows the calculation methods and sources of these parameters. High salinity and high sodium concentration in irrigation water are the main causes of soil salinization, which affects the growth of plants and crops [66]. The sodium adsorption ratio (SAR), recommended by Richards [43], is one of the important indices used to calculate the harm of sodium in irrigation water. In the Bangong Co Lake Watershed, the sodium adsorption ratio (SAR) value ranged from 0.03 to 16.72, with a mean value of 3.05 ± 4.92. Based on the classification of sodium adsorption ratio (SAR) values [43,67], water samples in the Bangong Co Lake Watershed are classified into good (n = 51, 85%), poor (n = 1, 2%), and unsuitable (n = 8, 13%). The water points with sodium adsorption ratio (SAR) values exceeding the limits are located in the west and middle parts of Bangong Co Lake. Electrical conductivity (EC), the salinity parameter of irrigation water, ranges from 23 to 7861 µS/cm with an average value of 1493 ± 2179 µS/cm. According to electrical conductivity (EC) value classification, water samples in the Bangong Co Lake Watershed can be classified as good (n = 35, 58%), poor (n = 17, 28%), and unsuitable (n = 8, 13%). Similarly, the water points with severe EC values are located in the west and middle parts of Bangong Co Lake, and points with poor electrical conductivity (EC) values are mainly located in the eastern part of Bangong Co Lake and groundwater near the lake's shoreline and lower sections of streams.
Doneen pointed out that poor soil permeability is usually caused by irrigation water with high Na + and HCO 3 − concentrations [63]. Thus, the permeability index (PI) is often used to evaluate the harm of a high Na ion concentration in soil in irrigation water. In this study, the permeability index (PI) ranged from 36 to 181, averaging 71 ± 26. According to Wilcox's classification, watershed water samples were classified as good (n = 49, 32%) or suitable (n = 41, 68%) [67,68].
Ayers and Westcot noticed that irrigation water with high Mg 2+ could lead to soil Ca 2+ deficiency and crop yield reduction [64]. In this study, the magnesium hazardous ratio (MHR) ranged from 10 to 98, with a mean of 54 ± 23. The results show that 55% of water samples are suitable for irrigation (MHR < 50), and the rest of the water samples (n = 27, 45%) are unsuitable (MHR ≥ 50). The unsuitable points for irrigation of MHR are mainly the lake's water points, some river water samples, and some groundwater samples.
The USSL classification suggested by US Salinity Laboratory Staff (1954) is a practical way to assess the suitability of irrigation water [43]. The USSL diagram best explains the combined effect of the salinity hazard and sodium hazard. In this paper, the USSL diagram ( Figure 11) shows that 13% of the total samples are in the C1-S1 field, including three glacier samples, two groundwater samples, and three river samples (distributed upstream) which indicates relatively low alkalinity and salinity hazards. In the C2-S1 field, there are 10 groundwater samples, 19 river samples, and 1 lake sample with a medium of alkalinity hazard and low salinity. A total of seven groundwater samples and four lake samples (mainly in the eastern lake) are classified as part of the C3-S1 region. These samples with low sodium (S1) and low, medium, or high alkalinity (C1, C2, C3) apply to all soil types of irrigation. Only one lake sample was plotted on the C4-S2 region, and western lake samples were almost all plotted on the C5-S4 region, which is unsuitable for crops and soil irrigation. In conclusion, glacial meltwater and river water are more suitable for irrigation than lake water and some groundwater. Watershed can be classified as good (n = 35, 58%), poor (n = 17, 28%), and unsuitable (n = 8, 13%). Similarly, the water points with severe EC values are located in the west and middle parts of Bangong Co Lake, and points with poor electrical conductivity (EC) values are mainly located in the eastern part of Bangong Co Lake and groundwater near the lake's shoreline and lower sections of streams. Doneen pointed out that poor soil permeability is usually caused by irrigation water with high Na + and HCO3 − concentrations [63]. Thus, the permeability index (PI) is often used to evaluate the harm of a high Na ion concentration in soil in irrigation water. In this study, the permeability index (PI) ranged from 36 to 181, averaging 71 ± 26. According to Wilcox's classification, watershed water samples were classified as good (n = 49, 32%) or suitable (n = 41, 68%) [67,68].
Ayers and Westcot noticed that irrigation water with high Mg 2+ could lead to soil Ca 2+ deficiency and crop yield reduction [64]. In this study, the magnesium hazardous ratio (MHR) ranged from 10 to 98, with a mean of 54 ± 23. The results show that 55% of water samples are suitable for irrigation (MHR < 50), and the rest of the water samples (n = 27, 45%) are unsuitable (MHR ≥ 50). The unsuitable points for irrigation of MHR are mainly the lake's water points, some river water samples, and some groundwater samples.
The USSL classification suggested by US Salinity Laboratory Staff (1954) is a practical way to assess the suitability of irrigation water [43]. The USSL diagram best explains the combined effect of the salinity hazard and sodium hazard. In this paper, the USSL diagram ( Figure 11) shows that 13% of the total samples are in the C1-S1 field, including three glacier samples, two groundwater samples, and three river samples (distributed upstream) which indicates relatively low alkalinity and salinity hazards. In the C2-S1 field, there are 10 groundwater samples, 19 river samples, and 1 lake sample with a medium of alkalinity hazard and low salinity. A total of seven groundwater samples and four lake samples (mainly in the eastern lake) are classified as part of the C3-S1 region. These samples with low sodium (S1) and low, medium, or high alkalinity (C1, C2, C3) apply to all soil types of irrigation. Only one lake sample was plotted on the C4-S2 region, and western lake samples were almost all plotted on the C5-S4 region, which is unsuitable for crops and soil irrigation. In conclusion, glacial meltwater and river water are more suitable for irrigation than lake water and some groundwater.

Conclusions
The hydrochemical characteristics of surface water and lakeshore groundwater in Bangong Lake Watershed, Northwest Tibet, China, were analyzed using a multivariate statistical method, Piper diagrams, Gibbs diagrams, the ion ratio method, and principal component analysis (PCA). At the same time, the drinking water quality index (DWQI), irrigation water indexes, and spatial analysis were used to evaluate the quality of irrigation water and drinking water. The main conclusions are summarized below.
Most of the water samples, such as lake and river samples in the study area, were found to be slightly alkaline. Glacier water was categorized as soft-fresh, and lake water was classified as hard-brackish. Groundwater and river water were categorized as soft-fresh and hard-fresh. Na-Cl was the primary water type for lake water, while glacier, river, and groundwater were dominated by the Ca-HCO 3 type. Only 78% of groundwater samples were Ca-Cl types. The ionic concentrations of river waters showed an increased trend from upstream to downstream in terms of their spatial distribution, and high concentrations were found in the lake's coastal region. Regarding the spatial patterns of water in the Bangong Co Lake, ionic concentrations increased from east to west. There was no obvious spatial distribution of groundwater ion concentrations in the lakeshore zone in the east-west direction.
Precipitation dominance was observed in rivers and groundwater, and evaporation and crystallization dominance was observed in lakes, according to the Gibbs diagram. The ion ratio showed that the dissolution of carbonates and evaporites plays a leading role in the regional hydrochemical characteristics, corresponding to the limestone and dolomite formations widely exposed in the region. The calculation of the saturation index showed that calcite and dolomite were supersaturated in the water body of the basin, except for glaciers and some upstream rivers. Although all the lake samples from the study area were undersaturated regarding halite, their saturation index showed an increasing trend from east to west, indicating strong evaporation and continuous halite concentration characteristics. The collected data revealed salt lakes supersaturated with halite in this area. Thus, the main mineralization is related to the dissolution of calcite, dolomite, and halite in the strata. PCA showed that the regional hydrochemistry was mainly affected by natural factors, and only a few river sections and groundwater points had an anthropogenic input of NO 3 − and Ca 2+ plasma, which may be related to animal husbandry activities in the valley and human activities in the river sector.
Regarding drinking water quality, the drinking water quality index (DWQI) value classified the groundwater as excellent (67%), good (5%), and poor (19%). The map of the suitability assessment for drinking water demonstrated that the closer the shallow groundwater was to the lake shore, the worse the water quality. With regard to the irrigation water quality, the irrigation water indexes indicated that glaciers, rivers, most groundwater, and a small part of the east lake are suitable for irrigation, while these indicators also indicated that the west and middle lake (n = 8, 13%) are not suitable for irrigation. The permeability index (PI), electrical conductivity (EC), and USSL classification showed that some shallow groundwater near the lake shoreline and the lake water on the east side have moderate hazards in terms of salinity and permeability, and it is necessary to consider rational use for different types of land irrigation. Overall, the primary contribution of this research will help the Bangong Co Lake area sustainably utilize water resources.