Water Quality Assessment by Measuring and Using Landsat 7 ETM + Images for the Current and Previous Trend Perspective : Lake Tana Ethiopia

Recently there are signs of water quality impairment in Lake Tana, the largest fresh water in Ethiopia. The lake is the growth corridor of the government and supports millions of livelihood around. In order to sustain the benefit and maintain the ecosystem of the lake, the lake health has to be kept safe. Therefore monitoring and evaluation of the water quality of lake is very vital. This study focuses on current and previous trends water quality of the lake through measurements and Landsat Images near entry of Gumera River. Statistical analysis of the physical (Turbidity and STD and biological (Cha-a,) and chemical (DPC) water quality parameters were done. Linear and non-linear regression models between water quality parameter and reflectance of Landsat 7 ETM+ images were fitted based on band combinations. Pervious trend in turbidity was analyzed based on the regression models. The results showed that reflectance and turbidity satisfactorily result with an R ranging from 0.61 0.68. Form 1999-2014 the turbidity of the lake has indicated an increasing trend. Delta development near the entry of Gumera River has been enlarged by 48% because of an increase sediment inflow. The sign in the decreasing water quality of the lake was attributed to the non-point source sediment and nutrient inflow to the lake with high erosion rate from the watersheds. Measures to reduce the non-point source sediment and nutrient inflow by targeting the source areas (hot spots) in the agricultural watersheds need to How to cite this paper: Moges, M.A., Schmitter, P., Tilahun, S.A., Ayana, E.K., Ketema, A.A., Nigussie, T.E. and Steenhuis, T.S. (2017) Water Quality Assessment by Measuring and Using Landsat 7 ETM+ Images for the Current and Previous Trend Perspective: Lake Tana Ethiopia. Journal of Water Resource and Protection, 9, 1564-1585. https://doi.org/10.4236/jwarp.2017.912099 Received: September 17, 2017 Accepted: November 27, 2017 Published: November 30, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Water quality of surface waters across the globe is impaired by several factors such as the non-point source contaminants through runoff [1].This affects the impairment of water quality and the aquatic life in the surface waters such as lakes [2] [3].Especially in developing country where the starting and existing water quality management is not satisfactory [4] becomes worse.One of the challenges in evaluating the effectiveness of management strategy is to monitor the water quality and analyses the trends in water quality data [4].Which could be aggravated by the lack of local knowledge, and the measurements with a number of water quality parameters to be monitored are too expensive to monitor.In countries like Ethiopia, water quality challenges are emerging and in addition to hardly enough data for monitoring the water quality in surface waters such as Lake Tana, the largest fresh water in Ethiopia.
In Lake Tana, previous studies have indicated that the lake water quality showing signs of degradation [5] [6].The rapid loss of vegetative cover and land use change, due to recent agricultural intensification, is a major driver of recent increases in sedimentation and biodiversity degradation in the lake [7] [8].The degradation of water quality has reduced the lake's capability to support aquatic life, and thus, impacts the sustainable utilization of the lake water resources.The incidental occurrence and spread of water hyacinths since 2011 [9] is attributed to the increased sediment and nutrient inflow into the lake.Studies indicate that eutrophication from anthropogenic drivers is the main cause for rapid spreading of water hyacinth over reservoirs [10].
The existing water quality studies such as [4] [11] in Lake Tana used grab samples by considering the smaller portion of the lake near Bahir Dar town and near the entrance of one river.The studies revealed water quality of lake is getting enriched with pollutants such as higher nutrient inflow reduced water quality of the lake water.However, these studies only cover a portion of the lake and hence, overall assessment of agricultural non-point sources and its effect on sediment and nutrient inputs remains challenging and requires a robust spatio-temporal assessment.Basing the assessment purely on grab samples will be costly as it requires a large number of samples both in space and time.However satellite imagery in combination with a robust sampling scheme could reduce the sample size and provide an alternative for water quality monitoring at a range of temporal and spatial scales [12].
Worldwide remote sensing imagery has been used to evaluate water quality in lakes by [13] [14] [15] [16] Various satellite images from Landsat 5, 7, 8, MODIS, and AVHRR have been used with Landsat images being one of most frequent ones [17].The water quality parameters evaluation using images (with different type of satellite and resolution) include, study of sediment transport [18], Chlorophyll-a for detecting over surface water [19] [20] [21] [22] [23] cyanobacteria and microphytes [24], suspended solids [25] and water clarity [11] [26] [27] [28] [29] [30].All of this study has been carried out of Ethiopia where fresh water quality is becoming a concern and monitoring has not been advanced as required.Specifically use of satellite imagery to assess changes in water quality in Ethiopian freshwater systems has been limited.The only existing satellite image based study in Ethiopia used MODIS/Terra images [31] to evaluate turbidity and total suspended solids.While using MODIS images gives the advantage of higher temporal resolution (i.e.daily), the coarser spatial resolution of the images renders them less suitable to map spatial distribution of the water quality parameters.The use of Landsat 7 ETM+ images such as with 30 m by 30 m resolution will help to overcome such shortcomings.
This study develops a method to continuously monitor the water quality (i.e.turbidity, Secchi disc Transparency Depth (STD), Dissolved Phosphorus Concentration (DPC), and Chl-a) in the lake using Landsat 7 ETM+ imagery.Furthermore, it uses the established protocol to evaluate the current situation in the lake as well as the historical changes in sediment load at the river deltas feeding into the lake.This assessment is crucial to understand the sustainability of the lake's ecosystem as it is an important growth corridor for the country's economy as it fulfills services such as transportation, hydropower, irrigation, fishery and ecotourism.

Description of Study Area
Lake Tana (Figure 1) is located northwest of the Ethiopian highlands.Geographically the lake lies between 278,946 meter E, 1,280,138 meter N and 349,403 meter E, 1,360,638 meter N. It is the largest fresh water resource in Ethiopia and contains half of the county's fresh water reserve [32] covering 20% of the drainage area of the 15,096 km 2 Lake Tana basin [33].The rural populations as well as those from Bahir Dar town are highly dependent on revenues from the lake economy as it fulfills services such as transportation, hydropower, irrigation, fishery and ecotourism [34].
It is a shallow Lake with a mean depth of 9 m and maximum depth of 15 m [35].Usually shallow lakes are susceptible for water quality impairment [36] and carrying out this study in Lake Tana vital.The Lake has an average temperature of 22˚C and mean annual rainfall of 1450 mm/year.The lake is located in a wide depression of the Ethiopian basaltic plateau and bordered by flood plains that northwest [32].The Lake is fed by more than 40 small streams [37] [38] [39] and receives 93% of its volume from the major rivers of Gilgel Abay, Gumara, Ribb and Megech [37].The point sources from Bahir Dar town and non-point sources from the agricultural watersheds draining to the lake, increasing the deterioration of the lake as it decreases water quality due to the increase of sediment and nutrient loads [40].

Methods
The general methodology of this study includes measurements of water quality parameters (Turbidity, Secchi disc Transparency Depth (STD), Dissolved Phosphorus Concentration (DPC), and Chl-a) over the surface of the Lake Tana from 26 locations near the entry of Gumara River at the moment that the satellite passed over the lake.Landsat 7 ETM+ images were preprocessed and image reflectance was extracted from single bands and band ratios (Table 1).The samples and images were used to establish a regression model between the measured Water Quality Parameters (WQP) and reflectance observed by the Landsat 7 ETM+ images.Historical assessment of turbidity changes since 1999 were derived by using the developed relationship and available archived Landsat 7 ETM+ images.A change in specific water quality parameter (turbidity) was furthermore assessed using the temporal range of satellite images and the developed regression model.

Water Quality Sampling
Sampling was carried out in three campaigns 10 th of Dec. 2014, 19 th of May 2015 and 7 th of Aug. 2015).The sampling dates were spread within a year to capture seasonal variations in sediment and nutrient inflow.Grab samples of one liter from the top 20 cm of the water surface were collected from 26 locations (Figure 2) near the entry of Gumara river (draining from 1354 km 2 watershed).The inter-sample distance was approximately 200 m to avoid repeated sampling over a single pixel (30 m × 30 m) and to ensure that the sampling area covered the spread of a sediment plume.The coordinates of each sampling location were recorded and samples were analyzed for turbidity, STD, DPC and Chl-a (Table 1).

Turbidity and Water Clarity
Turbidity was measured using the portable Hach 2100Q turbidimeter.The instrument was calibrated using the 2, 10, 20, 100, 200 and 800 NTU standards obtained with the instrument.Water clarity was measured as another indicator for lake turbidity using a sechi disk.The disk was submerged using a metal gaged rode or rope and depths where the disk could no long be seen from the surface were recorded.The STD depth was correlated with the turbidity measured at the same position in order to establish a transfer function.This would allow predicting turbidity using STD for further decreasing the costs of water quality  monitoring.

Dissolved Phosphorus Concentration (DPC)
Samples for dissolved phosphorus were taken in 10 th Dec 2014, 19 th may 2015.
The preserved water samples of 100 ml were prepared for water quality analysis within 24 hours after the sample has been collected.The collected samples were filtered using 0.25 μm filter paper and analyzed for dissolved phosphorus using wagtech plainest spectrophotometer using the automatic wavelength (wave length of 680 nm for P) with a detection limit of 0.001 mg•l −1 .

Chlorophyll-a (Chl-a)
For measuring Chl-a, the samples of 100 ml were first filtered using GF Whatman 0.2 μm pore diameter filter paper.The material remaining on the filter paper was used to extract Chl-a pigment.The filter paper was grinded manually and extraction of Chl-a pigment was done using 5 ml acetone (98 percent aqueous solution) upon which the sample was centrifuged for 20 minutes at 4000 rpm.The centrifuged samples were placed in a cool dark place to settle the suspended material and prevent growth.For the atmospheric correction a simple Dark Object Subtraction (DOS) was applied, which is a common and popular procedure in various water quality monitoring using related remote sensing, due to its simplicity [42] [43].In addition, DOS avoids the use of high demand input data as the case with complex atmospheric correction algorithms.Digital Number (DN) values from the image was converted to reflectance for each of the selected bands (band1-band5) and the ArcGIS 10.1 spatial analyst tool was used to extract the reflectance for each sampling point for the various bands and images.An example on image preprocessing is given in Figure S1.

Non-Linear Regression Model and Evaluate Trends in Water Quality
A regression model between measured water quality parameters (turbidity and  In order to assess the temporal variation of water quality parameters among the 3 campaigns Analysis Of Variance (ANOVA).Specifically, to evaluate differences between sampling periods.Afterwards, the trend of WQP was analyzed by using the Mann-Kendall (M-K) trend test [44] [45] using images Landsat from the archived image from 1999-2014 specifically images archived in December for each sample point.The sample points were clustered into three cluster zones as indicated in Figure 3. Clustering was based on the proximity of water quality parameter (WQP) values and WQP values were averaged to get one regional representative value for each cluster.For instance, the turbidity value for cluster 3 will be based on the weighted average of the turbidity measurements of each of the sample points in cluster 3.
In addition the trends in turbidity estimated from the satellite images will also be evaluated by using 1) analyzing with the trends in annual sediment inflow from Gumara watershed to the lake and 2) estimating identifying the enlargement of delta development near the entry of Gumara River.
The annual inflow of sediments from Gumara River was analyzed by generating the data with the newly developed rating curve for a monsoonal climate by [44].To evaluate how the trends of sediment inflow were in line with the water quality trend, the predicted sediment inflow from 1999-2014 was also analyzed using the Mann and Kendall [45] [46] trend test.
A change in delta development was assessed using the image from 1999 as a baseline and compared against images for similar month (December) in 2014.
Similar periods were taken to ensure similar water levels in the lake as those influence reflectance values and estimation of the delta expansion.Furthermore, a specific line of reference was selected within the delta to compare and estimate the delta expansion from 1999 till 2014.

Spatio-Temporal Variation in Water Quality Parameters
For the sampled months, turbidity has an average value of 17.1, 13.9 and 24.9 NTU in December, May and August respectively.The maximum turbidity values of 37 NTU were observed in August near the shore and the minimum of 5 NTU in May, further from the shoreline.This could be due to that the lake is shallow depth and re-suspension of sediments could also occur with the movement of lake water as the river joins and due to tides in the lake.The results also indicated the turbidity varies significantly temporally among the sampling months within a 95 percent confidence interval (Table S1).The values observed in August, during the rainy season were higher compared to those observed in the months May and December falling within the dry period.This was related to the higher sediment load coming from non-point sources in the intensively cultivated agricultural watershed.
The analysis of the STD data supported the turbidity findings.The average value of STD obtained from the observation was resulted 70.4,84.2 and 50.6 centimeters in December, May, and August, respectively.This result clearly indicated that in dry season the lake water has higher clarity than in the rainy season, which is likely more turbid due to erosion from runoff contributing areas at the upland.The maximum STD of 121 cm was observed in May, 2015 faraway from the shore and the minimum was 36 centimeters observed in August near the shore of the lake.The STD values showed an inverse linear relationship with turbidity both temporally and spatially over the surface of the lake (Figure 4).
Temporally, the ANOVA analysis of STD (p = 7.69 × 10 −12 ) showed that it varies significantly between the months where observation has been carried out (Table S2).
The DPC measured in December, 2014 and May, 2015 showed an average value of 0.17 and 0.23 mg•l −1 , respectively with maxima of 0.47 mg•l −1 (May) and 0.33 mg•l −1 (December).The minimum DPC observations were 0.06 mg•l −1 in December and 0.04 mg•l −1 in May.There is a significant difference between both months of minimum, average, and maximum values.Spatially the maximum DPC was found from the samples collected closer to shore and minimum DPC from the samples further from the shore (from Gumara river entrance to the lake).Spatially the amount of DPC and turbidity was inverse pattern indicating Figure 4. Relationship between turbidity and STD over the surface of Lake Tana near entrance of Gumara River to Lake Tana from December/2014 observations.the concentration of DPC is controlled by the available sediment.In dry period such as May, phosphorous might be released from the bottom deposited sediment of the lake [11] and in wet period, the phosphorous is absorbed by sediment and metallic elements [47].The results obtained were consistent with the observations made by the Tana sub Basin Organization [48].
The results of Chl-a indicated an average value varying from 3.76 µg•l −1 in December to 1.56 µg•l −1 in May.The maximum and minimum values of Chl-a observations were 8.5 µg•l −1 and 0.57 µg•l −1 in October and 3.1 µg•l −1 and 0.26 in May, respectively.The t-test indicated significant variation of Chl-a from wet to dry season within a 95 percent confidence interval.
Based on the Trophic State Index (TSI) classification [49] turbidity and chl-awere averaged for each month.The results indicated that the lake water quality ranges from oligo-trophic to meso-trophic (Table 3).The TSI of Lake Tana could be an indicator of the lake water quality deterioration (Table S3).

Estimating Water Quality through Remote Sensing
The relationship which provided an overall best fit was based on Equation (1) for turbidity and Equation (3) for Chl-a (Table 4 and Table 5).The non-linear regression equation was selected based on the R²between observed and predicted water quality parameter during calibration.A regression for STD was made given that results are expected to be inverse (R 2 = 0.7, Figure 4) to those found for turbidity.
The results shown in the Table 4 and Table 5 indicate that the relationship between the natural logarithm of turbidity with a single band had indicated relatively poor correlation (R 2 of 0.36) for all months except august 2015.In August, the best relationship was based on single band i.e. band 5 with R 2 = 0.66 and found to be the highest (Table 4).The natural logarithm of turbidity with normal band ratios (i.e.band 1 and band 2) provided a more reasonable relationship compared with the single bands for December 2014 (R 2 = 0.68) and May 2015 (R 2 = 0.48).The observed vs. predicted values showed a good correlation for the month of December (Figure 5).The results obtained in this study generally showed lower R 2 values than those obtained by [31] using MODIS.Even though the resolution of the Landsat images was higher than the MODIS images, the lesser R 2 likely due to the following two reasons, 1) MODIS data used by [31] study was acquired daily to develop a relationship between the average daily turbidity and reflectance values for the entire month and 2) the gap correction for the Landsat 7 ETM + images might influence the true reflectance values in parts of the image where the strip occurred.
For Chl-a the relationships using single bands were almost insignificant and nearly zero (data not shown).Using band 1 and band 2 as indicated (Equation (1)) resulted in the best regression model from Landsat 7 ETM+ image reflectance for Chl-a (Table 5).However, the R 2 was much lower as those obtained for turbidity which might be related to the fact that Chl-a absorbs light at lower wavelengths than turbidity.Hence, detecting Chl-a may require lower range bands than those used for band ratios in this study i.e. bands 1 -5.The regression model performed slightly better in December for Chl-a, as this is related to the dry season and the higher production of Chl-a due to increased sunlight hours.

Historical Changes in Turbidity
The Landsat ETM+ images from 1999 till 2013 were obtained for the month of December and used for trend analysis (Table S4).The dates of the downloaded images fell within a 10-daytime period around the 10 th of December, the date of the image used for the regression function (Equation ( 1)).Turbidity was selected to understand potential trends in water quality given the low correlation of the satellite images for Chl-a.An increasing trend in turbidity of lake water with Figure 5. Predicted versus observed turbidity obtained for December/2014.slope of 0.64 (based on the Mann-Kendall trend test) clearly showed that the average turbidity lake water was increased (Figure 6).The average sediment inflow of 38-ton ha −1 •yr −1 into the lake was estimated from [44] which were 46% higher than the country average of 26-ton ha −1 •yr −1 .The calculated annual sediment inflow showed the same trend as those observed in the 3 clusters for the same time period (Figure 6).

Delta Development
Results indicated that the in the delta area expanded by nearly 1.2 km 2 from 1999 to 2014 (Figure 7).This indicates a 48% increase in area in comparison to a study conducted by [50] estimated the area of delta development increase by 2.1 km 2 from 1959-2006.The difference between this result and the latter study was due to the time period difference for analyzing the results in delta development expansion.
Both the annual sediment load inflow increasing (Figure 8) trends and delta enlargement with increasing trend of turbidity from landsat 7 ETM+ image strongly illustrates the decreasing water quality trend.Therefore, to safeguard Predicted Turbidity (NTU) Observed Turbidity (NTU) Figure 6.Historical turbidity trends as predicted from Landsat 7 ETM+ image reflectance over surface of Lake Tana.
the lake water quality best management practices have to be implemented in agricultural watersheds, especially by targeting the hydrologic sensitive areas (HSAs), which are the main sources of non-point sources of nutrients such as phosphorus.In addition, establishing a buffer zone in the surrounding of the lake to sustain the existing wetlands is very important for trapping sediment and nutrient before it reaches the lake to cause eutrophication.

Conclusion
The corrected Landsat 7 ETM+ image could be used for monitoring previous trends current turbidity and suspended sediment and SDT but is less suitable for Chl-a for Lake Tana.The findings of this study indicated that the current water quality of Lake Tana is deteriorating.This could be caused by increased erosion which has resulted in an incremental increase of sediment inflow to the lake.The variation of water quality parameters over the lake was associated with the seasonal variation of regional climate.For example, turbidity is lower in the dry season and is greatest at the beginning of a rainy season.The continuous inflow

Figure 1 .
Figure 1.Relative location and geographical positioning of the study area.

Figure 3 .( 3 )
Figure 3. Cluster zones and sampling positions used in testing the previous trends in WQ.

Figure 7 .
Figure 7. Delta developments near the entry of Gumara for (a) 1999 and (b) 2013 (the blue line indicting the reference line with known ground control point and red line area of delta in both time periods).

Figure S1 .
Figure S1.The preprocessed image of Landsat 7 ETM+ (a) the downloaded image subset by ROI, (b) the gap filled image using PANCHROMA 4.3 software (DN values), (c) atmospherically corrected image (DN values) and (d) Reflectance over Lake Tana area (the gap filled and atmospherically corrected and converted to reflectance).

Table 1 .
Landsat 7 satellite over pass and water quality parameter measurement dates.

Landsat 7 ETM+ Image Pre-Processing and Analysis Landsat
Samples were analyzed four hours after storage using the emission spectrophotometer instrument at 663,664,665 and 750 nm.Chl-a was estimated using Environmental Protection Authority protocol (EPA 150.1 standard for Chl-a measurement [41].Chlorophyll-a was only measured in December 2014 and May 2015.7 ETM+ images with Path/row of 170/52 for 10 th of Dec. 2014, 19 th of May 2015 and 7 th of Aug. 2015 were downloaded from USGS website (Table 2).Sensor scanning and atmospheric correction was applied.Since 31 st of May 2003 strip/gaps in the Landsat 7 ETM+ images occurred (due to failure of scan line off corrector) which are corrected by using PANCROMA 6.42 scan line off corrector tool.The reason to use gap filled corrected Landsat 7 ETM+ images was to analyze the trend in turbidity from the previous Landsat archives. 2.4.

Table 2 .
Descriptive statistics for the water quality parameters measured over the surface of Lake Tana.

Table 3 .
Trophic level index and class for Lake Tana based on the measurement of water quality near the entry of Lake Tana.

Table 4 .
R 2 values for relationship between measured turbidity with Landsat 7 ETM + single band reflectance over the surface of Lake Tana.

Table 5 .
R 2 values for the relationship between measured turbidity and Chl-a with Landsat 7 ETM+ reflectance band ratios over the surface of Lake Tana.