Abstract
As the water source for the middle route of the South-to-North Water Transfer Project, the Han River in China plays a role of the world’s largest inter-basin water transfer project. However, this human-interfered area has suffered from over-standard pollution emission and water blooms in recent years, which necessitates urgent awareness at both national and provincial scales. To perform a comprehensive analysis of the water quality condition of this study area, we apply both the water quality index (WQI) and minimal WQI (WQImin) methods to investigate the spatiotemporal variation characteristics of water quality. The results show that 8 parameters consisting of permanganate index (PI), chemical oxygen demand (COD), total phosphorus (TP), fluoride (F-), arsenic (As), plumbum (Pb), copper (Cu), and zinc (Zn) have significant discrepancy in spatial scales, and the study basin also has a seasonal variation pattern with the lowest WQI values in summer and autumn. Moreover, compared to the traditional WQI, the WQImin model, with the assistance of stepwise linear regression analysis, could exhibit more accurate explanation with the coefficient of determination (R2) and percentage error (PE) values being 0.895 and 5.515%, respectively. The proposed framework is of great importance to improve the spatiotemporal recognition of water quality patterns and further helps develop efficient water management strategies at a reduced cost.
Similar content being viewed by others
Introduction
High-quality water resources are of crucial importance to maintain ecological integrity and promote sustainable socio-economic development1,2. However, water quality issues have been more intricate than ever before due to different ecosystem pressures from rapid urbanization and population explosion. Water deterioration has also become a public urgent concern worldwide and a serious threat to people’s livelihoods3,4. Taking China for illustration, there were 874 records of local water contamination during 2006 ~ 2016, resulting in great troubles for domestic life and making millions of economic loss5. Reasonable evaluation of water quality variation has been demonstrated as a practical tool for water quality warning and protection6,7,8.
Indeed, lots of efforts have been made to investigate reliable information of water quality variation9,10,11. The water quality assessment based on hydro-chemical monitoring methods and data sampling underlies the geological heterogeneity of water quality, the understanding of human activities and the control of water contamination12,13. There have been considerable research for water quality evaluation in recent years, including the single-factor evaluation method14, Namerow pollution method15 and water quality index (WQI) method16,17. Of them, the recent WQI method takes full advantage of water quality parameter information and converts all parameters into a clear normative status value of water quality. It has been a prevailing approach for water quality assessment in a wide range of studies. The calculation of WQI has also experienced modification and has been developed in different ways. For instance, the National Sanitation Foundation WQI (NSFWQI) that comprises 9 parameters has been applied in different regions across the world18. Bascaron WQI19 and Canadian WQI20 that resemble the NSFWQI also had wide applications in the various background.
However, the downside of inflexible water quality parameters limits the application of aforesaid different WQIs. Some studies are devoted to identifying the key water quality parameters and developing more efficient WQI methods at a low-cost level. The minimum WQI (WQImin) model consisting of the key water quality parameters that can deal with the information redundancy and high cost has been demonstrated effective. A highly linear correlation relationship between WQI and WQImin has been documented21,22, indicating that WQImin has a strong potential to reflect the variation of water quality economically. However, the results of WQImin models may show substantial differences using various parameters and assigned weights17,23. The parameters selection for a WQImin model requires a deliberate response in terms of the specific scientific issue. The stepwise linear regression analysis has been verified robust to identify critical parameters for reduction of data redundancy14,24,25, here, it is employed in this study.
As the largest tributary of the Yangtze River in China, the Han River has been subject to water blooms on several occasions. Since the first diatom bloom of the Han River in 1992, there have been 9 water bloom episodes of varying degrees during the period 1998 ~ 2016, all of which concentratedly occurred between February and April. The influenced river reached about 500 km above the estuary in severe cases, and it usually lasted as long as 20 days26. Some literatures have reported that water blooms are closely associated with water quality27,28. Meanwhile, the implementation of China’s Middle Route of South-to-North Water Transfer Project will enhance the water supply capacity in the north while reducing the downstream flow, which may further exacerbate the water quality deterioration of the downstream environment29,30. The unqualified river flow has a direct crash on the drinking water safety in riverine cities such as Wuhan, which has received intensive attention from the local government. Most of the previous studies regarding the middle and lower reaches of the Han River basin typically focused on the optimization of cascade reservoirs operation31 or water resources allocation32, the water quality has come to stand out in recent years. These studies aimed at analyzing the correlation between water bloom and water quality within a short period, while the temporal and spatial variations in water quality in different times and river sections were less investigated. Furthermore, only a few studies targeting assessing the water quality status but the single factor evaluation method or the Nemerow's pollution index was adopted as the evaluation method33. The WQI method had not been applied to Han River, let alone the minimal WQImin method. Consequently, we still lack a thorough and comprehensive understanding of the water quality in this high-profile area. The temporal and spatial variations in water quality, the application of WQI method and the development of the WQImin model consisting of the key water quality parameters need further investigation.
To bridge these gaps, we examined 15 representative water quality parameters collected from 11 monitoring stations during 2015–2017 for investigation of the water quality level. The main goals of this study are (1) to analyze the spatiotemporal variation of each water quality parameter in the study area; (2) to comprehensively evaluate the water quality status as well as the spatial profiles and seasonal patterns using the water quality index method; and (3) to identify the critical parameters and develop a WQImin model for more efficient and cost-effective water quality evaluation.
Materials and methods
Study area
Located in the middle of the Yangtze River economical belt, Han River is the largest tributary of the Yangtze River, China. It is usually segmented into three divisions, i.e., upstream from its source to the Danjiangkou Reservoir, middle stream from the Danjiangkou Reservoir to the Jingzhou City, and downstream from Jingzhou to the Wuhan City. The total length of middle and low streams is 676 km, with a drainage area of 64,000 km2. Characterized by the subtropical monsoon, the river basin has abundant annual average precipitation from 700 to 1800 mm which concentrates in the summer and autumn seasons.
With huge water resources potential, the area plays a significant role in the provincial granary, industry and national “one belt, one road” construction. Eleven water quality stations are located along the mainstream to monitor the water quality changes. More detailed information is presented in Fig. 1 and further provided in the supplementary material (Table S1).
Sample measurement and data collection
All archived data in this study were provided by the Hubei Provincial Academy of Eco-environmental Sciences. Water samples were collected on both sunny and cloudy days to eliminate the effect of precipitation. They were sampled monthly from January 2015 to December 2017 in 11 stations, spanning temporal and spatial variation. The Standard Methods for the Examination of Water and Wastewater34 were used for sample chemistry analyses. There were 15 water quality parameters in total, including pH, dissolved oxygen (DO, mg/L), permanganate index (PI, mg/L), chemical oxygen demand (COD, mg/L), five-day biochemical oxygen demand (BOD5, mg/L), ammonia nitrogen (NH3-N, mg/L), total phosphorus (TP, mg/L), fluoride (F-, mg/L), selenium (Se, μg/L), arsenic (As, μg/L), sulfide (mg/L), plumbum (Pb, μg/L), copper (Cu, μg/L), zinc(Zn, μg/L) and mercury (Hg, μg/L). The specific approach for collecting water samples can be found in the standard of the Guidance on Sampling Techniques35. All the samples were labeled with detailed information using waterproof markers on the bottles to prevent misdiagnosis. Additionally, procedure blank was also used at all the stations to control the accuracy of analyses. With respect to the measurement method, the Hydrolab Datasonde 5 Sensor (USA) was calibrated prior to sampling to measure pH and DO. Except that, titration assembly, UV spectrophotometer (UV 2450), Ion chromatograph system (ICS 2000) and other instruments were also employed for different water quality parameters. More detailed information concerning the instrumental and chemical analysis method could be found on the website of the Ministry of Ecology and Environment of China (http://www.mee.gov.cn/ywgz/fgbz/bz/).
Water quality index
We used a weighted sum of all fifteen water quality parameters to calculate the WQI, which can be expressed as follows.
where \(n\) is the total number of water quality parameters, \(C_{i}\) and \(P_{i}\) are the normalized value and assigned weight of ith parameter, respectively. All weights range between 1 (least impact) and 4 (highest impact on water quality), and the assigned weights listed in Table S2 for relevant water quality parameters have been suggested by previous literature36,37,38. Note that we referred to The Environmental Quality for Surface Water39 to obtain a normalized value of \(C_{i}\) for more accurate evaluation, which was shown as follows:
where \(T_{i}\) is the measured concentration of ith parameter; \(S_{i,k}\) and \(S_{i,k + n}\) are the standard thresholds of the ith parameter at level \(k\) and level \((k + n)\), respectively; \(I_{i,k}\) is the standard normalization value of the parameter level, i.e., 20, 40, 60, 80 and 100; \(m\) is the number of identical values of the threshold, and \(m\) is equal to 1 if there is no same threshold. It is worth noting that pH is a special parameter because it has no specific standard threshold, then \(C\) value is set to 100 when \(6 \le {pH} \le 9\), otherwise it is 0.
WQI value ranges from 0 to 100 and can be classified into five different types as follows: excellent (91–100), good (71–90), medium (51–70), and bad (26–50) and very bad (0–25)19. A larger WQI value indicates a better water quality condition. Particularly, the annual period in this study is divided into spring (March to May), summer (June to August), autumn (September to November), and winter seasons (December to the following February). However, the traditional WQI model involves too many parameters and has much uncertainty17,40. An improved WQImin model by identifying the key parameters is developed with the benchmark of WQI, and both weighted and non-weighted WQImin models are considered for comparison. The weighted WQImin (WQImin-w) model can be calculated by Eq. (1), while the non-weighted WQImin (WQImin-nw) model is calculated by Eq. (3).
With reference to the data split procedure in previous studies in Nong et al.24, Wu et al.14 and Uddin et al.41, the WQImin models in this study were established using the following steps: The WQI and Ci for each station in 2015 and 2016 were used as “training data” to calibrate the WQImin model while a test period in 2017 was used to verify the model performance.
Data analysis
The Mann–Kendall (M–K) test has been widely applied to analyze water quality trends in previous studies24,42. The calculation processes of the M–K test were shown in detail in Güçlü (2020). The results of the M–K test have illustrated the trends of water quality parameters in this basin and are shown in Fig. S1. The one-way analysis of variance (ANOVA) was applied to verify the spatial differences of parameters. The WQImin models in this study were established using the following steps: (1) The WQI and Ci for each station in 2015 and 2016 were used as “training data” to select the key parameters for the WQImin model; (2) The coefficient of determination (R2) was taken as the goodness-of-fit criterion and the Percentage Error (PE) was used to evaluate the forecasting precision of the WQImin models based on the “testing data” (i.e., the WQI and the Ci for each station in 2017)44. The data used in the stepwise multiple linear regression method were pre-processed by the log transformation (i.e., lg(x + 1)).
Transparency
The authors confirm that the manuscript is an honest, accurate, and transparent account of the study was reported; that no vital features of the study have been omitted; and that any discrepancy from the study as planned have been explained.
Results
Water quality characteristics
Table 1 presents the statistical summary of all water quality parameters in the middle and lower reaches of the Han River basin and Fig. 2 displays their station concentrations.
Biochemical and physicochemical parameters (pH, DO, PI, COD, BOD5, and sulfide)
The annual mean pH values were greater than 7.80, and both the maximum and minimum measured pH values occurred in 2016, which were 8.80 and 6.60, respectively. For the annual mean pH, the highest value occurred in BJW station (pH = 8.22), closely followed by LHZ (pH = 8.21). The ZD, HZ, SD and HNC stations obviously observed much lower values than other stations, particularly, the HNC occupied the lowest annual mean pH value of 7.48. The M–K test showed that ZD, HZ, SD and HNC stations with relatively lower pH values had significant upward trends.
The annual mean DO concentrations were higher than the Class I standard (7.5 mg/L). The highest DO concentration of 13.40 mg/L was observed in 2015, while the lowest was in 2016 with the value of 5.9 mg/L. The annual mean DO concentrations increased monotonically from 2015 to 2017. For all stations, LHZ ranked first in terms of the annual mean DO concentration with the value of 9.73 mg/L, while SD and HNC were at the bottom. Surprisingly, XG and ZF observed relatively higher DO concentration compared with the upstream stations. This might be attributed to the influence of the water temperature. Lower river water temperature contributed to the higher DO concentration, vice versa. However, the one-way ANOVA (p < 0.01) indicated that there was no spatial difference of DO. The results of the M–K test showed there was only one station, i.e., LHZ, having a downward trend of DO while the remaining (about 91%) exhibited unchanged. Hence, the DO concentration kept relatively stable in all stations during the periods.
PI, COD, and BOD5 are all the key parameters for measuring water pollution levels arising from organic compounds. The statistical summary showed the annual mean PI concentration of each year was greater than 2 mg/L, which implied that PI observation couldn’t meet the standard of Class I. The years from 2015 to 2017 observed the maximal PI concentration of 4.6 mg/L, 4.2 mg/L, and 4.4 mg/L, respectively. Figure 2c showed that SW had the lowest mean PI of 1.83 mg/L and the downstream stations had higher PI concentrations than the upstream in general. The last 4 stations all exhibited the mean PI concentration over 2.5 mg/L while the mean concentration of other stations was lower than that value. With respect to COD, the maximum COD of 21.30 mg/L was in 2016 and the minimum was observed the same in each year. The annual mean COD concentration in each year was observed below 11.0 mg/L and less than Class I. Nine consecutive stations showed annual mean COD concentration larger than 7.5 mg/L and the downstream concentration were generally higher than the upstream except ZK station. The M–K test indicated that there were significant upward trends for HNC, XG and ZG stations, which were to be controlled in the following years. In terms of BOD5, all the monthly BOD5 in all stations were less than 3.0 mg/L and satisfied the threshold value of Class I. The maximum BOD5 in each year were almost the same. On spatial scale, SD, HNC and ZK stations showed three highest annual mean BOD5 concentrations, respectively, while the lowest occurred in ZD, which implied that SD, HNC and ZK stations exhibited relatively higher organic pollution.
The annual mean sulfide concentrations were far below the threshold of Class I. They were only about 10% of the threshold value (0.05 mg/L). The maximum sulfide occurred in 2017 with a value of 0.09 mg/L. The area had significant spatial differences in sulfide. Four consecutive stations from LHZ to HNC had higher sulfide content than the remaining 7 stations. About 72.7% of stations were analyzed by the M–K test to have significant downward trends for sulfide.
Nutrients and soluble ions (NH3-N, TP, F-)
The annual mean NH3-N concentration in each year exceeded the threshold value of Class I-0.15 mg/L. The maximum content was observed in 2016 and the minimum was in 2015, with the values of 0.60 mg/L and 0.02 mg/L, respectively. Figure 2f showed that the consecutive 9 stations exhibited the annual mean NH3-N concentration over 0.15 mg/L and only 2 downstream stations met the water quality requirement of Class I. The one-way ANOVA test revealed that no spatial differences were shown from upstream to downstream. But the study area was exposed to heavy pollution caused by TP. The statistical results demonstrated that the TP content of each station was far over the Class I threshold. The maximum TP content was 0.18 mg/L in 2017, followed by 0.15 mg/L in 2015. Figure 2g showed that the area has spatial heterogeneity of TP (one-way ANOVA, P < 0.01) and the TP content increased gradually from upstream to downstream. The spatial distribution of TP might be closely correlated with the phosphate industry in Hubei Province and the accumulation of the domestic sewage discharge, the use of fertilizer and pesticide45. Furthermore, the TP concentration was closely related to the algal blooms, which should receive much attention.
The soluble ion, F-, showed a small variation ranging from 0.246 mg/L to 0.263 mg/L. Compared to TP, the water quality condition for F- was much better. The maximum measured F- concentration in each year from 2015 to 2017 was 0.45 mg/L, 0.48 mg/L and 0.57 mg/L, respectively. Although there was spatial heterogeneity with an increasing trend, all the stations were excellent in terms of F-. The M–K test (Fig. S1) indicated that BJW and HNC stations had a significant downward trend.
Heavy metal parameters (Se, As, Pb, Cu, Zn and Hg)
The maximum annual mean concentration occurred in 2016 for Se, As, Cu, Zn and Hg with the values of 0.458 μg/L, 2.602 μg/L, 2.096 μg/L, 10.154 μg/L, 25.308 μg/L and 0.038 μg/L, respectively. The remaining Pb parameter had its maximum in 2015. Among the 6 parameters, Se, As and Hg performed best, followed by Pb, Zn and Cu in order. Table 1 showed the maximum value of Pb in 2017, Zn in 2015 and 2016, and Cu in 2015, 2016 and 2017, which were all over the threshold of Class I. The mean concentrations of Se and Hg had some fluctuation yet no obvious spatial changes (one-way ANOVA, P < 0.01). The mean concentrations of As increased from upstream to downstream in general except for higher values of ZD and HZ. As for Zn, LHZ had the minimum annual mean concentration and its downstream had higher concentration than upstream. The M–K test showed that there were 6, 2, 2, 6, 5, and 8 stations showing upward trends for Se, As, Pb, Cu, Zn and Hg, respectively.
Water quality assessment using the WQI method
The seasonal and spatial patterns of the water quality variations were presented in Fig. 3. In general, the water quality in our study area can be classified into “good” or “excellent” status in most cases, with all average WQI values greater than 87. Figure 3a displayed the seasonal variation of WQI during the monitoring period. The 4 highest seasonal average WQIs occurred in the winter of 2014, the spring and summer of 2015, and the winter of 2017, with values of 90.33, 90.28, 90.13, and 90.21, respectively. The lowest seasonal WQI was in the summer of 2016, indicating the worst water quality condition compared to other seasons. However, it could still be categorized to be “good” water quality status. The seasonal variation of the WQI in 2016 behaved more dramatically than that in 2015 and 2017. In 2015 and 2017, both the seasonal WQI values decreased from spring to summer, followed by a drop to the bottom in autumn and then rose again. In 2016, the lowest WQI values occurred in the summer. The WQI values descended from spring to summer to reach the lowest and then rose to 88.33 in autumn and continued to rise in winter. The lowest WQI in 2016 summer could be tightly associated with the flood during that period46. The persistent heavy rain caused severe floods in Han River basin and masses of non-point pollutants entered the channel and contaminated the river. Therefore, a lower WQI occurred in that time.
The spatial profile was displayed in Fig. 3b. The maximum average WQI was in SW and the minimum was in SD and HNC. WQI had certain spatial differences and the upstream stations occupied higher WQI than downstream. It can be attributed to the contaminant accumulation from domesticity, production, and other sources. In addition, YJH owed the most stable water quality condition yet HZ and ZG experienced larger fluctuations.
The training and test performance of WQImin models
The results of the stepwise multiple linear regression were shown in Table S3. It showed Zn contributed most to WQI according to the training performance, i.e., Model 1, R2 = 0.411, P < 0.001; PI, NH3-N, TP, DO were recognized in sequence and the R2 values increased monotonically until up to 0.857. For models with more parameters, such as models 6–8, R2 values further improved. To comprehensively evaluate the performance of the WQImin models for assessing water quality, all the models in Table S3 were considered for further investigation. R2 and PE were the two selected evaluation criteria for the model performance evaluation. The results indicated that except model 3, other R2 values improved with the increasing number of parameters. Compared model 2 with model 3, an extra NH3-N would decrease R2 values but increase PE. Table 2 also showed that the WQImin-nw models had higher R2 values and lower PE values than the WQImin-w models, indicating that the former could better explain the WQI variations. Moreover, significant differences between WQI and WQImin existed in the training period as shown in Fig. 4, indicating that the first models were not suitable for WQI simulation, while the remaining models exposed their potential. The average and median of WQImin-nw7 model were close to that of the WQI model, which accounted for 90% of the WQI variation with the lowest PE of 2.64%. It was identified as the most suitable WQImin model in this study.
Regarding the performance of WQImin models in the test period (in Fig. 5), the WQImin-nw7 model still behaved well. It had an R2 value of 0.895 and a low PE of 5.515%. WQImin-w5/WQImin-nw5 model performed poorly with lower R2 yet higher PE; WQImin-nw6 model behaved well than WQImin-w6, WQImin-w7 or WQImin-nw7 models in terms of R2 but had a similar PE situation with WQImin-nw5; WQImin-w8 and WQImin-nw8 models showed slightly larger R2 values than WQImin-nw7 model whereas PE also increased a lot, nearly 46.84% (from 2.515% to 3.693%) and 30.93% (From 2.515% to 3.293%), respectively. Overall, WQImin-nw7 model could explain nearly 90% of the variations of WQI. The result demonstrated that WQImin-nw7 model had a powerful prediction ability, which was consistent with the training performance.
Discussion
Key parameter selection for WQImin model
According to the stepwise multiple linear regression, the WQImin model in this study consists of seven main parameters: Zn, PI, NH3-N, TP, DO, Pb and Cu. The WQImin explains the variations of WQI well, and is more efficient for water quality assessment. Zn was the first parameter chosen for the regression, and it contributed most to the WQI evaluation in the training dataset (R2 > 0.40, P < 0.001). PI was the second significant sign for the variation of WQI, which represented the organic pollution of the water body. NH3-N and TP were the third and fourth parameters, respectively, on behalf of the nutritional parameters of water quality. Figure 2 also depicted that the Han River basin was threatened by a significant spatial variance of NH3-N and TP. In fact, previous literature has declared the impact of NH3-N and TP on water quality25,47. DO was the fifth parameter that could substantially improve R2. It is the main force of shaping the aquatic environment and biochemical process. It was also a crucial parameter appearing in the WQImin model in other works. Pb and Cu were introduced into the model last, both reflected the heavy metals concentration and affected water quality. Heavy metals in the water body could also harm human health. More attention should be put on these two substances since they sometimes fail to reach the standard of Class I of water quality.
Except for the universal multiple linear regression, other methods can be employed and are found in some relevant studies. Hou et al.48 used principal component analysis (PCA) to identify the explanatory parameters for WQI variation. Additionally, a linear correlation analysis between parameters was also meaningful. It can be employed for reducing the number of parameters and selecting the key parameters for analysis. The Pearson’s linear correlation of physicochemical parameters was shown in Fig. 6. The colors signified the positive or negative correlations between parameters and the areas occupied by the clockwise rotation showed the strong or weak correlations. The results illustrated that pH had a strong correlation with many parameters except for Se, Cu, and Hg. DO is positively correlated with NH3-N but negatively with PI, Zn, and Hg. Positive correlations were identified between PI and COD, BOD5, TP, F-, As, sulfide, Cu, Zn. BOD5 exhibited positive correlations with sulfide, Cu and Zn but negative correlations with Se, and Hg. More elaborate results can be found in Fig. 6, which can make sense for the selection of key parameters.
Impact sources and variation of water quality
Water quality evaluation can provide insights for the water resources exploitation and contamination control. In this study, we mainly utilized the measured data to unclose the water quality condition and developed a new method for evaluating the water quality. However, human activities, such as agricultural non-point source discharge, industrial and domestic sewage discharge, and anthropogenic intervention, could have a direct impact on the WQI performance. Some previous studies have found and reported this phenomenon. For example, Liu et al.45 found that farmland and urbanization could deteriorate water quality in the Han River basin. Tian et al.49 demonstrated that Lihu Lake suffer from the worst pollution during 1997–2003 due to excessive TN and TP. Although these parameters have improved significantly after artificial management, the last two years have witnessed a slight rebound. The overall evaluation revealed that the water quality could be taken as “good” or “excellent” levels in most cases and the upstream stations presented better water quality than downstream with some parameters like PI and NH3-N exceeded the thresholds, which implied that more effort can be made for the downstream water quality protection and the crucial water parameters needs to be constrained in these years. Moreover, the South-to-North Water Transfer Project and Yangtze-to-Han Water Transfer Project in recent decades also altered the hydrodynamic condition of the Han River to some extent and further affect water quality parameters. Han River has been confronted with emerging issues in these years and what exactly drives the change of the water quality might alter both in space and time still needs further investigation. For instance, how the operation of complex water transfer projects will impact the variation of the discharge and contamination concentration in middle and lower Han River is a hotspot issue50. Building a coupled hydrodynamic and water quality simulation model could contribute to the prediction of the contamination concentration and unveil the spatio-temporal variation. Furthermore, the simulation results that provide data sources under different hydrology or hydrodynamic condition can be coupled with the WQI and WQImin model constructed in this study for achieving comprehensive water quality evaluation. The model coupled the physics-based hydrodynamic-water quality model and the statistical WQI model could shed light on the source of contamination and provide precise preventions and control measures in the changing environments. What exactly drives the change of the water quality in different sections and in different times still needs investigation. Due to data scarcity of wastewater discharge, urban expansion, and hydrodynamics, more results regarding impact factors remains explored to assist water quality management. The coupled model that can unveil the source and variation of the contamination and assess the comprehensive water quality under complex conditions is to be built in the further study.
Uncertainty in water quality evaluation
Water quality evaluation is disturbed by multiple uncertainty sources, such as water sampling, measurement variability, water quality standard and water quality parameters. For example, both the weather and sampling time determine the DO content, and DO in aquatic ecosystems occupied the highest weight in water quality evaluation25. In this study, the water quality samplings were carried out on sunny or cloudy days in most cases to eliminate precipitation disturbance. Nevertheless, it was difficult to require all the samples following a strict schedule in such a huge study area. The inevitable laboratory measurement uncertainty might influence the confidence in WQI evaluations. Furthermore, the uncertainty was also induced by the water quality standards17,51. The classification standards vary in different districts and in different protection objectives and thus introduce some uncertainties in the assessment. In this study, the normalization values reconcile with the Surface Water Quality Standards GB3838-200239, which conformed to the official guidance in China However, this standard system might not work for other regions or other protection objectives. The WQI evaluation results could be altered with the change of the standards. Another source of uncertainty is the assigned weights of water quality parameters17. In terms of the previous studies concerning the water quality evaluation, the weight allocated to each parameter showed large variation52,53. In some WQI aggregations, it is impractical to use a unit weight for each parameter due to their different levels of impact on the water status. The high concentration of a parameter with high weight could translate to a low WQI value resulting in a misunderstanding of the water quality condition. Therefore, a proper approach for assigning the weights is essential. The weights used in this study were adopted and revised from Pesce et al.19 and Sun et al.54, which has also been verified in other studies. Whereas the uncertainty of WQI weights has not been investigated in this study, which is out of the main scope of the research. To further consider the overall uncertainty in WQI evaluation in the next steps, the statistical uncertainties with respect to the water quality parameters can be investigated using Monte-Carlo simulation17.
Conclusions
The water quality condition of the middle and lower Han River basin in China from 2015 to 2017 and the spatio-temporal variation of 15 water quality parameters were investigated and analyzed in this study. The water quality from seasonal and spatial scales was comprehensively evaluated with the WQI method and the key parameters were selected to develop the WQImin model. The main conclusions are as follows:
-
(1)
Eight parameters, i.e., PI, COD, TP, F-, As, Pb, Cu and Zn, performed obvious spatial discrepancy. The physicochemical and nutrient parameters, except for PI and NH3-N, satisfied the standard of Class I during the period 2015–2017. The heavy metal parameters in the middle and lower reaches of the Han River basin stayed at a low level, while Cu, Zn, and Pb exceeded the threshold.
-
(2)
The water qualities of this study basin were evaluated as “good” and “excellent” in most cases, with the mean WQI values of stations and seasons varying from 87.07 to 92.33 and from 87.43 to 90.33, respectively. On a spatial scale, the upstream stations presented better water quality than downstream; on a temporal scale, the autumn season was found the appearance time of worst water quality in 2015 and 2017 and the summer season was replaced in 2016.
-
(3)
The WQImin model developed in this study consisted of seven key parameters, i.e., Zn, PI, NH3-N, TP, DO, Pb and Cu. It has more physical explanatory and better evaluation performance than WQI for water quality evaluation.
Given the study area to be the strategic water source of China and the core of the Han River ecological economic belt, more efforts on linkage between water quality and United Nations-Water Sustainable Development Goal 6 are necessary and recommended in the future.
Data availability
The datasets used and/or analysed during the current study available from the corresponding authors on reasonable request.
Code availability
The code applications used in this work can be obtained by contacting the corresponding authors.
References
Ghotbizadeh, M. et al. Spatiotemporal variations of total and dissolved trace elements and their distributions amongst major colloidal forms along and across the lower Athabasca River. J. Hydrol. Reg. Stud. 40, 101029 (2022).
Sin, Y. & Lee, H. Changes in hydrology, water quality, and algal blooms in a freshwater system impounded with engineered structures in a temperate monsoon river estuary. J. Hydrol. Reg. Stud. 32, 100744 (2020).
Brkić, Ž, Kuhta, M., Larva, O. & Gottstein, S. Groundwater and connected ecosystems: An overview of groundwater body status assessment in Croatia. Environ. Sci. Eur. 31, 75 (2019).
Stehle, S. & Schulz, R. Agricultural insecticides threaten surface waters at the global scale. Proc. Natl. Acad. Sci. 112, 5750–5755 (2015).
Xu, J. et al. Spatial-temporal distribution and evolutionary characteristics of water environment sudden pollution incidents in China from 2006 to 2018. Sci. Total Environ. 801, 149677 (2021).
De Andrade Costa, D., Soares de Azevedo, J. P., dos Santos, M. A. & dos Santos Facchetti Vinhaes Assumpção, R. Water quality assessment based on multivariate statistics and water quality index of a strategic river in the Brazilian Atlantic forest. Sci. Rep. 10, 22038 (2020).
Saber, M. et al. An integrated assessment approach for fossil groundwater quality and crop water requirements in the El-Kharga Oasis, Western Desert, Egypt. J. Hydrol. Reg. Stud. 40, 101016 (2022).
Zhang, M., Wang, L., Mu, C. & Huang, X. Water quality change and pollution source accounting of Licun River under long-term governance. Sci. Rep. 12, 2779 (2022).
Chen, W. et al. Assessment of extrinsic and intrinsic influences on water quality variation in subtropical agricultural multipond systems. Environ. Pollut. 276, 116689 (2021).
Kämäri, M., Tattari, S., Lotsari, E., Koskiaho, J. & Lloyd, C. E. M. High-frequency monitoring reveals seasonal and event-scale water quality variation in a temporally frozen river. J. Hydrol. 564, 619–639 (2018).
Wang, J. et al. Response of bacterial communities to variation in water quality and physicochemical conditions in a river-reservoir system. Glob. Ecol. Conserv. 27, e01541 (2021).
Liu, Y. & Mao, D. Integrated assessment of water quality characteristics and ecological compensation in the Xiangjiang River, south-central China. Ecol. Indic. 110, 105922 (2020).
Singh, K. P., Malik, A. & Sinha, S. Water quality assessment and apportionment of pollution sources of Gomti river (India) using multivariate statistical techniques—a case study. Anal. Chim. Acta 538, 355–374 (2005).
Wu, T., Wang, S., Su, B., Wu, H. & Wang, G. Understanding the water quality change of the Yilong Lake based on comprehensive assessment methods. Ecol. Indic. 126, 107714 (2021).
Tenodi, S. et al. Assessment of the environmental impact of sanitary and unsanitary parts of a municipal solid waste landfill. J. Environ. Manage. 258, 110019 (2020).
dos Santos Simoes, F., Moreira, A. B., Bisinoti, M. C., Gimenez, S. M. N. & Yabe, M. J. S. Water quality index as a simple indicator of aquaculture effects on aquatic bodies. Ecol. Indic. 8(5), 476–484 (2008).
Seifi, A., Dehghani, M. & Singh, V. P. Uncertainty analysis of water quality index (WQI) for groundwater quality evaluation: Application of Monte-Carlo method for weight allocation. Ecol. Indic. 117, 106653 (2020).
Hoseinzadeh, E., Khorsandi, H., Wei, C. & Alipour, M. Evaluation of Aydughmush river water quality using the national sanitation foundation water quality index (NSFWQI), river pollution index (RPI), and forestry water quality index (FWQI). Desalin. Water Treat. 54, 2994–3002 (2015).
Pesce, S. F. & Wunderlin, D. A. Use of water quality indices to verify the impact of Córdoba city (Argentina) on Suquı́a River. Water Res. 34, 2915–2926 (2000).
Ramos, M. A. G., Bueno de Oliveira, E. S., Pião, A. C. S., Nalinde Oliveira Leite, D. A. & de Franceschi de Angelis, D. Water quality index (WQI) of Jaguari and Atibaia rivers in the region of Paulínia, São Paulo, Brazil. Environ. Monit. Assess. 188, 263 (2016).
Akkoyunlu, A. & Akiner, M. E. Pollution evaluation in streams using water quality indices: A case study from Turkey’s Sapanca lake basin. Ecol. Indic. 18, 501–511 (2012).
Sánchez, E. et al. Use of the water quality index and dissolved oxygen deficit as simple indicators of watersheds pollution. Ecol. Indic. 7, 315–328 (2007).
Monferrán, M. V., Galanti, L. N., Bonansea, R. I., Amé, M. V. & Wunderlin, D. A. Integrated survey of water pollution in the Suquía River basin (Córdoba, Argentina). J. Environ. Monit. 13, 398–409 (2011).
Nong, X., Shao, D., Zhong, H. & Liang, J. Evaluation of water quality in the south-to-north water diversion project of China using the water quality index (WQI) method. Water Res. 178, 115781 (2020).
Wu, Z., Wang, X., Chen, Y., Cai, Y. & Deng, J. Assessing river water quality using water quality index in lake Taihu basin China. Sci. Total Environ. 612, 914–922 (2018).
Cheng, B., Zhang, Y., Xia, R., Zhang, N. & Zhang, X. Spatio-temporal variations of water quality in Hanjiang River and its influencing factors in recent years. Environ. Sci. 1–15 (2021) (in Chinese).
Wang, J.-H. et al. Meteorological factors and water quality changes of Plateau Lake Dianchi in China (1990–2015) and their joint influences on cyanobacterial blooms. Sci. Total Environ. 665, 406–418 (2019).
Xin, X. et al. Algal blooms in the middle and lower Han River: Characteristics, early warning and prevention. Sci. Total Environ. 706, 135293 (2020).
Kuo, Y.-M., Liu, W., Zhao, E., Li, R. & Muñoz-Carpena, R. Water quality variability in the middle and down streams of Han River under the influence of the middle route of south-north water diversion project, China. J. Hydrol. 569, 218–229 (2019).
Zhou, Y., Guo, S., Hong, X. & Chang, F.-J. Systematic impact assessment on inter-basin water transfer projects of the Hanjiang River basin in China. J. Hydrol. 553, 584–595 (2017).
He, S. et al. Multi-objective operation of cascade reservoirs based on short-term ensemble streamflow prediction. J. Hydrol. 610, 127936 (2022).
Liu, D. et al. Assessing the effects of adaptation measures on optimal water resources allocation under varied water availability conditions. J. Hydrol. 556, 759–774 (2018).
Song, F. et al. Pollution risk assessment of heavy metals in the sediments of upstream Hanjiang River, China. Earth Sci. Inform. 14, 655–668 (2021).
Ministry of Environmental Protection, China. Standard Methods for the Examination of Water and Wastewater (Version 4). (China Environmental Science Press, 2002) (in Chinese).
Ministry of Ecology and Environment, China. Water quality-Guidance on sampling techniques. (2009) (in Chinese).
Debels, P., Figueroa, R., Urrutia, R., Barra, R. & Niell, X. Evaluation of water quality in the Chillán River (Central Chile) using physicochemical parameters and a modified water quality index. Environ. Monit. Assess. 110, 301–322 (2005).
Koçer, M. A. T. & Sevgili, H. Parameters selection for water quality index in the assessment of the environmental impacts of land-based trout farms. Ecol. Indic. 36, 672–681 (2014).
Zhao, P., Tang, X., Tang, J. & Wang, C. Assessing water quality of three Gorges Reservoir, China, over a five-year period from 2006 to 2011. Water Resour. Manag. 27, 4545–4558 (2013).
Ministry of Environmental Protection, China. Environmental Quality Standards for Surface Water (GB3838–2002). (2002) (in Chinese).
Asadollah, S. B. H. S., Sharafati, A., Motta, D. & Yaseen, Z. M. River water quality index prediction and uncertainty analysis: A comparative study of machine learning models. J. Environ. Chem. Eng. 9, 104599 (2021).
Uddin, M. G., Nash, S., Rahman, A. & Olbert, A. I. A comprehensive method for improvement of water quality index (WQI) models for coastal water quality assessment. Water Res. 219, 118532 (2022).
Kisi, O. & Ay, M. Comparison of Mann-Kendall and innovative trend method for water quality parameters of the Kizilirmak River, Turkey. J. Hydrol. 513, 362–375 (2014).
Güçlü, Y. S. Improved visualization for trend analysis by comparing with classical Mann-Kendall test and ITA. J. Hydrol. 584, 124674 (2020).
Castrillo, M. & García, Á. L. Estimation of high frequency nutrient concentrations from water quality surrogates using machine learning methods. Water Res. 172, 115490 (2020).
Liu, C. et al. Spatial and temporal characteristics of water quality and response identification of land use patterns in middle and downstream of Hanjiang River. Res. Environ. Sci. 34, 910–919 (2021) ((in Chinese)).
Zhang, W. et al. Unraveling El Niño’s impact on the east Asian monsoon and Yangtze River summer flooding. Geophys. Res. Lett. 43, 11375–11382 (2016).
Xu, H. et al. Anthropogenic impact on surface water quality in Taihu Lake region, China. Pedosphere 19, 765–778 (2009).
Hou, W. et al. Assessing water quality of five typical reservoirs in lower reaches of Yellow River, China: Using a water quality index method. Ecol. Indic. 61, 309–316 (2016).
Tian, W. et al. Effect of a comprehensive improvement project on water quality in urban lakes: A case study of water quality variation in Lihu lake over the past 30 years. Environ. Sci. 41, 183–193 (2020) ((in Chinese)).
Lu, X., Wang, X., Ban, X. & Singh, V. P. Considering ecological flow in multi-objective operation of cascade reservoir systems under climate variability with different hydrological periods. J. Environ. Manage. 309, 114690 (2022).
Da Luz, N. & Kumpel, E. Evaluating the impact of sampling design on drinking water quality monitoring program outcomes. Water Res. 185, 116217 (2020).
Ewaid, S. H., Abed, S. A. & Kadhum, S. A. Predicting the Tigris River water quality within Baghdad, Iraq by using water quality index and regression analysis. Environ. Technol. Innov. 11, 390–398 (2018).
Noori, R., Berndtsson, R., Hosseinzadeh, M., Adamowski, J. F. & Abyaneh, M. R. A critical review on the application of the national sanitation foundation water quality index. Environ. Pollut. 244, 575–587 (2019).
Sun, W., Xia, C., Xu, M., Guo, J. & Sun, G. Application of modified water quality indices as indicators to assess the spatial and temporal trends of water quality in the Dongjiang River. Ecol. Indic. 66, 306–312 (2016).
Acknowledgements
This work is financially supported by the National Natural Science Foundation of China (No. 51909112) & the China Postdoctoral Science Foundation (2022M711493).
Author information
Authors and Affiliations
Contributions
Lele Deng: Conceptualization, Methodology, Software, Investigation, Visualization, Writing- original draft. Kebing Chen: Conceptualization & Writing-review. Zhangjun Liu: Data curation & funding acquisition. Boyang Wu: Data curation. Zekun Chen: Data curation. Shaokun He: Conceptualization, Supervision, Writing- review & editing, funding acquisition.
Corresponding authors
Ethics declarations
Competing interests
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Deng, L., Chen, K., Liu, Z. et al. Spatiotemporal variation evaluation of water quality in middle and lower Han River, China. Sci Rep 12, 14125 (2022). https://doi.org/10.1038/s41598-022-16808-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-16808-w
This article is cited by
-
Unraveling the response of water quality and microbial community to lake water backflowing in one typical estuary of Lake Taihu, China
Environmental Monitoring and Assessment (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.