Temporal water quality response in an urban river: a case study in peninsular Malaysia

Ambient water quality is a prerequisite for the health and self-purification capacity of riverine ecosystems. To understand the general water quality situation, the time series data of selected water quality parameters were analyzed in an urban river in Peninsular Malaysia. In this regard, the stations were selected from the main stem of the river as well as from the side channel. The stations located at the main stem of the river are less polluted than that in the side channel. Water Quality Index scores indicated that the side channel station is the most polluted, breaching the Class IV water quality criteria threshold during the monitoring period, followed by stations at the river mouth and the main channel. The effect of immediate anthropogenic waste input is also evident at the side channel station. The Organic Pollution Index of side channel station is (14.99) ~3 times higher than at stations at river mouth (4.11) and ~6 times higher than at the main channel (2.57). The two-way ANOVA showed significant difference among different stations. Further, the factor analysis on water quality parameters yielded two significant factors. They discriminated the stations into two groups. The land-use land cover classification of the study area shows that the region near the sampling sites is dominated by urban settlements (33.23 %) and this can contribute significantly to the deterioration of ambient river water quality. The present study estimated the water quality condition and response in the river and the study can be an immediate yardstick for base lining river water quality, and a basis for future water quality modeling studies in the region.


Introduction
Rivers are lifelines for human societies around the globe, embodying immense influence in shaping civilizations. The catchment area usually supports a wide variety of flora and fauna, creating a very diverse ecosystem composed of ecologically delicate and inter related, physical, chemical and biological entities. However, rivers are also the subsequent waste disposal arena for anthropogenic activities and are pathways of waste materials to coastal regions. Although the global water crisis tends to be viewed as a water quantity problem, water quality is increasingly being acknowledged as a central factor in the water crisis (Belayneh and Bhallamudi 2012). The quality of river water is a deterministic factor for the healthy, sustainable survival of the riverine ecosystem, which primarily depends on the Waste Assimilative Capacity (WAC) of the water. WAC is the natural ability of the river to withstand or assimilate a certain amount of pollutants without impairing ambient water quality Electronic supplementary material The online version of this article (doi:10.1007/s13201-015-0303-1) contains supplementary material, which is available to authorized users. conditions (Krom 1986;Tett et al. 2011). Increasing anthropogenic contaminants affect WAC and water quality is deteriorated beyond WAC of a water body (VishnuRadhan et al. 2014(VishnuRadhan et al. , 2015. This can disrupt the prevailing ecological homeostasis, ultimately affecting the riverine health. In water quality management, the determination of each water quality variable is important to obtain collective information on water quality, as it can provide concise information on overall environmental conditions (Chen et al. 2007). Practically, it is difficult to assess water quality based on each parameter/variable. The water quality indices of significant and influential parameters aim at giving a single value to the water quality of a source on the basis of a system which translates their existing concentrations in a sample into a single value. These values are used as communication tools by regulatory agencies to describe the quality or health of a specific environmental system (Abbasi and Abbasi 2012). The WQI is frequently utilized as a mathematical tool for evaluating water quality status around the globe (Shirodkar et al. 2010;Lumb et al. 2011;Gazzaz et al. 2012;Dede et al. 2013).
In urban areas, streams are often degraded as they are diverted through storm water runoff systems, removal of riparian vegetation, and the construction of roads, parking lots and buildings (Buffers 2000). Riparian zones have diversified functions that include preserving bank stability, functioning as habitats for streamside living organisms and also playing a critical role in preserving the water quality of rivers by filtering out pollutants from runoff (Zainudin et al. 2013). Water quality is a major factor impacted by anthropogenic action at landscape scales which is a principal threat to the ecological integrity of river ecosystems (Allan 2004). Land-use changes often affect the water quality over a long historical period (Garnier et al. 2013) and future land-use changes will exacerbate the water quality problems (Whitehead et al. 2013). The changes in ecosystem goods and services that result from landuse change revert on the drivers of land-use change (Lambin et al. 2003). Many studies have quantified the effect of population increase on land use/land cover (LULC) (Meyers and Turner 1992;Wu et al. 2013;Meyfroidt et al. 2013) and the associated anthropogenic activities can ultimately reflect on the water quality of natural waters. Thus, LULC can give a generalized impression on the state of a river's water quality in an urban area.
To understand the water quality situation and contribution of anthropogenic activities on the water quality degradation in an urban river, an investigation was performed using time series data on water quality parameters. In this regard, a suite of selected referred water quality indices, GIS and statistical techniques were attempted in the present study. The present study is first of its kind in Sg. Sri Melaka and will contribute to the baseline information for future water quality studies in the region.

Study area
Sg. Sri Melaka/Sg. Malim is a small urban river, located in the state of Melaka in southern Peninsular Malaysia. The climate of the region is generally characterized as humid tropical and peak rainfall commences in September and ends abruptly around November. In March, the rainfall amount rises again and a smaller peak occurs in April. The rainfall amount remains virtually constant until September (Asry et al. 2012). The river, which drains into the Straits of Malacca (Fig. 1a), is frequently used for recreational fishing. Near the river mouth, the waste water discharge is more domestic than industrial. Floodgates located on the channel trap pollutants from being flushed out to sea, as they are closed most of the time. This results in a grotesque condition on this particular stretch. Debris and floatables are observable, indicative of its polluted physico-chemical state. Moreover, there is an urban settlement that may be exposed to foul odor emanating from the channel due to anoxic biodegradation. Land reclamation works are also currently going on downstream (mouth region), paving the way for future development. The locations of the sampling stations are shown in Fig. 1b, whereas Table 1 describes each sampling station and its relevance to the study. Three stations were selected for the present study: one at the side channel, the second at the main stem of the river and the third one at the river mouth.

Water quality data
Time series data of water quality parameters from the three sampling stations were collected from the basin. In this regard, samples were collected at 3-h intervals for 3 days from 23 rd to 26 th of August 2010. The water quality parameters measured are pH, temperature, Dissolved Oxygen (DO), Biochemical Oxygen Demand (BOD), Chemical Oxygen Demand (COD), Total Suspended Solids (TSS), oil & grease, Phosphorous (P), Ammoniacal nitrogen (NH 3 -N), and E. coli. The methodologies adopted for the sample analysis are presented in Table 2.

Water Quality Index
The primary method employed to classify the river water quality was the Water Quality Index (WQI) and the National Water Quality Standards (NWQS), a set of standards derived based on beneficial uses of water in Malaysia. The NWQS defined classes I-V, referred to classification of rivers or river segments based on the descending order of water quality: Class I being the best and Class V being the worst (Zainudin 2010). A WQI ascribes quality value to an aggregate set of measured parameters. It usually consists of sub-index values assigned to each pre-identified parameter by comparing its measurement with a parameter-specific rating curve, optionally weighted, and combined into the final index. The purpose of a WQI is to summarize large amounts of water quality data for a specific river into simple values (i.e., one number and a statement such as ''good'') (Saffran et al. 2001).
The WQI primarily used in Malaysia, also referred to as the Malaysian Department of Environment-Water Quality Index (DOE-WQI), is an opinion-poll formula where a panel of experts is consulted on the choice of parameters and on the weightage to each parameter (DOE 1985). The WQI is calculated using six parameters WQI: DO, BOD, COD, TSS, NH 3 -N and pH with the inclusion of intermediate sub-indices. Calculations are performed on the water quality parameters to find out their respective sub-indices. The sub-indices are named SIDO, SIBOD, SICOD, SIAN, SISS and SIPH. The best fit equations used for the estimation of the six sub-indices are shown below (DOE 2007).
(a) Sub-index for DO (in % saturation): SIDO   where x is the concentration in mg/L for all parameters except pH.
Once the respective sub-indices have been calculated, the WQI can then be calculated using Eq. (1).
The summation of the weightages for all the sub-indices must have a value of unity. The respective class designation for the WQI scores is shown in Tables 3 and 4.

Organic pollution index
Organic Pollution Index (Wei et al. 2009) is an immediate and reliable measure of river water quality and pollution. The equation is as follows: where A is the OPI, while BOD i , COD i , NH 3 -N i and DO i are the monitored pollution concentrations in different segments and BOD 0 , COD 0 , NH 3 -N 0 and DO 0 are the guidelines set for the maximal amount of permitted pollution content. If A C 2, the river water begins to be contaminated with organic matter. We have used the following reference values for the calculation of organic pollution index (Class II).

Statistical analysis
The descriptive statistics and two-way ANOVA on the time series data of water quality parameters were performed to identify significant differences between stations and time. PROC MEANS procedure of Statistical Analytical Systems (SAS) 9.3 (SAS 2012) was used to estimate the descriptive statistics, viz. minimum value, maximum value, mean, standard error and coefficient of variation for various water quality parameters. The significant source of variation (station, time) was detected by analysis of variance adopting the two-way ANOVA using the PROC GLM procedure of SAS 9.3 (SAS 2012). The ANOVA was followed by Tukey's HSD test for analyzing the grouping among the factors (station or time) using the 'MEANS' statement in PROC GLM procedure of SAS 9.3 (SAS 2012). Further, the ten water quality parameters were subjected to factor analysis using PROC FACTOR procedure of SAS 9.3 (SAS 2012) to test whether the water quality parameters are effective in discriminating different stations. Factor analysis is a statistical method used to describe variability among observed, correlated variables in terms of a potentially lower number of unobserved variables called factors. For example, it is possible that variations in four observed variables mainly reflect the variations in two unobserved variables. Factor analysis searches for such joint variations in response to unobserved latent variables. The observed variables are modeled as linear combinations of the potential factors, plus ''error'' terms. The information gained about the interdependencies between observed variables can be used later to reduce the set of variables in a dataset. Computationally, this technique is equivalent to low rank approximation of the matrix of observed variables. In the factor analysis, the water quality parameters which loaded heavily on the first and second factors were identified. The variables loaded on different factors were selected based on Hatcher's scratching procedure (Hatcher 2003). The statistical procedures were carried out using software platforms of STATISTICA (Hill and Lewicki 2007) and SAS 9.3 (SAS 2012).

Land-use land cover classification
An area of 321 sq. km surrounding the sampling site is analyzed for the land-use land cover classification using ArcGIS 10.2 and ERDAS Imagine 2013 to explore the general contribution by anthropogenic activities towards water quality of the study area. Landsat 7 ETM ? satellite imagery having 30 meter resolution from USGS earth explorer is downloaded for the year 2010. This image is used to perform geo referencing and a final RGB band true color and false color imagery is prepared using layer stack tool in ERDAS Imagine 2013. Imagery resolution is increased to 15 m by pan merging technique to increase the accuracy of the results. The final image obtained after preprocessing the satellite imagery is used for classification. Supervised classification is carried out using ERDAS Imagine 2013 for land-use land cover assessment of the study area. The image is classified into four major classes, i.e., vegetation, water bodies, barren land and urban settlements. Area of classes is calculated using histogram value of different bands in true color imagery.

Results and discussion
The site SG1 is located at the side channel (adjacent to the residential area) of the river; SG2 is on the main channel and SG3 on the river mouth. Near stagnant water flow condition was observed at SG1 in comparison with other stations during the monitoring period (SI 1). The overall descriptive statistics (Table 5) for the water quality parameters showed that there is high variation in NH 3 -N (coefficient of variance (CV)-177.11 %), P (CV-60.67 %), E. coli (CV-49.27 %) and TSS (49.12 %). NH 3 -N ranged from 0.005 to 9.86 mg/L with a mean value of 1.52 ± 0.32 mg/L. DO ranged from 0.62 to 6.20 mg/L with mean of 2.84 ± 0.15 mg/L and temperature ranged from 25.80 to 31.80 with a mean value of 30.03 ± 0.10°C. The acidic pH, existing in the system, ranged from 4.13 to 6.23 with a mean value of 5.05 ± 0.06. COD values ranged from 4 to 71 mg/L with a mean value of 32.08 ± 1.99 mg/ L, while BOD values ranged from 2 to 8 mg/L with a mean  Appl Water Sci value of 3.64 ± 0.18 mg/L. P ranged from 0.10 to 1.14 mg/L with a mean value of 0.36 ± 0.03 mg/L, while the TSS ranged from 18 to 160 mg/L with a mean value of 61.17 ± 3.54 mg/L. Temperature values at all stations have shown a decreasing trend and acidic pH prevailed throughout the observation period. At SG2 and SG3, DO followed the decreasing trend, similar to temperature, and at SG1 low DO is observed as the flow is stagnant most of the time. High BOD and COD were observed at SG1 and SG2 in comparison with the station at the river mouth, SG3. The lower COD and BOD values at SG3 may be probably due to the fast pollutant flushing towards the sea. Similarly, Low TSS values observed at SG1 may be due to the stagnant nature of the side channel in comparison with the main channel. Almost similar trends of high phosphorous at all the stations indicated the anthropogenic addition from an adjacent residential area. High NH 3 -N and E. coli were observed at SG1 compared to other stations.
Interaction plots between time and stations were generated for different water quality parameters. The time wise distribution of water quality parameters was not significantly different except for temperature. On the other hand, station wise difference was observed for NH 3 -N, COD, DO, E. coli, Phosphorous and TSS. In addition, the temperature and pH did not show significant difference between the stations. The results are clearly presented in the interaction plots (Fig. 2). Two-way ANOVA (Table 6) showed that there is significant difference between the stations for NH 3 -N, COD, DO, E. coli, P and TSS. However, there was no significant difference between the stations for BOD, temperature and pH. Time wise significant difference was observed only in case of temperature. The Tukey's HSD test results and classifications for the different water quality parameters based on stations are shown in Table 7. SG1 is found to have higher values in comparison with SG2 and SG3 for NH 3 -N, COD, E. coli and P, while SG1 has lower values in the case of DO and TSS and it gives an impression that the pollution load would be more in SG1 in comparison with SG2 and SG3. However, there was no significant difference between stations for BOD, temperature and pH.
The factor analysis (Table 8) revealed that the variables with high loadings on the first two factors were found useful in distinguishing the stations. The first two factors (Factor 1-61.05 % and Factor 2-36.25 %) together explained about 97.3 % of the variation in the data. The plot (Fig. 3) between first and second factor scores indicated that these factors were efficient in discriminating the stations SG1, SG2 and SG3. On the basis of the first and  second factors, SG1 is found to be separated from SG2 and SG3. The variables loaded on the first factor were E. coli, NH 3 -N and COD with positive loadings. The second factor was loaded with TSS (positive loading) and P (negative loading). SG1 was separated from SG2 and SG3 on the basis of higher values of E .coli, NH 3 -N and COD in SG1. Similarly, SG2 and SG3 were separated from SG1 on the basis of higher values of TSS (Since P negatively loaded on the second factor). Looking at the water quality trend from a macro perspective, the WQI scores (Fig. 4) indicated that SG1 is the most polluted water body breaching the Class IV threshold during the monitoring period, followed by SG3 and SG2. BOD, COD and TSS at SG1, however, were between Class II and Class III of the NWQS. On the other hand, NH 3 -N was relatively high between Classes IV and V of the NWQS, which ultimately contributed to the WQI, breaching Class IV threshold. This was hardly surprising, considering that most of the in-stream flow is augmented, which propels biodegradation producing NH 3 -N. Flow retention also caused the in-stream DO levels to be between 1 and 2 mg/L (Class IV), which may be due to the lack of natural re-aeration (Haider et al. 2013). Flow retention inhibits pollutant washout, resulting in Fig. 3 The scatter plot of the first and second factor scores obtained from the factor analysis of water quality parameters across the sampling locations Appl Water Sci particulates settling on the underlying benthos, at the same time prolonging the presence of organic and inorganic constituents within the water body. The sedimentary oxidation further depletes DO levels in the water column. Although the water quality condition at SG1 can already be considered as polluted, further volumetric addition from additional discharge, for example from storm water, would likely result in severe pollution. In addition, the continuous and large volumetric discharge would have nowhere to go as the manual floodgate in the channel is closed most of the time. This presents the possibility of odor problems due to anaerobic decomposition, which may be a nuisance to the residents.
SG2 is located on the upstream segment of the river and water quality conditions here can be considered to be moderate, with a consistent Class III WQI rating. The relatively low pollutant levels indicated that the upstream land-use activities appear to exert a marginal impact on the in-stream water quality. The water quality status at SG2 (main stem) is observed to be between Classes II and III of the NWQS for majority of the constituents measured. DO levels, however, are still relatively low as reflected in the low in situ DO readings. SG3 is located near the mouth and reflects the impact of the side channel towards the water quality (SG1), postconfluence. Despite this fact, there is a very clear concentration increment post-confluence with the channel. The overall WQI score depletes to a lower value most of the time while still being within the Class III denotation. Some increment in organic levels (BOD and COD) was observable, though the most significant increase in SG2 was for NH 3 -N, P and E. coli.
The OPI (Fig. 5) can be relied on as an immediate measure of the overall health of a river ecosystem. The mean OPI at SG3 (4.11) lies between that of SG1 (14.99) and SG2 (2.57). SG3 is most susceptible to flushing by sea water and is normally expected to be least polluted. But, SG3 also receives the pollutant load from the side channel. SG2 is located upstream so that the effect of pollutant from the side channel is not seen at SG2. Finally, the land-use land cover classification of the study area shows that the region near the sampling sites is dominated by urban settlements compared to other classes. The image is classified into four major classes, i.e., vegetation, water bodies, barren land and urban regions (Fig. 6). Urban regions (Table 9) constitute about 106.70 sq km in the total studied area of 321 sq km, followed by vegetation (112.27 sq km), barren land (87.60 sq km) and water bodies (14.43 sq km). The urban area comprises 33.23 % of the total studied area. This is a major sign of anthropogenic activities which can be responsible for water quality deterioration in the river. Once the urbanization process commences, it will expand into nearby areas also. Urbanizing streams pose particular  challenges for management given an inherent changing nature (Chin 2006). If unchecked, the present pollution problem will evolve into a serious environmental problem in the near future and will affect the delicate environmental balances and interactions.

Conclusion
The short-term water quality trend is analyzed for an urban river. The low DO levels will lead to the mortality of aquatic flora/fauna, the decomposition of which further decreases the availability of already depleted DO levels. This can eventually drive the riverine ecosystem from oxic to hypoxic and then to anoxic conditions. Understanding the water quality trend is a prerequisite in sustainable management of the river water and adjacent ecosystem. This also governs the self-purification capacity of the river and will aid in following the national water quality standards and keep the ambient conditions of the river. The present study is also a base for future water quality modeling studies for predicting long-term changes in the era of climate change and increased anthropogenic pressures.