New Patterns of Temporal and Spatial Variation in Water Quality of a Highly Artificialized Urban River-Course—A Case Study in the Tongzhou Section of the Beiyun River

This paper aims to gain a better understanding of urban river pollution through evaluation of water quality. Data for 10 parameters at eight sites of the Tongzhou Section of the Beiyun River (TSBR) are analyzed. Hierarchical cluster analysis, fuzzy comprehensive assessment, discriminant analysis and Spearman’s correlation analysis were used to estimate the water situation of each cluster and analyze its spatial-temporal variations. Principal component analysis/factor analysis were applied to extract and recognize the sources responsible for water-quality variations. The results showed that temporal variation is greater than spatial and sewage discharge is the dominant factor of the seasonal distribution. Moreover, during the rapid-flow period, water quality is polluted by a combination of organic matter, phosphorus, bio-chemical pollutants and nitrogen; during the gentle-flow period, water quality is influenced by domestic and industrial waste, the activities of algae, aquatic plants and phosphorus pollution. In regard to future improvement of water quality in TSBR, the control of reclaimed wastewater from adjacent factories should first be put in place, as well as other techniques, for example, an increase of the impervious area, low-impact development, and integrated management practices should also be proposed in managing storm water runoff.


Introduction
Services provided by river-course ecological systems, such as rivers and lakes, contribute to water habitat and human welfare directly and indirectly [1].Rivers play an important role in assimilating and carrying point source pollution (mainly from domestic and industrial wastewater) and non-point source pollution (runoff from agricultural land) [2][3][4][5].As the most vulnerable water bodies to pollution, it is imperative to recognize river water quality and interpret spatial and temporal variation patterns [6][7][8][9].For over a century, large scales of modern engineering technologies have been employed in river renovation.Until now, about 60% of rivers around the world have gone through artificialization such as damming, banking and channelization to meet human needs [10].This may result in a different

Study Areas
Beiyun River is the main branch of the Hai River.Upstream of the Beiyun River (viz.Wenyu River) has a length of 47.5 km and a catchment area of 2403 km 2 .Downstream of the Beiyun River has a length of 41.9 km and catchment area of 1945 km 2 .The climate belongs to continental monsoon climate zone, with a mean annual precipitation of 620 mm, and average temperature of 11.3 • C. May-September is the rainy season, accounting for 85% of total precipitation of the year.Average evaporation of the surface water is between 400-500 mm.The study area has undergone a sharp reduction of precipitation in recent years.evaporation of the surface water is between 400-500 mm.The study area has undergone a sharp reduction of precipitation in recent years.

Sampling and Parameters
Water samples were collected from eight sites along the main channel of the Beiyun River in Spring (March, April, May), Summer (June, July, August), and Autumn (October) of 2014 and 2015 (Figure 1).Duplicate water samples were collected with a plastic container at each selected sampling site.The samples were placed in 1000 mL plastic bottles that had been soaked in dilute nitric acid for 24 h, and rinsed thoroughly with deionized water.These samples were taken to the laboratory and stored in a refrigerator at −4 °C for further analysis.Ten parameters were measured including DO, Chl-a, SS, TN, NH4-N, NO3-N, COD, BOD5, TP, PO4-P.Among them, DO, Chl-a, SS were measured on-site by Hydrolab Datasonde5.TN was determined using Multi N/C UV, and PO4-P was determined by using ammonium molybdate tetrahydrate spectrophotometry.Other determinations measurements were based on the "Environmental Quality Standards for Surface Water" of China (GB 3838-2002).

Fuzzy Comprehensive Assessment (FCA)
For water quality assessment, single index/parameter may not be enough in illustrating water conditions, and FCA emerged to overcome the disadvantages of the single index method.Fuzzy evaluation methods comprehensively evaluate the contributions of various pollutants according to predetermined weights, and decrease the fuzziness by using membership functions.Therefore, the sensitivity is quite high compared to other index evaluation techniques [25].FCA has proved effective in solving problems of fuzzy boundaries and in controlling the effect of monitoring errors on assessment results [26][27][28].

Multivariate Statistical Methods
Hierarchical cluster analysis (HCA), discriminant analysis (DA) and principal component analysis/factor analysis (PCA/FA) have been applied on the standardized data set of surface water quality in TSBR.All the calculations are implemented by Office Excel 2010 (Microsoft Corporation Ten parameters were measured including DO, Chl-a, SS, TN, NH 4 -N, NO 3 -N, COD, BOD 5 , TP, PO 4 -P.Among them, DO, Chl-a, SS were measured on-site by Hydrolab Datasonde5.TN was determined using Multi N/C UV, and PO 4 -P was determined by using ammonium molybdate tetrahydrate spectrophotometry.Other determinations measurements were based on the "Environmental Quality Standards for Surface Water" of China (GB 3838-2002).

Fuzzy Comprehensive Assessment (FCA)
For water quality assessment, single index/parameter may not be enough in illustrating water conditions, and FCA emerged to overcome the disadvantages of the single index method.Fuzzy evaluation methods comprehensively evaluate the contributions of various pollutants according to predetermined weights, and decrease the fuzziness by using membership functions.Therefore, the sensitivity is quite high compared to other index evaluation techniques [25].FCA has proved effective in solving problems of fuzzy boundaries and in controlling the effect of monitoring errors on assessment results [26][27][28].

Multivariate Statistical Methods
Hierarchical cluster analysis (HCA), discriminant analysis (DA) and principal component analysis/factor analysis (PCA/FA) have been applied on the standardized data set of surface water quality in TSBR.All the calculations are implemented by Office Excel 2010 (Microsoft Corporation of Washington, USA) and SPSS 21.0 (Norman H. Nie, C. Hadlai (Tex) Hull, and Dale H. Bent, the Headquartered is in Chicago, USA.).
HCA, the most widely used, is an exploration method developing from CA.Its primary function is to assemble objects based on characteristics' closeness.Based on their similarity, objects are clustered step by step until only one cluster remains [29].The similarity can be expressed in two ways: distance between samples or classes [30].In this paper, HCA is used to interpret the data and indicate patterns [3].Ward's method with squared Euclidean distances is used for similarity measurement.The results are typically illustrated by a dendrogram (tree diagram), which provide a visual summary of the clustering processes [31,32].DA is used to confirm the cluster found by HCA, while at the same time determining the most significant parameters responsible for spatial and temporal variations.Standard and stepwise modes of DA have been used in this paper.The clusters which HCA show is the grouping variables, and the measured parameters are the independent variables; the accuracy of DA results is tested by cross-validation.
PCA is a pattern recognition technique that attempts to explain the variance of a large set of inter-correlated variables by transforming them into a smaller set of independent variables [2,6].PCA is executed to extract significant principal components (PCs), then FA further reduces the contribution of less significant variables obtained from PCA and the new group of variables known as variance factors (VFs) is extracted through rotating the axis defined by PCA.A VF can include unobservable, hypothetical, and latent variables, while a PC is a linear combination of observable water-quality variables [3,24].When applies to water quality analysis, it is mainly used to extract pollution factors and recognize pollution sources [33].

The Overall Status of Water-Quality Characteristics
Summary characteristic of each parameter and guideline of national standards for surface water are shown in Table 1.By comparison between each parameter's mean value and associated national standard, it can be seen that the whole water quality in TSBR is poor.Seven out of 10 parameters are beyond the national standard V. N and P pollutions are the most severe.The mean value of TN reaches 22.25 mg•L −1 , five times higher than standard V. TP is four times higher with an average of 1.85 mg•L −1 .In view of variable coefficients, Chl-a (139.15%) and DO (103.33%)rank the top, followed by NO 3 -N (65.22%),SS (56.30%),BOD 5 (47.65%), which means a higher degree of dispersion and an uneven distribution in different periods and sites among these parameters.Similar situations were also found in other studies, which showed a poor water quality in urban areas.For example, Surindra [34] suggested that ranges of key water quality parameters such as, DO, BOD, COD and TDS, were several times higher than the prescribed standards, and that urban runoff and industrial wastewater discharges were the major threats to river water quality at Ghaziabad.Overall, the water quality of Hindon was relatively poor with respect to its use for domestic purposes.Alam, Hossain and Khan's [35] study on deterioration of water quality of the Surma River indicated that the water of the Surma River was polluted from various sources like sewage and garbage disposal.The river was found to be highly turbid in the monsoon season.Concentration of BOD and fecal coliform was found to be higher in the dry season.The water was unfit for drinking purposes without any form of treatment, but for various purposes of other surface water usage, it still could be considered quite acceptable.

Temporal and Spatial Distribution Characteristics of Water Quality
HCA assembles the objects into different clusters based on their similarity, and is used for recognizing the temporal and spatial distribution recognition.Then, standard and stepwise modes of DA are used to confirm the cluster found by HCA and to determine the most significant parameters determining spatial and temporal variations.In combination with the DA and Spearman's correlation analysis, the temporal and spatial distribution characteristics of water environmental element can be obtained.

Temporal Variations
Both of seasons and points are non-numeric parameters, so an assignment is necessary before the statistical analysis.The dendrogram of temporal HCA is shown in Figure 2a.It is observed that all the monitoring months can be obviously classified into two statistically meaningful clusters at (D link /D max )*100 < 15 (the quotient between the linkage distance for a particular case divided by the maximal distance).Cluster 1 contains May-October, and Cluster 2 contains March-April.As indicated before, 85% of the precipitation is concentrated in May-September; therefore, Cluster 1 and 2 basically correspond to the rapid-flow period and gentle-flow period, respectively.

Temporal and Spatial Distribution Characteristics of Water Quality
HCA assembles the objects into different clusters based on their similarity, and is used for recognizing the temporal and spatial distribution recognition.Then, standard and stepwise modes of DA are used to confirm the cluster found by HCA and to determine the most significant parameters determining spatial and temporal variations.In combination with the DA and Spearman's correlation analysis, the temporal and spatial distribution characteristics of water environmental element can be obtained.

Temporal Variations
Both of seasons and points are non-numeric parameters, so an assignment is necessary before the statistical analysis.The dendrogram of temporal HCA is shown in Figure 2a.It is observed that all the monitoring months can be obviously classified into two statistically meaningful clusters at (Dlink/Dmax)*100 < 15 (the quotient between the linkage distance for a particular case divided by the maximal distance).Cluster 1 contains May-October, and Cluster 2 contains March-April.As indicated before, 85% of the precipitation is concentrated in May-September; therefore, Cluster 1 and 2 basically correspond to the rapid-flow period and gentle-flow period, respectively.Integrated evaluation of water quality can be derived through FCA, which demonstrates in Table 2 that water quality are both beyond level V in the selected two periods.The results are largely affected by TN and TP's heavy weighting, which means severe N and P pollution in TSBR.By comparison between the pollution levels between two periods, it is observed that the pollution level in the rapid-flow period (0.976) is higher than gentle-flow period (0.942), which reveals the scouring movement by rainfall in the rapid-flow period.In addition, runoff carries large amounts of land and atmosphere pollutants into the river-course, so non-point source pollution becomes more serious [36].On the other hand, due to the runoff erosion, endogenous load within the river course (which is mainly from point source pollutants) can be released, which in turn increases pollutant concentration in rapid-flow period [37][38][39].Integrated evaluation of water quality can be derived through FCA, which demonstrates in Table 2 that water quality are both beyond level V in the selected two periods.The results are largely affected by TN and TP's heavy weighting, which means severe N and P pollution in TSBR.By comparison between the pollution levels between two periods, it is observed that the pollution level in the rapid-flow period (0.976) is higher than gentle-flow period (0.942), which reveals the scouring movement by rainfall in the rapid-flow period.In addition, runoff carries large amounts of land and atmosphere pollutants into the river-course, so non-point source pollution becomes more serious [36].On the other hand, due to the Water 2018, 10, 1446 6 of 14 runoff erosion, endogenous load within the river course (which is mainly from point source pollutants) can be released, which in turn increases pollutant concentration in rapid-flow period [37][38][39].Tables 3 and 4 show classification functions and verification obtained from the standard, stepwise modes of DA.For temporal DA, standard mode assigning 97.8% of cases means the selected 10 parameters can explain temporal variation of water quality successfully.While for stepwise mode, correct assignations reach 93.3% using only three parameters, which means Chl-a, TN, PO 4 -P should be enough to illustrate temporal variation.This suggests that the monitoring intensity of these three parameters should be appropriately increased in future works.
The results of SCA are shown in Table 5. NH 4 -N, TN, TP, PO 4 -P, BOD 5 , Chl-a are found significantly (sig.< 0.01) correlated with the periods, suggesting that these parameters can be easily affected by temporal factors like water stage or flow rate.Through the results from SCA and DA, three parameters selected by SCA cannot be found through DA.NH 4 -N, TP, BOD 5 which correlates with the periods.The temporal variation is explained by.NH 4 -N, TP, whose loss of discriminant ability is mainly caused by flexible water discharge from factories along the main river course in Tongzhou section.In addition, the degradation ability of microbial factors varies along with the scale and frequency of the water discharge, which explains BOD 5 fails to show seasonal variations.The results also suggest that NH 4 -N, TP, BOD 5 are highly vulnerable to human activities.
As identified by DA, box and whisker plots of the selected parameters showing temporal trends are demonstrated in Figure 3. Significant change occurs in all the 3 parameters.NH 4 -N varies from 3.14 to 19.20 mg/L in the rapid-flow period and 1.87-16.30mg/L in the gentle-flow period.Average value in the gentle-flow period (11.31 mg/L) is higher than that in the rapid-flow period (7.56 mg/L), which can be explained as in the gentle-flow period.The microbial degradation rate of nitrogen is limited by the lower temperature.TP varies from 0.71 to 2.87 mg/L in the rapid-flow period and 0.57-1.63mg/L in the gentle-flow period.Average value in the rapid-flow period (1.69 mg/L) is higher than that in the gentle-flow period (1.15 mg/L).Different kinds of non-point source pollution (including phosphorous pollutant) flow into the river course along with rainfall runoff in the rapid-flow period, and a higher temperature in the rapid-flow period may enhance the sediment phosphorus release [40].For BOD 5 , the average value in the gentle-flow period (12.14 mg/L) is higher than that in the rapid-flow period (7.60 mg/L).Studies have shown a significant correlation between BOD 5 and Chl-a.While BOD 5 peaks, algal blooms at the greatest level at the same time.Scarce water volume, insufficient water power, and higher nutrient concentrations can all enhance the respiration of phytoplankton, finally resulting in higher BOD 5 in the gentle-flow period.Affected by natural conditions and human activities, the variation of water quality in urban river courses may have a unique pattern of temporal variation; this paper shows that different study objects may lead to different seasonal variation ranges.Pejman [11] indicated that the natural parameters (temperature and discharge), the inorganic parameter (total solid) and the organic nutrients (nitrate) were the most significant parameters contributing to water-quality variations for all seasons.Prakash Raj Kannel's partition was based on precipitation duration (pre-monsoon season and post-monsoon seasons), high level of BOD5, COD and low DO indicated that the level of organics was higher in the pre-monsoon season, while during post-monsoon season, the main pollutants were nutrients [41].

Spatial Variations
The dendrogram of spatial HCA is shown in Figure 2b; eight sampling stations can be grouped into three statistically meaningful clusters at (Dlink/Dmax)*100 < 18.From the site location and distribution of supplement sites in Figure 1, these three clusters correspond to up-, middle and down stream of the reach which is filled with reclaimed water.Three sites from Cluster 1 (site 1, site 6 and site 7) comprise the closer sites to drainage channels.Cluster 2 (site 2 and site 3) corresponds to moderate distance sites.Cluster 3 (site 4, site 5 and site 8) makes up the farthest ones.Affected by natural conditions and human activities, the variation of water quality in urban river courses may have a unique pattern of temporal variation; this paper shows that different study objects may lead to different seasonal variation ranges.Pejman [11] indicated that the natural parameters (temperature and discharge), the inorganic parameter (total solid) and the organic nutrients (nitrate) were the most significant parameters contributing to water-quality variations for all seasons.Prakash Raj Kannel's partition was based on precipitation duration (pre-monsoon season and post-monsoon seasons), high level of BOD 5 , COD and low DO indicated that the level of organics was higher in the pre-monsoon season, while during post-monsoon season, the main pollutants were nutrients [41].

Spatial Variations
The dendrogram of spatial HCA is shown in Figure 2b; eight sampling stations can be grouped into three statistically meaningful clusters at (D link /D max )*100 < 18.From the site location and distribution of supplement sites in Figure 1, these three clusters correspond to up-, middle and down stream of the reach which is filled with reclaimed water.Three sites from Cluster 1 (site 1, site 6 and site 7) comprise the closer sites to drainage channels.Cluster 2 (site 2 and site 3) corresponds to moderate distance sites.Cluster 3 (site 4, site 5 and site 8) makes up the farthest ones.
Water 2018, 10, 1446 9 of 14 Spatial FCA is shown in Table 2, upstream of the reach has the largest membership degree to level V (0.982), followed by middle-stream (0.969) and downstream (0.960).The results show that water pollution has a slightly decrease trend when it travels from the upper stream to downstream.
As shown in Table 6, standard mode DA yields 64.4% accuracy, which means spatial variation exists in TSBR, but is not as obvious as the temporal variation.Only one parameter (NO 3 -N) is selected by stepwise DA with the accuracy of less than 52.0%, which indicates the fact that only one parameter cannot be enough in describing spatial variation of water quality.Therefore, the "F" value is changed to improve the determination accuracy.With "F = 2" in the arithmetic [42], three parameters are identified, and the determination accuracy rises to 64.2% (69.6%, 64.3% and 56.3% for up, middle and downstream, respectively).NO 3 -N (Sig.< 0.01) and DO (Sig.< 0.05) are correlated to spatial variations through SCA.Relatively fewer parameters in spatial SCA indicate that there is lack of dissimilarity in spatial scale; all the three river-courses may have a similar pollution condition.According to the results from SCA and DA, TN and Chl-a can be determined by DA but not SCA, which means the variation of these two parameters only exist in specific parts of the river course.DO is significant in SCA but not in SCA, and loss of discriminant ability may contribute to the unset distribution of pollutants, which are caused by point source pollution.
Box and whisker plots of the selected parameters of spatial trends are illustrated in Figure 4. From the point of variation amplitude, NO 3 -N varies between 1.33 and 15.18 mg/L at upstream, 4.62-12.30mg/L at middle-stream, and 1.20-9.65 mg/L at downstream, the great variation means an uneven distribution in different parts of the river course.TN has the greatest variation at up-stream.For Chl-a, the greatest variation appears at down-stream.From comparison between average values in the three parts of the river course, NO 3 -N does not show a dramatic variation, as well as TN and Chl-a, which further explains the low accuracy in spatial DA.As the main sewage drainage river course, TSBR is largely affected by domestic wastewater.Automatically, the high NO 3 -N, which is the degradative product of nitrogen organic matter by microbial action, is mainly from domestic wastewater [43].Besides, TN, Chl-a are infected by the fact that factories along the river course do not have obvious geographical characteristics or discharge frequency.

Source Identification
PCA/FA is applied in recognizing main factors of water pollution.As pointed out in Section 3.2.2,temporal variation is greater than spatial variation, thus the paper focuses on recognizing pollution factors in different periods, as delineated by HCA.Eigenvalues greater than one are taken as the criterion for extraction of the factors required to explain the sources of variances in the data.Each extraction factor contains comprehensive information of the original variables, so it may be hard to explain.Therefore, a rotation must be made upon the extraction factors in order to make clear the meaning of each factor.In this paper, varimax rotation is adopted to minimize the variables with the highest load in each factor.The new group of variables is known as variance factors (VFs).Classification of factor loading is thus "strong", "moderate" and "weak", corresponding to absolute loading values of >0.75, 0.75-0.50,and 0.50-0.30,respectively [44].Corresponding VFs, variables loading and variance in two periods are presented in Table 7; strong and moderate loading values have been highlighted by underline and bold.For the dataset from the rapid-flow period, three VFs accounts for 71.82% of total variance.VF1 explaining 31.30% of total variance, which has strong positive loadings on TP, PO4-P, SS, TN and COD.Phosphorus indicators imply that part of this pollution is from sewage discharge; on the other hand, organic pollutants have been intensely decreasing in supplementary (reclaimed water from treatment plants); a positive correlation of organic pollutant (COD and SS) means that these are mainly from non-point source pollution like agricultural runoff and atmospheric sedimentation but not from sewage discharges.Thus, VF1 represents a combination of point source and non-point

Source Identification
PCA/FA is applied in recognizing main factors of water pollution.As pointed out in Section 3.2.2,temporal variation is greater than spatial variation, thus the paper focuses on recognizing pollution factors in different periods, as delineated by HCA.Eigenvalues greater than one are taken as the criterion for extraction of the factors required to explain the sources of variances in the data.Each extraction factor contains comprehensive information of the original variables, so it may be hard to explain.Therefore, a rotation must be made upon the extraction factors in order to make clear the meaning of each factor.In this paper, varimax rotation is adopted to minimize the variables with the highest load in each factor.The new group of variables is known as variance factors (VFs).Classification of factor loading is thus "strong", "moderate" and "weak", corresponding to absolute loading values of >0.75, 0.75-0.50,and 0.50-0.30,respectively [44].Corresponding VFs, variables loading and variance in two periods are presented in Table 7; strong and moderate loading values have been highlighted by underline and bold.For the dataset from the rapid-flow period, three VFs accounts for 71.82% of total variance.VF1 explaining 31.30% of total variance, which has strong positive loadings on TP, PO 4 -P, SS, TN and COD.Phosphorus indicators imply that part of this pollution is from sewage discharge; on the other hand, organic pollutants have been intensely decreasing in supplementary (reclaimed water from treatment plants); a positive correlation of organic pollutant (COD and SS) means that these are mainly from non-point source pollution like agricultural runoff and atmospheric sedimentation but not from sewage discharges.Thus, VF1 represents a combination of point source and non-point source pollution [38].VF2, explaining 22.01% of total variance, has positive loading on BOD 5 , Chl-a and DO, which can be interpreted as bio-chemical pollution and eutrophication caused by sewage discharge [6,10].VF3, explaining 18.51% of total variance, has a positive loading on NH 4 -N and negative loading on NO 3 -N.Nitrification in the water body is the main reason for NO 3 -N decrease and NH 4 -N increase [45].Thus, VF3 represent nitrogen pollution caused by rainfall and surface runoff.
Three VFs accounts for 75.67% of total variance in the gentle-flow period.Among the three VFs, VF1 explains 29.59% of total variance, has strong positive loadings on BOD 5 , COD, NH 4 -N, and negative loading on NO 3 -N.As for local field situations, VF1 basically represent domestic and industrial sewage.VF2, explaining 24.58% of total variance, has positive loadings on Chl-a, DO and negative loading on NO 3 -N.Previous studies show that microbial activities can increase denitrification and consume the inorganic nitrogen in water [46][47][48].Therefore, VF3 can be interpreted as the effects of algae and other aquatic plants.VF3 has positive loadings on TP and PO 4 -P, both of which can be regarded as phosphorus pollution caused by domestic discharge.
In summary, for an urban river course, the constitution of pollutants differs from natural water bodies, but both of them are obviously disturbed by human activities.Point-source pollution, like domestic sewage and industrial wastewater, is the main pollution source of both urban and natural river water bodies, which is consistent with previous research [34,[49][50][51].

Conclusions and Discussions
As a significant part of an urban ecosystem, an urban river course plays an important role in carrying municipal industrial wastewater and agricultural runoff [52].Research on urban river water-quality status as well as pollution distribution may serve as important guidance for urban water environment improvement.
This paper reveals that water environment has been seriously degrading in TSBR influenced by both internal and external source pollution.Most of the monitoring parameters are far beyond the national standard V. N and P are the most severe ones, while Chl-a, DO make up the uneven distribution ones.This indicates that water eutrophication is serious and the self-cleaning ability of the water is poor in the study area.With the gradual decrease of rainfall and the increase of wastewater discharge in recent years, and the rapid development of urbanization, this has led to an increase in sewage discharge of the Beiyun River.In addition, due to the drought in successive years, the amount of clear water resources is insufficient and the proportion of sewage is increasing, resulting in a decrease in the dilution and natural purification of rivers.The external reason is that the Tongzhou section of the Beiyun River is a gate-controlled river, and the river channel is silted.The water flow is relatively small, which leads to the reduction of the water's self-purification capacity.The deterioration of water quality in the urban river course of the Beiyun River is greatly affected by the inflow of sewage into the river.The interception of sewage is a prerequisite for the improvement of the water quality of the Tongzhou section of the Beiyun River in the future.
In this paper, HCA is applied to the river water quality, and the whole data set can be divided into two parts: rapid-flow period (May-October), and gentle-flow period (March-April).From the dendrogram in Figure 2a, it can be seen that May is included in dry season along with March and April, but the month also shows similar characteristics with wet seasons like June and August.This means precipitation is not the key factor resulting in the change of pollution condition in the urban river course; sewage discharge is the dominant factor of the seasonal distribution of water quality.Therefore, it can be concluded that highly artificialization characteristics have resulted in the changes of nature feature of the urban river course.The urban river course suffers from point and non-point source pollution at the same time, presenting a far more serious water condition.For different periods, water pollution in the rapid-flow period is more severe than the gentle-flow period as interpreted by FCA, which means the infusion of runoffs during the rapid-flow period does not dilute the pollutant concentration within water body.On the contrary, the infusion of runoff becomes the secondary pollution source.As stepwise mode assigns 93.3% of cases correctly, three parameters (Chl-a, TN, PO 4 -P) demonstrated their good performance in explaining temporal variation of water quality successfully; in monitoring the frequency of these three parameters should be appropriately increased in future works.The results of SCA and DA comprehensively suggest that NH 4 -N, TP, BOD 5 are highly vulnerable to human activities.
TSBR is the typical river course protection built by laying bricks or stones which are carried out on both sides of the river bank and its bottom; this highly artificialized action leads to a great variation for hydrochemical and hydrological characteristics compared to a natural river course.Because of the limited catchment area compared to natural rivers, the paper fails to find a dramatic spatial variation of water quality in different reaches.Thus, the paper only focuses on pollution identification in different periods.The results show that influences by multivariate factors like precipitation, temperature, sources/constituents are different between the two periods.During the rapid-flow period, water is polluted by a combination of organic matter (non-point source pollutant) and phosphorus (point source pollutant), bio-chemical pollutants and nitrogen (point source pollutant).For the gentle-flow period, water quality is influenced by domestic and industrial waste (point source pollutant), the activities of algae and aquatic plants, and phosphorus pollution (point source pollutant).Point-source pollution is blamed for water deterioration in both periods.Due to a lack of precipitation and runoff, agricultural runoff is no longer responsible for water deteriorating during the gentle-flow period.Thus, future works should focus on the restraint of domestic and industrial releases for the gentle-flow period, while during the rapid-flow period, various protection measures should be put into effect for reducing non-point source pollution.
Strategies on storm water management have been explored and put into practice in some areas around the world, but there is still increasing concern about how to better utilize storm water in the future, especially for countries that lack water resources.Furthermore, some techniques should be adopted, such as the increase of impervious area and the application of low-impact development or integrated management practices, in managing storm water runoff in the future.

Figure 1 .
Figure 1.Map of the study area and surface water sampling sites in the Tongzhou Section of the Beiyun River (TSBR).

Figure 1 .
Figure 1.Map of the study area and surface water sampling sites in the Tongzhou Section of the Beiyun River (TSBR).

Figure 2 .
Figure 2. Dendrogram of cluster analysis for spatial and temporal variation in TSBR.(a) Temporal cluster analysis.(b) Spatial cluster analysis.

Figure 2 .
Figure 2. Dendrogram of cluster analysis for spatial and temporal variation in TSBR.(a) Temporal cluster analysis.(b) Spatial cluster analysis.

Figure 3 .
Figure 3.The trend of temporal variations of significant water quality parameters in TSBR.

Figure 3 .
Figure 3.The trend of temporal variations of significant water quality parameters in TSBR.

Figure 4 .
Figure 4.The trend of spatial variations of significant water quality parameters in TSBR.

Figure 4 .
Figure 4.The trend of spatial variations of significant water quality parameters in TSBR.

Table 1 .
Statistics of measured parameters and the national standards for surface water quality.

Table 2 .
Integrated evaluation form of water quality.

Table 2 .
Integrated evaluation form of water quality.

Table 3 .
Classification functions for discriminant analysis (DA) of temporal and spatial variations.
3RF/GF stand for rapid-flow period/gentle-flow period, US/MS/DS stand for up-stream/middle-stream/down-stream above.

Table 4 .
Verification and classification functions for DA of temporal variations.

Table 5 .
Correlation coefficient of Spearman's correlation analysis for water quality.Sig.< 0.01 correspond to extremely significant of correlation, Sig.< 0.05 correspond to significant of correlation.

Table 4 .
Verification and classification functions for DA of temporal variations.

Table 5 .
Correlation coefficient of Spearman's correlation analysis for water quality.

Table 6 .
Verification and classification functions for DA of spatial variations.

Table 7 .
Factor loadings value and explained variance of water quality parameters in two periods.

Table 7 .
Factor loadings value and explained variance of water quality parameters in two periods.