Analysis of physicochemical water quality parameters for streams under agricultural, urban and forest land-use types: in the case of gilgel Gibe catchment, Southwest Ethiopia

Streams draining to Gilgel Gibe catchment cross agricultural and urban land uses receiving a different pollutant that challenges water quality. A total of 21 sampling sites were selected from seven streams of agricultural (n = 3), urban (n = 3) and forest (n = 1) land-use types. Composite samples were collected from upstream, middle and downstream of all land-use types. Twenty-three physicochemical parameters were measured from each sampling site. Temperature, DO, pH, EC, turbidity, width, depth, current velocity and discharge were measured onsite. Two milliliters of unfiltered water samples was collected from every site for laboratory analysis. Mean of NO3-N was highest in agricultural streams than forested and urban streams. In contrast, mean of SRP, NH4-N and COD and BOD5 concentration was greater in urban streams followed by agricultural streams, whereas forest streams are lowest. Concentrations of nutrients, EC and turbidity were recorded in increasing manner across land-use gradient from forested to agricultural and urban streams. The analysis of one-way ANOVA showed that all physicochemical parameters were significantly different among all sites with different land-use types (P < 0.05), except for water temperature (one-way ANOVA: F = 0.987, P = 0.494). NMDS and cluster analyses have discriminated the sites into three groups of land-use types. Then, we conclude that water quality of urban stream is highly impaired than agricultural streams, whereas forested streams have better water quality. Therefore, stream restoration projects, reforestation, conservation of riparian vegetation appropriate waste disposal need to be encouraged in the study area for sustainable management of freshwater resources.


Introduction
Freshwaters are precious and very scarce resources (Gardner and Engelman 1997;Rock 1998;Ridoutt and Pfister 2010; Arjen and Hoekstra 2010) on which our daily life is based in order to sustain on this planet. Nowadays, the fundamental challenge of freshwater quality is derived from anthropogenic activities such as industrialization (Ebenstein 2012;Teng et al. 2011;Temesgen and Seyoum 2018), urbanization (Luo et al. 2017;Ren et al. 2003;Pires et al. 2015) and intensive agricultural activities (Buda and DeWalle 2009;Hall and Leavit 1999;Lenat 1984;Ren et al. 2003). Conversion of land from one land-use type to other land use can adversely affect the quality of freshwater (Bernot et al. 2010;Yu et al. 2016;Likens and Bormann 1974).
Modification of land surface for human use from forest to agriculture and then from agriculture to urban is common practices in developing countries (Pires et al. 2015;Buda and DeWalle 2009;Yu et al. 2016;Ren et al. 2003;Wallace et al. 1997;Mulholland 1992) like Ethiopia. These have drastically inflicted water quality parameters (Kail et al. 2012;Hall and Leavitt 1999). As streams pass through different land-use types, they receive various chemical contaminants (Ding et al. 2016;Bernot et al. 2010;Likens and Bormann 1974). The pollutant types that worsen water quality may vary for different land-use types (Yu et al. 2016;Zhou et al. 2016;Bu et al. 2014). These depend on the purpose 1 3 234 Page 2 of 8 for which land is used like crop cultivation, pastoral, healthcare centers, residential places and institutions that discharge wastes without enough treatment (Chang 2008;Waziri and Ogugbuaja 2010;Mereta et al. 2012).
Overgrazing and crop cultivation lead to poor water quality (Yu et al. 2016;Mereta et al. 2012;McGrane 2016;Luo et al. 2017). Chemical residues such as manure, fertilizers and pesticides (De Gerónimo et al. 2014;Mekonen et al. 2016) from intensive cultivation of crops percolate to groundwater and move to streams posing health risk. Fertilizers have adverse impact on water resources as increased nutrient concentration can cause pathogens and vectors to disseminate (Wondim et al. 2016;Buda and DeWalle 2009). In addition, increased concentration of nutrients such as phosphate and nitrate in stream directs to excessive algal growth, which deteriorates water quality, poses health risk, depletes dissolved oxygen and increases water temperature (Lenat 1984;Li et al. 2008).
As urbanization is expanding with population growth, different organic and inorganic chemical compounds are discharged to rivers (Temesgen and Seyoum 2018;Buda and DeWalle 2009;Bernot et al. 2010;Yu et al. 2016;McGrane 2016). Runoff takes oil, grease, household chemicals and domestic wastes to the streams. Due to these streams are threatened by excessive and high loads of organic waste (Wondim et al. 2016;Buda and DeWalle 2009;Zhou et al. 2016). River pollution first influences chemical quality, then followed by disquieting aquatic organisms and delicate food web. Because of these, the physicochemical qualities of streams crossing through different land-use types become major concerning issue for today's authors.
The catchment of Gilgel Gibe River is characterized by fast urban expansion and intensive agricultural activities. Thus, physicochemical qualities of streams crossing these land uses were threatened by human activities at the catchment area (Mereta et al. 2012). In order to assess the extent to which different land-use types affect water quality, we have included samples from urban and agricultural land uses as study site and forest streams as reference sites.
There was no research done on physicochemical water quality of streams draining across agricultural, urban and forest streams. Therefore, the aim of this study was to analyze the physicochemical water quality parameters of streams under agricultural, urban and forest land uses at Gilgel Gibe catchment.

Study area
This study was conducted in Gilgel Gibe catchment streams located in Jimma zone of Southwest Ethiopia (Fig. 1). Data were collected from three agricultural streams (Merewa, Gibe and Gulufa), three urban streams (Kito, Kochi and Awetu) and one forest stream. The study area was characterized by fast-growing urban and intensive agricultural activities. As it was reported by Mereta et al 2012 overgrazing, brick preparation, vegetation removal and land conversion to cropland, drainage, municipal waste discharge and crop cultivation are the major threats of human activities around these study streams. Maize (Zea mays), teff, sorghum, barley and wheat cultivation is commonly practiced at the Gilgel Gibe watershed. Conversely, forest stream was characterized by the absence of human activities. All sites of forest stream have no any proportion of agricultural land use and human disturbance; rather, natural vegetation is restricted to state protection in upper stream. Therefore, overall hills were covered by vegetation and protected by government.

Water sampling sites
A total of 21 sampling sites were taken from seven streams that drain agricultural (n = 3), urban (n = 3) and forest (n = 1) catchments. All steams are permanent and they flow throughout the year. Samples were taken from upstream, middle stream and downstream of all land-use types. From urban streams, Awetu (A1, A2, and A3 sampling sites), Kochi (KO1, KO2 and KO3 sampling sites) and Kito (K1, K2, and K3 sampling sites) were selected. Correspondingly, from agricultural streams, Merewa (M1, M2 and M3 sampling sites), Gibe (G1, G2 and G3 sampling sites) and Gulufa (GU1, GU2 and GU3 sampling sites) were selected. For forest stream, F1, F2 and F3 were chosen. In these symbols, letters refer to the name of streams and numbers 1, 2 and 3 illustrate upstream, middle stream and downstream, respectively.

Water sampling and measurements
Water samples were collected during the wet season of the region from May to June 2017. Onsite measurement was suggested for parameters that will be changed overtime due to chemical reactions or biological changes (APHA 1996). Thus, dissolved oxygen (DO), temperature, pH and electrical conductivity (EC) were measured with multiparameter probe meter (HATCH hd401). Turbidity was also measured onsite using turbid meter (Wag-WT3020). The physical features measured include stream width, depth and adjacent landuse pattern. Current velocity was measured with Vale port's Flow Meter Model 001, and depth was recorded at three to ten evenly spaced points across the wetted width of each site. Discharge was calculated as the wetted width of each site multiplied by its average depth and velocity. Moreover, for incorporation of GIS (Geographic Information System) database, altitude, longitude and latitude sampling sites were recorded using GPS (global positioning system) reading.

Water sampling and laboratory analysis
Composite sampling technique was used to take water samples. Two liters of unfiltered water samples was collected using acid-washed polyethylene bottles and transported to Jimma University Environmental Health Laboratory with ice box (below 4 °C) for analysis.

Statistical analysis
Nonmetric multidimensional scaling (NMDS) was used to summarize variation in physicochemical parameters among sampling sites. Dissimilarity matrices based on the Bray-Curtis coefficients (Bray and Curtis 1957) were employed. Goodness of fit of the ordination was assessed by the magnitude of the associated stress value; a value of < 0.2 corresponds to a good ordination. Differences in water chemistry were compared among sites using one-way analysis of variance (ANOVA) followed by Tukey's post hoc multiple comparison tests. Cluster analysis for all land-use types was made using PAST software. All statistical analyses were performed with Statistica7 (Version7, 2004, Stat Soft, Tulsa, Oklahoma), unless otherwise indicated.

Results and discussion
The concentrations of stream inorganic nutrients were generally increased across the land-use gradient from forested (Table 3) to agricultural (Table 1) and urban (Table 2) streams. Mean of NO 3 -N was highest in agricultural streams (ranging from 4.8 mg/l to 7.58 mg/l) than forested and urban streams (lowest with range of 0.37 mg/l to 2.71 mg/l). These may be because of landscape change leading to nutrient load in streams sediment (Ding et al. 2016;Bu et al. 2014;Turner and Rabalia 2003;Mouri and Oki 2011).
Among agricultural streams, Gulufa stream sites have the highest NO 3 -N concentration when compared with Merewa and Gibe streams. According to Ding et al. 2016 andZhou et al. 2016, these reveal that intensive agricultural activities have caused high nutrient measures in the streams. The concentration of nitrate decreased at downstream sites when compared with upstream site for agricultural streams. As it was reported by Zhou et al. 2016, there can be high nutrient concentration at upper stream but the lower concentration recorded in our result at the downstream may be due to denitrification process according to Bohlke and Denver 1995. From urban streams, the mean of nitrate was highest for Kito followed by Awetu, while the lowest concentration was in Kochi stream. As it was reported by Bu et al. 2014, these have direct relationship with human activities in the city that discharge organic and inorganic waste in Awetu and Mean of NH 4 -N is greater in agricultural and urban streams than forest stream. From urban stream site (Table 2), Kochi was the highest and Kito was the lowest in NH 4 -N concentration. Agricultural streams were greater than forest stream sites and lower than urban streams sites in terms of NH 4 -N concentration. As reported by other authors such as Ding et al. 2016, Bu et al. 2014, Bohlke and Denver 1995and Mulholland 1992, increased human activities lead to  increase in the concentration of ammonium nitrogen in the freshwater, indicating that quality of water under different land-use types is different. From agricultural streams (Table 1), Gibe stream has the highest NH 4 -N concentration when compared with Merewa and Gulufa streams. Upstream of Gibe site (G1 and G2), there is intensive cultivation of crops and there is fertilizer use that can be moved to streams in line with report of other authors (Ranalli and Macady 2010;Mouri and Oki 2011). In contrast, mean SRP concentration was greater in urban streams followed by agricultural streams, while in forest stream lowest concentration was absorbed. According to Zhou et al. 2016;Iseyemi et al. 2018;Jones et al. 2001;and Ding et al. 2016, these is the indicator of land-use effect on water quality.
DO saturation was high in forest stream sites than agricultural streams and urban streams. DO concentration was highest in forest stream sites when compared with agricultural streams that range from 6.57 mg/l to 5.53 mg/l. In line with other authors like Ding et al. 2016 andBernot et al. 2010, this is due to riparian land cover change. Thus, forest streams have high dissolved oxygen and lower water temperature because of riparian vegetation cover. For urban streams, DO ranges from 7.11 mg/l to 0.73 mg/l higher in Kito stream, medium in Kochi stream and lower in Awetu streams. It shows decrement downstream in all sites of the land uses. As reported by Bernot et al. 2010, dissolved oxygen is different from land use to forested streams and from upstream to downstream responding to organic pollution and water temperature.
Measures of pollution indicators, COD and BOD 5 , were greater at urban stream than forest and agricultural streams. Correspondingly, the value of COD in agricultural stream was greater than forest stream (Tables 1, 2 and 3). Among agricultural streams, Gulufa stream sites have the highest value of COD when compared with Merewa stream and Kochi stream has the highest value (nearly same to Gulufa stream). BOD 5 was higher in all urban streams and lower in forest and agricultural streams. Other authors such as Mouri andOki 2011, andDing et al. 2016 have reported that wastes discharged from households, hotels and institutions are the main causes of organic waste leading to COD and BOD loads. In line with Chang 2008 and Ranalli and Macady 2010, our result reveals the more severely impaired urban and agricultural streams and then the more organic loads.
Electrical conductivity and turbidity were higher in urban and agricultural streams than forest stream sites as it was reported by Yu et al. 2016;Turner and Rabalias 2003;Ding et al. 2016. Water temperature varies among all sites of land uses, and the highest value was observed in agricultural stream and urban, while reference stream was similar to agricultural streams in line with Bu et al. 2014;Zhou et al. 2016;Likens and Bormann 1974;Bohlke and Denver 1995;and Mulholland 1992. The analysis of one-way ANOVA showed that all of the parameters were significantly different among all sites with different land-use types (P < 0.05). However, water temperature did not differ statistically among sites (one-way ANOVA: F = 0.987, P = 0.494).
NMDS ordination analysis discriminated the studied streams into three groups (Fig. 2). All sites of agricultural and forest streams were similar, while sites of urban streams were categorized into two groups: A2, A3 and KO1, KO2, KO3. Physicochemical parameters such as BOD 5 , SRP, current and turbidity were directly correlated with A2 and A3, while ambient temperature, ammonium, COD, EC and discharge were with KO1, KO2 and KO3. As reported by Chang 2008;Bernot et al. 2010;Bohlke and Denver 1995;and Mulholland 1992, this shows that physicochemical parameters such as DO, nitrate, water temperature and depth and wet width were discriminating variation among agricultural streams from urban sites.
The cluster analysis was performed by using major physicochemical parameters, which shaped three main clusters (Fig. 3). Clusters clearly discriminated water quality among the dominant land-use types running along the study catchment. The forest stream and agricultural streams, namely Merewa, Gibe and Gulufa, were clustered together. The second cluster comprises only urban streams (Kito and Kochi), whereas the third cluster separately contains Awetu which is a heavily impaired urban stream as compared to reference sites. It was reported by Roy 2018 that agricultural and urban streams have similar physicochemical parameters. But here agricultural streams resemble more forested streams (reference site) than urban streams, indicating that urbanization severely impaired water quality than agriculture. Therefore, cluster analyses have definitely illustrated the difference in study streams in physicochemical water quality due to landuse types.

Conclusion
The physicochemical water quality of study streams decreased from forest to agricultural and urban streams, respectively. The highest concentration of NO 3 -N was measured in heavily polluted agricultural stream sites (Gulufa stream site). In contrast to agricultural stream sites, concentration of NO 3 -N for urban stream is greater for slightly affected sites (Kito), while heavily polluted sites (Awetu) have moderate concentration. NO 3 -N concentration has no significant change from upstream to downstream of forest and agricultural streams, but for urban steams it has shown increment from upstream to downstream. Greater mean NH 4 -N was measured at agricultural and urban stream than forest streams. Concentration of SRP is greater for urban stream sites followed by agricultural streams; conversely, dissolved oxygen is higher for forest followed by agricultural streams. The higher measurement for COD and BOD 5 was recorded at urban sites, while the least was recorded for forest stream. NMDS has discriminated the sites into three groups of land-use types. The performed cluster analyses also categorized streams into three main groups based on physicochemical water quality. Then, we conclude that water quality of urban stream is highly impaired than agriculture streams, whereas forested streams have improved water quality. Based on our result, we recommend that inappropriate disposal of wastes in urban area and intensive agricultural activities need to be monitored by concerned body, since they are degrading water quality of freshwater resources.
Stream restoration projects, reforestation and conservation of riparian vegetations should be encouraged in the study area for sustainable management of freshwater resources.