Evaluation and Prediction of Water Quality of Typical Wetlands in the Source Region of the Yangtze River

: Wetlands play an important role in water storage and water conservation, but with global climate change, the degradation of wetland ecosystems is accelerating. In this study, we conducted research on the current situation and future prediction of water quality in typical wetlands in the source region of the Yangtze River to provide a scientiﬁc basis for the protection and restoration of wetlands in the source region of the Yangtze River. The Bayesian water quality assessment method and Yao Zhiqi evaluation method were used to evaluate the water quality of typical wetlands in the source region of the Yangtze River from 2016 to 2021 and based on the climate change scenarios of three RCPs (Representative Concentration Pathways) under the CMIP5 (Coupled Model Intercomparison Project Phase 5) global climate model and SWAT (soil and water assessment tool) hydrological model, the wetland water quality in the source region of the Yangtze River from 2022 to 2100 was predicted. The results show that the inter-annual changes in COD Mn , NH 3 -N, and TN in a typical wetland show a downward trend, while the temperature and DO concentration show an upward trend from 2016–2021. The changes in COD Mn , temperature, and conductivity within the year are abundant season > ﬂat season > dry season; and DO, NH 3 -A, TN, and TP concentrations within the year are opposite. The water quality of typical wetlands in the source region of the Yangtze River has reached Class II and above. From 2022 to 2100, under climate change in the future, TN, TP, COD Mn , NH 3 -N, and temperature in the wetland water in the source region of the Yangtze River will continue to rise, and the concentration of DO will continue to decline. Therefore, the pressure on water resources in the source region of the Yangtze River is further aggravated, so it is urgent to strengthen water resources protection.


Introduction
Wetlands play a major role in regulating and storing water resources and conserving water sources [1]. The wetland ecosystem is the most sensitive to climate change [2,3]. Under the background of global change, the wetland ecosystem is accelerating degradation all over the world [4]. The wetland ecosystem has been part of the most threatened ecosystems in the world [5]. The wetland region of the Qinghai Tibet Plateau, about 13.19 × 10 4 km 2 and accounting for 20% of China's wetland region, is a sensitive region to global climate change and an important ecological security barrier in China and even Asia [6]. The quality of the wetland water body and evaluation of wetland water quality have traditionally been the focus of attention, and the work has developed accordingly. The moisture content is significantly reduced, swampy meadows begin to evolve into alpine meadows, water content of the soil layer is reduced, shrinkage of wetlands is intensified, vegetation cover is reduced, and desertification is intensified. Therefore, it is urgent to protect the wetlands in the source region of the Yangtze River.
At present, there are few studies on the current situation of the water quality of typical wetlands in the source region of the Yangtze River, and most of them are short-term studies. There are some problems, such as one-sided water quality indicators and single evaluation methods; at present, the research on the SWAT hydrological model in the source region of the Yangtze River has achieved some results, but there are some problems such as short research space and time scale [30], and the research on future water quality prediction by climate change coupled with the SWAT hydrological model has not been reported on. It is of great significance to study the present situation and future prediction of the water quality of typical wetlands in the source region of the Yangtze River for the protection and restoration of wetlands in the source region of the Yangtze River and formulate the water environment protection and restoration plan for the development and utilization of water resources in the source region of the Yangtze River. This study analyzes the current situation of seven water quality indexes (permanganate index, ammonia nitrogen, total nitrogen, total phosphorus, temperature, dissolved oxygen, and conductivity) of typical wetland water quality in the source region of the Yangtze River from 2016 to 2021. The Bayesian water quality evaluation method and comprehensive pollution index method are used to evaluate the current situation of typical wetland water quality in the source region of the Yangtze River, as well as three RCP (RCP2.6, RCP4.5, and RCP8.5) scenarios coupled with the SWAT hydrological model to establish the change in sequence of water quality in the source region of the Yangtze River from 2022 to 2100 to provide a scientific basis for the protection of typical wetlands and the protection of water resources in the source region of the Yangtze River and the safety of water resources in the Yangtze River Basin under future climate change.

Overview of the Study Region
The source region of the Yangtze River is situated above the Zhimenda hydrological station, with an average altitude of about 4500 m, geographical coordinates of 32 • Table 1). The western part of the region is the Inland Lake region of Qiangtang Mountains in northern Tibet. The northwest is the eastern section of the Kunlun Mountains, and the southwest is the middle of Tanggula Mountain and the eastern section. The average altitude of the two mountains is 5500-6000 m. The top of the mountain is the subject with perennial snow or glaciers, and between the two mountains is a vast plateau [32]. It is a wetland ecosystem with a large region, high altitude, and the most significant biodiversity in the world. The source region of the Yangtze River is rich in biological resources and is the habitat of some terrestrial organisms, such as Tibetan antelope, rock sheep, rabbit, snow leopard, Tibetan fox, black can, etc. The aquatic resources include Herzensteinia microcephalus, Ptychobarbus kaznakovi, Triplophysa stoliczkae, Triplophysa stenura, and Triplophysa leptosoma, etc. It is likewise the main source of China's freshwater resources and an important barrier for ecological security [33]. The region has a high proportion of endemic species, a large number of Grade I and II protected animals, and many rare plant species and is one of the areas with a rich distribution of endemic plant and animal resources in China, as well as an important treasure house of biological resources, species resources, and biological genes in China. The source region of the Yangtze River is a very important ecological security barrier in China's "two screens and three belts" ecological security strategy pattern (The "two screens" include the Qinghai-Tibet Plateau Ecological Barrier and Loess Plateau-Sichuan-Yunnan Ecological Barrier; the "three belts" refer to the northeastern forest belt, southern hilly mountain belt, and northern sand control belt.), and it is also a key ecological function area in China's key ecological function areas, carrying important ecological functions such as water conservation and storage, maintenance of biodiversity, soil and water conservation, and wind and sand control. The ecological quality of the source region of Yangtze River and its water recharge capacity are not only directly related to the socio-economic development of the source area, but also have a direct or indirect impact on many areas in the middle and lower reaches of each basin. The perennial permafrost in the source region of the Yangtze River is widespread, forming a stable water barrier that is not restricted by geological features. The regional groundwater includes three special groundwater types of water above the frozen layer, water below the frozen layer, and groundwater in the thaw zone. The river and lake melt area in the source area is mostly the convergence center of groundwater and surface water, and the lithology of the aquifer is mostly quaternary alluvial sand and gravel layer, with large thickness, high permeability, and a large volume of water, which is the main object of groundwater mining.
well as an important treasure house of biological resources, species resources, and biological genes in China. The source region of the Yangtze River is a very important ecological security barrier in China's "two screens and three belts" ecological security strategy pattern (The "two screens" include the Qinghai-Tibet Plateau Ecological Barrier and Loess Plateau-Sichuan-Yunnan Ecological Barrier; the "three belts" refer to the northeastern forest belt, southern hilly mountain belt, and northern sand control belt.), and it is also a key ecological function area in China's key ecological function areas, carrying important ecological functions such as water conservation and storage, maintenance of biodiversity, soil and water conservation, and wind and sand control. The ecological quality of the source region of Yangtze River and its water recharge capacity are not only directly related to the socio-economic development of the source area, but also have a direct or indirect impact on many areas in the middle and lower reaches of each basin. The perennial permafrost in the source region of the Yangtze River is widespread, forming a stable water barrier that is not restricted by geological features. The regional groundwater includes three special groundwater types of water above the frozen layer, water below the frozen layer, and groundwater in the thaw zone. The river and lake melt area in the source area is mostly the convergence center of groundwater and surface water, and the lithology of the aquifer is mostly quaternary alluvial sand and gravel layer, with large thickness, high permeability, and a large volume of water, which is the main object of groundwater mining.

Data Sources
Water samples were collected from 9 wetland monitoring points in the source region of the Yangtze River. The geographical location of the source region of the Yangtze River is shown in Figure 1, and the information of each monitoring point is shown in Table 1.

Data Sources
Water samples were collected from 9 wetland monitoring points in the source region of the Yangtze River. The geographical location of the source region of the Yangtze River is shown in Figure 1, and the information of each monitoring point is shown in Table 1. season (April, May, September), and wet season (June, August). Each sampling time is from 10 a.m. to 2 p.m. Sampling shall be conducted in accordance with technical guidance for water quality sampling (HJ494-2009) [34]. Three water indices were measured directly by a Merck portable multi-parameter tester (WTM) in the field: water temperature, conductivity, and dissolved oxygen. Water samples from the monitoring sections were transported to the laboratory for analyses within 12 h in polyethylene bottles.TN was determined by potassium persulfate oxidation ultraviolet spectrophotometry (GB11893-89) [34], ammonia nitrogen (NH 3 -N) was determined by the flow analyzer method, total phosphorus (TP) was determined by the potassium persulfate digestion method (GB11893-89) [34], and the permanganate index (COD Mn ) was determined by the sulfuric acid method (The sample was added with a known amount of potassium permanganate and sulfuric acid, heated in boiling water for 30 min to oxidize some organic and inorganic reducing substances, and after the reaction, excess sodium oxalate was added to reduce the remaining potassium permanganate, and then the excess sodium oxalate was back-dropped with the standard solution of potassium permanganate.) (GB11892-89) [34]. Table 2 shows the I-V standard limits of dissolved oxygen, ammonia nitrogen, total nitrogen, total phosphorus, and permanganate and the environmental quality standard for typical wetlands in China (GB3838-2002) [34].

Bayesian Water Quality Evaluation Method
Bayesian formula [35]: where i is the standard type, i = 1, 2, . . . , s; i is the indicator, j = 1, 2, . . . , s, j for each monitoring point indicator observation value; y ji is the water quality-type standard value. P(y ji ) represents the a priori probability, on behalf of the monitoring point water quality may belong to some level of water quality. P(x j /y ji ) is a likelihood probability, representing the likelihood of a particular water quality indicator corresponding to different levels of water quality at the monitoring site.
Derive the posterior probability that the water quality level at the monitoring point is i.
where w j is for the first j indicator of the impact of the weight of water quality. P i for m water quality indicators after the integrated post-test probability; in this study, w j is uniformly taken as 1/5, that is, the various water quality indicators accounted for the same weight.
where P h is the final level of water quality at the monitoring site; a low value of P h means the better the water quality.

Yao Zhiqi Evaluation Method
Yao Zhiqi Evaluation Formula [36]: where X is the highest score index, that is, the highest value of the ratio of each measured value to the standard value. Y is the average score index, that is, the mean value of the ratio of each measured value to the standard value. If p value is smaller, the water quality is better. The classification bounds of Yao Zhiqi's evaluation method are: I represents excellent, P ≤ 0.90. II represents good, P ≤ 1.51. III represents light pollution, P ≤ 2.63. IV represents moderate pollution, P ≤ 3.55. V represents heavy pollution, P > 3.55 [36].

Water Quality Prediction Method
The typical wetland water quality prediction flow chart of the source region of Yangtze River is shown in Figure 2.
water quality indicators after the integrated post-test probability; in this study, is uniformly taken as 1/5, that is, the various water quality indicators accounted for the same weight.

= max
where is the final level of water quality at the monitoring site; a low value of means the better the water quality.

Yao Zhiqi Evaluation Method
Yao Zhiqi Evaluation Formula [36]: where X is the highest score index, that is, the highest value of the ratio of each measured value to the standard value. Y is the average score index, that is, the mean value of the ratio of each measured value to the standard value. If p value is smaller, the water quality is better. The classification bounds of Yao Zhiqi's evaluation method are: Ⅰ represents excellent, P ≤ 0.90. II represents good, P ≤ 1.51. III represents light pollution, P ≤ 2.63. IV represents moderate pollution, P ≤ 3.55. V represents heavy pollution, P > 3.55 [36].

Water Quality Prediction Method
The typical wetland water quality prediction flow chart of the source region of Yangtze River is shown in Figure 2.

Establishment of SWAT Model Database of Typical Wetlands in the Source Region of the Yangtze River
The pre-treated high-precision DEM map was used to extract the wetland river network in the source region of the Yangtze River, and the outlet locations were obtained according to the calculation. The typical wetland water system in the source region of the

Establishment of SWAT Model Database of Typical Wetlands in the Source Region of the Yangtze River
The pre-treated high-precision DEM map was used to extract the wetland river network in the source region of the Yangtze River, and the outlet locations were obtained according to the calculation. The typical wetland water system in the source region of the Yangtze River is divided into several sub-basins, and then each sub-basin is divided into different hydrological response units (HRUs) based on the minimum threshold of land use type/soil type area, soil type, and wetland runoff data.

Parameter Sensitivity Analysis and Calibration of the SWAT Model
The built-in module of the SWAT model was used to analyze the parameters of the runoff, and the method was LH-OAT. The LH-OAT method combines the Latin square (LH) sampling method with one-way (OAT) sensitivity analysis, providing both the robustness of LH sampling and accuracy of the OAT algorithm [37]. Sensitivity analysis was performed on the parameters of the established SWAT model. Sensitivity parameters include base flow attenuation coefficient (ALPHA-BF), shallow aquifer water storage threshold (GWQMN), soil evaporation compensation factor (ESCO), saturated water conductivity of soil layer (SOL-K), shallow groundwater REVAPMN, effective water conductivity of the main riverbed (CH-K2), and effective water capacity of the soil layer (SO) L-AWC), SCS runoff curve number (CN2) ( Table 3). Sensitivity analysis was performed to obtain the range of optimal parameters by SWAT-CUP calculation, and then according to the value range of the best parameter, the parameter range is continuously reduced by iterative analysis, and finally the rated value of the model (optimal parameter) is obtained (Table 4). According to the results of sensitivity analysis, the final rate values of each data are within the obtained value range, and the simulation effect is good. The optimal parameter values are brought into the SWAT model to verify and compare the actual measured values of the model.

Validation of the SWAT Model
The model has a rate period of 1961 to 2016 and a validation period of 2017 to 2019 ( Figure 3). The rate period R 2 is 0.84, indicating that the measured and simulated values in the rate period fit well, and the SWAT model simulates the net flow in the source region of the Yangtze River in a realistic way. The R 2 of the validation period was 0.89, indicating that the SWAT model has strong applicability and high reliability in the source region of the Yangtze River.
The model has a rate period of 1961 to 2016 and a validation period of 2017 to 2019 ( Figure 3). The rate period R 2 is 0.84, indicating that the measured and simulated values in the rate period fit well, and the SWAT model simulates the net flow in the source region of the Yangtze River in a realistic way. The R 2 of the validation period was 0.89, indicating that the SWAT model has strong applicability and high reliability in the source region of the Yangtze River.

Coupling of CMIP5 Global Climate Model and SWAT Hydrological Model
Three RCP (RCP2.6, RCP4.5, RCP8.5) greenhouse gas emission scenarios released by IPCC are adopted, in which RCP8.5 refers to the climate scenario with high CO2 emission, while RCP2.6 and RCP4.5, respectively, represent the climate scenario with low CO2 emission and medium and low CO2 emission [38]. Combined with the statistical downscaling climate data of CMIP5 global climate model from 2021 to 2100 released by the IPCC data Center, the simulation results of 21 CMIP5 global climate models (Table 5) were interpolated, and the calculated results were uniformly downscaled to the same resolution. Then, the simple average method was used to conduct multi-model aggregation. The monthly average data under the emission scenarios of RCP2.6, RCP4.5, and RCP8.5 were made. The weather data is then fed into the calibrated and validated SWAT model.
The statistical downscaling method was adopted to downscale the output results. The processed meteorological data and collected spatial, soil, hydrological, and water quality data of the source region are taken as the initial data of the determined SWAT hydrological model and imported into the model for calculation. The responses of CODMn, NH3-N, conductivity, TN, TP, temperature, and DO to CMIP5 global climate change under RCP2.6, RCP4.5 and RCP8.5 scenarios were further simulated to obtain future data of wetland water quality changes in the source region of the Yangtze River under different climate scenarios.  Three RCP (RCP2.6, RCP4.5, RCP8.5) greenhouse gas emission scenarios released by IPCC are adopted, in which RCP8.5 refers to the climate scenario with high CO 2 emission, while RCP2.6 and RCP4.5, respectively, represent the climate scenario with low CO 2 emission and medium and low CO 2 emission [38]. Combined with the statistical downscaling climate data of CMIP5 global climate model from 2021 to 2100 released by the IPCC data Center, the simulation results of 21 CMIP5 global climate models (Table 5) were interpolated, and the calculated results were uniformly downscaled to the same resolution. Then, the simple average method was used to conduct multi-model aggregation. The monthly average data under the emission scenarios of RCP2.6, RCP4.5, and RCP8.5 were made. The weather data is then fed into the calibrated and validated SWAT model. The statistical downscaling method was adopted to downscale the output results. The processed meteorological data and collected spatial, soil, hydrological, and water quality data of the source region are taken as the initial data of the determined SWAT hydrological model and imported into the model for calculation. The responses of CODMn, NH3-N, conductivity, TN, TP, temperature, and DO to CMIP5 global climate change under RCP2.6, RCP4.5 and RCP8.5 scenarios were further simulated to obtain future data of wetland water quality changes in the source region of the Yangtze River under different climate scenarios.

Data Processing
The changing trend chart of water quality data of nine typical wetlands in the source region of the Yangtze River from 2016 to 2021 was obtained by Origin 2021. Data such as COD Mn , NH 3 -N, TN, TP, conductivity, temperature, and DO in the SWAT hydrological model were imported into Excel. The trends of CODMn, NH 3 -N, conductivity, TN, TP, temperature, and DO in the forthcoming period (2022-2100) were plotted under three different climate scenarios (RCP2.6, RCP4.5, and RCP8.5).

Analysis of Water Quality of Typical Wetlands in the Source Region of the Yangtze River from 2016 to 2021
The COD Mn , NH 3 -And, and TN in typical wetland water bodies in the source region of the Yangtze River show an overall downward trend. The highest value of all parameters appears in 2016 and the lowest value of all parameters arise in 2021. Temperature and dissolved oxygen (DO) in the overall rise, its peak appear in 2021 (Figure 4). The TP concentration and water conductivity do not vary significantly. In general, water quality of typical wetlands in the source region of the Yangtze River has improved year by year. Figure 5 shows, among the changes in the year, the changes of DO, NH 3 -N, and TN concentration in the year are dry season > flat season > abundant season, while the COD Mn , temperature, and conductivity change in the opposite trend during the year, and the TP fluctuated slightly during the monitoring period. dissolved oxygen (DO) in the overall rise, its peak appear in 2021 (Figure 4). The TP con-centration and water conductivity do not vary significantly. In general, water quality of typical wetlands in the source region of the Yangtze River has improved year by year. Figure 5 shows, among the changes in the year, the changes of DO, NH3-N, and TN concentration in the year are dry season> flat season > abundant season, while the CODMn, temperature, and conductivity change in the opposite trend during the year, and the TP fluctuated slightly during the monitoring period.   centration and water conductivity do not vary significantly. In general, water quality of typical wetlands in the source region of the Yangtze River has improved year by year. Figure 5 shows, among the changes in the year, the changes of DO, NH3-N, and TN concentration in the year are dry season> flat season > abundant season, while the CODMn, temperature, and conductivity change in the opposite trend during the year, and the TP fluctuated slightly during the monitoring period.

Bayesian Water Quality Evaluation Results
The Bayesian evaluation method was used to evaluate the water quality indicators, and the evaluation results are shown in Table 6. From 2016 to 2021, the water quality at each monitoring point is above the Class II water standard. The water quality in 2016 and 2019 was mostly Class II water. The evaluation grade tends to improve in 2016 and 2017, but significantly declines in 2019. In 2019, except for the monitoring point 9, the water quality grade reached Class I; the remaining monitoring points were all listed as Class II. In 2020 and 2021, the water quality of most of the monitoring points had improved, and the water quality grade of most of the monitoring points had reached the Class I water quality standard. However, at the end of 2020, monitoring points 4, 6, and 7 were still classified as Class II water. 2021 is the best year for water quality evaluation, and all monitoring points are of first-class water quality. Being dependent on the above results, the water quality of typical wetlands in the source region of the Yangtze River fluctuates and rises during the six years from 2016 to 2021.

Yao Zhiqi Method Water Quality Evaluation Results
Depending on Yao Zhiqi's evaluation method, the higher the p value, the worse the water quality. As shown in Table 7, water quality is above class II, but P 2019 > P 2016 > P 2017 > P 2018 . The p value obtained in 2019 was larger than that in other years. The water quality in 2019 was the worst. In other years, the water quality of typical wetlands is class I, with P2021 > P2020. The water quality value in 2021 is the best. The water quality of typical wetlands from 2016 to 2021 was generally good, and the water quality rating from 2016 to 2019 was Class II, indicating good water quality. 2020 is the turning point. The water quality in 2019 is the worst, but from 2020 onwards, the water quality gradually improves and reaches Class I. In 2021, the water quality reaches the best, P 2021 = 0.62.  Figure 6a and Table 8, line (a) show the inter-annual variation trend of COD Mn in wetland water in the source region of the Yangtze River under three RCP climate scenarios from 2022 to 2100 and the linear relationship between different climate scenarios. The fluctuation range of COD Mn in wetland water is the smallest, which is 3.960 mg/L under RCP2.6, and the largest, which is 6.117 mg/L under RCP8.5, and the fluctuation range is RCP8.5 > RCP4.5 > RCP2. 6. In the future, COD Mn in wetland water will show an obvious linear growth trend. In the RCP2.6 scenario, COD Mn in wetland water is mainly characterized by inter-annual oscillation. In the RCP4.5 scenario, COD Mn in wetland water showed an upward trend at first, and a downward trend in the medium term. Under RCP8.5, the increase in COD Mn in wetland water was gentle and shows an overall upward trend. Under the three situations, the overall rising trend was RCP8.5 > RCP4.5 > RCP2.6, significant changes in COD Mn concentrations under three climate scenarios. Under the RCP2.6 scenario, COD Mn reaches the lowest value in 2023, and the maximum value appears for the first time in 2096; under the scenarios of RCP4.5 and RCP8.5, the minimum value of COD Mn appears in 2021 and the maximum value appears for the first time in 2097. COD Mn in the wetland water body in the source region will continue to increase, and the water quality will gradually deteriorate. der the three situations, the overall rising trend was RCP8.5 > RCP4.5 > RCP2.6, significant changes in CODMn concentrations under three climate scenarios. Under the RCP2.6 scenario, CODMn reaches the lowest value in 2023, and the maximum value appears for the first time in 2096; under the scenarios of RCP4.5 and RCP8.5, the minimum value of CODMn appears in 2021 and the maximum value appears for the first time in 2097. CODMn in the wetland water body in the source region will continue to increase, and the water quality will gradually deteriorate.    Note: RCP2.6 represents the lowest greenhouse gas emissions and emission concentrations; RCP4.5 represents the decline in the use of some non-renewable fossil energy sources, such as coal, and the significant reduction in greenhouse gas emissions due to the use of clean energy; RCP8.5 is the highest greenhouse gas emission concentration without climate policy intervention. Figure 6b and Table 8, line (b) show the climate models of NH 3 -N in the wetland water body in the source region of the Yangtze River and the linear relationship between different scenarios under three RCP climate model scenarios from 2022 to 2100. Under the three scenarios, the fluctuation range of NH 3 -N content in the wetland water body is the smallest, which is 0.462 mg/L under RCP2.6 and the largest, which is 0.886 mg/L under RCP4.5, that is, the fluctuation range is RCP4.5 > RCP8.5 > RCP2. 6.

Change Trend Analysis of NH 3 -N under Future Climate Scenarios
It can be observed in Figure 6b that NH 3 -N in wetland water in the source region of the Yangtze River shows an obvious linear growth trend in the near future. Under the RCP2.6 scenario, NH 3 -N in wetland water shows an upward trend in the early, middle and late stages, and the rising rate in the initial stage is greater than that in the middle and late stage, showing an overall upward trend. In the RCP4.5 scenario, NH 3 -N in wetland water showed an upward trend at the beginning of this century and a downward trend in the middle and late stage. Under the RCP8.5 scenario, the growth rate of NH 3 -N at the beginning of this century is slightly higher than that in the middle and late of this century, showing an overall upward trend. Under the RCP4.5 scenario, the upward trend of NH 3 -N is more obvious, while under the RCP2.6 scenario, the upward trend is more stable, that is, the overall upward trend is RCP4.5 > RCP8.5 > RCP2.6, with a significant difference (p < 0.05). In the RCP2.6 scenario, NH 3 -N will reach the lowest value in 2025; the highest value will appear for the first time in 2088 and again in 2099. In the RCP4.5 scenario, NH3-N will reach the lowest value in 2022 and the maximum value first will appear in 2099. In the RCP8.5 scenario, NH 3 -N will reach the lowest value in 2022, and the maximum value also will appear in 2099 for the first time. NH 3 -N in wetland water in the source region of the Yangtze River will continue to increase. Figure 6c and Table 8, line (c) indicate the inter-annual variation trend of DO in the wetland water body in the source region of the Yangtze River and the linear relationship between different scenarios under three RCP climate model scenarios from 2022 to 2100. The fluctuation ranges of the three climate scenarios are: RCP8.5 > RCP2.6 > RCP4. 5. In the future, the content of DO in wetland water in the source region of the Yangtze River will show an obvious linear downward trend. Under the RCP2.6 scenario, DO in the wetland water body shows a downward trend in the initial stage; the decline rate in the initial stage is greater than that in the middle stage, and there is an upward trend in the late stage, which shows a downward trend as a whole. Under the RCP4.5 scenario, DO in the wetland water showed a downward trend at the beginning, middle and late of this century. Under the RCP8.5 scenario, DO in the wetland water body increased at the beginning of this century and decreased in the middle and late of this century, while the decline in the middle is slightly larger, showing a downward trend as a whole. Under the three scenarios, the downward tendency of DO was RCP8.5 > RCP2.6 > RCP4.5, with a significant difference (p < 0.05). Under the RCP2.6 scenario, DO reaches the highest value in 2025 and the lowest value appears for the first time in 2077. Under the RCP4.5 scenario, DO reaches the highest value in 2022 and the lowest value appears for the first time in 2096. Under the RCP8.5 scenario, DO reaches the highest value in 2040, and the lowest value first will appear in 2098. DO in the wetland water body in the source region will still decline.

Analysis of the Change Trend of TN under Future Climate Scenarios
The inter-annual variation trend of TN in the wetland water body in the source region of the Yangtze River under the three RCPs climate model scenarios from 2022 to 2100 and the linear relationship between different scenarios are shown in Figure 6d and Table 8, line (d). The lowest fluctuation range of TN content in wetland water is 0.261 mg/L under RCP4.5, and the largest fluctuation range is 0.441 mg/L under RCP8.5.
In the future, TN in wetland water will show an obvious linear growth trend. Under the RCP2.6 scenario, TN in wetland the water body shows an upward trend in the early stage, a downward trend in the middle and late stage, and the decline range in the late stage is greater than that in the middle stage, and the overall trend is upward. Under the scenario of RCP4.5, TN in wetland water show an upward trend in the early and late stages of this century, and the upward range in the late stage is greater than that in the early stage. Under the RCP8.5 scenario, the growth rate of TN at the beginning of this century is slightly higher than that in the middle and late of this century, showing an overall upward trend. In the RCP8.5 scenario, the rising trend of TN is the most obvious. While in the RCP2.6 scenario, the rising trend is the gentlest. The overall rising trend is RCP8.5 > RCP4.5 > RCP2.6, with a significant difference (p < 0.05). Under the RCP2.6 scenario, TN reaches the lowest value in 2039, and the maximum value appears for the first time in 2096. Under the scenarios of RCP4.5 and RCP8.5, TN reached the lowest value in 2022, and the maximum value appeared for the first time in 2097. In the future, TN in the wetland water body in the source region will continue to rise.  between different scenarios under three RCP climate model scenarios from 2022 to 2100. The response of TP in the water body to climate change is complex. Under the three scenarios, the inter-annual fluctuation of TP in the wetland water body in the source region is different. The fluctuation range of TP content in the wetland water body is the smallest, which is 0.033 mg/L under RCP4.5 and the largest, which is 0.063 mg/L under RCP8.5, that is, the fluctuation range is RCP8.5 > RCP2.6 > RCP4.5.
It can be seen from Figure 6e that TP in wetland water showed an obvious linear growth trend at the beginning of the twenty-second century. Under the scenario of RCP2.6, TP in the wetland water body has an upward prospect in the initial stage, a downward prospect in the middle and late stage, and a greater decline in the late stage than in the middle stage, which is an upward prospect as a whole. Under the scenario of RCP4.5, TP in wetland water has the prospect of rising in the early and late stage of this century, and the rising range is large in the late stage. Under the RCP8.5 scenario, the growth rate of TP at the beginning of this century is slightly higher than that in the middle and late of this century. Under the RCP8.5 scenario, the rising trend of TP is obvious. While under the RCP2.6 scenario, the rising trend is relatively stable, that is, RCP8.5 > RCP4.5 > RCP2.6; the difference is very significant. In the RCP2.6 scenario, TP reaches the lowest value in 2038, and the maximum value appears for the first time in 2050. In the RCP4.5 and RCP8.5 scenarios, TP reaches the lowest value in 2022, and the maximum value appears for the first time in 2097. However, in the RCP4.5 scenario, the maximum value appears again in 2099 and 2100. In the future, TP in the wetland water body in the source region of the Yangtze River will continue to rise.

Analysis of Temperature Trends in Future Climate Scenarios
According to Figure 6f and the Table 8, line (f), the lowest fluctuation range of wetland water temperature is 0.35 • C under RCP8.5, and the largest range of fluctuation is about 4 • C according to RCP2.6.
In the future, the wetland water temperature will show a linear growth trend. Under the RCP2.6 scenario, the water temperature of the wetland shows an upward trend in the early and late stages, and a downward trend in the middle stage, and the increase in the early stage is slightly larger than that in the late stage, which is an upward trend on the whole. Under the RCP4.5 scenario, the water body temperature of wetland in the early, middle and late stages of this century shows an upward trend, with a large increase in the middle stage, and an overall upward trend. Under the scenario of RCP8.5, the rising range of wetland water body temperature in the middle of this century is slightly larger than that in the later and early stages, and generally shows an upward trend. The rising trend of wetland water temperature under RCP2.6 is the most obvious, and the rising trend of wetland water temperature under RCP8.5 is relatively stable, with a significant difference (p < 0.05). Under the RCP2.6 scenario, the wetland water body temperature reaches the lowest value in 2054, and the maximum value first appears in 2096. Under the RCP4.5 scenario, the wetland water body temperature reaches the lowest value in 2023, and the maximum value first appears in 2090. Under the RCP8.5 scenario, the wetland water body temperature reaches the lowest value in 2022 and the highest value in 2100. The water temperature of wetland will continue to increase, and the water quality may gradually deteriorate. In the future, the conductivity of wetland water will show a linear growth trend, but the growth trend is not obvious (p > 0.05). Under the RCP2.6 scenario, the conductivity of the wetland water body shows a downward trend in the initial and middle stages; the downward range in the initial stage is greater than that in the middle stage, and the upward trend in the late stage, which is an upward trend as a whole. Under the RCP4.5 scenario, the conductivity of the wetland water body showed an upward trend at the beginning, middle and late stages of this century, with a large increase in the middle term. Under the RCP8.5 scenario, the increase in wetland water conductivity in the middle of this century is slightly greater than that in the later period. But the increase is very small. In the early period, it is a downward trend and an overall upward trend. In the RCP 8.5 scenario, the conductivity of wetland water increased significantly, while in the RCP2.6 scenario, the upward trend was relatively stable, and the difference was significant (p < 0.05). Under the RCP 2.6 scenario, the conductivity of the wetland water body reaches the lowest value in 2022 and 2091, and the maximum value first appears in 2098. Under the RCP4.5 scenario, the conductivity of the wetland water body reaches the lowest value in 2026, and the maximum value first appears in 2062. This causes a jump in conductivity in 2060-2065 due to the rapid increase in water temperature in 2060-2065. Under the RCP8.5 scenario, the conductivity of the wetland water body reaches the lowest value in 2032 and maximum value in 2091. The water conductivity of the wetland will show an upward trend, but the growth trend is not obvious (p > 0.05), which has little impact on water quality.

Discussion
This study shows that the water quality of typical wetlands in the source region of the Yangtze River from 2016 to 2021 is excellent and basically not disturbed by human activities, which is consistent with the research results of Liu et al. [38] and Zhao et al. [39]. The source region of the Yangtze River is sparsely populated, with less agricultural land and a relatively low degree of industrialization. Relatively few pollutants enter the water body through precipitation impact, atmospheric sedimentation, or sewage discharge.
From 2016 to 2021, the annual variation of NH 3 -N, DO, TN, and TP concentrations in the wetland water body of the source region of the Yangtze River in the dry season (It refers to the period when the surface water flow in the basin is depleted and relies mainly on groundwater to recharge the water source.) > flat season (It refers to the period when the river is at normal water level.) > abundant season (It refers to the period when the river water flow is mainly recharged by rainfall or snowmelt.), which is basically consistent with the research results of Su [40] on the runoff variation characteristics and variation trend of the source region of the Yangtze River in nearly 60 years (1956-2012) and Yan [41] on the variation characteristics of wetlands in different source regions of the three rivers from 1975 to 2007. And the results are basically consistent with the consequences of Zhao [42] and others in the study of surface water quality in the source regions of the three rivers in Qinghai Province from 2005 to 2013. This may be due to the large rainfall in the wet season. Due to the influence of precipitation, there are relatively few pollutants discharged into the water by atmospheric sedimentation and sewage, and the water itself has a good self-purification capacity [43]. From 2016 to 2021, the water quality of the wetland is relatively good, which is related to the grassland enclosure, reduction of over grazing and strengthening of ecological environment protection in the source region of the Yangtze River in recent years.
The water quality of different monitoring points is evaluated by the Bayesian evaluation method, and the overall water quality of the source region of the Yangtze River is comprehensively evaluated by the Yao Zhiqi evaluation method [36]. The evaluation results intuitively and comprehensively show the change in water quality of typical wetlands in the source region of the Yangtze River and have received excellent results. This study also found that monitoring point 3 had most of the class II water quality from 2016 to 2021. Through analysis, it was found that monitoring point 3 had a low altitude (3500 m), was situated in a densely populated region, and was close to the highway. There were animals grazing around monitoring point 3. At the same time, there were wild animal activities in this region, which was greatly affected by human activities, resulting in slightly poor wetland water quality. The water quality of monitoring point 9 is mostly Class I water during the monitoring period, which was related to the fact that the monitoring point is located at the edge of Cocosili Nature Reserve. The monitoring point has a high altitude (4692 m), less interference from human activities, and less pollution from surrounding point sources and non-point sources, leading to better water quality of the wetland.
Hassanjabbar, et al. [44] used GCMS climate models and artificial neural networks (ANNs) to study the response of the flow regime and water quality indicators to climate change. The results show that global climate change will have a negative impact on water quality in the future. Liu [45] and Fang [46] simulated and analyzed the response of hydrology and water quality of the Changlejiang basin and Dahuofang basin to climate by using A2 and B2 emission scenarios and the SWAT hydrological model. Under the GCMS climate model, the concentration of TP increased and decreased, but overall, it showed an upward trend.
In this study, RCP2.6, RCP4.5, and RCP8.5 scenarios under the CMIP5 climate model are coupled with the SWAT hydrological model to predict the change trend in TN and TP concentrations in wetland water in the source region of the Yangtze River in the future (2022-2100). The results show that the concentrations of TN and TP in wetland water in the source region of the Yangtze River are increasing. It shows that future climate change will have a negative impact on water quality.
At present, relevant studies have used multi-mode sets to simulate and evaluate the future temperature or precipitation in the source region of the Yangtze River but there are generally 4-14 model sets, which are analyzed and compared under single or two future emission scenarios. At the same time, most of the existing studies focus on the middle of the 21st century, and the time scales of temperature and precipitation series are short, which makes the data simulation accuracy not high, resulting in a significant increase in the uncertainty of temperature or precipitation change prediction [47]. In this study, the simulation results of 21 models in the CMIP5 model were used to predict the change of wetland water quality in the source region of the Yangtze River, and more models were adopted to overcome the uncertainty of prediction caused by a single model or fewer models [48]; and the prediction results were more reliable, providing scientific support for the protection of wetlands in the source region of the Yangtze River under climate change.

Conclusions
(1) From 2016 to 2021, the concentrations of COD Mn , NH 3 -And, and TN in typical wetland water bodies in the source region of the Yangtze River showed a downward trend; temperature and dissolved oxygen concentration generally showed an upward trend. The inter-annual fluctuation range of water conductivity and total phosphorus (TP) is not obvious. The changes in COD Mn , temperature, and conductivity in the year are abundant season > flat season > dry season, and the changes of DO, NH 3 -And, and TN concentration in the year are dry season > flat season > abundant season. Overall, the water quality of typical wetlands in the source region of the Yangtze River is better in the wet season.
(2) According to the Bayesian evaluation method, the water quality of the nine monitoring stations in the source region of the Yangtze River is Class II in 2016 and 2019, and most of the stations in other years are Class I. According to Yao Zhiqi's evaluation method, the water quality is Class II from 2016 to 2019. The two evaluation methods show that the overall water quality has obviously improved since 2019.
(3) In the future period, TN, TP, COD Mn , NH 3 -N, and temperature in the wetland water in the source region of the Yangtze River will continue to rise; the changes of DO will continue to decline. In the forthcoming period (2022-2100), the wetland water quality in the source region of the Yangtze River will gradually deteriorate, and the wetland protection work is an urgent necessity.
(4) In the future period, the water quality of wetlands in the source region of the Yangtze River is gradually deteriorating, and wetland management becomes an urgent task. The source region of the Yangtze River should continuously improve the strategy to cope with climate change, strengthen the research capacity of ecological and environmental resources in the source region of the Yangtze River, and increase the control of greenhouse gas emissions; further, improve the information sharing mechanism and enhance the capacity of water resources optimization.