Next Article in Journal
Adsorption of Sunscreen Compounds from Wastewater Using Commercial Activated Carbon: Detailed Kinetic and Thermodynamic Analyses
Previous Article in Journal
Three-Dimensional Static and Dynamic Analyses of an Embedded Concrete-Face Rockfill Dam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Chemical Weathering of Rocks and Its Carbon Sink Effect in the Naqu River Basin of the Nujiang River Source Area, Southwest China

1
School of Water Conservancy, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
2
Henan Key Laboratory of Water Resources Conservation and Intensive Utilization in the Yellow River Basin, Zhengzhou 450046, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(23), 4191; https://doi.org/10.3390/w15234191
Submission received: 3 November 2023 / Revised: 29 November 2023 / Accepted: 30 November 2023 / Published: 4 December 2023

Abstract

:
Carbon plays an important role in global climate change. The mechanisms of carbon sources and carbon sinks have also received wide attention from society, and the physical and chemical characteristics of riverine ions can reflect the chemical weathering of rocks and carbon sink capacity of river basins. Based on the data on river, rainwater, and rock samples from 2019, this study used various methods, such as ion ratio diagrams and ternary diagrams, to analyze the chemical characteristics of water; the chemical weathering and carbon sink effects of rocks were also calculated while assuming three scenarios based on the main sources of ions in the Naqu River. The results showed that for the whole catchment, the main ion sources in the river were: carbonate rock chemical weathering > silicate rock chemical weathering > evaporite dissolution > atmospheric precipitation input. According to the calculations, in the three scenarios, the carbonate weathering rates were 16.84, 11.32, and 14.08 t/km2/yr, and the carbon sink capacities were 66.47, 121.13, and 93.80 mol/km2/yr, respectively; the evaporite weathering rates were 2.20, 9.63, and 5.92 t/km2/yr, respectively. The silicate chemical weathering rate and carbon sink capacity did not change significantly in either scenario, with 6.82 t/km2/yr and 248.6 mol/km2/yr, respectively. This study quantified the ion sources in the Naqu River basin and accurately analyzed their chemical genesis, which helps in understanding the role of the rivers of the Qinghai–Tibet Plateau in the global carbon cycle and global climate change, in addition to providing a reference for the scientific development of the Nujing River.

1. Introduction

Chemical weathering plays an important role in the mechanism of carbon sinks. River water can be used as a carrier of chemical weathering products, such as solutes and suspended particles [1]. Solutes in rivers may originate from various factors, including rock weathering, soil erosion, atmospheric precipitation inputs, and anthropogenic influences [2]. Analyzing the chemical characteristics of river water can reflect the chemical weathering capacity and carbon sink capacity of a basin to a certain extent [3,4].
The chemical weathering of rocks is impacted by hydrological processes, and it produces carbon monoxide and carbon dioxide; the ensuing carbon sinks have the potential to address a portion of the global carbon imbalance. The escalation of carbon dioxide concentrations resulting from rock weathering contributes to global warming, and heightened temperatures further stimulate rock weathering and the absorption of more carbon dioxide [5]. The succession of climatic extremes and the phenomena of the drying of land stemming from global warming unavoidably exert both positive or negative impacts on the weathering processes of rocks [6]. There is international recognition of the existence of missing carbon sinks in the global carbon cycle, some of which are stored in solid silicate minerals resulting from silicate weathering. The carbon sinks formed through the chemical weathering of carbonates are also considered as parts of these missing sinks [7]. Beyond accelerating the rate of chemical weathering in rocks, warming has the additional effect of reshaping the total ion content produced by the chemical weathering of rocks in river water [2]. Carbonates demonstrate a quicker pace of weathering kinetics when contrasted with silicates, and in virtually all watersheds, the predominant sources of dissolved inorganic carbon are carbonate rocks [8]. Through photosynthesis, dissolved inorganic carbon undergoes a process with aquatic organisms to generate organic carbon, enhancing the carbon sink effect by converting inorganic carbon into organic forms [9].
At present, there is abundant research about chemical weathering and carbon sequestration in basins. Due to factors such as its high altitude and harsh climatic conditions, it is difficult to study the Qinghai–Tibet Plateau. Yu et al. [10] meticulously gauged the vigor of chemical weathering, the velocity of CO2 consumption, and the nuanced influences at play in the upstream and midstream sections of the Yarlung Tsangpo. Their discerning results unveiled a narrative of subdued chemical weathering intensity and restrained CO2 sequestration, standing in stark contrast to the grandiose demeanor of most expansive rivers coursing down the Tibetan Plateau. Yu et al. [11,12,13] studied the central basin of the Qinghai–Tibet Plateau and found that the chemical compositions of the rivers in the Niyaqu and Qugaqie basins mainly included Ca2+ and HCO3−. Seasonal changes had little effect on the solute sources in the Niyaqu basin, but they had a significant impact on the Qugaqie basin. In addition, glaciation enhanced the chemical weathering of the central Qinghai–Tibet Plateau. Zhang et al. [14] studied the oxidative weathering process of pyrite in the eastern part of the Qinghai–Tibet Plateau, and the results showed that in 2012–2013, the CO2 emissions due to carbonate weathering coupled with sulfuric acid exceeded the CO2 consumption by the silicate weathering coupled with carbonic acid, but it did not exceed that in 2013–2014. Kang et al. found that the weathering process of rocks in an alpine basin was predominantly shaped by the chemical weathering of carbonate rocks, with additional influences from climatic factors such as temperature and precipitation [15]. The study of chemical weathering in large watersheds can allow one to understand the overall situation of a region, but the study of small watersheds can allow one to accurately analyze the chemical causes and influencing factors, so this is also very necessary. At present, some studies have been carried out on the source region of the Nujiang River on the Qinghai–Tibet Plateau. Huang et al. [4] found that the tributary of the Salween River had a pH between 8.1 and 8.7, the EC was between 33.7 and 54.0 mS/m, and the conductivity of Yuqu, a tributary of the Salween River, was relatively low. Wu et al. [16] found that, compared with the chemical weathering rates of other rocks in Nujiang, carbonate rock had the highest CO2 consumption flux due to chemical weathering, reaching 5.9 × 105 mol/km2/yr, and the total consumption was 6.5 × 1010 mol/yr in their study. By analyzing chemical weathering and carbon sink effects, Tao et al. [17] found that the hydrochemical type of the Nujiang River was mainly that of Ca-HCO3, and the dissolution of evaporite made little contribution to the Nujiang River. The chemical weathering rate of carbonate rock was higher (33.54 t/km2/yr), and its CO2 consumption flux was 2.22 × 105 mol/km2/yr in the Nujiang River. Previous research mainly focused on the middle and lower reaches of the Nujiang River, and studies of the carbon sink of the Naqu River watershed, which is the source area of the Nujiang River, are relatively scarce.
Based on hydrochemical data from 2016 to 2018, the predominance of rock weathering in the Naqu River basin was found in our previous study through cluster analysis, triangular plots, correlation analysis, and other methods [18]. But our previous study did not quantitatively analyze the ionic sources in the Naqu River. Therefore, through an examination of the hydrochemical traits within the Naqu River basin, this study postulates three scenarios grounded in the principal origins of ions. Subsequently, both the rate of chemical weathering of rocks and the carbon sink capacity are computed. The overarching objective is to delve into the hydrochemical attributes of the Naqu River, unraveling insights into the intricate processes of rock weathering and the carbon sequestration potential within a prototypical plateau river basin. Such investigations bear profound significance, as they contribute to a comprehensive understanding of the roles played by rivers on the Qinghai–Tibetan Plateau within the larger context of the global carbon cycle and climate dynamics.

2. Materials and Methods

2.1. Study Area

The study area was located in Naqu City, North Tibet Plateau, Tibet Autonomous Region, Southwest China [19], as shown in Figure 1. The Naqu River is the source of the Nujiang River, the geographical position of which is from 30°54′ to 32°43′ N and from 91°12′ to 92°54′ E [18]. The Naqu River, which is 460 km long, covers 16,350 km2 of total drainage area [20]. The Naqu River is composed of four main tributaries, which are, Chengqu, Mugequ, Gongqu, and Luoqu from north to south [21]. Most of the Naqu River basin is located in a plateau subfrigid monsoon semi-humid climate zone, with intensive solar radiation and thin air. According to the analysis of the precipitation data, it can be seen that the average precipitation in the Naqu River basin is 478 mm, indicating the total rainfall [22]. The rainy season occurs from July to August and gradually decreases from east to west [23]. From 1980 to 2000, it was shown that the mean annual evaporation was approximately 1058 mm [24]. In addition, the average temperature from 2015 to 2019 ranged from 0.28 to 1.15 °C.
The upper reaches of the Naqu River basin are relatively gentle, and the terrain is flat and open; the lower reaches are relatively steep, with valley-shaped terrain [16]. The carboniferous to Permian, Jurassic, Tertiary, and Quaternary rocks all have different degrees of exposure in the basin. The main exposed rocks are slate, shale, sandstone, graywacke, granite and other magmatic rocks, a small number of metamorphic rocks, and a fourth system of clay, sand, pebbles, gravel, and so on. Both silicate and carbonate rocks are distributed in the upper and lower reaches of the Naqu River, and the middle reaches of the Naqu River are dominated by silicate rocks. The tributaries of the Sangqu, Basuoqu, Mumuqu, Zongqingqu Mugequ, and Gongqu all flow through carbonate-dominated areas [23,25].
The Naqu River basin is an important region of Naqu City, which has a total of 13 townships. The population density is low, and the economic development is dominated by animal husbandry [26]. The main plants in the domain are temperate subtropical alpine meadows, and the main vegetation is alpine tarragon, alpine dwarf tarragon, alpine mat vegetation, and alpine sparse vegetation.

2.2. Sample Collection and Analysis Methods

2.2.1. Sample Collection

The sample for this study was collected in August 2019, and the sampling points are shown in Figure 2. The sampling frequency for atmospheric precipitation was once per month, and 200 mL polyester plastic bottles were used as sampling bottles for the water quality samples. When collecting river samples, the sample bottle was rinsed with river water 2 to 3 times before sampling, and then it was filled with river samples, sealed, and marked. Rock samples were collected near the river sampling sites. A Horiba U50 multi-parameter water quality analyzer (HORIBA Ltd., Kyoto, Japan) was used on site to measure the temperature, dissolved oxygen (DO), electrical conductivity (EC), oxidation–reduction potential (ORP), and total dissolved solids (TDS) of the river. These indicators are commonly used for water quality monitoring. Among them, DO is an important limiting factor for aquatic organisms, as it depends on the concentration of water-soluble salt, biological action, the water temperature, and the partial pressure of atmospheric gases on the water surface. EC is used to determine water quality and detect the presence of chemical substances. The higher the concentration of dissolved minerals in water, the higher the EC value. ORP indicates the cleanliness of water by measuring the oxidation of water, and the lower the ORP value is, the higher the likelihood of contamination will be. All river water samples and rock samples were determined in the laboratory of the First Geological Environmental Investigation Institute of Henan Bureau of Geology and Mineral Exploration and Development. Inductively coupled plasma atomic emission spectroscopy (ICP-AES) (ICAP 6300, Thermo Fisher Scientific, Waltham, MA, USA) was used to measure the concentrations of K+, Na+, Ca2+, and Mg2+ in the river, and the concentrations of Cl, HCO3−, and CO32− were measured through titration with a 50 mL acid burette [27]. The concentrations of SO42−, NO3, and SiO2 in the river were measured with a T9CS dual-beam UV-VIS spectrophotometer (Beijing General Instruments Co., Ltd., Beijing, China) [28], and an inductively coupled plasma mass spectrometer (ICP-MS) (Thermo Fisher Scientific, USA) was used to measure the Sr concentration in the river water samples [29].
If the normalized inorganic charge balance (NICB; NICB = (TZ+ − TZ)/TZ+) was close to 0, this indicated that the positive and negative charges in the river water were closer to the conservation of charge. If the NICB was large, this indicated that there were ions in river that had not been detected. Figure 3 depicts the relationship between the EC and TZ+ in the Naqu River basin. The results show that EC was significantly positive correlated with TZ+ in 2019. This good correlation proved that the hydrochemical data obtained in the field and laboratory measurements in this study were reliable because the sum of the cations in the solution was proportional to the conductivity.
Based on the sample data, a Piper diagram was used in this study to determine the hydrochemical types in the river.

2.2.2. Qualitative Methods for Analyzing the Sources of River Solutes

The combined chemical weathering of rocks is one of the main sources of river water ions. According to differences in their chemical properties, weathered surface rocks can be divided into three major categories: carbonate, silicate, and evaporite [30]. The chemical weathering of rocks requires the participation of acidic media. The reaction equations related to the chemical weathering of rocks are as follows.
Hydrolysis of carbonate rocks:
( C a 1 x M g x ) C O 3 + C O 2 + H 2 O 1 x C a 2 + + x M g 2 + + 2 H C O 3
Oxidation of pyrite:
F e S 2 + 7 O 2 + 2 H 2 O 2 F e S O 4 + 2 H 2 S O 4
Reaction between sulfuric acid formed by the oxidation weathering of pyrite and carbonate rocks:
2 ( C a 1 x M g x ) C O 3 + H 2 S O 4 2 1 x C a 2 + + 2 x M g 2 + + 2 H C O 3 + S O 4 2
Weathering reaction of anhydrite and halite:
C a S O 4 = C a 2 + + S O 4 2
N a C l = N a + + C l
The classical weathering reactions of silicate mineral hydrolysis:
2 N a A l S i 3 O 8 + 2 C O 2 + 11 H 2 O A l 2 S i 2 O 5 O H 4 + 2 H C O 3 + 2 N a + + 4 H 4 S i O 4
2 K A l S i 3 O 8 + 2 C O 2 + 6 H 2 O A l 2 S i 4 O 10 O H 2 + 2 H C O 3 + 2 K + + 2 H 4 S i O 4
C a A l 2 S i 2 O 8 + 2 C O 2 + 3 H 2 O A l 2 S i 2 O 5 O H 4 + 2 H C O 3 + C a 2 +
M g 2 S i O 4 + 4 C O 2 + 4 H 2 O 4 H C O 3 + 2 M g 2 + + H 4 S i O 4

2.2.3. Quantification Methods for Analyzing the Sources of River Solutes

According to the data on the concentrations of dissolved ions in river water and atmospheric precipitation, the forward model method was used to analyze the sources of ions [31]. The equation is expressed as follows:
X r i v = X c y c l + X e v a p + X c a r b + X s i l + X a n t h
In the formula, [X]riv represents the dissolved ions in river water, [X]cycl represents the input of atmospheric precipitation, [X]evap represents the contribution of dissolved minerals in evaporite rocks, [X]carb represents the contribution of the chemical weathering of carbonate rocks, [X]sil represents the input of chemical weathering of silicate minerals in the basin, and [X]anth represents the contribution of human activities in the basin to dissolved ions in river water; the unit of ion X is mmol/L. The application of the forward evolution method was based on the following assumptions: (1) The Cl in the river only came from evaporite and precipitation; (2) Na+ and K+ were not affected by carbonate rocks. In this study, the rates of contribution of precipitation, carbonate, silicate, and evaporite to the cations dissolved in river water were quantitatively analyzed based on forward modeling.
Atmospheric precipitation is one of the main sources of river ions. The calculation methods for the input of atmospheric precipitation mainly include the water balance method, rainwater ratio method, and mixing method [30]. According to the available data, the rainwater ratio method was adopted in this study to correct the input ions of atmospheric precipitation. Because the minimum Cl concentration in the collected river samples was 12.96 μmol/L [18], this value was selected as the Cl concentration for the input of atmospheric precipitation in the Naqu River basin. The equation is as follows [32]:
X c y c l = C l c y c l × X / C l a t m
In the formula, X = SO42−, Na+, Ca2+, Mg2, and K+ (mmol/L), and (X/Cl)atm is the ratio of an ion X in rainwater to the concentration of Cl.
The rate of contribution of atmospheric precipitation input to the river water can be calculated with the following formula [32]:
T Z c y c l + T Z + = N a c y c l + K c y c l + 2 × C a c y c l + 2 × M g c y c l T Z +
T Z + = N a r i v + K r i v + 2 C a r i v + 2 M g r i v
In the formula, the unit of ion X is mmol/L.
The input from the chemical weathering of rocks in the basin is generally composed of three parts: evaporite chemical weathering, silicate chemical weathering, and carbonate chemical weathering.
The output ion concentration and ion contribution rate of the chemical weathering of the evaporite in the Naqu River basin can be calculated according to the following formula:
C l e v a p = C l r i v C l c y c l
N a e v a p = C l e v a p  
[ S O 4 ] e v a p = α e v a p × [ S O 4 ] r i v [ S O 4 ] c y c l
[ S O 4 ] s u l f u r i c   a c i d = 1 α e v a p × [ S O 4 ] r i v [ S O 4 ] c y c l
C a e v a p = [ S O 4 ] e v a p  
T Z e v a p + T Z + = C l e v a p + 2 × [ S O 4 ] e v a p T Z +
In the formula, αevap is the proportion of the input of SO42− from gypsum to the remaining SO42− ions after deducting the inputs of human activities and atmospheric precipitation, while [SO4]sulfuric acid is the input of SO42− into the river from the H2SO4 erosion of carbonate rocks. Then, according to the relevant research, it was assumed in this study that, after deducting the inputs of human activities and atmospheric precipitation, the calculation of the proportion of SO42− generated through the hydrolysis of gypsum and the action of sulfuric acid in the chemical weathering of rocks was conducted in three different scenarios.
Scenario 1: The SO42− ions in the river water of the Naqu River basin were all derived from sulfuric acid after deducting the input from atmospheric precipitation; αevap = 0.
Scenario 2: The SO42− ions in the river water in the Naqu River Basin were all derived from gypsum after deducting the input from atmospheric precipitation; αevap = 1.
Scenario 3: Gypsum dissolution and H2SO4 erosion of carbonate rocks contributed in a ratio of 1:1 with respect to SO42− in the water of the Naqu River.
In all three of these scenarios, the rate of chemical weathering of rocks was calculated, and the carbon sink effect of the watershed was analyzed.
According to the forward method, after deducting the inputs of atmospheric precipitation and evaporite rock, the Na+ and K+ in the river mainly came from the chemical weathering of silicate rock. Therefore, the input of Na+ and K+ contents due to the chemical weathering of silicate rocks can be calculated according to the following formula [32]:
N a s i l = N a r i v N a c y c l N a e v a p
K s i l = K r i v K c y c l
According to relevant studies, the concentrations of the input of Ca2+ and Mg2+ due to the chemical weathering of silicate rocks (mainly anorthite) in rivers can be calculated according to the following formula [33]:
C a s i l = N a s i l × C a / N a s i l
M g s i l = N a s i l × M g / N a s i l
The values of (Ca/Na)sil and (Mg/Na)sil were approximately 0.08 and 0.03 in river sand samples from the Naqu river. The molar concentration of the input of Si into the river due to the chemical weathering of silicate rock and the rate of contribution of the ion input into the river can be calculated according to the following formula:
S i s i l = S i r i v
T Z s i l + T Z + = N a s i l + K s i l + 2 × C a s i l + 2 × M g s i l T Z +
The concentration of the input of Ca2+ and Mg2+ due to the chemical weathering of carbonate rocks can be calculated according to the following formula:
T Z c a r b + T Z + = 2 × C a r i v C a c y c l C a e v a p C a s i l + 2 × M g r i v M g c y c l M g s i l T Z +

2.2.4. Methods for Analyzing Chemical Weathering and Carbon Sink Effects in Rocks

The chemical weathering rates of rock minerals (t/km2/yr) can be measured by using the concentration of dissolved ions that are fed into a river by the weathering products of rock minerals per unit area [34,35]. The contribution of HCO3 to the chemical weathering of silicate and carbonate rocks can be calculated according to the following formula [36]:
[ H C O 3 ] s i l = K s i l + N a s i l + 2 C a s i l + 2 M g s i l
[ H C O 3 ] c a r b = [ H C O 3 ] r i v [ H C O 3 ] s i l
In this formula, [HCO3]sil and [HCO3]carb are HCO3 ions produced by the chemical weathering of silicate rock and carbonate rock, respectively. The units of [HCO3]sil and [HCO3]carb are both mmol/L, and the other symbols have the same meaning as those mentioned above.
The HCO3 ions produced through the H2SO4 erosion of carbonate rocks can be calculated as follows [36]:
[ H C O 3 ] s u l f u r i c   a c i d = 2 × [ S O 4 ] r i v [ S O 4 ] c y c l [ S O 4 ] e v a p
In the formula, [HCO3]sulfuric acid is the HCO3 ion produced through the erosion of carbonate rock by H2SO4, and the unit of each element in the formula is mmol/L.
The chemical weathering rate of evaporative salt rock (EWR), the chemical weathering rate of silicate rock (SWR), and the chemical weathering rate of carbonate rock (CWR) in the Naqu River basin can be calculated according to the following formula [32]:
E W R = C a e v a p + N a e v a p + [ S O 4 ] e v a p + C l e v a p × F l o w A
S W R = K s i l + N a s i l + C a s i l + M g s i l + S i O 2 × F l o w A
C W R = C a c a r b + M g c a r b + 0.5 [ H C O 3 ] c a r b + 0.5 [ H C O 3 ] s u l × F l o w A
In the formula, Flow represents the flow rate, A represents the watershed area, and the unit of Xevap, Xsil, and Xcarb is mg/L.
The carbon sink effect of the chemical weathering of rocks in the watershed was mainly reflected by the CO2 in the atmosphere or soil through the HCO3 ions generated by the chemical weathering of rocks. The CO2 consumption flux of the chemical weathering of silicate and carbonate in the Naqu River basin can be calculated according to the following formula [32]:
C O 2 s i l = [ H C O 3 ] s i l × F l o w A
C O 2 c a r b = 0.5 × ( [ H C O 3 ] c a r b [ H C O 3 ] s u l f u r i c   a c i d ) × F l o w A
In the formula, [CO2]sil and [CO2]carb are the CO2 consumption flux during the chemical weathering of silicate rocks and carbonate rocks, respectively.

3. Results

3.1. Runoff

The precipitation and annual average temperature were 433.2 mm and 0.88 °C in 2019. About 73% of the precipitation occurred in July, August, and September, with the largest precipitation occurring in July, reaching 169.8 mm. The winter temperature in the Naqu River basin was below 0 °C, the average temperature was −9.2 °C, and the lowest temperature reached −10.6 °C.
As can be seen from Table 1, in 2019, the areas of the upper-, middle-, and downstream parts of the Naqu River were 6412.6 km2, 3962 km2, and 4338.8 km2, respectively, and the annual runoff was approximately 0.15 × 108 m³, 0.16 × 108 m³, and 0.32 × 108 m³, respectively. The four main tributaries of the Naqu River, Mugequ, Chengqu, Gongqu, and Luoquhave, had drainage areas of 2103 km2, 1090 km2, 1232 km2, and 1479 km2, respectively. Their annual runoff was approximately 0.019 × 108 m³, 0.005 × 108 m³, 0.012 × 108 m³, and 0.024 × 108 m³, respectively. Due to the relatively small annual runoff for Mugequ, Chengqu, and Gongqu, this study mainly focused on a runoff analysis of the upper-, middle-, and downstream parts of the Naqu River and its main tributary, Luoqu. The monthly precipitation, temperature, and runoff depth in the downstream area of the main Naqu stream are shown in Figure 4.
The discharge pattern of the Naqu River basin exhibited a singular mode that was characterized by a peak manifesting in the month of August, which coincided with the zenith of air temperatures and the heightened ablation activity experienced by the glaciers. The total runoff of Naqu River from July to September accounted for 70.90% of the water flux. The minimum, which was barely 0.9% of the runoff in August, happened in February. The peak occurred in August and was slightly lower in July and September in the middle of Luoqu. The total runoff from July to September accounted for 77.78% of the water flux. January had the least runoff, which made up just 0.6% of the August runoff. Overall, the seasonal fluctuation in runoff matched the changes in precipitation that the meteorological station recorded.

3.2. Physicochemical Parameters

The statistics of the major ions in the Naqu River basin are shown in Table 2. The data were obtained by clustering and analyzing the data from the samples collected from the Naqu River Basin in 2019; since the water samples were from the same river, they received basically the same classification. As seen in Table 2, the SD values of HCO3, TDS, and EC were large, which indicated that the values of HCO3, TDS, and EC in the datasets were far away from their mean values, which may have been related to the sampling time. By analyzing the sample data, it was seen that the HCO3 concentration of the river water in the upper reaches of the main stream of the Naqu River was the maximum value in the whole basin for four consecutive years, and the TDS and EC values of the river water in the Naqu River Basin were the maximum values in the whole basin in each sampling year.
It also can be seen in Table 2 that the measured temperature was 10.55–21.7 °C, and the highest temperature of the river appeared in the upstream area of the Chengqu river. The Naqu River basin was found to be alkaline with an average pH value of 8.93, which varied from 7.65 to 10.08 in August 2019. The coefficient of variation (CV) of the pH in the Naqu River basin was relatively small (7.38%), indicating that the pH value of the river may have changed less as a function of space. The turbidity of the Naqu River basin was between 4.36 and 172 NTU, with an average of 52.17 NTU. Its CV was the highest, as it exceeded 1, indicating that the turbidity may have had a significant variation as a function of space. The EC ranged from 53 to 586 μS/cm, with an average of 220.57 μS/cm. The variation trends of EC and TDS in the Naqu River basin were basically consistent. The value of DO was between 8.28 and 31.27 mg/L, with an average of 17.26 mg/L.
As calculated by using the values in Table 2, the total cationic charge (TZ+, TZ+ = 2Ca2+ + 2Mg2+ + K+ + Na+) of the river samples from the Naqu River basin was between 0.87 meq/L and 10.10 meq/L, with an average value of 4.06 meq/L. The total anionic charge (TZ, TZ = HCO3 + Cl + 2CO32− + 2SO42−) of the anions was between 1.06 meq/L and 10.02 meq/L, with an average value of 4.08 meq/L. Most of the NICB in the river was between −3% and 3% (Table 3), indicating that the anions and cations may have changed their relative balance and showing that the test results were credible. According to the mean value of the ionic concentration, the order of cation concentration of the Naqu River basin was Ca2+ > Na+ > Mg+ > K+, and the order of anion concentration was HCO3 > SO42− > Cl > NO3. HCO3 and SO42− were the main anions in the river. The mean concentration of HCO3 was lower than that in the middle and upper reaches of the Yarlung Tsangpo river (63.44 mg/L) [37] and about 3.34 times the average concentration in rivers around the world (51.8 mg/L) [38]. However, the CV of SO42− was large (93.95%), indicating that the content of SO42− in the main stream and tributaries of the Naqu River may have varied greatly as a function of space.
The spatial variation in the TDS content in the Naqu River basin is shown in Figure 5. It can be seen in Table 2 and Figure 5 that, in the whole Naqu River basin, the TDS values were between 53 mg/L and 586 mg/L, with an average value of 220.57 mg/L in August 2019. Combined with the analysis of the sampling data from 2016 to 2018, it can be seen that during these three years, the TDS values in the Naqu catchment were between 85.67 and 548.89 mg/L, with an average value of 241.36 mg/L. The TDS values were between 85.67 and 548.89 mg/L in 2016, while they were between 115 and 455 mg/L in 2017, and they were between 120 and 512 mg/L in 2018. It can be seen that from 2016 to 2019, the TDS values in the Naqu River basin had little change. From 2016 to 2018, the minimum value of the TDS content had an upward trend, while the maximum value decreased, which may have been related to the runoff of the Naqu River basin in each sampling year. The maximum TDS values appeared in the upper reaches of the Naqu River, while the TDS showed a decreasing trend from upstream to downstream, which may have been caused by the increasing runoff.
In this study, the hydrochemical types were used to study the hydrochemical characteristics of the Naqu River basin [33,39]. Through the analysis of river sample data collected in August 2019, the chemical characteristics of the main stream and a tributary stream in the Naqu River basin were determined and are shown in Table 3. The milli-equivalent proportions (meq%) of predominant ions are graphically delineated on a Piper diagram [40], which was subsequently extrapolated onto a central diamond field to discern and appraise the hydrogeochemical facies and typologies inherent in these aqueous specimens, as shown in Figure 6.
Considering these results in combination with those of our previous study [22], from 2016 to 2018, the most common type of water chemistry in the Naqu River basin was that with Ca-HCO3, which accounted for about 79.24%. The hydrochemical type was that of Ca-HCO3 in the upper reaches of Naqu, Sangqu, Basuoqu, and Mumuqu over the four years of the study (2016–2019); the hydrochemical type in the middle reaches of the Naqu River was that of Mg-HCO3, which indicated that the weathering of silicate may have affected the sources of ions in this section of the river. Except for 2019, the main stream of the Gongqu River showed the hydrochemical type of Ca-HCO3; for the other three years, the hydrochemical type was Ca-SO4, indicating the contribution of evaporative or sulfide minerals in those three years [41].

4. Discussion

4.1. Provenance of Solutes in River Water

4.1.1. Qualitative Analysis of Different Sources of Riverine Solutes

Our previous research used ion ratio diagrams and the Pearson correlation coefficient to qualitatively analyze the sources of ions in the Naqu River basin from 2016 to 2018 [22]. Two different analysis methods can be mutually corroborated. Combining the formulas and the Pearson correlation coefficient, it was found that during 2016–2018, Ca2+ had the strongest correlation with Mg2+ among the three cations of K+, Na+, and Mg2+, indicating that Ca2+ and Mg2+ may have had a common source. Na+ and K+ had the strongest correlation, indicating the chemical weathering of silicate rocks containing Na+ and K+. Through the analysis of the ions sampled in the Naqu River basin in 2019, the same results were found. The difference was that in 2019, the correlation coefficient of SiO2 in the river was analyzed, and it was found that SiO2 had the strongest correlation with K+ among the four cations of K+, Na+, Ca2+, and Mg2+, followed by Na+, thus further indicating the chemical weathering of silicate rocks containing Na+ and K+. SiO2 and HCO3 had the strongest correlation among the anions, which further indicated that HCO3 may have been mainly involved in the chemical weathering of silicate rocks in the Naqu River basin [42].
The analysis based on the Pearson correlation coefficient was also confirmed by ion ratio diagrams. Combining our results with those of previous studies, when analyzing ions in the Naqu River basin using ion ratio diagrams, it was found that during 2016–2018, carbonic-acid-mediated reactions in the Naqu River basin played a leading role in rivers other than the main stream of Gongqu. The ion analysis results of all rivers in the Naqu River basin in 2019 showed that the carbonic-acid-mediated reaction was mainly the result of the chemical weathering of rocks. The analysis of water samples from the Naqu River basin’s middle and upper regions from 2016 to 2018 showed that the chemical weathering of gypsum (CaSO4·2H2O) or anhydrite (CaSO4) had a certain effect on the chemical composition of the Naqu River basin [22]. It was also speculated that the chemical weathering of carbonate rocks had a great influence on the chemical composition of the water of Basuoqu, Mumuqu, Zongqingqu, and Luoqu. The chemical weathering of silicate rock may have had a great impact on the source of ions in the main stream of Mugequ. The water chemistry characteristics of Chengqu, the upstream area of the main Naqu stream, and the downstream area of the main Naqu stream were predominantly governed by the processes of carbonate and silicate weathering.

4.1.2. Quantification of Different Sources of Riverine Solutes

The average elevation of the Naqu River basin is as high as 4500 m, with little human activity. Moreover, the concentration of NO3 in the rivers of the Naqu River basin was low, with an average value of 1.58 mg/L in 2019. Therefore, the contribution of human activities to dissolved ions in the Naqu River was not considered.
The calculations of the atmospheric precipitation inputs are shown in Table 4. The rate of contribution of inputs of atmospheric precipitation to the ions in the rivers of the Naqu River basin was between 0.25% and 2.89%, and the average contribution rate was 0.88%.
This study examined the chemical weathering of rocks under the action of gypsum hydrolysis and sulfuric acid in three different scenarios. The corresponding αevap values were 0, 1, and 0.5, respectively. The calculations are shown in Table 5; the average rate of contribution of the chemical weathering of evaporite to the ions in the rivers of the Naqu River basin was the greatest in scenario 2. In scenario 3, the chemical weathering of evaporite played an important role in the upper regions of the Naqu River basin, and its rate of contribution was 27.25%. This also confirmed the qualitative analysis of the chemical weathering in the upstream area of the main Naqu stream. In addition, the chemical weathering of evaporite had a relatively high rate of contribution (21.7%) to the ion composition of the Luoqu River’s primary tributary.
According to the formula that was used for the calculations, the rate of contribution of the chemical weathering of silicate in the Naqu River basin ranged from 3.55% to 39.78%, with an average rate of contribution of 21.73%. The silicate rocks’ chemical weathering had the highest rate of contribution (39.78%) to the ionic composition of the primary tributary of the Mugequ River. This observation aligns seamlessly with the antecedent qualitative scrutiny of the impact arising from chemical weathering [22].
In concordance with the three input scenarios involving evaporite chemical weathering, the percentage of major ions in the Naqu River attributed to the chemical weathering of carbonate rock markedly exceeded that originating from the chemical weathering of evaporite rock and silicate rock. In scenario 3, a preeminence in the Basuoqu River was discerned, with the chemical weathering of carbonate registering the highest rate of contribution of 78.03%. Furthermore, the rates of contribution of the weathering of carbonate rock for Mumuqu, Zongqingqu, the main Gongqu stream, and the main Luoqu stream all exceeded 70%, which was in agreement with the qualitative analysis expounding the impact of chemical weathering on these riverine systems. In general, carbonates weather more easily than silicates do [43]. Consequently, a quantitative analysis of the contributions of different river-water ion sources revealed that carbonates were the most weathered. In the whole Naqu River basin, the sources contributing the main ions were: carbonate rock chemical weathering > silicate rock chemical weathering > evaporite rock chemical weathering > atmospheric precipitation input.
Combining this with the results of our previous qualitative research on the origin of ions in the Naqu River basin, atmospheric precipitation had a small contribution to the Naqu River basin, while the chemical weathering of carbonate rock contributed more to the sources of ions, and the chemical weathering of silicate rock and evaporite rock was also an important ion source [22]. It can be seen that the results of this quantitative analysis are basically consistent with the results of the qualitative analysis, confirming the reliability of this research.
The Nujiang–Salween River originates on the Tibetan Plateau. The maximum NO3 concentration in the Naqu River basin was 3.74 mg/L, which was higher than the maximum NO3 concentration of 2.82 mg/L in the Nujiang–Salwen River according to previous studies [44,45]. The maximum SiO2 concentration in the Nujiang–Salwen River was 11.5 mg/L, and the CV was 23.82%. Similarly, the CV of SiO2 in the Nujiang–Salwen River was 21.62% [46], indicating that the SiO2 concentration may change relatively smoothly. Moreover, the concentration of SiO2 in the samples of the Nujiang–Salween River in Hpa’an, Myanmar was mostly higher than that in the primary tributary of the Nujiang–Salween River in China [47], which may have been because the weathering of silicate rocks occurred more intensely in the Nujiang–Salween River of Myanmar [44].

4.2. Analysis of the Chemical Weathering Rates of Rocks and Carbon Sink Effects in Watersheds

Table 6 shows the calculation results for the chemical weathering rates of rocks in the Naqu River basin’s mid-upper regions.
In the three scenarios, the rate of chemical alteration of carbonate rocks was the highest. In scenario 1, the maximum rate of chemical weathering of carbonate rocks reached 16.84 t/km2/yr, and the minimum rate of chemical weathering of evaporative rock salt was 2.20 t/km2/yr. In scenario 3, as shown in Table 7, the order of rates of chemical weathering of different lithologies in the Naqu River basin’s mid-upper regions was: carbonate rock > silicate rock > evaporite rock.
The rate of chemical metamorphism of silicate rocks in the middle and upper regions of the main Naqu stream was 6.82 t/km2/yr. As shown in other studies, the rates of chemical weathering of silicate rocks in the Nujiang River basin were 5.90 t/km2/yr according to Wu [35], 4.27 t/km2/yr according to Tao [17], and 5.10 t/km2/yr according Zhang [44]. It can be seen that the rate of chemical weathering of silicate rocks in the middle and upper regions of the main Naqu stream was higher than that in the middle and lower reaches of the Nujiang River.
The carbon sink capacity of chemical weathering In the Naqu River basin’s mid-upper regions is shown in Table 7. Under the three conditions, the CO2 consumption fluxes of the chemical weathering of silicate and carbonate in the upper and middle reaches of the Naqu River were 315.07 × 103 mol/km2/yr, 369.72 × 103 mol/km2/yr, and 342.40 × 103 mol/km2/yr. According to the analysis, the carbon sink capacity of chemical weathering in the mid-upper regions of the main Naqu stream showed that the CO2 consumption flux of the chemical weathering of silicate was more intense than that of the chemical weathering of carbonate. Zhang’s study [45] showed that under the alkaline conditions of the Qinghai–Tibetan Plateau, the weathering of siliciclastic evolutionary rocks and alkalization of lakes are mutually reinforcing, resulting in the precipitation of large quantities of carbonate rocks. However, the IPCC’s Fifth Assessment Report emphasized that the chemical weathering of carbonate rocks not only produces carbon sinks, but the resulting sinks are also stable [46]. Thus, inland aquatic systems assume a pivotal role in the intricate dynamics of the global carbon cycle, and their strategic position in combating climate change is gradually coming to the fore, while the strong chemical erosion of the waterways of the Tibetan Plateau and the resulting high rate of carbon dioxide depletion play an important role in slowing down climate warming [47].
The calculation methods for the CO2 consumption rates of other rivers with corresponding lithologies and climate conditions were analyzed. The CO2 consumption flux of the chemical weathering of silicate in the middle and upper reaches of the Naqu River was 248.60 × 103 mol/km2/yr, which was 2.26 times that in the Nujiang River basin (110 × 103 mol/km2/yr) [34], 5.43 times that in the Tuotuo River (45.8 × 103 mol/km2/yr) [48], and about 55.41% of that in the Ganges River (448.611 × 103 mol/km2/yr) [34].
Table 7. Comparison of the carbon sink capacities of different catchments.
Table 7. Comparison of the carbon sink capacities of different catchments.
RiverDrainage Area[CO2]sil[CO2]carbSumReference
106 km2108
mol/yr
103
mol/km2/yr
108
mol/yr
103
mol/km2/yr
108
mol/yr
103
mol/km2/yr
The upstream and middle areas of Naqu River basin0.0125.79248.606.9066.4732.69315.07Scenario 1
0.0125.79248.6012.57121.1338.36369.72Scenario 2
0.0125.79248.609.7393.8035.52342.40Scenario 3
Nujiang River0.11120110650590890810Wu et al. [35]
Tuotuo River0.0167.345.821.4135--Li et al. [49]
Ganges River1.054710448.6112360224.8-692Gaillardet et al. [50]

5. Conclusions

Based on river samples collected in 2019, in this study, the primary ion fluxes in the main stream and tributaries of the Naqu River were estimated in August of each year from 2016 to 2019. Based on the forward model, the contributions of atmospheric precipitation, the chemical weathering of evaporated saline rocks, the chemical weathering of silicate rocks, and the chemical weathering of carbonate rocks to the dissolved ions in river water were quantitatively estimated, and the rate of chemical weathering of rocks in the mid-upper watersheds of the Naqu River and the carbon sink effect were analyzed; finally, the following conclusions were drawn:
a
The changes in terms of space of the TDS and EC values in the Naqu River watershed were the same; they decreased in concentration from upstream to downstream. The pH changed little as a function of space, but the difference was that the turbidity changed greatly. Ca-HCO3 was the main type of water chemistry in the Naqu River basin.
b
Based on the forward model, the contributions of different sources of river water ions were quantitatively analyzed, and they are listed in descending order: chemical weathering of carbonate rocks > chemical weathering of silicate rocks > chemical weathering of evaporite saline rocks > atmospheric precipitation input.
c
By analyzing the chemical weathering and carbon sink effects in the watershed, three scenarios were hypothesized based on the percentage of SO42− production from rock weathering, with carbonate chemical weathering rates being the largest. In the upstream and middle of the Naqu River basin, the chemical weathering of silicate rocks consumed more CO2 than that of carbonate rocks did.

Author Contributions

All authors contributed to the study conception and design. Data curation, H.Z. (Heng Zhao), H.Z. (Honglu Zhang), F.W., and P.K.; Supervision, H.Z. (Heng Zhao), H.Z. (Honglu Zhang), F.W., and P.K.; Writing—original draft, S.R. and X.H.; Writing—review and editing, H.Z. (Heng Zhao), H.Z. (Honglu Zhang), F.W., and P.K. All authors commented on previous versions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This project was financially supported by the National Natural Science Foundation of the People’s Republic of China (52279014), the Henan Province Key Research and Development and Promotion Project (Science and Technology) (232102320257), and the Key Research and Development Program of Ningxia Hui Autonomous Region (2021BEG02012).

Data Availability Statement

The original contributions presented in this study are included in the article; further inquiries can be directed to the corresponding author.

Acknowledgments

The authors would like to express their sincere gratitude to the anonymous reviewers for their constructive comments and useful suggestions that helped us improve this study.

Conflicts of Interest

The authors declare that this research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.

References

  1. Yang, R.; Sun, H.; Chen, B.; Yang, M.; Zeng, Q.; Zeng, C.; Huang, J.; Luo, H.; Lin, D. Temporal variations in riverine hydrochemistry and estimation of the carbon sink produced by coupled carbonate weathering with aquatic photosynthesis on land: An example from the Xijiang River, a large subtropical karst-dominated river in China. Environ. Sci. Pollut. Res. Int. 2020, 27, 13142–13154. [Google Scholar] [CrossRef] [PubMed]
  2. Gong, S.; Bai, X.; Luo, G.; Li, C.; Wu, L.; Chen, F.; Ran, C.; Xi, H.; Zhang, S. Climate change has enhanced the positive contribution of rock weathering to the major ions in riverine transport. Glob. Planet. Chang. 2023, 228, 104203. [Google Scholar] [CrossRef]
  3. Bastia, F.; Equeenuddin, S.M. Chemical weathering and associated CO2 consumption in the Mahanadi river basin, India. J. Asian Earth Sci. 2019, 174, 218–231. [Google Scholar] [CrossRef]
  4. Huang, X.; Sillanpää, M.; Duo, B.; Gjessing, E.T. Water quality in the Tibetan Plateau: Metal contents of four selected rivers. Environ. Pollut. 2008, 156, 270–277. [Google Scholar] [CrossRef] [PubMed]
  5. Zhang, S.; Bai, X.; Zhao, C.; Tan, Q.; Luo, G.; Wang, J.; Li, Q.; Wu, L.; Chen, F.; Li, C.; et al. Global CO2 consumption by silicate rock chemical weathering: Its past and future. Earth’s Future 2021, 9, e2020EF001938. [Google Scholar] [CrossRef]
  6. Xi, H.; Wang, S.; Bai, X.; Tang, H.; Luo, G.; Li, H.; Wu, L.; Li, C.; Chen, H.; Ran, C.; et al. The responses of weathering carbon sink to eco-hydrological processes in global rocks. Sci. Total Environ. 2021, 788, 147706. [Google Scholar] [CrossRef]
  7. Jia, B.; Zhou, G. Estimation of global karst carbon sink from 1950s to 2050s using response surface methodology. Geospat. Inf. Sci. 2023, 1–18. [Google Scholar] [CrossRef]
  8. Liu, Z.; Dreybrodt, W. Significance of the carbon sink produced by H2O–carbonate–CO2–aquatic phototroph interaction on land. Sci. Bull. 2015, 60, 182–191. [Google Scholar] [CrossRef]
  9. Liu, Z.; Macpherson, G.L.; Groves, C.; Martin, J.B.; Yuan, D.; Zeng, S. Large and active CO2 uptake by coupled carbonate weathering. Earth-Sci. Rev. 2018, 182, 42–49. [Google Scholar] [CrossRef]
  10. Yu, Z.L.; Yan, N.; Wu, G.J.; Xu, T.L.; Li, F. Chemical weathering in the upstream and midstream reaches of the Yarlung Tsangpo basin, southern Tibetan Plateau. Chem. Geol. 2020, 559, 119906. [Google Scholar] [CrossRef]
  11. Yu, Z.; Wu, G.; Li, F.; Chen, M.; Tran, T.V.; Liu, X.; Gao, S. Glaciation enhanced chemical weathering in a cold glacial catchment, western Nyaingêntanglha Mountains, central Tibetan Plateau. J. Hydrol. 2021, 597, 126197. [Google Scholar] [CrossRef]
  12. Yu, Z.L.; Wu, G.J.; Li, F.; Huang, J.; Xiao, X.; Liu, K.S. Small-catchment perspective on chemical weathering and its controlling factors in the Nam Co basin, central Tibetan Plateau. J. Hydrol. 2021, 598, 126315. [Google Scholar] [CrossRef]
  13. Yu, Z.; Wu, G.; Keys, L.; Li, F.; Yan, N.; Qu, D.; Liu, X. Seasonal variation of chemical weathering and its controlling factors in two alpine catchments, Nam Co basin, central Tibetan Plateau. J. Hydrol. 2019, 576, 381–395. [Google Scholar] [CrossRef]
  14. Zhang, D.; Zhao, Z.Q.; Li, X.D.; Li, X.D.; Zhang, L.L.; Chen, A.C. Assessing the oxidative weathering of pyrite and its role in controlling atmospheric CO2 release in the eastern Qinghai-Tibet Plateau. Chem. Geol. 2020, 543, 119605. [Google Scholar] [CrossRef]
  15. Kang, X.; Li, Y. Chemical weathering of small watersheds in the upper reaches of Jinsha River on the eastern edge of the Qinghai–Tibet Plateau. J. Hydrol. 2023, 623, 129844. [Google Scholar] [CrossRef]
  16. Wu, W.F.; Xu, S.J.; Yang, J.D.; Yin, H.W. Silicate Weathering and CO2 Consumption Deduced from the Seven Chinese Rivers Originating in the Qinghai-Tibet Plateau. Chem. Geol. 2008, 249, 307–320. [Google Scholar] [CrossRef]
  17. Tao, Z.H.; Zhang, Z.Q.; Chen, D.; Li, X.D.; Liu, C.Q. Chemical Weathering in the Three Rivers (Jinsha River, Lancang River and Nujiang River) Watershed, Southwest China. Chin. J. Ecol. 2015, 34, 2297–2308. (In Chinese) [Google Scholar] [CrossRef]
  18. Zhao, Y.; Zhao, H.; Wang, F.Q.; Chen, X.; Ci, R.Q.J. Multivariate Statistical Analysis on Hydrochemical Characteristics of Naqu River in the Source Region of Nujiang River. J. North China Univ. Water Resour. Electr. Power Nat. Sci. Ed. 2020, 41, 63–71. (In Chinese) [Google Scholar] [CrossRef]
  19. Cui, B.L.; Li, X.Y. Runoff processes in the Qinghai Lake Basin, Northeast Qinghai-Tibet Plateau, China: Insights from stable isotope and hydrochemistry. Quat. Int. 2015, 380–381, 123–132. [Google Scholar] [CrossRef]
  20. Yan, D.; Liu, S.; Qin, T.; Weng, B.S.; Liu, J.J. Evaluation of TRMM precipitation and its application to distributed hydrological model in Naqu River Basin of the Tibetan Plateau. Hydrol. Res. 2017, 48, 822–839. [Google Scholar] [CrossRef]
  21. Chen, X.; Wang, G.L.; Wang, F.Q.; Yan, D.H.; Zhao, H. Characteristics of Water Isotopes and Water Source Identification During the Wet Season in Naqu River Basin, Qinghai-Tibet Plateau. Water 2019, 11, 2418. [Google Scholar] [CrossRef]
  22. Wang, F.Q.; Zhao, Y.; Chen, X.; Zhao, H. Hydrochemistry and Its Controlling Factors of Rivers in the Source Region of the Nujiang River on the Tibetan Plateau. Water 2019, 11, 2166. [Google Scholar] [CrossRef]
  23. Yang, Y.; Wu, Q.; Hou, Y.; Zhang, Z.; Zhan, J.; Gao, S.; Jin, H. Unraveling of permafrost hydrological variabilities on Central Qinghai-Tibet Plateau using stable isotopic technique. Sci. Total Environ. 2017, 605–606, 199–210. [Google Scholar] [CrossRef]
  24. Liu, S.; Yan, D.; Qin, T.; Weng, B.; Lu, Y.; Dong, G.; Gong, B. Precipitation phase separation schemes in the Naqu River basin, eastern Tibetan plateau. Theor. Appl. Climatol. 2018, 131, 399–411. [Google Scholar] [CrossRef]
  25. Deshmukh, R.Y.; Sharma, P.; Dhote, P.R. Review on hydrochemical analysis of rivers. Concern 2020, 7. [Google Scholar]
  26. Kai, J.; Wang, J.; Ju, J.; Huang, L.; Ma, Q.; Daut, G.; Zhu, L. Spatio-temporal variations of hydrochemistry and modern sedimentation processes in the Nam Co basin, Tibetan Plateau: Implications for carbonate precipitation. J. Great Lakes Res. 2020, 46, 961–975. [Google Scholar] [CrossRef]
  27. Ministry of Geology and Mineral Resources of the People’s Republic of China. Testing Methods of Underground Water Quality-Determination of Carbonate, Bicarbonate and Hydroxide-Titrimetric; China Standards Press: Beijing, China, 1993. (In Chinese) [Google Scholar]
  28. Ministry of Geology and Mineral Resources of the People’s Republic of China. Testing Methods of Underground Water Quality-Determination of Sulfate-Turbidimetry; China Standards Press: Beijing, China, 1993. (In Chinese) [Google Scholar]
  29. Ministry of Environmental Protection of the People’s Republic of China. Water Quality-Determination of 32 Elements-Inductively Coupled Plasma Optical Emission Spectrometry; China Environmental Sciences Press: Beijing, China, 2015. (In Chinese) [Google Scholar]
  30. Moquet, J.S.; Crave, A.; Viers, J.; Seyler, P.; Guyot, J.L. Chemical Weathering and Atmospheric/Soil CO2 Uptake in the Andean and Foreland Amazon Basins. Chem. Geol. 2011, 287, 1–26. [Google Scholar] [CrossRef]
  31. Xie, C.J.; Gao, Q.Z.; Tao, Z. Review and Perspectives of the Study on Chemical Weathering and Hydrochemistry in River Basin. Trop. Geogr. 2012, 32, 331–337. [Google Scholar] [CrossRef]
  32. Jiang, P.; Yu, G.; Zhang, Q.; Zou, Y.; Tang, Q.; Kang, Z.; Sytharith, P.; Xiao, H. Chemical weathering and CO2 consumption rates of rocks in the Bishuiyan subterranean basin of Guangxi, China. Sci. Rep. 2020, 10, 11677. [Google Scholar] [CrossRef]
  33. Vespasiano, G.; Apollaro, C.; De Rosa, R.; Muto, F.; Larosa, S.; Fiebig, J.; Mulch, A.; Marini, L. The Small Spring Method (SSM) for the definition of stable isotope–elevation relationships in Northern Calabria (Southern Italy). Appl. Geochem. 2015, 63, 333–346. [Google Scholar] [CrossRef]
  34. Wu, W.; Yang, J.; Xu, S.; Yin, H. Geochemistry of the headwaters of the Yangtze River, Tongtian He and Jinsha Jiang: Silicate weathering and CO2 consumption. Appl. Geochem. 2008, 23, 3712–3727. [Google Scholar] [CrossRef]
  35. Wu, W. Hydrochemistry of inland rivers in the north Tibetan Plateau: Constraints and weathering rate estimation. Sci. Total Environ. 2016, 541, 468–482. [Google Scholar] [CrossRef] [PubMed]
  36. Li, Q.; Wu, J.; Shen, B.; Zeng, H.; Li, Y. Water Chemistry and Stable Isotopes of Different Water Types in Ta-jikistan. Environ. Process. 2018, 5, 127–137. [Google Scholar] [CrossRef]
  37. Jiang, L.G.; Yao, Z.J.; Wang, R.; Liu, Z.; Wang, L.; Wu, S. Hydrochemistry of the middle and upper reaches of the Yarlung Tsangpo River system: Weathering processes and CO2 consumption. Environ. Earth Sci. 2015, 74, 2369–2379. [Google Scholar] [CrossRef]
  38. Meybeck, M. Global occurrence of major elements in rivers. Treatise Geochem. 2003, 5, 207–223. [Google Scholar] [CrossRef]
  39. Wu, L.; Huh, Y.; Qin, J.; Gu, D.; Lee, S.V.D. Chemical weathering in the Upper Huang He (Yellow River) draining the eastern Qinghai-Tibet Plateau. Geochim. Cosmochim. Acta 2005, 69, 5279–5294. [Google Scholar] [CrossRef]
  40. Piper, A.M. A graphic procedure in the geochemical interpretation of water-analyses. Eos Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar]
  41. Dlamini, A.E.; Demlie, M. Integrated hydrogeological, hydrochemical and environmental isotope investi-gation of the area around the Kusile Power Station, Mpumalanga, South Africa. J. Afr. Earth Sci. 2020, 172, 103958. [Google Scholar] [CrossRef]
  42. Kobus, S.; Glińska-Lewczuk, K.; Sidoruk, M.; Skwierawski, A.; Obolewski, K.; Timofte, C.M.; Sowiński, P. Effect of hy-drological connectivity on physicochemical properties of bottom sediments of floodplain lakes-Acase study of the Lyna River, Northern Poland. Environ. Eng. Manag. J. 2016, 15, 1237–1246. [Google Scholar] [CrossRef]
  43. Meybeck, M. Global chemical weathering of surficial rocks estimated from river dissolved loads. Am. J. Sci. 1987, 287, 401–428. [Google Scholar] [CrossRef]
  44. Zhang, L.L.; Zhao, Z.Q.; Zhang, W.; Tao, Z.H.; Huang, L.; Yang, J.X.; Wu, Q.X.; Liu, C.Q. Characteristics of water chemistry and its indication of chemical weathering in Jinshajiang, Lancangjiang and Nujiang drainage basins. Environ. Earth Sci. 2016, 75, 506. [Google Scholar] [CrossRef]
  45. Zhang, S.; Shao, M.; Wang, T.; Pei, Y.; Chen, B. Geochemistry of lacustrine carbonate rocks in southwestern Qaidam: Implications of silicate weathering and carbon burial triggered by the uplift of the Tibetan Plateau. Int. J. Coal Geol. 2023, 265, 104167. [Google Scholar] [CrossRef]
  46. Bai, X.; Zhang, S.; Li, C.; Xiong, L.; Song, F.; Du, C.; Li, M.; Luo, Q.; Xue, Y.; Wang, S. A carbon-neutrality-capactiy index for evaluating carbon sink contributions. Environ. Sci. Ecotechnol. 2023, 15, 100237. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, H.; Ju, P.; Zhu, Q.; Xu, X.; Wu, N.; Gao, Y.; Feng, X.; Tian, J.; Niu, S.; Zhang, Y.; et al. Carbon and nitrogen cycling on the Qinghai–Tibetan Plateau. Nat. Rev. Earth Environ. 2022, 3, 701–716. [Google Scholar] [CrossRef]
  48. Chapman, H.; Bickle, M.; Thaw, S.H.; Thiam, H.N. Chemical fluxes from time series sampling of the Irrawaddy and Salween Rivers, Myanmar. Chem. Geol. 2015, 401, 15–27. [Google Scholar] [CrossRef]
  49. Li, Z.J.; Li, Z.X.; Song, L.L.; Ma, J.Z.; Song, Y. Environment significance and hydrochemical characteristics of supra-permafrost water in the source region of the Yangtze River. Sci. Total Environ. 2018, 644, 1141–1151. [Google Scholar] [CrossRef]
  50. Gaillardet, J.; Dupré, B.; Louvat, P.; Allègre, C.J. Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers. Chem. Geol. 1999, 159, 3–30. [Google Scholar] [CrossRef]
Figure 1. A map showing the location of Naqu River basin.
Figure 1. A map showing the location of Naqu River basin.
Water 15 04191 g001
Figure 2. A map of the Naqu River basin showing the sampling sites, DEM, and water system.
Figure 2. A map of the Naqu River basin showing the sampling sites, DEM, and water system.
Water 15 04191 g002
Figure 3. The correlation between EC and TZ+ in the Naqu River basin.
Figure 3. The correlation between EC and TZ+ in the Naqu River basin.
Water 15 04191 g003
Figure 4. Precipitation, air temperature, and runoff in the downstream area of the main Naqu stream.
Figure 4. Precipitation, air temperature, and runoff in the downstream area of the main Naqu stream.
Water 15 04191 g004
Figure 5. The variation in TDS values in the Naqu River basin.
Figure 5. The variation in TDS values in the Naqu River basin.
Water 15 04191 g005
Figure 6. Piper diagram showing the concentrations of cations and anions in the stream water of the Naqu River basin (meq/L).
Figure 6. Piper diagram showing the concentrations of cations and anions in the stream water of the Naqu River basin (meq/L).
Water 15 04191 g006
Table 1. Drainage and runoff of the main rivers in the Naqu River basin.
Table 1. Drainage and runoff of the main rivers in the Naqu River basin.
RiverDrainage Area
/km2
Runoff
/108 m3
Runoff Depth
/mm
Upstream of the main Naqu stream6412.60.152.29
Middle of the main Naqu stream39620.164.01
Downstream of the main Naqu stream4338.80.327.33
Mugequ 21030.0190.92
Chengqu 10900.0050.47
Gongqu12320.0120.94
Luoqu14790.0241.60
Table 2. Statistical results of the hydrochemical characteristics of the Naqu River.
Table 2. Statistical results of the hydrochemical characteristics of the Naqu River.
ParameterMinMaxMeanSDCV/(%)
Ca2+/(mg/L)12.181.2537.7719.2350.92
Mg2+/(mg/L)1.9637.8313.599.7671.83
K+/(mg/L)0.395.542.601.8169.45
Na+/(mg/L)2.1963.122.4218.1881.09
Cl/(mg/L)2.4828.687.146.7494.43
HCO3/(mg/L)43.63321.15173.2487.6550.60
SO42−/(mg/L)11.72188.8547.2144.3693.95
NO3/(mg/L)0.583.741.580.8352.71
SiO2/(mg/L)5.0411.57.421.9926.82
Sr/(mg/L)0.0481.690.270.421.57
TDS/(mg/L)53586220.57137.8462.49
pH7.6510.088.930.667.38
Turbidity/NTU4.3617252.1753.46102.48
DO/(mg/L)8.2831.2717.266.9040.00
ORP/(mV)32181121.1435.4129.23
EC/(μS/cm)81915341.43213.8562.63
Temperature/(°C)10.5521.715.243.5223.10
Note: SD means standard deviation; CV means the coefficient of variation.
Table 3. Hydrochemical types of the main rivers in the Naqu River basin.
Table 3. Hydrochemical types of the main rivers in the Naqu River basin.
RiverElevationSampling TimeHydrochemical TypeNICB
(m)(yyyy-mm-dd)
Upstream of the main Naqu stream47562019-08-15Ca-HCO30.0081
Middle of the main Naqu stream45392019-08-15Mg-HCO3−0.0117
Downstream of the main Naqu stream44552019-08-16Na-HCO30.0382
Sangqu46282019-08-15Ca-HCO30.0055
Basuoqu47102019-08-15Ca-HCO3−0.0106
Mumuqu46252019-08-15Ca-HCO3−0.0018
Upstream of Chengqu stream45222019-08-15Ca-HCO30.0036
Downstream of Chengqu stream44992019-08-15Ca-HCO3−0.0503
Zongqingqu45652019-08-15Ca-HCO3−0.1216
Mugequ main stream45902019-08-14Na-HCO30.0058
Mugequ basin tributary46972019-08-14Ca-HCO3−0.0314
Gongqu main stream45742019-08-16Ca-HCO30.2221
Gongqu basin tributary44892019-08-16Ca-HCO3−0.0584
Luoqu main stream49862019-08-16Ca-HCO3−0.2119
Table 4. Ionic concentration due to atmospheric precipitation.
Table 4. Ionic concentration due to atmospheric precipitation.
BasinCa2+Mg2+K+Na+SO42−Average Contribution Rate
μmol/L%
Naqu River4.127.950.161.009.390.88
(0.25–2.89)
Table 5. Contributions of solutes in each reservoir of the Naqu River basin.
Table 5. Contributions of solutes in each reservoir of the Naqu River basin.
RiverElevation
(m)
Precipitation
(%)
Evaporite (%)Silicate
(%)
Carbonate (%)
Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3
Upstream of the main Naqu stream47560.257.8746.6427.2525.0166.8728.1147.49
Middle of the main Naqu stream45390.505.0530.6517.8533.1461.3235.7248.52
Downstream of the main Naqu stream44550.515.6021.8913.7537.0356.8540.5748.71
Sangqu46280.413.3728.7316.0526.8269.4044.0456.72
Basuoqu47100.492.4714.318.3913.0983.9572.1178.03
Mumuqu46250.992.2312.987.6114.4282.3671.6176.99
Upstream of Chengqu stream45220.502.7216.549.6337.1959.5845.7652.67
Downstream of Chengqu stream44990.523.4221.6812.5532.5563.5145.2554.38
Zongqingqu45651.213.6832.3017.996.0289.0860.4774.78
Mugequ main stream45900.847.5320.5514.0439.7851.8538.8345.34
Mugequ basin tributary46971.012.6733.5318.1014.4381.8951.0466.46
Gongqu main stream45741.244.2431.0617.657.3887.1460.3273.73
Gongqu basin tributary44891.023.1040.1721.6313.8482.0544.9863.51
Luoqu main stream49862.898.8034.6021.703.5584.7558.9571.85
Table 6. Chemical weathering rates of rocks in the upstream and middle parts of Naqu River basin in 2019.
Table 6. Chemical weathering rates of rocks in the upstream and middle parts of Naqu River basin in 2019.
EWRSWRCWRSum
Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3Scenario 1Scenario 2Scenario 3
t/km2/yr
2.209.635.926.8216.8411.3214.0825.8727.7826.82
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ren, S.; Zhao, H.; Hou, X.; Zhang, H.; Wang, F.; Kang, P. The Chemical Weathering of Rocks and Its Carbon Sink Effect in the Naqu River Basin of the Nujiang River Source Area, Southwest China. Water 2023, 15, 4191. https://doi.org/10.3390/w15234191

AMA Style

Ren S, Zhao H, Hou X, Zhang H, Wang F, Kang P. The Chemical Weathering of Rocks and Its Carbon Sink Effect in the Naqu River Basin of the Nujiang River Source Area, Southwest China. Water. 2023; 15(23):4191. https://doi.org/10.3390/w15234191

Chicago/Turabian Style

Ren, Suming, Heng Zhao, Xinli Hou, Honglu Zhang, Fuqiang Wang, and Pingping Kang. 2023. "The Chemical Weathering of Rocks and Its Carbon Sink Effect in the Naqu River Basin of the Nujiang River Source Area, Southwest China" Water 15, no. 23: 4191. https://doi.org/10.3390/w15234191

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