A Comprehensive Method for Water Environment Assessment considering Trends of Water Quality

Changjiang River Scientific Research Institute, Wuhan, Hubei 430010, China Key Lab of Basin Water Resource and Eco-Environmental Science, Hubei Province, Wuhan, Hubei 430010, China Hunan Provincial Key Laboratory of Key Technology on Hydropower Development, Changsha 410000, China South-to-North Water Diversion Project Midline Water Source Co.,Ltd., Shiyan 442700, China School of Civil and Hydraulic Engineering, Huazhong University of Science and Technology, Wuhan, Hubei 430010, China


Introduction
Sufficient access to freshwater resources is essential to the survival of human and natural ecosystems, and its role is inextricably coupled in economic sectors [1]. e actual state of water resources is a sensitive and critical issue in many countries, and the deterioration of water quality has become a global problem. An appraisal of available water resources and their quality is the first step to ensuring water security; thus, research on water quality evaluation has been an area of increasing interest and necessity [2].
Traditional methods for water quality assessment, which include measuring chemical parameters in water and normalized measurements against existing standards, are not sufficiently effective [3]. With the advantage of fully using water quality parameter information under a normative value, the water quality index (WQI), proposed by Horton [4] and Brown et al. [5], was widely applied to classify surface quality [1,[6][7][8][9]. WQI pairs concentrations of several water quality parameters such as dissolved oxygen (DO), pH, nitrate, phosphate, ammonia, chloride, hardness, and metals with associated weights [10]. e lower the WQI, the better the water quality, while a higher value indicates degraded water quality. e index value depends on monitored constituents and threshold values of specific domestic or local standards. In recent years, the WQI concept has evolved into other formulations including the National Sanitation Foundation WQI (NSFWQI) and the integrated water quality index (IWQI) [11,12].
With the advent of environmental monitoring networks and chemometrics, there have been associated increases in the application of multivariate statistical techniques, such as cluster analysis (CA), principal component analysis (PCA), and factor analysis (FA) [13][14][15][16]. ese techniques can reduce the dimensionality of a multivariate dataset while still maintaining the original structure, such as principal components or factors [14]. With the advantage of effective dimension reduction, these techniques have been frequently employed to classify water quality data and detect similarities among variables in many research studies [17], which help with the interpretation of results and in turn proffer a less subjective and more objective process.
Previous methods for assessing water quality in the environment were almost always based on historical data. When, in fact, the conveyance and storage systems were dynamic, water quality in the cleanest rivers could be in the process of deteriorating, while the opposite could be true for the most polluted rivers. A dynamic change in water quality could lead to inconsistencies between results of a one-time evaluation and reality.
In this paper, we proposed a water quality assessment method which takes into account water quality trends. First, a comprehensive pollution index (CPI) was calculated at several monitoring locations within a larger field site. Next, Spearman's rank correlation analysis was used to determine CPI trends at each location. Finally, principal component analysis (PCA) was deployed to assess water quality at the field scale using CPI trends as one of the evaluation parameters. e process of the method is illustrated in Section 2.4.

Comprehensive Pollution Index (CPI) Method.
e CPI method is one of the most common approaches used in China to qualitatively evaluate the water quality [18][19][20][21]. e following equation calculates the index P as where C i is the monitored value of water quality parameter i; n is the number of water quality parameters monitored; and S i is the standard permissible limit on water quality parameter i according to the environmental quality standards for surface waters (GB3838-2002) in China (Table 1) [22]. With this approach, P increases proportionally to the concentration increase of each constituent normalized by its water quality standard; hence, the value of P could be used to classify the degree of water contamination ( Table 2).

Spearman's Rank Correlation Analysis.
Spearman's rank correlation coefficient is defined as the Pearson correlation coefficient between two ranked random variables. It is used when data are not normally distributed between two variables. Spearman's rank correlation coefficient, rs, can be calculated with the following formula [23]: where d i � X i − Y i is the difference between the two ranked variables and n is the total number of observations. In its formulation, values of rs are constrained within −1 ≤ rs ≤ +1.
If rs � +1, it indicates a perfect positive correlation, while rs � −1 indicates a perfect negative correlation. In this approach, X i is the ranking sequence according to listed P values, and Y i is the ranking sequence according to the year of sampling. As rs approaches +1, the more obvious the uptrend in the value of P over time, hence indicating a progressive deterioration of water quality. Conversely, as rs approaches −1, the downtrend in P over time is more obvious as is the progressive improvement in water quality over time.

Principal Component Analysis.
In this framework, the main pollution sources were determined by principal component analysis (PCA).
is method is widely used because it can extract data and reduce system dimensionality. In PCA, the dataset contains P uncorrelated variables each of which is called a principal component (PC). e most important parameters are extracted in the first component (PC 1 ), while other parameters with relatively less importance fall within subsequent components (PC 2 , . . ., PC n ). Algebraically, for n original variables, x 1 , x 2 ,. . .,x n , the expression of the principal component PC i is determined as where a ii are the eigenvectors extracted from the covariance or correlation matrix and the number of components n is equal to the number of variables. As stated above, most of the information found within the n original variables is captured within the first few principal components, thus reducing the practical dimensionality of the dataset. Further details can be found elsewhere [14,24]. e score of each principal component can be obtained from where F i is the score of PC i ; Fac i is the value of factor i; and λ i is the eigenvalue of PC i , and i is the index of the PCs. e comprehensive score can be calculated using the following equation: where α i is the contribution rate of PC i and n is the number of PCs. e relative level of pollution at any given location can be determined by ranking the comprehensive scores of those sampled sections. e higher the comprehensive value, the lower the water quality.

e Process of the Method.
e process of the comprehensive method can be described as follows (Figure 1).
Step 1. Calculate the annual composite pollution index P for each monitoring location using equation (1) and at least 4 years of monitoring data to generate 4 yearly averages to be used in the next step.
Step 2. Use Spearman's rank correlation coefficient method to calculate CPI trends (rs) at each monitoring location using equation (2).
Step 3. Using PCA and the original water quality measurements of pH, DO, COD Mn , COD Cr , BOD 5 , NH 3 -N, TN, TP, and the rs coefficient at each sampled section, calculate and rank the scores of each section.
Step 4. e levels of pollution in all rivers are identified according to the score ranking.

Study Area.
e Danjiangkou Reservoir (32°20′-33°45′N, 110°40′-111°50′E) is located at the border of Hubei and Henan provinces, China, and has a surface area of 1050 km 2 when water reaches its normal level of 170 m [25,26]. It is the main drinking water source for the Middle Route of China's South-to-North Water Transfer Project (SNWTP). e success of the SNWDP depends on the water quality in the Danjiangkou Reservoir. ere are approximately 200 tributaries in the Danjiangkou Reservoir Catchment, of which 16 main tributaries make up 90% of the whole drainage area [25]. Since July 2016, water quality monitoring sections have been set up near the mouths of the 16 main tributaries, and the water quality samples were collected monthly. e location map showing the Danjiangkou Reservoir and monitoring sites (sections) is described in Figure 2 and Table 3.

3.2.
Annual P Value for Each River Section. Monitoring data from July 2016 to June 2020 were divided into four periods to calculate the average annual water quality concentrations at each river section. Water quality measurements of COD Mn , COD Cr , BOD 5 , NH 3 -N, TN, and TP were used in the calculation of annual P values for each period and at every river section using equation (1). Results are shown in Table 4. Table 4, Spearman's rank correlation coefficient, rs, was calculated using equation (2), which equates the slope of the CPI trendline listed in Table 5. Of the 16 major tributaries of the Danjiangkou Reservoir, 14 showed CPI trends of water quality improvement; moreover, the Shengding River and the Quyuan River showed continuous improvement trends. e two exceptions were the Jiangjun River and the Jianhe River.

Evaluation of Water Quality with Trends Considered.
ere are 9 original variables which were used in PCA. For each sampled tributary, the first eight variables were water quality indicators which included pH, DO, COD Mn , COD Cr , BOD 5 , NH 3 -N, TN, and TP, while the last variable was the CPI trend. Using PCA, the comprehensive scores of each tributary were obtained. Again, a higher comprehensive score equates to a tributary with lower water quality. PCA results were obtained using the commercial computer program SYSTAT version 19 from SPSS Inc. e 9 original variables were standardized (Table 6) and then subjected to PCA. e factor correlation coefficient considered significant is the one that is greater than 85.0%.
PCA on 9 parameters yielded two principal components explaining 88.3% of the sample variance. Table 7 shows component 1 describing 76.24% of the sample variance and component 2 capturing 12.09%. Figure 3 represents the 2D plot of both principal components. COD Mn , COD Cr , BOD 5 , NH 3 -N, TP, and TN had a positive influence on PC1, while pH and DO had a negative impact. As for PC2, it was dependent on rs coefficients, representing CPI trends.
According to the above analysis, scores from the first two principal components and the comprehensive score are obtained and shown in Table 8. As seen from this table, the pollution conditions in the tributaries were confirmed according to the comprehensive scores. When the CPI trends were considered, the three most polluted rivers were the Sihe River, Shending River, and Jianghe River, and the three least polluted were the Taohe River, Baihe River, and Qihe River, respectively.

Advances in Civil Engineering
Using PCA again, we repeated our evaluation without CPI trends and compared results to those presented in Table 8 in Figure 4. As shown in the figure, pollution rankings for most rivers were different with the exception of the two most polluted rivers. For example, the Jianghe River changed its pollution rating from 4th to 3rd when its rs   Advances in Civil Engineering coefficient was considered. e rs coefficient is greater than 0 which indicates that water quality was deteriorating. erefore, particular attention should be paid to this tributary when the local government plans to carry out pollution mitigation efforts.

Conclusions
In this paper, a comprehensive method for water environment assessment considering water quality trends was proposed and employed to analyze pollution in tributaries of the Danjiangkou Reservoir. A preliminary assessment was conducted using the CPI method and annual averaged water quality data. Using the annual P values from each tributary, CPI trends were quantified under Spearman's rank correlation analysis. CPI trends were shown to be important parameters in PCA and in the determination of comprehensive scores. By employing the comprehensive method developed herein, polluted tributaries were identified as candidates for future environmental management and pollution control.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors' Contributions
Sheng Bi, the corresponding author of this article, contributed to the main contents for the abstract, Sections Figure 4: e comparison results of pollution rankings considering rs coefficient or not.
Advances in Civil Engineering 7