Water Quality Evaluation and Variation Trend Analysis of Rivers Upstream of the Dahuofang Reservoir in China

: The Dahuofang Reservoir is one of the most important water sources in Liaoning Province, China, so it is critical to identify the status and evolution characteristics of its water quality. Six monitoring indicators were selected to analyze water quality differences and variation trends of each indicator in three inlet sections of the reservoir during different hydrological periods from 2003 to 2021, and an improved comprehensive pollution index method was proposed to study the pollution variation trends. The results showed three ﬁndings. (1) The water quality of the three rivers is better in high water periods than that in low water periods. (2) In terms of the spatial state of the water environment, water quality of the Hun River is the worst, the Suzi River is poor, and the She River is better. The worst indicator of the three rivers, TN (total nitrogen), has exceeded the standard for many years (Grade IV–Inferior Grade V). TP (total phosphorus) in the Hun River, which has deteriorated severely since 2013 and is positively correlated with rainfall; it is mainly inﬂuenced by pollution from agricultural activities. (3) The P value obtained by the improved method is lower than that of original method, which is mainly because TN is relatively stable, and the exceeding standard of TP is not as serious as TN. The improved method takes into account the interactions and ﬂuctuations of indicators, so that it can reﬂect the pollution situation more scientiﬁcally. These results are helpful to evaluate the pollution status of surface water. It is suggested that water be transferred appropriately to improve water quality and take necessary management measures to reduce TN and TP in the Hun River.


Introduction
Water is an important resource for human survival and development [1,2]. The "13th Five-Year Plan" has greatly improved China's poor water environment, with 83.4% of grade I~III surface water improving in 2020; however, this still falls short of the 85% target of the "14th Five-Year Plan". Therefore, it is still necessary to strengthen the protection of water quality in China. Water quality evaluation and variation trend analysis are helpful to propose targeted management measures [3][4][5]. As a water supply reservoir for the seven cities of Fushun, Shenyang, Liaoyang, Anshan, Yingkou, Panjin, and Dalian, in Liaoning Province, China, the quality of water in Dahuofang is crucial to the health of residents [6][7][8][9]. There is a correlation between the water quality of the reservoir and that of the rivers entering the reservoir [10,11]. Therefore, analyzing variation trends in the upstream rivers of the reservoir and scientifically evaluating their pollution are crucial to improving water quality.
At present, water quality analysis of rivers and reservoirs is usually carried out with the single factor evaluation method, the comprehensive pollution index method, the principal component analysis method, the grey relational analysis method, or the The water quality data of three inlet sections (the Beizamu section of the Hun River, the Gulou section of the Suzi River, and the Taigou section of the She River), from 23 April 2003 to 30 December 2021, were obtained from the Dahuofang Reservoir Administration and the Fushun Sub-Center of the Liaoning Water Environment Monitoring Center. From 2003 to 2009, the Dahuofang Reservoir Administration evenly set six sampling points on the sections, collected six bottles of 2000 mL three times a month, and sent them to the laboratory to detect indicators, such as biochemical oxygen demand (BOD5), total phosphorus (TP), total nitrogen (TN), and so on. Since 2009, unmanned monitoring stations have been set up in the sections, equipped with automatic samplers, online monitoring, and other equipment, which can automatically sample and transmit data to the Fushun Sub-Center of the Liaoning Water Environment Monitoring Center every four hours (6 times a day) through a public network VPN. In order to ensure the stability of the data and avoid fluctuations in water quality caused by inter-annual variations in rainfall, which may affect the trend analysis of the data, six monitoring indicators during the low (May) and high (September) water periods were selected for this study. Analytical indicators included dissolved oxygen (DO) and biochemical oxygen demand (BOD5), which are important indicators of the self-purification capacity of water bodies [36,37]; the permanganate index (CODMn), which is a common indicator of the pollution of organic and inorganic oxidizable substances in water bodies [38,39]; total phosphorus (TP), which is the main indicator for controlling eutrophication in water bodies [40,41]; total nitrogen (TN), which is an indicator of the degree of nutrient pollution in water bodies [42]; and ammonia nitrogen (NH3-H), which is an indicator of the level of inorganic nutrients in the water body [43,44]. The water quality data of three inlet sections (the Beizamu section of the Hun River, the Gulou section of the Suzi River, and the Taigou section of the She River), from 23 April 2003 to 30 December 2021, were obtained from the Dahuofang Reservoir Administration and the Fushun Sub-Center of the Liaoning Water Environment Monitoring Center. From 2003 to 2009, the Dahuofang Reservoir Administration evenly set six sampling points on the sections, collected six bottles of 2000 mL three times a month, and sent them to the laboratory to detect indicators, such as biochemical oxygen demand (BOD 5 ), total phosphorus (TP), total nitrogen (TN), and so on. Since 2009, unmanned monitoring stations have been set up in the sections, equipped with automatic samplers, online monitoring, and other equipment, which can automatically sample and transmit data to the Fushun Sub-Center of the Liaoning Water Environment Monitoring Center every four hours (6 times a day) through a public network VPN. In order to ensure the stability of the data and avoid fluctuations in water quality caused by inter-annual variations in rainfall, which may affect the trend analysis of the data, six monitoring indicators during the low (May) and high (September) water periods were selected for this study. Analytical indicators included dissolved oxygen (DO) and biochemical oxygen demand (BOD 5 ), which are important indicators of the self-purification capacity of water bodies [36,37]; the permanganate index (COD Mn ), which is a common indicator of the pollution of organic and inorganic oxidizable substances in water bodies [38,39]; total phosphorus (TP), which is the main indicator for controlling eutrophication in water bodies [40,41]; total nitrogen (TN), which is an indicator of the degree of nutrient pollution in water bodies [42]; and ammonia nitrogen (NH 3 -H), which is an indicator of the level of inorganic nutrients in the water body [43,44].

Water Quality Evaluation Standards
The surface water quality classification standards of China (GB3838-2002), shown in Table 1, were used in this paper to evaluate the classification of each index in different hydrological periods of the three inlet sections shown in Figure 1 [45].

Comprehensive Pollution Index Method
The comprehensive pollution index method took the ratio of the actual monitoring concentration of the index to the evaluation standard as the pollution index (P i ), and then weighted it to obtain the comprehensive pollution index (P). As a source of drinking water, the water quality of Dahuofang Reservoir should be maintained above grade II. Therefore, this paper used grade II in Table 1 as the evaluation standard, and adopted the comprehensive pollution level of its upstream sections. The larger the P value, the more serious the pollution, while the smaller the P value, the better the water quality. The DO was the opposite. The specific pollution evaluation criteria are shown in Table 2 [18,46]. The equation to calculate P is presented as follows: where P is the comprehensive pollution index, P i is the pollution index of water quality factor i, C i is the measured concentration of water quality factor i, S i is the standard limit of water quality factor i in the "classification standards of China" (GB 3838-2002), and n is the total number of indexes.

Improved Method
Exceeding multiples and concentration fluctuations are the main factors that affect water quality. In the comprehensive pollution index method, P i considers the index exceeding multiples, but each indicator is weighted equally, which ignores the importance of differences between pollutants. This paper redistributed the weights of each index according to the dimensionless coefficient of variation, which meets the unit requirements. In addition, this is a physical quantity that measures the dispersion of data, and indicators with greater variability have a greater impact on water quality assessment; thus, weighting based on the coefficient of variation can be considered a more scientific reflection of the state of water pollution.
The calculation formula is as follows: where C vi is the coefficient of variation of water quality factor i, σ is the standard deviation of water quality factor i, X is the mean value of water quality factor i, W i is the weight of water quality factor i, and the other letters have the same meanings as defined previously.

Statistical Test Methods
In order to test whether the variances of water quality indicators were equal, the normal distribution test and the homogeneity of variance test were carried out. In this paper, the number of samples was less than 2000; thus, Shapiro-Wilk was used for the distribution test. At present, there are mainly Cochran, Hartley, Bartellet, and Levene tests for homogeneity of variance in statistics [47]; among them, only the Levene test could be used for both normally distributed populations, non-normally distributed populations, or unknown populations [48], and there were two indicators with skewed distributions; thus, the Levene test was used in this paper. The steps of the Levene test are as follows: (i) Test the hypothesis.
(ii) Convert the original data to the new variable values.
where Y ij is the original data, and Y i· is the median of the i-th sample in the original data.
(iii) Calculate the test statistic value W.
where k is the number of sample groups, N i is the content of the i-th sample, N is the sum of the contents of each sample, and Z ·· is the total mean of all the data.
where the degrees of freedom v 1 and v 2 are k − 1 and N − k, respectively.
(v) Principle of judgment. The Levene test statistic obeys the F distribution with degrees of freedom v 1 and v 2 . When W ≥ (Fα, k − 1, N − k), then P ≤ α, reject H 0 and accept H 1 ; it can be considered that the variance of each sample is uneven; when W < F(α, k − 1, N − k), then P > α, H 0 is not rejected, and the variance can be considered to be homogenous.

Statistical Analysis of Water Quality in Different Hydrological Periods
According to the test results of Shapiro-Wilk in Figure 2, based on the yearly mean values of the six indicators in dry and wet periods, the fitting curves of BOD 5 (sig = 0.065 > 0.05), COD Mn (sig = 0.055 > 0.05), TN (sig = 0.2 > 0.05), and DO (sig = 0.2 > 0.05) were completely symmetrical, which fitted the normal distribution. However, frequency peaks of NH 3 -H and TP were all shifted to the left, showing a positively skewed distribution. Therefore, the Levene test, which had weaker requirements for normality than others, was used to determine whether the variances were homogeneous. The results showed that the significance was less than 0.05, thus the nonparametric test (Mann-Whitney test) was used to test the correlation.
where the degrees of freedom and are − 1 and − , respectively.
(v) Principle of judgment.
The Levene test statistic obeys the F distribution with degrees of freedom and . . When ， − 1， − , then , reject and accept ; it can be considered that the variance of each sample is uneven; when is not rejected, and the variance can be considered to be homogenous.

Statistical Analysis of Water Quality in Different Hydrological Periods
According to the test results of Shapiro-Wilk in Figure 2, based on the yearly mean values of the six indicators in dry and wet periods, the fitting curves of BOD5 (sig = 0.065 > 0.05), CODMn (sig = 0.055 > 0.05), TN (sig = 0.2 > 0.05), and DO (sig = 0.2 > 0.05) were completely symmetrical, which fitted the normal distribution. However, frequency peaks of NH3-H and TP were all shifted to the left, showing a positively skewed distribution. Therefore, the Levene test, which had weaker requirements for normality than others, was used to determine whether the variances were homogeneous. The results showed that the significance was less than 0.05, thus the nonparametric test (Mann-Whitney test) was used to test the correlation.  The results of the Mann-Whitney test and the data differences between the six monitoring indicators, of the three river sections during the low and high water periods, are shown in Tables 3-5. Accordingly, the P values of the DO in Beizamu, the DO in Gulou, and the TN in Taigou section were less than 0.05, which had statistical significance in different hydrological periods. The following was a statistical description of the raw data for each indicator. The maximum value, minimum value, mean value, and standard deviation of DO in the Beizamu section during the high water period (September) were significantly larger compared to the low water period (May); the mean values of TP and TN were larger than those in the low water period, with the mean value of TN in the high water period found to be 18 times larger than that in the low water period, while BOD 5 and COD Mn showed the opposite, and NH 3 -H was less different ( Table 3). The water level in the Beizamu section during the high water period is usually about 1.5 m higher than that during the low water period. The high water level and turbulence of the river are conducive to atmospheric oxygen dissolving in the water, which may increase the DO content, improve the self-purification ability of the water body, and reduce BOD 5 and COD Mn . The area of arable land in Qingyuan County upstream of the Beizamu section is 39,400 hectares, with an annual fertilizer application of about 18,000 t. Rainfall runoff during the abundant water period may wash excess waste into the river, resulting in increased TN and TP concentrations. This statistic showed that the DO value increases during the high water period, which can improve the self-purification ability of the water body to a certain extent, but the increases in TN and TP during the abundant water period are also significant, which may lead to the deterioration of the water body if measures to address this are not taken in the future.
The maximum value, minimum value, and standard deviation of DO in the Gulou section during the high water period were significantly increased compared to the low water period, the mean values of BOD 5 , COD Mn , and TN were reduced compared to the dry water period, and the differences in the remaining two indicators were slight ( Table 4). The upstream confluence flow of the Gulou section is much smaller than that of the Beizamu, and the pollution caused by arable land is limited. The flow velocity does not increase significantly during the high water period, thus the oxygen in the water flow has less opportunity to contact the water. Additionally, the rainfall is diluted, making the mean DO value slightly smaller during the low water period. Table 5 shows that the maximum value, mean value, and standard deviation of DO and TN during the high water period were slightly smaller than those during the low water period, while the differences in the mean values of the other four indicators were not significant. The cultivated area of Taigou is the smallest of the three sections at only 4776.5 hm 2 , and the flow rate is less than one third of that of the Gulou section, which may make the variation of each indicator insignificant.       There was a clear trend of interannual variation in TP, and the extent of its exceedance can be divided into two stages. The pre-study period, from 2003 to 2012, was basically satisfactory, and the post-study stage, from 2013 to 2021, exceeded the standard. The maximum exceedance multiple for TP was 7.12 mg/L in 2013, which may have been caused by heavy rainfall in that year. Multi-year exceedance multiples and exceedance rates were as follows: high water period (1.54 times, 53%) > low water period (1.8 times, 64%) (Figure 4). Figure 5 shows the relationship between rainfall and TP exceedance. It can be seen that TP exceedance was serious in the years when rainfall was high, which may have been caused by large-scale rainfall; in turn, this may have enhanced phosphorus substance accumulation in the soil due to factory discharge, agricultural fertilization, and farming to enter the river, all of which led to TP exceedance. In addition, the exceedance of the TP standard was more serious during the low water period than during the high water period, which showed that the TP pollution in the Hun River is heavily influenced by surface pollution.

Trends in the Various Water Quality Indicator Categories
substance accumulation in the soil due to factory discharge, agricultural fertil farming to enter the river, all of which led to TP exceedance. In addition, the of the TP standard was more serious during the low water period than duri water period, which showed that the TP pollution in the Hun River is heavily by surface pollution. During the low water period, it can be seen from Figure 6a that TN of the Gulou section was the worst; it was Grade V to Inferior Grade V for many years. TP showed a trend of first getting worse and then becoming better, and all indicators reached the standard after 2018. BOD 5  During the low water period, it can be seen from Figure 7a that TN was most severe in the Taigou section, which was Grade V to Inferior Grade V. TP showed a trend of exceedance in 2016, and in the last two years. NH 3 -H and COD Mn basically met the standard for many years, and the interannual variation characteristics were not significant. BOD 5 was the best indicator, maintaining Grade I for many years. From Figure 7b, DO, NH 3 -H, COD Mn , and BOD 5 showed no significant trends and were largely compliant during the high water period. TN deteriorated from Grade V at the beginning of the study to Grade V to Inferior Grade V by its end. TP had the same repeated fluctuations as during the low water period, and could be divided into three stages overall: poor and exceeding the standard from 2003 to 2008, better and maintained as Grade I from 2009 to 2015, deteriorating since 2016, and becoming Grade V and Inferior Grade V in 2020 and 2021, respectively.
Water 2022, 14, x FOR PEER REVIEW CODMn met the standard except for when it was Grade IV in 2020. NH3-H, BOD all reached the standard over the years, without significant changes. During the low water period, it can be seen from Figure 7a that TN was m in the Taigou section, which was Grade V to Inferior Grade V. TP showed a tr ceedance in 2016, and in the last two years. NH3-H and CODMn basically met the for many years, and the interannual variation characteristics were not signific was the best indicator, maintaining Grade I for many years. From Figure 7b, DO CODMn, and BOD5 showed no significant trends and were largely compliant d high water period. TN deteriorated from Grade V at the beginning of the study V to Inferior Grade V by its end. TP had the same repeated fluctuations as durin water period, and could be divided into three stages overall: poor and exce standard from 2003 to 2008, better and maintained as Grade I from 2009 to 2015 rating since 2016, and becoming Grade V and Inferior Grade V in 2020 and 202 tively.

Analysis of Comprehensive Pollution Index Results
The results of p for the three river sections under the original and improved methods are shown in Figures 8-10. The original weights and improved weig three sections are shown in Table 6. The mean p-values for the three sections fro 2021 under the two methods are shown in Table 7.

Analysis of Comprehensive Pollution Index Results
The results of P for the three river sections under the original and improved analysis methods are shown in Figures 8-10. The original weights and improved weights of the three sections are shown in Table 6. The mean P-values for the three sections from 2003 to 2021 under the two methods are shown in Table 7.

Analysis of Comprehensive Pollution Index Results
The results of p for the three river sections under the original and improved analysis methods are shown in Figures 8-10. The original weights and improved weights of the three sections are shown in Table 6. The mean p-values for the three sections from 2003 to 2021 under the two methods are shown in Table 7.      It can be seen from Figure 8a that the P value of the Beizamu section in the low water period was "multi-peak", and the water quality fluctuates greatly; this can be divided into three stages. The first stage was from 2003 to 2009, and the water quality was basically qualified; the second stage was from 2010 to 2017 when the P value was mostly greater than 1 and the water quality was affected by serious pollution; and the third stage was from 2018 to 2021 when the water quality was basically qualified except for minor pollution in 2020. The original method showed that 68% of the sections of Beizamu were polluted and 32% were heavily polluted during the low water period. However, the improved method showed that 42% were basically qualified and 58% were polluted. According to Figure 8b, the P value in the high water period showed a "double peak", which was larger in 2013 and 2020, of 1.242 > 1 and 1.458 > 1, respectively, all of which were polluted. The original method showed that 84% of the sections were polluted and 16% were heavily polluted; the improved method showed that 84% of the sections were basically qualified while 16% were polluted.
From Figure 9a, the original method showed 84% of the river sections were polluted and 16% heavily polluted during the low water period in the Gulou section, while the improved method showed 63% were basically qualified, 32% polluted, and 5% heavily polluted. Figure 9b shows that the P value of the Gulou section fluctuated upwards during the high water period, with a maximum of 2.041 > 2 in 2020, which was heavily polluted, and that the pollution status could be divided into two stages. The water quality in 2003 to 2012 was in the basic qualified state, except for 2007 when it was lightly polluted, and in 2013 to 2021 when it was basically qualified and heavily polluted. The original method showed that the Gulou section was 95% polluted and 5% heavily polluted during the high water period, while the improved method showed 53% basically qualified, 42% polluted, and 5% heavily polluted.
As can be seen from Figure 10a, the water quality of the Taigou section during the low water period of the study was more stable in the early stage, and more volatile in the later stage, with the maximum integrated pollution index occurring in 2021 at 3.336 > 2, indicating that it was heavily polluted. Water quality pollution could be divided into two phases, with better water quality from 2003 to 2015, mostly in the basically qualified and lightly polluted states, and from 2016 to 2021 when water quality became worse, and was mainly polluted to heavily polluted. The original method showed that 95% of the Taigou section was polluted and 5% was heavily polluted during the low water period, while the improved method showed 73% basically qualified, 16% polluted, and 11% heavily polluted. According to Figure 10b, the P value fluctuated greatly during the high water period, and the maximum value of 4.017 > 2 occurred in 2021, when it was heavily polluted. The original method showed 16% basically qualified, 73% polluted, and 11% heavily polluted during the high water period in the Taigou section, while the improved method showed 67% basically qualified, 21% polluted, and 11% heavily polluted.
From Table 6, it can be seen that the fluctuations of the indicators for the three sections are as follows: NH 3 -H > TP > BOD 5 > COD Mn > TN > DO. The minimum weights of NH 3 -H and TP in the improved method were 0.278 and greater than 0.167, respectively; BOD 5 was similar to the original method, and the weights of COD Mn , TN, and DO were lower than those of the original method. The original method results in Table 7 showed that the multi-year average comprehensive pollution index of the three sections was between 0.902 and 1.759, indicating polluted and basically qualified, respectively. However, in the improved method, the Gulou and Taigou in the wet season were 0.91 and 0.92, respectively, which were less than 1 and basically qualified.

Discussion
The average P values of the three sections during the high water period under both methods were less than those in the low water period. For instance, the average P value of the improved method in the high water period (1.08) was less than that in the low water period (1.119) in the Beizamu section ( Table 7). The water quality in the high water level season was better than that in the low water level season, which may be due to the dilution effect of large flow during the high water level season. This finding is consistent with the study of Honghu Lake (2021) by Chen et al. [49], the seasonal variation characteristics analysis of water quality in Huntai Basin explored (2021) by Li Yao et al. [50], and a study on the water quality of the Segamat River in Malaysia (2018) by Faridah et al. [51]. However, it is not very consistent with the conclusions of Choi et al. [52] on the water quality evaluation of the Jinwi River watershed in South Korea in 2021. This may be because the upper Jinwei River is an urban area with relatively high organic matter and nutrient load compared with the Dahuofang Reservoir. Therefore, the large flow during the wet season will increase pollutant loads and worsen the water quality. To sum up, increasing the flow volume of the upstream river can improve the inflow water quality to a certain extent. Therefore, strengthening the hydrological connectivity of the river, in addition to appropriate water diversion and transfer, are important measures to improve the water quality of the river.
In the comprehensive pollution index method, the exceeding multiple of each pollutant is reflected by the pollution index P i , and the weight of each water quality factor is the same; The Levene test showed that the fluctuations of each index were different, and the improved comprehensive pollution index method determines the corresponding weight according to the variation coefficient of each water quality factor. The greater the coefficient of variation, the greater the weight. For the most serious pollution index TN and the second serious index TP, the weight changed significantly (see Table 6) before and after improvement of the analysis method. TN index exceeds the standard greatly, but the fluctuation is small; thus, the calculated weight (0.16) is less than the original weight (0.167). The TP index is just the opposite, with small exceedance but large fluctuation; thus, the calculated weight (0.43) is greater than the original weight (0.167), making the final calculated P value lower than the original method. In addition, the P value trend calculated by the improved method is consistent with the original method, and can reflect that the water quality in the wet season is better than that in the dry season. To sum up, it is more scientific and reasonable to determine the weight based on the coefficient of variation.
Compared with the original method, the improved method can reflect the pollution level of rivers in different periods in detail. For example, the comprehensive pollution index of the Gulou section in the low water period (1.002, polluted) is greater than that in the high water period (0.930, basically qualified), which shows that measures such as water diversion and transfer can improve water quality. Compared with the principal component analysis method, grey correlation analysis method, neural network method and other methods, the improved method has a simpler calculation process and stronger operability. Therefore, the improved method can more specifically analyze the local water quality, help to provide more targeted environmental management policies and suggestions, and play a positive role in the restoration of the water environment.
The calculation results of the comprehensive pollution index method show that the pollution index in the wet season is ranked from large to small: Beizamu (1.708) > Gulou (1.483) > Taigou (1.421), and the pollution index in the dry season is Beizamu (1.759) > Gulou (1.540) > Taigou (1.507). The improved method shows that the pollution index during the high water period is Beimagi (1.08) > Gulou (0.93) > Taigou (0.902), while during the low water period the index is Beizamu (1.119) > Gulou (1.002) > Taigou (0.999). The three rivers are ranked from high to low in pollution degree: Hun River > Suzi River > She River, which may be closely related to the natural conditions and industrial structure of the basin. This finding is similar to the water quality evolution trend of the Dahuofang Reservoir analyzed by Yang Jinwei et al. [32] in 2016, and the water quality evolution results of inflow rivers analyzed by Wang Qiong et al. [53] in 2019. Hun River Basin has a population of 345,000, while Suzi River Basin and She River Basin have populations of 308,000 and 38,381 respectively. Therefore, Hun River Basin has the largest domestic sewage discharge. The cultivated land area of Suzi River is 51,000 hectares, while that of Hun River and She River are 39,400 hectares and 4777 hectares, respectively. Therefore, the pollution of Suzi River may come from chemical fertilizers and pesticides to a great extent. There are some industrial enterprises in the upper reaches of Hun River that may cause point source pollution, such as Fushun Hongtoushan Mining Co., Ltd., Fushun Longye Chemical Nanzamu Co., Ltd., etc., which are involved in metal smelting and chemical discharge. Although a series of environmental remediation measures have been taken, and enterprises are required to discharge wastewater after reaching the standard, the water quality is still not optimistic. There are paper mills and food factories in the Suzi River Basin, but their scale is far smaller than that of the mining plants in the Hun River Basin; therefore, the pollution of the Hun River is the most serious. She River Basin has the smallest population, the smallest cultivated land area, and no polluting enterprises; thus, the water quality is the best, with basically closed-circuit circulation and no sewage outflow.
TN concentration (Grade IV to Inferior Grade V) in Hun River is mainly influenced by point source pollution, and TP concentration (1.055-1.800 times higher than the standard) is mainly affected by agricultural non-point source pollution, which is positively correlated with rainfall. This is similar to the research results of Pang et al. [54] in 2015 and Beizamu et al. [55] in 2018 on the relationship between agricultural land and rivers.

Conclusions
This paper proposed, for the first time, a comprehensive pollution index method based on the coefficient of variation, to analyze the variation characteristics and pollution state of water quality in the upper reaches of the Dahuofang Reservoir. Our study resulted in four findings. (1) The improved method based on the coefficient of variation is reasonable. The improved method comprehensively considers the influence of pollution index fluctuations and the interactions between evaluation indexes, so that the weight coefficient can more scientifically reflect the degree of influence of each index on water quality. (2) The water quality of the three sections in the upper reaches of the Dahuofang Reservoir during the high water period is better than that in the low water period; thus, water quality supervision during the low water period should be strengthened in the future. (3) The pollution level was as follows: Hun River > Suzi River > She River. Strengthening the treatment of domestic sewage and industrial sewage from mines in the Hun River, and controlling the application of chemical fertilizers on arable land in the Suzi River basin would be useful. (4) TP has deteriorated since 2013 in the Hun River and continues to exceed the standard from Grade III to Grade IV. This is mainly a result of fertilization and pesticide application in farmland; thus, the use of chemical fertilizers and pesticide should be limited in the future else the water quality of the Hun River will become worse.
This study can provide reference and guidance for water quality management of the Hun River Basin, and has broad application prospects. However, the research also has some shortcomings. Future research could include the following two aspects: (1) adding more indicators to better reflect pollution status, especially heavy metal indicators such as iron and copper; and (2) formulating a variety of management measures, and estimating the treatment effect of each measure, with particular attention to TN and TP.