Surface and groundwater quality assessment based on multivariate statistical techniques in the vicinity of Mohanpur, Bangladesh
Md. Mahtab Ali Molla1, Narottam Saha2, Sayed Mohiuddin Abdus Salam1, Md. Rakib-uz-Zaman1
1 Department of Applied Chemistry and Chemical Engineering, Environmental and Tracer Studies Laboratory, Faculty of Engineering, University of Rajshahi, Rajshahi 6205, Bangladesh
2 Department of Applied Chemistry and Chemical Engineering, Environmental and Tracer Studies Laboratory, Faculty of Engineering, University of Rajshahi, Rajshahi 6205, Bangladesh; School of Environmental Sciences, Faculty of Science, University of East Anglia, Norwich, Norfolk NR4 7TJ, United Kingdom
|Date of Web Publication||27-May-2015|
School of Environmental Sciences, Faculty of Science, University of East Anglia, Norwich, Norfolk NR4 7TJ, United Kingdom
Source of Support: University of Rajshahi, Conflict of Interest: None
Aims: This work evaluated the surface and groundwater quality of Mohanpur area, Rajshahi district, Bangladesh. Multivariate statistical techniques were also applied to determine the possible sources of water contamination.
Materials and Methods: Water samples were collected from randomly selected ten different sampling sites for analyzing the chemical parameters including pH, electrical conductivity, total dissolved solids, total hardness, total alkalinity, Cl− , NO3− and some heavy metals such as Mn, Pb, Cd, and As concentrations. Concentrations of heavy metals were determined using atomic absorption spectrometer (AAS).
Results: Based on hydrochemical characteristics, surface and groundwater in the study area were, in general, fresh, hard, and alkaline in nature. All chemical parameters were within the WHO water quality guidelines. Whereas, among four analyzed heavy metals Pb, and Cd concentrations exceeded the WHO recommended values. Pearson correlation matrix showed a number of statistically significant associations (P < 0.01 and P < 0.05) among the examined water quality parameters. Moreover, principal component (PC) analysis (PCA) and cluster analysis (CA) were used to analyze the water quality dataset. PCA analysis identified two PCs as responsible for the data structure explaining 72.53% of the total variance in water quality. PCA indicated that the water quality variations were mainly of anthropogenic origin through agricultural and municipal discharges. Results of CA revealed three significant groups of similarity among the 10 sampling sites.
Conclusions: It could be deduced from the present results that water contamination was occurred to some extent throughout the area, and is likely to be continued in the near future. Improvement of local sanitation system along with frequent training and awareness programs can help in developing water quality in the studied area.
Keywords: Groundwater, hierarchical cluster analysis, multivariate statistics, principal component analysis, surface water
|How to cite this article:|
Molla MA, Saha N, Salam SA, Rakib-uz-Zaman M. Surface and groundwater quality assessment based on multivariate statistical techniques in the vicinity of Mohanpur, Bangladesh. Int J Env Health Eng 2015;4:18
|How to cite this URL:|
Molla MA, Saha N, Salam SA, Rakib-uz-Zaman M. Surface and groundwater quality assessment based on multivariate statistical techniques in the vicinity of Mohanpur, Bangladesh. Int J Env Health Eng [serial online] 2015 [cited 2022 Jan 22];4:18. Available from: https://www.ijehe.org/text.asp?2015/4/1/18/157717
| Introduction|| |
The surface and groundwater quality is a very sensitive issue, and a major factor affecting the human health as well as ecological systems.  The surface water quality is controlled by anthropogenic factors such as urbanizations, industrial and agricultural practices, and also by natural processes like erosion and the regional climatic conditions.  Likewise, the groundwater quality depends on many factors such as soil characteristics, manner of groundwater circulation through rock types, topography of the region and human activities on the ground.  However, water pollution is a critical problem in Bangladesh including Rajshahi city and the surrounding areas like Mohanpur. Most of the surface water bodies in the Mohanpur area are linked with drains, which are connected with the latrines. Moreover, the increased use of agrochemicals such as insecticides, pesticides and chemical fertilizers to accelerate the crop productions is intensifying the water pollution of the present study area.  According to the Rajshahi development authority report,  release of untreated urban effluent into the surrounding water bodies and the subsequent use of this polluted water in agricultural fields is immensely affecting the public health. Although the human health hazards caused by heavy metals associated with polluted water have been known for a long time, exposure to heavy metals is increasing in particular in less developed countries like Bangladesh. For example, manganese is a known mutagen. The chronic ingestion of Mn in drinking water is associated with neurologic damage.  Lead is considered as a possible human carcinogenic due to the inconclusive evidence of human and sufficient evidence of animal carcinogenicity.  Acute exposure to lead is known to cause renal failure and liver damage.  Moreover, recent research has shown that prolonged low-level exposure to lead may diminish the intellectual capacity of children.  Cadmium in our environment is a matter of concern since 1960s, when a painful bone disease "itai-itai" was reported in Japan. Moreover, it has recently been shown that Cd acts as an endocrine-disturbing substance and may lead to the development of prostate cancer and breast cancer.  Arsenic, a deadly poisonous metal, is unique among the metalloids and oxyanion-forming elements (e.g., Se, Mo) in its susceptibility to mobilization in various forms under the pH conditions of 6.5-8.5. The arsenic contamination of groundwater is among the most challenging environmental problems nowadays, threatening the well-being and livelihood of millions of people in South and Southeast Asia.  In Bangladesh, 30 million people drink arsenic-contaminated water without having alternative resources.  According to the estimate of the WHO, in the future, chronic consumption of such toxic water may lead to 1 in every 10 adult deaths caused by arsenic-related cancer.  Arsenic-contamination of natural water leads to the development of cancer, cardiovascular disease and inhibits the mental development of children. 
It is worth noting here that depending on the geological formation, the quality of water varies from place to place. Thus, a number of physicochemical parameters that results in a large data matrix need a complex data interpretation. Principal component (PC) analysis (PCA) and cluster analysis (CA) coupled with Pearson's product moment correlation analysis allow us to resolve this complex data matrices. These multivariate approaches help to identify possible factors/sources that influence the water systems and offers a valuable tool for reliable management of water resources, as well as a rapid solution on pollution problems. ,, In general, CA groups the samples based on inter-sample similarities and illustrates the overall similarity of variable data set. On the other hand, PCA identify the most important factors contributing to the data structure by transforming many original, interrelated variables into fewer, uncorrelated variables named PCs. Usually, the minor PCs are not taken into consideration in order to simplify the analysis because of their poor interpretation of data structure.  A number of reports emphasize on the importance of multivariate statistical analysis as a tool in the treatment of analytical and environmental data. ,
With these views, the present work was designed and carried out for the following purposes: To study the chemical properties and heavy metals concentrations of surface and groundwater resources in Mohanpur area of Rajshahi, Bangladesh; to compare the present status of water quality parameters with the recommended guidelines; and to extract information relating to the influence of possible sources (natural or anthropogenic) on water quality parameters using multivariate statistical analysis.
| Materials and methods|| |
The present study area, Mohanpur, belongs to Rajshahi district is situated on the northern bank of the Padma River, Bangladesh. The distance between Rajshahi city and Mohanpur town is about 30 km. It covers an area-extent of 163 km 2 and is distributed in 24,828 units of household. It experiences a tropical wet and dry climate and is generally marked with monsoon, high temperature, considerable humidity, and moderate rainfall. Ten villages were selected for the present study and the area map indicating the sampling sites is shown in [Figure 1].
|Figure 1: Sampling spots showing Mohanpur Upazila (sub-district), Rajshahi district, Bangladesh|
Click here to view
Twenty water samples (10 groundwater samples from hand tube-wells and ten surface water samples from different ponds or canals) were collected from predetermined sampling sites in the month of April to December, 2010. Samples were collected from Ghasigram, Dhurail, Borokuri, Mohanpur, Jahanabad, Raighati, Keshari Hat, Bakshimai, Dhopaghata Hat, and Maugachhi. The tube-wells water around the study area was mostly used for drinking purposes. Clean and dry 500-ml polyethylene bottles were used to collect surface and groundwater samples. Before the collection of the groundwater samples, tube-wells were pumped for 5 min to wash out the stagnant water inside the tube-wells and to get fresh groundwater. After rinsing with fresh water, the sample water was poured into bottles. Prior to collection of the surface water samples, the surface of the ponds or canals were cleaned, and the bottles were rinsed with cleaned surface water. The sample bottles were then labeled, tightly packed, transported immediately to the laboratory, and stored at 4°C for chemical analyses. A portion of collected water samples were acidified (2 ml concentrated HNO 3 per L) and used in the determination of metallic constituents.
Water samples were analyzed for the chemical parameters including pH, electrical conductivity (EC), total dissolved solids (TDS), total hardness (TH), total alkalinity (TAlk), chloride (Cl− ), nitrate (NO 3− ), and also for the heavy metals such as Mn, Pb, Cd, and As. pH of the samples were measured at the sites of collection using a portable pH meter (KRK, KP-5Z, Japan), previously standardized through pH 4, 7, and 10 standard buffer solutions (procured from Merck, India). EC of the water samples were measured as soon as possible after collection using a digital EC meter (HANNA instruments HI 9033, Singapore), previously calibrated through 0.01 M, 0.1 M and 1.0 M KCl solution. The TDS (EC × 0.7) of the samples were recorded by the TDS probe of the EC meter. After determining the pH, EC, and TDS on the sites, the samples were carried to the laboratory for remaining analyses. TH as CaCO 3 , TAlk, and Cl− were determined following standard methods recommended by American Public Health Association.  NO 3− contents were analyzed by a UV-Visible spectrophotometer (WPA S104, England). The concentrations of Mn, Pb, Cd, and As were determined by an atomic absorption spectrophotometer (AAS; Model AA-6800, Shimadzu Corporation, Japan) at their respective wavelength in the central science laboratory, Rajshahi university. Due to expected low concentrations of the metals in the water samples and limited instrument sensitivity, preconcentration of the water samples were done by evaporating 100 ml of water to 8 ml on a hot plate. The samples were then digested by adding 5 ml of 11.1 M HNO 3 and heated on the hot plate for 30 min. Thereafter, 10 ml of 16.3 M HCl was added, and digestion was continued until the solution turned into light brown or colorless, and the volume was then adjusted to 25 ml with double distilled water. In order to assure the precision of the measurement, reference standard solutions with a known concentration of each measured element were used as control samples. After every ten samples, the control sample was analyzed to check the accuracy of analysis. Each sample was measured at least two times in order to assess the reproducibility of the measurement. Samples were reanalyzed if the relative standard deviation of the measurement exceeded 10%. The reagents, including indicators and buffers, were of Analytical Reagent grade. The aqueous solutions were prepared using deionized water.
Analysis of variance (ANOVA) test was performed to evaluate the inter-sample locations variation of chemical parameters, as well as heavy metals concentrations. Moreover, to identify the relationship among the examined water quality parameters in the studied samples and to infer their sources (natural or anthropogenic), Pearson's correlation coefficient analysis, PCA, and CA were performed. In Pearson's correlation coefficient analysis, the values of correlation coefficients indicate the strength of inter-relationship between two chemical parameters. PCA and CA are the most common multivariate statistical methods for classification and interpretation of large datasets from environmental monitoring programs that allow the reduction of the dimensionality of the data and the extraction of information.  Data were processed using Statistical Package for Social Science (SPSS 17.0 for Windows; IBM, USA). Besides this, the STATISTICA 8.0 (StatSoft, Inc., USA) software package was employed for CA.
| Results|| |
The results obtained from the chemical and heavy metals analyses of ground and surface water samples in Mohanpur Upazila (sub-district); Bangladesh are presented in [Table 1]. One-way ANOVA analysis showed that the examined parameters in surface and groundwater samples differ significantly (at 95% confidence interval) except Mn (F = 0.045, P = 0.834) and As (F = 0.172, P = 0.683). Pearson's correlation coefficient analysis identified significant association between measured water quality parameters and the results are given in [Table 2]. Moreover, PCA extracted two PCs controlling the variability of examined parameters in this study [Table 3].
|Table 2: Pearson correlation between chemical parameters of water samples|
Click here to view
| Discussion|| |
The observed pH ranges of the surface and groundwater in the study area are 7.7-8.6 and 6.9-7.1, respectively [Table 1]. The higher pH value of surface water than that of groundwater could be due to increased photosynthetic assimilation of dissolved inorganic carbon by plankton. While pH usually has no direct impact on consumers, it is one of the most important operational water quality parameters. The recorded pH values of ground water samples were well within the WHO  prescribed standard for pH of 6.5-8.5. Whereas, surface water pH values were in line with the upper end of the WHO recommended standard limit for pH. According to environmental quality standards (EQS) for Bangladesh,  the maximum permissible value of pH in industrial water, fishing water, and drinking water are 6.0-9.5, 6.5-8.5, and 6.5-8.5, respectively. So the present values of pH indicated that the examined samples of groundwater were not objectionable for drinking, fish culture, irrigation, industrial and other purposes. However, 80% of surface water samples were in line with the permissible limits of pH recommended by EQS. Scatter plots [Figure 2] have been used for illustrating how the metals solubility varies with water pH and their concentrations as well. The solubility plots indicate a general decrease in heavy metals solubility with increasing pH, which is the usual trend with cationic metals.
|Figure 2: Solubility plots show the relationship between the pH and heavy metal concentrations in surface and groundwater samples|
Click here to view
In this study, the values of EC were ranged between 233 and 645 μS/cm for surface water and from 701 to 987 μS/cm for groundwater indicating high mineralization of groundwater compared to surface water [Table 1]. However, no studied samples were beyond the WHO  recommended maximum EC level of 1,500 μS/cm in drinking water.
Total dissolved solids values were varied between 161.1 and 690.9 mg/L [Table 1]. Low TDS (161.1-451.5 mg/L) was observed for surface water and high TDS (490.7-690.9 mg/L) for groundwater, could be due to the higher contact period of groundwater with rock than that of surface water. According to WHO,  the maximum acceptable concentration of TDS in natural water for domestic purpose is 500 mg/L and the highest permissible limit is 1,500 mg/L. All the studied samples were within the permissible limit of WHO.  The TDS indicates the salinity behaviors of water and based on TDS, natural water quality can be classified as fresh (if TDS <1,000 mg/L), brackish (if TDS = 1,000-10,000 mg/L), saline (if TDS = 10,000-1,000,000 mg/L) and brine (if TDS >1,000,000 mg/L).  According to this classification, the surface and groundwater of the studied area were in the fresh category. However, water with extremely low concentrations of TDS may also be unacceptable because of its flat, insipid taste.
The acceptable limit of TH (as CaCO 3 ) is 200 mg/L, which can be extended up to 600 mg/L in case of non-availability of any alternate water source.  The TH in the studied samples were ranged from 112 to 181 mg/L in surface water and from 222 to 400 mg/L in groundwater. Based on TH, Dufor and Becker  classified water as soft (if TH = 0-60 mg/L), moderately hard (if TH = 61-120 mg/L), hard (if TH = 121-180 mg/L), and very hard (if TH >180 mg/L). According to this classification, about 40% of surface water samples were belonged to moderately hard type, 50% were belonged to hard type and the rest (surface water sample collected from Jahanabad with TH = 181 mg/L) corresponded to very hard type water. Whereas, all the examined samples of groundwater were under the category of very hard water [Table 1].
The values of TAlk in the present study were measured from 398.2 to 539.7 mg/L in surface water and from 373.8 to 495.8 mg/L in groundwater [Table 1]. In natural water, most of the alkalinity is caused due to dissolution of CO2 in water. However, all the surface and groundwater samples contained TAlk more than TH, which is indicative of excess water alkalinity (EA).  It is also noted that the EA varied from 21 to 171 mg/L (not mentioned in the table), which could be due to the salt of NaHCO 3 .
The concentrations of chloride in the examined water samples were between 13.1 and 151.2 mg/L [Table 1]. Surface water showed higher chloride concentration (81.8-151.2 mg/L) than that of groundwater (13.1-57.4 mg/L). Chloride in drinking water originates from natural sources, sewage and industrial effluents, urban runoff and saline intrusion. No health-based guideline value has been proposed for chloride in drinking-water. However, chloride concentrations in excess of about 250 mg/L can give rise to detectable taste in water. 
The nitrate concentration in groundwater and surface water is normally low but can reach high levels because of agricultural runoff, refuse dump runoff, or contamination with human or animal wastes. The NO3− concentrations were varied from 2.10 to 5.20 mg/L (surface water) and from 0.08 to 2.80 mg/L (groundwater). In natural conditions, concentration of NO3− does not exceed 10 mg/L in the water.  All the analyzed samples contained nitrate concentrations within the tolerance limit of 45 mg/L. 
The manganese concentrations were varied from 0.90-2.86 to 0.72-3.66 μg/L in surface and groundwater accordingly. The average concentration (1.83 μg/L) of Mn in groundwater was slightly higher than the average Mn content (1.76 μg/L) of surface water. The maximum permissible limit of Mn concentration in drinking water is 500 μg/L set by WHO.  Potable water quality standard in terms of Mn is 100 μg/L set by Bangladesh Centre for Advanced Studies.  Mn concentrations of the present study were several folds lower than both WHO and BCAS suggested values.
The maximum permissible limits of Pb suggested by WHO  and USEPA  are 10 and 15 μg/L respectively. [Table 1] revealed that about 50% of surface water samples and almost all the groundwater samples contained higher Pb concentrations than the WHO recommended guideline value; clearly demonstrating anthropogenic impact.
The maximum permitted concentration of Cd in drinking water is 3 μg/L set by WHO.  Results showed that five surface water samples out of ten satisfied the WHO permissible level. On the other hand, all the groundwater samples were beyond the permissible limit.
In this study, arsenic-concentrations were varied from 1.72to 3.35 μg/L in surface water and from 2.08 to 3.16 μg/L in groundwater at different sampling locations [Table 1]. Arsenic enters into the water naturally from rocks and sediments by coupled biogeochemical and hydrologic processes, some of which are presently affected by human activity.  The WHO  and USEPA  recommended limit of arsenic in water is 10 μg/L. Considering WHO guideline, all the water samples in the present study area were free from arsenic-contamination.
Interrelations of water quality parameters
The interrelations among the chemical parameters were evaluated using Pearson correlation coefficient (r) model [Table 2]. Correlation coefficient is commonly used to measure and to establish the relationship between two variables. It is a simplified statistical tool to show the degree of dependency of one variable to the other.  The results showed that the parameters in the water samples were correlated to each other at P < 0.01 and P < 0.05 level. A significant positive correlations were found between pH and TAlk (r = 0.686), Cl− (r = 0.858), NO3− (r = 0.846) at P < 0.01 levels, whereas pH was negatively correlated with EC (r = −0.858), TDS (r = −0.858), TH (r = −0.857), Pb (r = −0.802) at P < 0.01 levels and with Cd (r = −0.484) at P < 0.05 levels. The negative correlation between heavy metals and pH suggested that release of heavy metals to water might be enhanced by acidic pH and oxidizing conditions. EC and TDS (r = 1.00, P < 0.05) showed highly positive correlation with TH (r = 0.785 and 0.784), Cl− (r = 0.774 and 0.774), NO3− (r = 0.764 and 0.764), and Pb (r = 0.828 and 0.828). The positive correlation indicated that the water was mainly controlled by TH (i.e., Ca 2+ and Mg 2+ ), Cl− , and NO3− ions, which depend upon the mineral dissolution, mineral solubility, ion exchange, evaporation, and anthropogenic activities. TAlk exhibited good positive correlation with Cl− (r = 0.730), NO3− (r = 0.602), and negative correlation with EC (r = −0.626), TDS (r = −0.626), TH (r = −0.587), Pb (r = 0.562) at P < 0.01 levels. Mn and As showed no significant correlation between them and with other chemical parameters. However, some correlations were also observed between Pb and other parameters.
Principal component analysis
Principal component analysis is a multivariate statistical technique used for deciphering patterns within large sets of data and to explain the variance of a large dataset of inter-correlated variables with a smaller set of independent variables (PCs) by orthogonal transformation.  PCA is a powerful tool for the characterization of anthropogenic and geogenic loads. In this study, each dataset was subjected to PCA using the correlation matrix in order to standardize each variable, meaning the analysis was not influenced by the differences in data magnitude and measurement scales. ,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, A varimax rotation was applied to aid interpretation of the results; this works by loading variables more strongly on fewer PCs. Significant PCs were then selected based on the Kaiser principal.  Under this principal, only component with eigenvalues ≥1 would be accepted as possible sources of variance in the data, with the highest priority ascribed to component that has the highest eigenvector sum. [Table 3] shows the results of the PC loadings with a varimax rotation, as well as the eigenvalues, % variance explained and % cumulative variance. The results indicated that two PCs were obtained with eigenvalues >1 summing almost 72.53% of the total variance in the water datasets [Table 3]. The scree plot [Figure 3] was used to identify the number of PCs to be retained in order to comprehend the underlying data structure. The scree plot indicated that first two components captured the most significant variation in the data. However, the elbow occurred at PC 2, therefore, up to PC 2 was accepted to account for the variation in the data. The first PC, accounting for 60.266% of the total variance was correlated with pH, EC, TDS, TH, TAlk, Cl− , NO3− , Pb, and Cd. PC 1 may be interpreted as representing influences of anthropogenic point and nonpoint sources such as land use change, municipal waste, agricultural runoff and atmospheric deposition. The second PC described 12.264% of the total variance, showed a high negative loading on Mn and moderate positive loading for As. Water samples in the study area were found to be below the WHO  recommended guideline values in terms of Mn and As. Thus, it could be stated that PC 2 represented natural soil leaching processes. The negative loading on Mn indicated that as the As concentration increases, Mn level decreases showing existence of the inverse relation with As in the water. However, [Figure 4] shows the score plot for the two PCs explaining 72.53% of the total variance. Two groups of water (groundwater and surface water), could be classified through the diagram of the scores for PC 1 versus PC 2, where groundwater samples showed more variations than the surface water samples.
|Figure 3: Scree plot showing point (principal component 2) at which the factor extraction is terminated|
Click here to view
|Figure 4: Principal component scores for the surface and ground water samples|
Click here to view
CA comprises a series of multivariate methods that are used to find true groups of data.  Hierarchical CA that has been used in this study is the most widely applied technique for analyzing water quality data. Hierarchical clustering joins the most similar observations and then successively the next most similar observations. The levels of similarity at which observations are merged are used to construct a dendrogram. In this study, hierarchical agglomerative CA was performed on the data set by means of City-block (Manhattan) distance as similarity measurement and Ward's method as the amalgamation rule. The Dendrogram of the clustering of 11 chemical parameters shows in [Figure 5]. It exhibited that the chemical parameters could be initially grouped into two main clusters, cluster A and cluster B. Nevertheless, in our opinion it would be possible to distinguish among three secondary clusters labeled as clusters I to III in [Figure 5]. Cluster I showed a close association between As and Mn. This finding corroborates the result of PCA and CA. Cluster II formed by Cd, TH, Pb, TDS and EC. TDS and EC exhibited very close association and completely accordance with the correlation analysis (r = 1.00, P < 0.05). Cluster III formed by TAlk, NO3− , Cl− and pH. Moreover, the CA was also used to identify the spatial similarity between 10 sampling sites based on the level of chemical concentration and three statistically significant clusters were obtained [Figure 6]. Cluster I consisted of Dhopaghat Hat to Dhurail; cluster II consisted of Bakshimail to Jahanabad; and cluster III from Mohanpur to Ghasigram. The clustering procedure revealed the groups of similar sites in a quite convincing way. These clusters included sampling sites with similar characteristics and natural background that were affected by similar types of sources. The CA carried out using our data indicates that this approach could be used to design future spatial sampling strategy in an optimal way in the Mohanpur area. For instance, the number of the sampling sites could be optimized in such a way that for rapid quality assessment studies only the representative sites from each cluster (not all monitoring sites) can be used. This will reduce the number of analysis and also the cost of the risk assessment procedure. However, two distinct clusters (cluster A and cluster B) emerged from the grouping of ground and surface water samples [Figure 7]. Cluster A consisted of groundwater samples whereas cluster B consisted of surface water samples, which is completely identical to the PC score plot for ground and surface water samples [Figure 4].
|Figure 5: Dendrogram of clustering of water quality parameters using City-block (Manhattan) distance and Ward's agglomerative method|
Click here to view
|Figure 6: Dendrogram of clustering of sampling sites using City-block (Manhattan) distance and Ward's agglomerative method|
Click here to view
|Figure 7: Dendrogram of clustering of samples using Cityblock (Manhattan) distance and Ward's agglomerative method; SW = Surface water and GW = Ground water|
Click here to view
| Conclusions|| |
In general, the surface and groundwater samples in Mohanpur area of Rajshahi district were alkaline in nature that could be due to the salt of NaHCO3. Groundwater samples were belonged to the category of very hard, whereas surface water samples were moderately hard to very hard. The pH values and the concentrations of EC, TDS, chloride, and nitrate in the examined water samples were less than the WHO recommended maximum values. The trace metals, Pb and Cd, in 50% surface water samples and all groundwater samples exceeded the recommended standards by WHO indicating the influence of anthropogenic activities. However, As and Mn concentrations were well within the recommended limits. Pearson correlation coefficients among the chemical parameters showed a number of strong associations. Multivariate statistical methods were applied to the water quality data set. The two PCs explaining 60.27 and 12.26% respectively of the total variance were obtained from PCA. The PC 1 demonstrated the anthropogenic impacts on water quality, whereas the PC 2 was related to natural activities. Hierarchical CA helped to group the 10 sampling sites into three clusters of similar water quality characteristics. Based on obtained information, it is possible to design a future optimal sampling strategy, which could reduce the number of sampling sites and associated cost. In addition, the dendrogram of 11 chemical parameters was plotted and grouped into three main clusters.
In conclusion, one can view that surface and groundwater contamination occurred throughout the area to some extent and is likely to continue in the future. We are still in a position to make choices on how best to use, manage, and protect the valuable resource, water. It is our obligation to try our utmost to check the misuses of water as every drop of it is highly precious. Prevention of water contamination is a more effective strategy than cure. Improvement of local sanitation system along with training and awareness programs (on hygiene, water handling and system maintenance) can help in improving surface and groundwater quality in the studied area.
| References|| |
Wang Y, Wang P, Bai Y, Tian Z, Li J, Shao X, et al
. Assessment of surface water quality via multivariate statistical techniques: A case study of the Shoghua River Harbin region, China. J Hydro Environ Res 2013;7:30-40.
Rahman MS, Saha N, Molla AH. Potential ecological risk assessment of heavy metal contamination in sediment and water body around Dhaka export processing zone, Bangladesh. Environ Earth Sci 2014;71:2293-308.
Arslan H. Application of multivariate statistical techniques in the assessment of groundwater quality in seawater intrusion area in Bafra Plain, Turkey. Environ Monit Assess 2013;185:2439-52.
Saha N, Zaman MR. Evaluation of possible health risks of heavy metals by consumption of foodstuffs available in the central market of Rajshahi City, Bangladesh. Environ Monit Assess 2013;185:3867-78.
RDA. Rajshahi Metropolitan Development Plan. Vol. 1. Structure Plan and Master Plan. Prepared by RDA under the Ministry of Housing and Public Works, Government of the Peoples's Republic of Bangladesh; 2004.
Kondakis XG, Makris N, Leotsinidis M, Prinou M, Papapetropoulos T. Possible health effects of high manganese concentration in drinking water. Arch Environ Health 1989;44:175-8.
Frisbie SH, Ortega R, Maynard DM, Sarkar B. The concentrations of arsenic and other toxic elements in Bangladesh's drinking water. Environ Health Perspect 2002;110:1147-53.
Rahman MS, Molla AH, Saha N, Rahman A. Study on heavy metals levels and its risk assessment in some edible fishes from Bangshi River, Savar, Dhaka, Bangladesh. Food Chem 2012;134:1847-54.
Järup L. Hazards of heavy metal contamination. Br Med Bull 2003;68:167-82.
Pan J, Plant JA, Voulvoulis N, Oates CJ, Ihlenfeld C. Cadmium levels in Europe: Implications for human health. Environ Geochem Health 2010;32:1-12.
Phuong NM, Kang Y, Sakurai K, Sugihara M, Kien CN, Bang ND, et al.
Arsenic contamination in groundwater and its possible sources in Hanam, Vietnam. Environ Monit Assess 2012;184:4501-15.
Nordstrom DK. Public health. Worldwide occurrences of arsenic in ground water. Science 2002;296:2143-5.
Chetia M, Chatterjee S, Banerjee S, Nath MJ, Singh L, Srivastava RB, et al.
Groundwater arsenic contamination in Brahmaputra river basin: A water quality assessment in Golaghat (Assam), India. Environ Monit Assess 2011;173:371-85.
Fendorf S, Michael HA, van Geen A. Spatial and temporal variations of groundwater arsenic in South and Southeast Asia. Science 2010;328:1123-7.
Morales MM, Martí P, Llopis A, Campos L, Sagrado S. An environmental study by factor analysis of surface seawaters in the Gulf of Valencia (Western Mediterranean). Anal Chim Acta 1999;394:109-17.
Reghunath R, Sreedhara Murthy TR, Raghavan BR. The utility of multivariate statistical techniques in hydrogeochemical studies: An example from Karnataka, India. Water Res 2002;36:2437-42.
Wunderlin DA, Del Pilar DM, Valeria AM, Fabiana PS, Cecilia HA, De Los Ángeles BM. Pattern recognition techniques for the evaluation of spatial and temporal variations in water quality. A Case Study: Suquía River basin (Córdoba-Argentina). Water Res 2001;35:2881-94.
Bagur MG, Morales S, López-Chicano M. Evaluation of the environmental contamination at an abandoned mining site using multivariate statistical techniques - the Rodalquilar (Southern Spain) mining district. Talanta 2009;80:377-84.
Chen K, Jiao JJ, Huang J, Huang R. Multivariate statistical evaluation of trace elements in groundwater in a coastal area in Shenzhen, China. Environ Pollut 2007;147:771-80.
APHA. Standard Methods for the Examination of Water and Wastewater. Washington DC, USA: American Public Health Association; 1992.
Massart DL, Kaufman L. The Interpretation of Analytical Chemical Data by the Use of Cluster Analysis. the University of Michigan: Wiley; 1983.
WHO. Guidelines for Drinking-Water Quality. Geneva: World Health Organization; 2004.
EQS. Environmental Quality Standards for Bangladesh. Department of Environment, Government of the People's Republic of Bangladesh; 1991.
Tiri A, Boudoukha A. Hydrochemical analysis and assessment of surface water quality in Koudiat Medouar reservoir, Algeria. Eur J Sci Res 2010;41:273-85.
WHO. Guidelines for Drinking-Water Quality. 2 nd
ed., Vol. 1. Geneva: World Health Organization; 1997.
Todd DK. Groundwater Hydrology. New York: Wiley; 1980.
Dufor CN, Becker E. Public Water Supplies of the 100 Largest Cities of the United States, 1962: U.S. Govt. Print. Off.; 1964.
Chow VT. Handbook of Applied Hydrology. London: McGraw-Hill; 1964.
Cushing EM, Kantrowirtz IH, Taylor KR. Water Resources of the Delmarva Peninsula: U.S. Govt.Print. Off.(Washington); 1973.
WHO. Guidelines for Drinking-Water Quality. 2 nd
ed., Vol. 1. Geneva: World Health Organization; 1993.
BCAS. Guide to the Environmental Conservation Act 1995 and rules 1997. Bangladesh Centre for Advanced Studies; 1998.
USEPA. Drinking Water Regulations and Health Advisories. Washington DC, USA: U.S. Environmental Protection Agency; 1996.
Belkhiri L, Boudoukha A, Mouni L. A multivariate statistical analysis of groundwater chemistry data. Int J Environ Res 2010;5:537-44.
Simeonov V, Stratis JA, Samara C, Zachariadis G, Voutsa D, Anthemidis A, et al.
Assessment of the surface water quality in Northern Greece. Water Res 2003;37:4119-24.
Delvalls TÁ, Forja JM, González-Mazo E, Gómez-Parra A, Blasco J. Determining contamination sources in marine sediments using multivariate analysis. TrAC Trends Anal Chem 1998;17:181-92.
Emmerson RH, O'Reilly-Wlese SB, Macleod CL, Lester JN. A multivariate assessment of metal distribution in inter-tidal sediments of the Blackwater Estuary, UK. Mar Pollut Bull 1997;34:960-68.
Webster R. Statistics to support soil research and their presentation. Eur J Soil Sci 2001;52:331-40.
Kuppusamy MR, Giridhar VV. Factor analysis of water quality characteristics including trace metal speciation in the coastal environmental system of Chennai Ennore. Environ Int 2006;32:174-9.
Danielsson A, Cato I, Carman R, Rahm L. Spatial clustering of metals in the sediments of the Skagerrak/Kattegat. Appl Geochem 1999;14:689-706.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7]
[Table 1], [Table 2], [Table 3]
|This article has been cited by|
||Is road-side fishpond water in Bangladesh safe for human use? An assessment using water quality indices
| ||M. Shahidul Islam, Kei Nakagawa, M. Abdullah-Al-Mamun, Md. Abu Bakar Siddique, Ronny Berndtsson |
| ||Environmental Challenges. 2021; : 100434 |
|[Pubmed] | [DOI]|
||Spatiotemporal mapping of groundwater salinity in Al-Batinah, Oman
| ||Ammar Abulibdeh,Talal Al-Awadhi,Noura Al Nasiri,Ali Al-Buloshi,Montasser Abdelghani |
| ||Groundwater for Sustainable Development. 2021; 12: 100551 |
|[Pubmed] | [DOI]|
||Heavy metal pollution assessment in the groundwater of the Meghna Ghat industrial area, Bangladesh, by using water pollution indices approach
| ||Mirza A. T. M. Tanvir Rahman,Moutushi Paul,Nikhil Bhoumik,Mahmud Hassan,Md. Khorshed Alam,Zakia Aktar |
| ||Applied Water Science. 2020; 10(8) |
|[Pubmed] | [DOI]|
||Evaluation of water quality near the Malanjhkhand copper mines, India, by use of multivariate analysis and a metal pollution index
| ||Bably Prasad,Abhay Kumar Soni,Anusha Vishwakarma,Ratnesh Trivedi,Krishna Kant Kumar Singh |
| ||Environmental Earth Sciences. 2020; 79(11) |
|[Pubmed] | [DOI]|
||Comprehensive radioecological monitoring of freshwater ecosystems in the vicinity of Rooppur NPP (Peopleæs Republic of Bangladesh)
| ||R. A. Mikailova,D. N. Kurbakov,E. V. Sidorova,I. V. Geshel,N. V. Andreeva,Yu. V. Sorokin,A. V. Panov |
| ||Marine Biological Journal. 2020; 5(3): 30 |
|[Pubmed] | [DOI]|
||Impact of Fly Ash Placement in an Abandoned Opencast Mine on Surface and Ground Water Quality: A Case Study
| ||Bably Prasad,Deblina Maiti,Krishna Kant Kumar Singh |
| ||Mine Water and the Environment. 2019; |
|[Pubmed] | [DOI]|
||Evaluation of seasonal variation of water quality using multivariate statistical analysis and irrigation parameter indices in Ajakanga area, Ibadan, Nigeria
| ||S. A. Ganiyu,B. S. Badmus,O. T. Olurin,Z. O. Ojekunle |
| ||Applied Water Science. 2018; 8(1) |
|[Pubmed] | [DOI]|
||Assessment of groundwater salinity using GIS and multivariate statistics in a coastal Mediterranean aquifer
| ||H. Ferchichi,M. F. Ben Hamouda,B. Farhat,A. Ben Mammou |
| ||International Journal of Environmental Science and Technology. 2018; |
|[Pubmed] | [DOI]|
||Assessing groundwater quality and its sustainability in Joypurhat district of Bangladesh using GIS and multivariate statistical approaches
| ||A. R. M. Towfiqul Islam,Shuanghe Shen,Muhammed Atikul Haque,Md. Bodrud-Doza,K. W. Maw,Md. Ahosan Habib |
| ||Environment, Development and Sustainability. 2017; |
|[Pubmed] | [DOI]|
||Spatial distribution and contamination assessment of six heavy metals in soils and their transfer into mature tobacco plants in Kushtia District, Bangladesh
| ||Narottam Saha,M. Safiur Rahman,Yeasmin Nahar Jolly,Atiqur Rahman,M. Abdus Sattar,M. Abdul Hai |
| ||Environmental Science and Pollution Research. 2016; 23(4): 3414 |
|[Pubmed] | [DOI]|
||Characterization of groundwater quality using water evaluation indices, multivariate statistics and geostatistics in central Bangladesh
| ||Md. Bodrud-Doza,A.R.M. Towfiqul Islam,Fahad Ahmed,Samiran Das,Narottam Saha,M. Safiur Rahman |
| ||Water Science. 2016; |
|[Pubmed] | [DOI]|