Appraisal of water quality and ecological sensitivity with reference to riverfront development along the River Gomti, India

The conflict between the vitality of natural ecosystem versus artificially developed systems has existed since decades. The ecological sensitivity and socio-economic aspects associated with riverfront development along rivers have attracted the attention of environmentalists and ecologists across the globe. The present study evaluates the impacts of channelization and riverfront development on the water quality of river Gomti through Water Pollution Index (WPI) and other statistical tools. Of the total studied sites, 75% were found to be in the ‘highly polluted’ category even after the development of riverfront. An approximate increase of 274.5% and 171.76% was witnessed in the WPI values at the midstream sites of Kudiaghat and Daliganj, respectively. This increase in the WPI values clearly stated the deteriorated water quality of river Gomti after the channelization. The major issue of domestic sewage discharge with partial or no treatment into the river seems to be unresolved even after a considerable period of riverfront development. This study can provide a reference database toward development of such projects across the globe.


Introduction
The rehabilitative development of urban streams through riverfront construction is a globally employed methodology (Bockelmann et al. 2004;Che et al. 2012;Lu et al. 2019;Mitsch 2014; Thompson et al. 2018). The development of riverfront projects is transforming the freshwater ecosystems through alterations in water flow, fluvial habitat, floodplains, and water quality. The development of riverfront projects has attracted contradictory views with researchers emphasizing on the vitality of the natural state of rivers, while some arguing toward the need of artificial modifications. The lack of balance between channel engineering and ecological perspectives has created complexities for the river systems in the Ganga Basin, India (Dutta et al. 2018). The reduction in the density of fish population consequent to the destruction of potential natural habitats post-channelization has been mentioned in various studies (Blake and Rhanor 2020;Dutta et al. 2018;Jurajda 1995). The fragmentation of habitats into sub-clusters makes the survival of aquatic flora and fauna difficult (Jurajda 1995;Khan et al. 2014). The recolonization of such clusters due to the species sensitivity with respect to habitats has even caused extinction (Collinge 1998). The interdisciplinary understanding of ecology and channel engineering is very vital in such projects. The whole river system is very dynamic in nature; hence, development of a small part of the river further modifies the channel limiting the possibilities of restoration.
The River Gomti, a unique groundwater-fed river in the Ganga Alluvial Plain (GAP), is an important source for drinking water in many Indian cities and towns along the river basin (Shukla and Saxena 2020a, b, c, d). Residents in these cities, especially in Lucknow (capital of Uttar Pradesh), are exposed to various pollutants from the point and nonpoint sources including domestic sewage, industrial effluents, agricultural, and livestock waste which are complex to monitor, assess, and control Tangri et al. 2018). The population of Lucknow city expanded by approximately 38% from 2001 to 2011, which had caused an excessive increase in the water demand for drinking purposes. The stress on the river is also highlighted by the fact that roughly 60% of the total water demand of 415 MLD (millions of liter per day) is fulfilled by river Gomti and remaining 170 MLD by subsurface water resources (Goel et al. 2018).
Hence, to augment towards the increasing demand and deteriorating water quality, the 'riverfront development' project for river Gomti was initiated in April 2015, which was completed in March 2017 in the Lucknow city. It included straightening, narrowing, and lining of the river channel, along with development of intercepting drains for efficient management of sewage. The beautification of the riverbank with parks was also included in the project.
Very limited studies discuss the impact of riverfront development and channelization on the water quality of the river Gomti. A study by Dutta et al. (2018) explained the need of maintaining environmental flows, ecological balance, and changes in water quality after the riverfront development. However, the study was more focused on morphological aspects of the riverfront development project, and no comparisons were made for water quality during preand post-channelization. Thus, the present study assesses the detailed impact of riverfront development on water quality of the river pre-and post-channelization. Because the riverfront development does not cover the entire stretch in the city of Lucknow, the present study discusses the water quality along various stretches of the river. An upstream site located at ~ 30 km from Lucknow city was also selected to present a comparative view of the river water quality, at sites with and without riverfront development. The primary aim of this study was to assess the water quality of river Gomti in both the pre-and post-riverfront development phases. Various pollution zones are also identified through Q-mode hierarchical cluster analysis and represented through spatial maps. Finally, the source and inter-dependability of various water quality parameters were also assessed through Pearson's correlation analysis.

Study area
The Gomti River Basin (GRB) exists between the 80°00′-83°10′ E longitude and 24°40′-28°40′ N latitude (Khan et al. 2021a, b, c). It forms the northwest part of the Indo-Gangetic alluvial plain having a catchment area of 30,437 km 2 . The GRB is categorized by two primary geologic units of Quaternary age (younger and older alluvium). The younger alluvial plain stretches along the river Gomti and forms a wide flood plain which supports agricultural activities throughout its stretch. The older alluvium occurs at higher elevation and is characterized by alluvial deposits comprising sand, clay, and kankar. The channel sediments of the river consist of quartz, feldspar, rock fragments, mica, and illite (major clay mineral). The origin of river Gomti is from 'Fulhar Jheel' at Madhotanda in Uttar Pradesh, and it meets the river Ganga after meandering through an approximate stretch of 960 km. It has two major tributaries, i.e., the Saryu and Sai rivers, which join the river in Mohamadi and Jaunpur, respectively.
River Gomti divides the city of Lucknow, into two parts (cis and trans), and water from the river is lifted at Aishbagh waterworks for water supply in the city. The climate of Lucknow city varies between semiarid to subtropical monsoon type with a mean annual rainfall of 963 mm (Goel et al. 2018). The river flows from N-W direction toward S-E direction across the city. The river water is a major drinking water source and also fulfills requirements of industrial and agricultural sectors in Lucknow city. However, increasing population and water demand have contributed towards the deterioration of river water. Furthermore, the untreated domestic and industrial wastewater is released through 26 major drains causing a further decline in the water quality (Khan et al. 2021c). The riverfront development project has potentially been unable to tap these drains and still a large quantity of untreated wastewater is released into the river, which poses human health risks and endangers the aquatic population in the river as well.

Sample collection and analysis
Twenty-four water samples have been collected from river Gomti across the study area, with three samples from each location (one at each bank, and one from the center of the river). The present study is in continuation with the research conducted by Goel et al. (2018) in the year 2015 (before riverfront development), to evaluate the current status and impacts of riverfront development project on river Gomti. The sampling was done in the post-monsoon period of year 2019 at sites S1 (Chandrika Devi), S2 (IIM Road), S3 (Kudiaghat), S4 (Daliganj), S5 (Lakshman Mela), S6 (Bhaisakund), S7 (Dilkusha), S8 (Shahid Path) (Fig. 1). The channelization of the riverbanks starts from S3 at one of its banks, and from S5, it continues at both the banks of the river. The water samples were analyzed for pH, dissolved oxygen (DO), biological oxygen demand (BOD), total coliform (TC), potassium (K + ), calcium (Ca 2+ ), phosphate (PO 4 2− ), nitrate (NO 3 − ), and fluoride (F − ) to compute the WPI. The water samples were analyzed in accordance with the standard methods of APHA (2012). The water quality parameters pH and DO were measured at the site using potable pH meter and digital DO Meter (Lab Junction, LJ-831), respectively. BOD was computed using the Winkler's method in the laboratory, while NO 3 − and F − were analyzed using advanced research grade benchtop ion selective electrode (Hanna, HI5522). K + , Ca 2+ , and PO 4 2− were measured using ion chromatography (Metrohm 850).

Quality assurance of results
Analytical grade chemicals (purity >97%, from Sigma-Aldrich) were used for complete analytical procedures, to maintain the quality control and assurance (QA/QC) standards. Further, all the glasswares were soaked for one hour in diluted nitric acid (a 1% nitric acid solution) and subsequently rinsed with distilled water. To remove the analytical inaccuracy, duplicate samples were collected and analyzed for each sampling location. The water samples were analyzed using calibrated equipment with acceptable uncertainties. An internal quality assurance system is applied regularly at the laboratory. During analyses, different concentrations of standard solutions were measured as an internal quality control. The relative standard deviation (RSD) values were less than 5%.

Water pollution index
Water quality is a quantification parameter for the nature of water concerning the anthropogenic requirements (Bempah and Ewusi 2016;Shukla and Saxena 2020a). Water quality assessment is primarily done using various physical, chemical, and biological parameters (Shukla and Saxena 2020b). The evaluation of the status of water pollution of rivers is a critical and emerging area of interest around the world, requiring data collection, assessment, and interpretation (Shukla et al. 2020c). There have been many approaches for assessing the overall water quality (Gorgij et al. 2019;Li et al. 2018;Su et al. 2019Su et al. , 2020Tian and Wu 2019). The efficacy of WPI is higher in comparison with that of conventional indexing methodologies. The use of varying weights and theoretical ideal value of any parameter can possibly affect the water quality indexing approach. The unsegregated approach of WPI is more accurate considering the conversion of input variables to a single index value. Hence, the slightest variation in the input parameter concentration can affect the WPI category of water quality.
In the current study, nine water quality parameters including pH, DO, BOD, Total Coliform, F − , Ca 2+ , K + , NO 3 − , PO 4 3− were selected for pollution load computation. The number of water quality parameters n can vary considering the flexible approach adopted in the WPI methodology.
The first step in pollution load computation (PLi) of i th parameter is done using equation below: where Ci is the observed concentration, and Si is the standard or highest permissible limit for the respective parameter. The equation for PLi computation is different for various pH values, i.e., pH value 7 is considered as neutral, but < 7 or > 7 pH values are supposedly detrimental. With this view following equations are recommended for different pH ranges: If pH value is < 7, then Eq. 1.1 is recommended where Si a is the minimum acceptable pH value, i.e., 6.5, and Eq. 1.1 is used: If pH is > 7, then Si b is the maximum acceptable pH value, i.e., 8.5, and Eq. 1.2 is used: Finally, the pollution status of a water sample or water pollution index (WPI) with n number of variables is calculated through aggregation of PLi values and finally dividing with n, Eq. 1.3 is used.

Statistical analysis
The temporal/seasonal variation in the dataset from pre-riverfront development to the dataset in post-riverfront development was evaluated with the nonparametric Mann-Whitney U test (Woldeab et al. 2018). The spatial variation in water quality parameters among different sampling sites was evaluated with the Kruskal-Wallis H test (Fatema et al. 2014). Furthermore, Pearson's correlation matrix and hierarchical cluster analysis (HCA) were performed using OriginPro 2020b software package (OriginLab Corporation, Northampton, USA).
Two correlation matrices were generated, each for the dataset from 2015 to 2019, respectively. The Pearson's correlation analysis is a widely used tool which estimates the linear dependence between various parameters (Batabyal and Chakraborty 2015;Khan et al. 2021b;Wu et al. 2014Wu et al. , 2020Li et al. 2019;Ren et al. 2021). The value of Pearson's correlation coefficient, 'r,' lies between ± 1, suggesting the positive or negative correlation, and there is no correlation between the parameters when 'r' is zero. Moreover, when 'r' lies between ± 0.9 and ± 1, a 'very strong' correlation exists between the parameters. Similarly, a 'strong' correlation exists if values of 'r' vary between ± 0.76 and ± 0.89, a 'good' correlation is there when the values of 'r' lie in the range of 0.51 to ± 0.75, and the correlation is called 'poor' for 'r' values of 0 to ± 0.50 (Batabyal and Chakraborty 2015).
The datasets of water quality parameters from the 2015 to 2019 were further subjected to Q-mode HCA (Q-HCA). HCA uses either similarities or dissimilarities between the parameters within the datasets and classifies the dataset into several clusters according to the distance of similarity/ dissimilarity between these clusters (Adimalla et al. 2020;Elumalai et al. 2020;Loh et al. 2020;Shukla and Saxena 2020d). Q-HCA helps in classification of the monitoring sites based on their similar chemical composition geochemistry, which suggests a probable origin of contaminants (Zhu et al. 2017). For Q-HCA, the Ward's method with squared Euclidian distances was used, which is considered to provide the best results. PLi

Variation in water pollution index post-channelization
The study by Goel et al. (2018) highlighted the impact of untreated domestic sewage discharged in the river and discussed the possible inefficiency of sewerage systems, suggesting possible remediation measures. The riverfront development project included straightening and shortening of river channel, affecting its width, shape, and riverbed, and comprising the construction of a diaphragm wall on both banks along a stretch of ~ 8 km (Dutta et al. 2018).
The water samples at all sites fell in the 'highly polluted' category during both pre-and post-channelization surveys. The increased WPI values after riverfront development clearly stated the deteriorated water quality of the river. A very steep increase in the WPI values in the post-riverfront development phase was witnessed at all sites in the order S3 > S4 > S5 > S8 > S6 > S7 > S2 > S1 as illustrated in Fig. 2. The presence of only one water sample in the 'good' category water quality after the development of the Gomti riverfront was seen in the study. The highest deterioration in the water quality was seen at site S3 (~ 275%) signifying the possible influence of small-scale local manufacturing units, dyeing of textile, washing clothes, etc., in the vicinity. The water quality at site S4 showed a significant deterioration, and WPI value increased by ~ 172%, signifying the direct discharge of untreated domestic sewage into the river. The WPI values at sites S5 and S6 showed an increase of ~ 70.1% and ~ 38.6%, respectively. Various cremation-related rituals performed in the vicinity of both sites and consequent anthropogenic influence emerge as possible causes of deteriorated water quality. The WPI of the water sample at site S7 showed an approximate increase of 24.03% after the channelization. The WPI value at downstream site S8 increased by 40.26% after the channelization. The upstream site S1 located ~ 30 km from the city and is not channelized, and it showed the minimum increase in the WPI value (~ 5.76%). The WPI value at this site was 0.71 which felt in the 'good' water quality category.

Cluster analysis
For a better representation and understanding of the impacts of channelization, a stretchwise assessment was required. Hence, the variation in water quality parameters is assessed as per the various clusters obtained through Q-HCA. The stretchwise summary of various water quality parameters and their statistical measures is presented in Table 1. The results of Q-HCA were also plotted as dendrograms for both assessment periods (pre-and post-channelization). The dendrograms present the extent of similarities between various sampling locations, and similar sites are kept in same cluster. The sampling locations were categorized grouped into two clusters for pre-channelization period (Fig. 3a) and three clusters for post-channelization period of assessment (Fig. 3b). For the pre-riverfront development period, Cluster  (Table 1). Similarly, cluster II grouped five sites, Chandrika Devi (S1), IIM Road (S2), Kudia Ghat (S3), Daliganj (S4), and Lakshman Mela (S5), because of similarities between their values of BOD, EC, NO 3 − , K + , and F − , were the reason behind grouping of these sites into a cluster. Cluster II can be considered as 'low pollution' zone based on the mean values of the water quality parameters (Table 1).
In the post-riverfront development period, cluster I had only one site, viz. Shaheed Path (S8), which was the most polluted site across the whole stretch considered in this study, with maximum values of all the parameters reported at this location (Fig. 3b). Cluster II had two sub-clusters, with Daliganj (S4) and Lakshman Mela (S5) in one group, and Kudia Ghat (S3), Bhaisakund (S6), and Dilkusha (S7) in another group. Cluster II can be categorized as representing 'moderate to high pollution' and has similar values of the water quality parameters at all sites. Further, cluster III can be categorized as the least polluted zone, having Chandrika Devi (S1) and IIM Road (S2) grouped together. Based on the results of stretchwise cluster analysis, it can be concluded that riverfront had a critical deteriorating impact on the water quality parameters. Sites S3, S4, and S5, which were in the low pollution zone in the pre-channelization period, were found to be in 'moderate to high pollution' zone after the channelization.

Pearson's correlation matrix
The significance of spatial and temporal differences between water quality parameters was verified using the Kruskal-Wallis H test and the Mann-Whitney U test, respectively, with respect to the temporal variation between the dataset from pre-and post-riverfront development, and the Mann-Whitney U tests indicated that the distribution cannot be considered as significantly different at the significance level of 0.05. Moreover, the Kruskal-Wallis H test indicated that the water quality parameters were significantly different among the sampling sites with p < 0.05 and a Chi-square value of ~ 70. It can be concluded from these results that various anthropogenic activities can be attributed behind the variation among the water quality parameters from pre-to post-channelization. Further, the correlation matrices for the water quality parameters are presented in Fig. 4a, b. It can be seen that pH was negatively correlated with all parameters except DO in both sampling periods. Similarly, DO had a 'very strong' negative correlation with BOD and TC pre-and post-riverfront development, whereas BOD and TC were very strongly positively correlated, suggesting that sewage is the primary source of pollution in River Gomti. All the cations had 'good' to 'very strong' correlation with EC, suggesting that variation of EC is controlled by these ions in the River Gomti. Nitrate exhibited a 'strong' Fig. 3 Dendrograms for Q-mode hierarchical cluster analysis (Q-HCA) for the dataset from a pre-riverfront development and b post-riverfront development 13 Page 8 of 12 correlation with K + in both sampling periods, suggesting that contribution from agricultural practices and subsequent runoff can be responsible for the occurrence of these ions in the river water. Moreover, Ca 2+ and F − also had a 'strong' correlation with K + , suggesting common anthropogenic origin of these contaminants.

Spatial distribution of WPI
The spatial variation of WPI across different stretches of River Gomti in Lucknow city was determined through inverse distance weighted (IDW) interpolation using Arc-Map 10.3. The WPI variation along the river clearly indicated the increasing pollution load and depleting water quality as illustrated in Fig. 5. The minimum pollution was witnessed at sites S1 and S2 with water sample showing WPI values 0.71 (good water quality) and 0.77 (moderately polluted quality), respectively. The upstream location of the sites emerged as the primary reason for better water quality. This finding was in coherence with the results from cluster analysis in which site S1 and S2 occurred in the least polluted zone, i.e., Cluster III (Fig. 3b). The midstream sites (S3, S4, S5, and S6) showed "highly polluted" water quality which was again in coherence with their grouping in Cluster II by the Q-HCA. Sites S7 and S8 had the worst water quality with maximum pollution load. The downstream location of the sites and sluggish flow of river water because of the development of riverfront prominently affects the water quality of sites S7 and S8. The placement of site S8 in Cluster I confirms our findings from the values of WPI. Thus, it can be clearly seen that the water quality at all sites throughout the stretch after the riverfront development has not witnessed any improvement. The potential causes of deteriorating water quality include the reduced flow in the river caused due to shortening of the river channel and a reduction in channel width through the construction of a diaphragm wall along both banks in a stretch ~ 8 km. The lining of the riverbank has minimized the possibility of adsorption of pollutants on sediments. The extensive development of riverfront infrastructure, including weirs and dams, has caused disruption to the natural river habitats. The reducing aquatic population because of the disruption of natural riffle-pool sequences, and excessive dredging during construction have been highlighted in a previous study (Dutta et al. 2018). The reduction of fish biomass in river Gomti has been reported after channelization highlighting the negative impact of excessive dredging and disturbing river banks. The reduced pace of reinstatement of fish population and even no recovery has been reported after channelization in various studies (Kennedy and Turner 2011;Sharan 2016;Vaughan and Ormerod 2010).

Riverfront development and its impacts on river restoration
Riverfront development and subsequent restoration projects across the globe aim at the maintenance and improvement of riverine ecosystem goods and services. But the maintenance of balance between environment, ecology, and channel engineering should be considered in such river restoration projects. The dominance of channelization engineering can cause major deterioration of the water quality of rivers, as seen in river Gomti. Modifications of ecology, floodplains, and other key fluvial characteristics as a result of channelization have been witnessed in the river Gomti (Goel et al. 2018). The channelization, covering floodplains with concrete, and filling of wetlands have contributed towards modification of habitats critically affecting the fish diversity. A huge-scale variation in the deposition of sediments in unpredictable ways can be witnessed at sites below channelized river sections (Kennedy and Turner 2011).
River Gomti supports various vegetation patches that perform the function of sustaining processes of channel erosion and deposition. The existence of pools and riffles is vital for the optimal survival of fishes, considering their use as areas for feeding, cover, and breeding. The tendency of pools to scour at high flow and fill at low flow and vice versa in riffles is very important in the maintenance of poolriffle sequence morphology. Such sequences, in turn, act as dwellings for some organisms, which consequently serve as food for other organisms (Keller 1978). The removal of pool-riffle sequences from the river Gomti has occurred as a result of the heavy channelization. The dredging action has led to degradation of natural habitats, remobilization of contaminants, increases in suspended sediment concentrations, and uneven sedimentation. Pool-riffle sequences are composed of various materials, which include inputs from benthonic species, including fishes (Jurajda 1995). The development of riverfront modification of channel width for river straightening led to the loss of natural habitats. The channelization of the river Gomti has also led to a scarcity of lentic zones and side arms with aquatic vegetation. The diaphragm wall across the stretch of River Gomti has removed some fish habitats considering the absence of required mild slope shorelines. The local fisherman, during an interaction at site visits (Fig. S1), expressed immense worry and trouble considering the reduction of the fish population. The role of in-stream habitats in the maintenance of heterogeneity of fish population is very vital. The absence of many primary fish habitats after the river channelization has also been reported (Dutta et al. 2018). The homogenization of various fish-fauna habitats at a regional scale can lead to a potential loss of future habitats causing depletion of fish diversity. The straightening and the reduction in width of the river have disrupted the natural state of the channel, modifying flow velocity, gradient, and depth of the river channel. The changed flow and constant release of untreated/ partially treated wastewater from various drains across the Lucknow city have directly affected the breeding and survival of aquatic flora and fauna populations. These changes have consequently led to a reduction in the fish population. A report by Ahmad (2013) had also mentioned about the need to strictly avoid reducing the width of the river channel (< 250 m) during the riverfront development and channelization project. However, the design of diaphragm walls seems to lack a coherent approach toward incorporating the suggestions and providing to the maintenance of river ecology. The balance between river ecology and channel development was vital in this riverfront development project, which seems not to have been maintained.

Conclusions
Riverfronts can prove to be a rare resource toward urban development and eco-environmental protection when developed with ecological sensitivity. However, in case of river Gomti the aim towards development of water-land and human-nature interaction zone could not be attained with complete efficacy. An average increase of ~ 81% in the WPI values across all the sites post-riverfront development was assessed in the study. This increase stated the deteriorated river water quality highlighting the unstopped discharge of sewage with partial or no treatment, raising concerns toward the operation and management of sewerage systems. The lowest water quality was witnessed at the downstream site of Shaheed Path with a WPI value of '7.04.' The lowest depletion in the water quality (0.71% and 0.77%) was witnessed at the two upstream sites (S1 and S2) which are notably not channelized. The artificial modification through straightening and decrease in the width of the river channel have affected the natural flow along with the associated selfcleaning capacity of River Gomti. Thus, the current study highlights the need of a balanced and coherent approach with ecologists, environmentalists, and civil engineers working together considering all aspects associated with the impacts of riverfront development on the natural state of surface water resources.