Sulfur Isotope and Stoichiometry–Based Source Identiﬁcation of Major Ions and Risk Assessment in Chishui River Basin, Southwest China

: Hydrochemistry and sulfur isotope ( δ 34 S–SO 42 − ) of Chishui River watershed in Southwest China were measured to identify the sources of riverine solutes, the potential impact of human activities, water quality, and health risk. The main ﬁndings indicated that the HCO 3 − (2.22 mmol/L) and Ca 2+ (1.54 mmol/L) were the major ions, with the cation order of Ca 2+ (71 ± 6%) > Mg 2+ (21 ± 6%) > Na + + K + (8 ± 3%) and the anion sequence of HCO 3 − (55 ± 9%) > SO 42 − (41 ± 9%) > Cl − (4 ± 3%). The riverine δ 34 S–SO 42 − values ﬂuctuated from − 7.79‰ to +22.13‰ (average +4.68‰). Overall, the water samples from Chishui River presented a hydrochemical type of Calcium–Bicarbonate. The stoichiometry and PCA analysis extracted three PCs that explained 79.67% of the total variances.


Introduction
Watershed-scale hydrochemical evolution and water environmental quality, the most important hot point for river water resource research, is the basis of effective hydrospheric environment planning and sustainable usage of freshwater resources, in particular, under the background of the imbalance between unevenly distributed water resources and water requirements. Both anthropogenic and natural processes can affect hydrochemistry and the major ion concentration level [1,2]. With the accelerated improvement of economy, the increased anthropogenic activities (e.g., industrial/domestic wastewater, agricultural emissions) have rapidly enhanced the riverine pollutants [3,4], which will further result in a series of environmental issues and risks [5,6], such as the destruction of soil aggregate structure by polluted river water irrigation, and the health risks via drinking water [7][8][9].
The major ions are the most important part of river dissolved loads. Exploring the origins of these ions could greatly benefit for the understanding of watershed-scale geochemical dynamics (e.g., weathering processes) and distribution/transformation regulation of pollutants. Moreover, the ion source identification is also benefit for the environmental supervision of the government [10]. Generally, in addition to regularly direct monitoring of cation and anion concentrations, statistic approaches, ion ratios, and isotopic methods are widely applied in source identification [11]. Statistical methods, such as principal component analysis and correlation analysis could explore the common sources of riverine solutes due to similar physicochemical characteristics and potential co-origins of some ionic species [12]. The potential dilution effect can be avoided by the fluvial ion ratios, reflecting the mixing processes of the sources [13]. The isotopic methods, such as sulfur (S) isotope of riverine SO 4 2− , could distinguish different sources via specific isotope compositions of these end-members, based on the inheritance of S isotopes of riverine sulfate from both natural and anthropogenic sources [2,14].
In the karst area of southwest China, the typical carbonate geomorphology is developed, such as the sinkhole, depression, cockpit, and cone, which further result in the barren and thin soil [15,16]. Thus, the karstic ecosystem, in particular river system is extremely sensitive and vulnerable due to the unique geological conditions and strong karstification [17,18]. In order to obtain more knowledge of karst fluvial hydrochemistry, this study presents the detailed investigation of river water stoichiometry and S isotope of Chishui River (a tributary of upper Yangtze River) watershed in southwest China karst region. The main aims are to: (i) Carify the ionic compositions and S isotope compositions of river water, (ii) identify the source of major ions, and (iii) explore the water quality as well as the potential irrigation and health risks.

Study Region
The Chishui River is one of the first-level tributaries in upper Yangtze River that has not been dammed in the main channel. With a total length of 445 km for main stream, the Chishui River begins in the Zhenxiong County, Yunnan Province and flows from southwest to northeast through Yunnan-Guizhou Plateau and Sichuan Basin in southwestern China, and finally enters into the Yangtze River in Hejiang County, Sichuan Province (Figure 1) [7]. In Chishui River watershed (CRW, 27 • 15 ~28 • 50 N, 104 • 44 ~107 • 1 E), there is a wide catchment area (~1.89 × 10 4 km 2 ) with various tributaries influenced by different degrees of human activities. As shown in Figure 1a, the lithology exposed in CRW is mainly composed of carbonate sedimentary rocks (~44.6%), siliciclastic sedimentary rocks (~24.5%), mixed sedimentary rocks (~30.1%), and basic volcanic rocks (~0.7%). In more detail, the strata distributed in upper stream are mainly dolomite of Dengying Formation of Ediacaran with the minerals resources of phosphorite, barite, and fluorite. The middle reaches are mainly composed of limestone and dolomite of Cambrian, Ordovician, Silurian, Permian, Triassic, and Jurassic, followed by mudstone, sand shale and coal-bearing rock group. The middle reaches also developed the minerals resources of coal, pyrite (iron sulfide), gypsum, and barite with the same geological age. The lower stream mainly distributed Jurassic-Cretaceous siltstone and mudstone with sporadic oil shale. The CRW covers various land use and a large elevation range (205 to 2237 m). The land use of CRW mainly includes water area, urban area, grass land, unused land, cropland, and forest land (Figure 1b). Forest land (~73.4%) and cropland (~19.6%) are the dominant land uses in CRW. The CRW is affected by the sub-tropical monsoon climate, with the average annual air temperature of 11~13 • C. The annual rainfall ranges from 800 to 1200 mm. The agriculture is distributed in upper-middle stream with a vulnerable eco-environment, while the industries are mainly located in the middle stream. The water quality of CWR is of great importance due the significant supporting role of water resources in the watershed.

Sample Collection and Analyses
A systematic watershed survey was conducted from November to December in 2012 (dry season was selected to avoid the dilution effect and biologic effect as much as possible). According to the natural features (landuse and lithology) and the population distribution, 38 sample sites were chosen in the mainstream and tributaries of CRW (2 sites were chosen in Changjiang River, Figure 1). In total, 38 river water samples were obtained and further filtered (0.45-µm membrane, Millipore) and saved in the clean polyethylene sample bottles. Each sample was divided into three parts for the measurement of sulfur isotope and the major ions. The storage conditions of all samples were lightless and refrigerated (4 • C). The river water parameters (pH, T, DO, and EC) were measured in the field using the multi-parameter meter (Multi Line 3320, WTW, Munich, Bavaria, Germany). HCO 3 − concentrations were detected by HCl-titrated method. The Na + , K + , Mg 2+ , Ca 2+ , F − , Cl − , NO 3 − , and SO 4 2− concentrations were measured by ion chromatograph (ICS-1100, Thermo Fisher, Waltham, MA, USA) at the Key Laboratory of Karst Geological Resources and Environment, Ministry of Education, Guizhou University. The analysis was conducted with replicate samples and procedural blanks to maintain the accuracy of measurement. The measured result of replicate samples suggested acceptable repeatability for all ions (relative standard deviations were within ±5%), and the procedural blanks for all ions were generally lower than the detection limit.
The sulfur isotope of riverine SO 4 2− was detected according to previous studies [19]. Briefly, 10% BaCl 2 solution were added to water samples to convert the dissolved SO 4 2− into BaSO 4 precipitation, and the mixture was filtered by 0.45-µm membrane filters after 48 h. Then, the BaSO 4 precipitation on membranes was further calcined (800 • C, 40 min) to gain the BaSO 4 solid. The S isotope compositions (in δ notation relative to Vienna-Canyon Diablo Troilite, VCDT; Equation (1)) were detected by Elemental Analyzer-IsoPrime MS (IsoPrime, GV Instruments, Manchester, UK) at the Institute of Geochemistry, Chinese Academy of Sciences [20].
The standard reference materials (NBS 127) and internal laboratory standard reference materials (pre-calibrated Lab-SO1, Lab-SO 4 , Lab-Sigma) for S isotope were also detected to ensure the analysis accuracy, which were within the recommended value (2SD smaller than 0.02‰).

Assessment Method
Chishui River is the most important water source of both drinking and agricultural irrigation water within the watershed, which is necessary to assess the suitability for drinking and irrigation purpose. The riverine parameters (pH) and ion concentrations of Chishui River water were appraised by Chinese (GB 5749-2006) and WHO drinking water quality guidelines. For the irrigation water quality, the soil quality attributes can be influenced by the different salinity and alkalinity level of irrigation water, and further changed the yields from farmland. To gain an overall evaluation of irrigation water salinity and alkalinity hazard, the commonly-applied indicators such as sodium adsorption ratio (SAR), soluble sodium percentage (Na%), and residual sodium carbonate (RSC) were calculated based on the ion equivalent concentrations (meq/L) of river water as follows [11]: 5 (2) Na% = 100% × Na + /(Na + + K + + Mg 2+ + Ca 2+ ) Regarding to the human health risk, the health threat could be occurred via ingestion and dermal exposure to river water with high-concentration ions [21]. In comparison, the health risk through ingestion intake is the most noteworthy. To assess this threat (noncarcinogenic health risk), the Hazard quotient (HQ) suggested by the U.S. Environmental Protection Agency was calculated for F − , NO 3 − , and NH 4 + [11]: where ADD ingestion , C, IR, EF, ED, BW, AT, and RfD represents the daily ingestion intake doses, ion concentrations (mg/L), daily ingestion rate (0.6 L/day for children, 1.0 L/day for adults), exposure frequency (365 days/year), exposure duration (12 and 25 years for children and adults), body weight (16 and 56 kg for children and adults), average time for non-carcinogenic effect (4380 and 10950 days for children and adults), reference dose of different ions (0.04, 1.6, and 0.97 ppm/day for F − , NO 3 − , and NH 4 + ), respectively. If the HQ of each ion (or total hazard quotient, HQ t ) is >1, the local people will be threatened by the corresponding ion-derived risk.

Software for Data Analyses
For the statistical analyses of riverine ion concentrations and other parameters, the principal component analysis (PCA, a common method for exploring the potential origins of ions) and Piper diagram analysis were carried out by SPSS 21.0. All the data-based figures were illustrated by Origin 2018.

Overview of Hydrochemical Compositions
The hydrochemical compositions and parameters of all river water samples are presented in Table S1, with the statistical results, including ranges, and mean values of these data summarized in Table 1. The pH values ranged from 7.63 to 8.87 (mean value 8.35), presenting the slightly alkaline river waters, further indicated that the inorganic carbon species was dominated by bicarbonate (HCO 3 − ) based on the carbonate equilibrium, while CO 3 2− and dissolved CO 2 were negligible. The EC value averaged 403 µS/cm with a range of 128~738 µS/cm. The ionic balance was also applied for the river water samples, that is, the total cation concentration (TZ + = Na + + K + + Mg 2+ + Ca 2+ , in meq/L) and total anion concentration (TZ − = F − + Cl − + NO 3 − + SO 4 2− + HCO 3 − , in meq/L) of river water, which presented well ion balances (R 2 = 0.987, p < 0.01). The normalized ionic charge balance ([TZ + −TZ − ]/TZ − ) is smaller than 11% (Table S1), suggesting the potential effect of the organic anionic species [22] like oxalate. Note: a the data of concentration is from [23]; b the unit of related values in Chinese guideline and WHO guideline are converted to mmol/L.
The major chemical ions of the Chishui River water are presented in Figure 2. It is very distinct that the principal cation in Chishui River is Ca 2+ , which accounts for 71 ± 6% (mean ± SD) of the cations, followed by Mg 2+ (21 ± 6%) and Na + + K + (8 ± 3%). The anions sequence is HCO 3 . HCO 3 − is therefore the predominant anion in Chishui River, which is a widely confirmed product (as well as riverine Ca 2+ and Mg 2+ ) of CO 2 -associated weathering process without the influence of anthropogenic pollution [13,24]. The massive carbonate sedimentary rock distribution in such a karst landscape river watershed creates the advantageous conditions for the chemical weathering of Ca/Mg-contented rock. Although the lower reaches flow through the area where the siliciclastic sedimentary rock developed, the inheritance of river water will make the lower reach reveal the characteristics of dissolved loads from upper and middle streams. Generally, the weathering rate of carbonate rock is much higher than other rocks in the same situations [25,26], implying that the HCO 3 − , Ca 2+ , and Mg 2+ , derived from carbonate sedimentary rock weathering and exported into Chishui River in uppermiddle stream, are far more than the potential Na + and K + originated from the weathering of siliciclastic sedimentary rock in the lower stream. Therefore, the hydrochemical type of Chishui River water is a Ca-Mg-HCO 3 − type controlled by carbonate weathering, which is highly in agreement with the studies of the whole Changjiang River (Yangtze River) basin ( Figure 2).

Figure 2.
Piper diagrams representing the percentages of hydrochemical species in Chishui River (including 2 water samples collected in Changjiang River before and after the Chishui River drained into, the solid gray circle) and Changjiang River (data reported by previous studies [10,[27][28][29][30]).
The potential sources of riverine SO 4 2− mainly include weathering process, atmospheric precipitation, industrial/mining activities, and sewage inputs [10]. Chemical weathering related studies have found that the accelerated carbonate dissolution occurred with relatively high SO 4 2− concentrations (as an agency of weathering processes) in the karst region [31,32]. Thus, the considerable fraction of SO 4 2− (41 ± 9%, close to HCO 3 − ) is non-negligible.

Stoichiometry-Revealed Sources of Solutes
In hydrochemistry studies, the high Cl − concentration could be applied as a significant index of anthropogenic inputs (mainly in urban area) and domestic sewage [33]. Generally, the river nitrate originated from agricultural synthetic fertilizers exhibits a high NO 3 − concentration and a high value of NO 3 − /Cl − ratio, while the domestic sewage presents low NO 3 − /Cl − ratio and high concentration of Cl − due to highly organic matter [34]. As shown in Figure 3a, the high NO 3 -/Clratio and low level of Cl − concentration (range 0.04~0.55 meq/L, mean value 0.18 ± 0.14 meq/L) were observed in all water samples of Chishui River, suggesting the major contribution by agricultural input (fertilizer), which can be supported by the large distribution area of cropland (~19.6%) in the whole watershed [7]. By contrast, the potential Cl − inputs from the anthropogenic domestic sewage in urban region may be assuaged by various landscape setting and the well-mixed processes of river water, further resulting in a low concentration of riverine Cl − . Moreover, although previous work have revealed that rainwater NO 3 − could be a significant contributor (even up to 71%) of riverine NO 3 − during storm-frequency season [35], however, the river water NO 3 − concentration (0.19 ± 0.13 meq/L) in the study period (dry season) was much higher than the karst rainwater NO 3 − concentration (typically < 0.04 meq/L) [36][37][38]. Therefore, the contribution of wet deposition to Chishui River water NO 3 − is very limited.
In addition, the river water SO 4 2− /Ca 2+ ratios were much higher than the NO 3 − /Ca 2+ ratios (Figure 3b), indicating an additional potential impact of industrial activities and mining (e.g., pyrite in middle reaches) on fluvial solutes [39,40]. The weathering-derived major ions always present typical ratios, thus the ratios of Ca 2+ /Na + and HCO 3 − /Na + are the useful index to track the weathering origins of ions [13,19]. The Ca 2+ /Na + and HCO 3 − /Na + ratios of Chishui River water were scattered between the silicate and carbonate end-members and more inclined to the carbonate source (Figure 3c), implying the leading role of carbonate weathering products on riverine ions with a slightly mixing process of silicate weathering products. It is noteworthy that the two samples form Changjiang River were more trended to silicate end-members (Figure 3c) Figure 3d, most of samples distributed between the two types of acid-involved weathering results, implying that both H 2 SO 4 and H 2 CO 3 weathering processes occurred. The samples were relatively closer to the H 2 SO 4 -involved weathering results, indicating the importance of sulfuric acid weathering processes, which may notably accelerate the rock dissolution. The generally low Mg 2+ /Ca 2+ ratios (<0.5) of most river water samples further distinguished the contributions of different types of carbonate minerals, that is, the riverine Ca 2+ and Mg 2+ were primarily controlled by calcite dissolution (the dolomite dissolution dominated river water Mg 2+ /Ca 2+ ≈ 1) [26]. In contrast, if the extensive silicate weathering occurs, the [Ca 2+ + Mg 2+ ]/HCO 3 − ratios will decreased and the sample points in Figure 3d will move downwards systematically due to the simultaneously production of Na + , K + , HCO 3 − , and SO 4 2− in silicate weathering [11,41].
According to low concentrations of riverine K + (0.05 mmol/L), F − (0.01 mmol/L), and NH 4 + (0.003 mmol/L) [23], these three ions related ratios were not presented in Figure 3. Even so, previous studies have concluded that K + is contributed both by rock weathering (silicate) and human emissions [10,41,43]. Thus, here we hold the opinion that K + is mainly controlled by rock source in such a low riverine K + level (anthropogenic input is negligible) in Chishui River. Riverine F − is also primarily originated from bedrocks [44]. As for NH 4 + , the agriculture-related processes, mainly ammonia-contained fertilizer, can be an interpretation for the samples with relatively high ammonia concentration [45].

PCA Analysis
The PCA was used to further explore the potential sources of riverine ions in Chishui River. Three principal components (PCs) with eigenvalues exceeding 1 were distinguished and illustrated in Figure 4, with more details presented in Table 2. These three PCs contained a total of 79.67% variance. The PC 1, PC 2, and PC 3 accounted for 33.8%, 23.9%, and 23.1% of total variances, respectively. Four ions (K + , Mg 2+ , F − , and HCO 3 − ) as well as EC presented notably positive loadings in PC 1, and the relatively strong loading of Ca 2+ (0.48) is also found in PC 1 (Figure 4 and Table 2). Based on the stoichiometry discussed in Section 3.2.1, here we infer that PC 1 was mainly contributed by natural rock weathering inputs in the Chishui River watershed due to the high loadings of feature ions of rock weathering (primarily carbonate) [46]. Na + and Cl − presented clear positive loadings in PC 2, but these two ions were weakly loaded in PC 1 (natural source) and PC 3 (SO 4 2− -dominated human inputs). Combined with the negligible contribution of silicate weathering and the limited anthropogenic Cl − inputs mentioned before, the atmospheric input/deposition could be a reasonable contributor, which is also supported by previous studies [42]. The strongest loadings of SO 4 2− and NO 3 − were observed in PC 3 (Figure 4), considering as typical anthropogenic inputs. It is noteworthy that the PC 3 also presented a significantly positive loading of Ca 2+ , further supporting the H 2 SO 4 -involved weathering processes. Additionally, the same positive loading of EC (reflect total riverine ion content) was found in PC 1 and PC 3 ( Table 2), indicating that the contribution of natural and anthropogenic sources to total fluvial ions was comparable.

Sulfur Isotope-Based Source Identification of Riverine Sulfate
Although limited impact of fluvial SO 4 2− for drinking and irrigation purpose is reported, the potential accelerated weathering may result in more solutes input to river system, which further causes other environmental issues (e.g., carbon source/sink variations) [11]. Therefore, to further constrain the SO 4 2− sources, the effective indicator δ 34 S-SO 4 2− was applied. Generally, the riverine SO 4 2− is derived from four sources: atmospheric inputs (mainly rainwater); sulfide oxidation (e.g., pyrite); human inputs (e.g., agricultural fertilizers, sewage, and industrial effluent); gypsum dissolution (a kind of evaporitic rock) [12,14,47]. These four sources can be well constrained by combining the distinguished δ 34 S-SO 4 2− values and SO 4 2− concentrations ( Figure 5). The sulfur isotope compositions of atmospheric depositions in south China were well explored in previous studies, such as the 950-day based study in Wuhan urban area (δ 34 S-SO 4 2− = +4.5 ± 1.3‰) [48], rainwater in Sichuan Basin (δ 34 S-SO 4 2− = +1.5+ 6.0‰) [49], and the rainwater sulfur isotope study carried out in a karst agricultural catchment (δ 34 S-SO 4 2− = +1.3 ± 6.2‰) [50]. Based on the similar geographical location, climatic condition, lithology and agriculture-dominated industrial structure, the rainwater δ 34 S-SO 4 2− values (+1.3 ± 6.2‰) of the karst agricultural catchment can be regarded as an atmospheric input end-member. The sulfide oxidation derived SO 4 2− usually inherit the sulfur isotope compositions of the initial sulfide due to limited isotope fractionation during oxidation processes [19]. The average δ 34 S-SO 4 2− value of pyrites in the study area (Sichuan, Guizhou, and Yunnan provinces) was −7.8‰ [51], which is within the typical range of sulfide (−13.0~−1.1‰) [19]. Moreover, although the sulfide oxidated SO 4 2− concentrations were rarely reported in the study area, previous work have shown that its concentration should be exceed 0.1 mol/L [47], even up to 2 mol/L, that is, 1/[SO 4 2− ] = 0.5~10. The anthropogenic inputs, mainly including agricultural fertilizers, domestic and industrial sewage, present a wide range of sulfur isotope compositions. Therefore, the typical δ 34 S-SO 4 2− values of water in farmland and livestock farm areas (+2.8~+13.9‰) were regarded as the end-member of human inputs [12,52], and the potentially high SO 4 2− concentration (1/[SO 4 2− ] below~17.5) [19] was selected to distinguish the overlapped part of the δ 34 S-SO 4 2− values with atmospheric input. Moreover, the gypsum dissolution source presented the highest δ 34 S-SO 4 2− value (up to~+30‰) [2] and high SO 4 2− concentration. As shown in Figure 5, the ternary diagram presented the relative contributions of different origins to fluvial SO 4 2− via the relationship between δ 34 S-SO 4 2− values and 1/[SO 4 2− ], this is, riverine SO 4 2− were controlled by the sources of human inputs and sulfide oxidation. The middle-stream-distributed pyrite deposits and the relatively high population density (282 persons per km 2 ) and wide urbanization area (1.2%) [7,23,40] also underpinned these results. It is noteworthy that four samples were distributed over the range of human inputs and tended into the gypsum dissolution sources ( Figure 5). This suggests the potential influence of gypsum dissolution on river water in some sampling sites to a certain extent, but it was negligible on a whole watershed scale, which was also supported by previous work [19]. Although the range of δ 34 S-SO 4 2− value of atmospheric inputs was partially overlapped with other sources, the much lower SO 4 2− concentrations of atmospheric source were significantly separated from human input and sulfide oxidation sources ( Figure 5), confirming again that the contribution of atmospheric inputs to riverine sulfate was very limited. This can also be supported by the average investigation results of global rivers (atmospheric input contributed to~3% sulfate) [13]. In comparison, the water samples obtained in larger-scale Zhujiang River watershed (originated in karst region of Guizhou province) in Figure 5 presented a distribution closer to the sulfide oxidation end-member, with some samples trended to the atmospheric source [11], indicating a complicated mixing processes and differential contributions of sulfate sources in different scaled catchments.

Irrigation and Guideline-Based Water Quality
Chishui River servers as primary water resources of the agriculture, industry, and local residents (population > 5 million). As summarized in Table 1, according to the Chinese and WHO drinking water quality guidelines, the desirable pH value for drinking is recommended between 6.5~8.5, thus most of the water samples (81.6%) had suitable pH values (Table S1). The high pH values (>8.5) of partial samples were mainly caused by carbonate weathering. Most of the guidelines-covered ions, including F − , Cl − , NO 3 − , and NH 4 + of all river water samples were below the recommended limits. However, although the average value of riverine SO 4 2− (0.85 mmol/L) was lower than that of guidelines limits (2.60 mmol/L), the maximum value of riverine SO 4 2− (2.64 mmol/L) exceeded the permissible limits, which needs more attention. The previous study on trace metal also assessed the 19 metals-based water quality index (WQI) of Chishui River, which showed the assessment results of excellent water quality (WQI < 50, only from the perspective of metals) without extensive trace metal contamination [7].
The Na% and SAR indicators can reflect the Na hazard to agricultural land by irrigation-influenced soil aggregate [19]. According to the calculated SAR and Na% values, all water samples in the Chishui River could be regarded as excellent/good quality (SAR < 0.61; Na% < 19.0%) with no samples out of the desirable limits (SAR = 1; Na% = 30.0%). For the residual sodium carbonate (RSC), there was no water samples present a RSC value exceeded 1.25 (Table S1), indicating the well suitability of irrigation water. Moreover, the United States Salinity Labortory (USSL) diagram and Wilcox diagram were plotted based on the EC, SAR, and Na% values (Figure 6), and all the Chishui River water samples were scattered in the C1S1 and C2S1 zone of USSL diagram and the "Excellent to Good' area of Wilcox diagram. Overall, the Chishui River water is desirable for agricultural irrigation and will not bring about the soil hazard. However, the continuous long-term monitoring is still a worthy mission due to the enhanced human activities.

Health Risk Assessment
Among the major ions in natural water, excessive exposure (ingestion intake) of F − , NO 3 − , and NH 4 + can cause typical non-carcinogenic hazards, while SO 4 2− does not result in health issue [11]. Thus, riverine F − , NO 3 − , and NH 4 + were involved in the health risk assessment in this study. The average concentration-based HQ values were calculated and presented in Figure 7. The HQ values of three assessed ions were in the order of HQ-NO 3 − (0.28 ± 0.19) > HQ-F − (0.13 ± 0.06) > HQ-NH 4 + (0.002 ± 0.000) for children and HQ-NO 3 − (0.11 ± 0.07) > HQ-F − (0.05 ± 0.02) > HQ-NH 4 + (0.001 ± 0.000) for adult. The results were similar to other researches [44], that is, NH 4 + occupied a tiny part of the total HQ, while NO 3 − and F − were the majority. For both children and adults, HQs of each ion and total HQ were all lower than the hazard level (<1), indicating a lower health risk of assessed ions. It is noteworthy that the HQ value for children was always higher than that for adults, implying that children were facing higher risks from riverine ions compared to adults. In addition, the potential harmful health effects can also occur if the HQ value is greater than 0.1 for children [21]. Therefore, the riverine ion-related health risk is not completely negligible, particularly for children.

Conclusions
The present study aimed to investigate hydrochemical compositions and sulfur isotope compositions of Chishui River water, southwest of China. Based on the isotope tracer, stoichiometry and PCA method, the sources of major ions were identified. The rock (mainly carbonate) weathering inputs were the primary sources of K + , Mg 2+ , Ca 2+ , F − , HCO 3 − , and atmospheric contribution was the source of Na + and Cl − , while the anthropogenic input was responsible for SO 4 2− and NO 3 − . Sulfur isotope compositions further reflected the relative contribution of different sources, that is, the contribution of human inputs (main) and sulfide oxidation controlled the riverine SO 4 2− , but the input of atmospheric deposition was limited. Both the water quality and hazard quotient assessment produced good results, indicating the suitable aims of river water for drinking and irrigation. This work and the possible continuing study would be benefit to the prevention of watershed water environment and further help planning the sustainable water resources in Chishui River watershed.