# Introduction n the current decades the heavy metal accumulation in the soils is a growing concern due to its potential health risks as well as its detrimental effects on soil ecosystems (McLaughlin et al., 1999;Qishlaqi & Moore, 2007). Heavy metals have characteristics including that they are non-biodegradable (Facchinelli et al. 2001) and they can be necessary or beneficial to plants at certain levels, but can be toxic when exceeding specific thresholds (Qishlaqi & Moore, 2007;Bilos et al., 2001) Sources of these elements in soils mainly include natural occurrence derived from parent materials and anthropogenic activities. Anthropogenic inputs are associated with industrialization and agricultural activates, deposition, such as atmospheric deposition, waste disposal, waste incineration, emissions from traffic, fertilizer application and long-term application of wastewater in agricultural land (Qishlaqi & Moore, 2007;Bilos et al., 2001;McLaughlin et al., 2001;Koch et al., 2001). For the present day's environmental researchers, knowledge of the heavy metal accumulation in soil, the potential source of heavy metals and their possible interactions with soil are one of the prime focuses. Different statistical analysis tools can provide such knowledge and can be very helpful for the interpretation of environmental data (Tuncer et al., 1993;Sena et al., 2002;Einax et. al 1999). In recent times, the statistical methods (univariate or multivariate) have been applied widely to investigate heavy metal concentration, accumulation and distribution in soils. (Vega et al. 1998;Wunderlin et al., 2001;Grande et al., 2003;Simeonov et al., 2003;Pekey et al., 2004;Singh et al., 2004;Astel et al., 2006;Kowalkowski et al., 2006;Shrestha & Kazama, 2007Salman et al., 1999). Once, Dhaka City, the capital of Bangladesh had excellent natural drainage system even 40 years ago. The city was interlaced with numerous natural channels/canals and wetlands. It is estimated that there were up to 45 natural drainage canals that span the city. In the course of rapid expansion of the city, most of the natural drainage canals as well as wetlands has been intervened and destroyed. Now only few canals exist but these have become contaminated wetland because of disposal of solid waste, toxic industrial waste which are the potential sources of heavy metal pollution (Subramanian, 2004;Karn & Harada, 2001). Now even after a medium size rainfall, the streets of the city get flooded for hours at a time because water has no way to drain out easily. Although some drainage structures have been built over the last two decades, they are woefully inadequate. In addition, due to haphazard design and construction of these drainage structures and lack of proper maintenance, over the years these have lost their carrying capacity due to severe clogging. The effects of water logging causes serious suffering of the city dwellers as well as damage the roads and thereby increasing the road maintenance cost. On the other hand, water supply and sanitation infrastructures have been being become ineffective due to unwanted water logging in the city. There is thus an immediate need for rehabilitation and development of the natural I drainage network and find ways to properly operate and maintain the already constructed drainage structures, such as box culverts and drain lines. For the restoration of these canals in Dhaka city, it is very important to explore the current pollution status of the sludge deposited at the bottom of canals over the year. The present study stems from the above concerns, with the primary focus on the current status of distribution of some selected toxic metals in the sludge sample collected from 10 major canals in the Dhaka city Corporation (DCC) area by multivariate statistical analysis. # II. Major Canals in Dcc Area and Their Present Conditions The major canals (Khals) are located in the various part of the city and their other particulars are shown in the Table 1. From our field survey, it was found that illegal encroachment and nearby pollution activities are the major concerns to maintain the natural conditions in most of the canals in Dhaka city. A variety of small industries are building up near the bank of the canals. Various Residential plots are now under construction on the different parts of the canals. The water of the canals water gets polluted and become blackish with lots of waste including construction debris, vegetations, chemical waste, polythene sheets and as a result there is partial drainage blocking. Some photographs, which were taken during our field survey, are presented in Fig. 1 # Materials and Methods For the present study, sediment (sludge) samples were collected from 10 canals (name of the canal mentioned in Table 2) on the months of March and April, 2011 and samples were analyzed for the metals, Cd, Cr, Cu, Mn, Fe and Pb. Above mentioned metals concentration in the collected sludge samples were determined by total extraction with Aqua-Regia. Shimadzu, Japan, AA6800). Standard QA/QC protocol was followed throughout, including replicate analysis (1 in every 5 samples), checking of method blanks (1 in every 10 analysis) and standards (1 in every 10 analysis). The estimated metal levels were compared with the permissible safe levels for the sediment sample proposed by USEPA. Multivariate statistical techniques (Lin et 1976) were adopted to assess the metal contamination in the sludge. For this purpose, the well founded techniques of Pearson correlation analysis, Principal Component Analysis (PCA), and Cluster Analysis (CA) were jointly used, the first affording a direct measure of interdependence of the set of variables under investigation while the latter two provides the visual grouping of the data to help understand the interrelated metal clusters produced (Hopke 1992). SPSS software (Version 16.0) was used to perform the multivariate statistical analyses. IV. # Results and Discussion The Sludge samples analysis results for different heavy metals are presented in Table 2. The comparison between the metals concentration present in the sludge samples of DCC canals with permissible metal concentration limit proposed by USEPA for the sediment sample shows that some canals in DCC area are facing heavy metal pollution. # Note: BDL-Below Detection Limit According to USEPA guideline, the sludge samples of Mohakhali Khal, Mirpur Housing Khal, Segunbagicha Khal, Jirani Khal are heavily polluted with Cu. Baunia Khal, Kalyanpur Shakha 'Gha" (Shewrapara), Segunbagicha Khal, Jirani Khal are heavily polluted with Cr. Pollution level of Pb in Mirpur Housing Khal is also exceed the USEPA heavily polluted criteria. Pollution level of Hazaribagh Khal, Kalyanpur "Kha" Khal , Kalyanpur main Khal ,Digun Khal are comparatively low for all tested heavy metals except Cr. By reviewing the pollutant limits of USEPA against the result from sludge samples analysis, it can be easily stated the pollutants concentration in the sludge of the selected canals are below the prescribed hazard limit for land application but some canals exceed the EPA guideline for heavily polluted sediments for some metal. Before forming a judgment on the observed distribution of metal levels and interrelationship among them, the metal data was first examined on the basis of linear correlation between metal pairs in terms of significant positive correlation coefficient. Strong positive correlations were observed for Mn -Fe (r = 0.860), Pb -Cu (r = 0.786), Pb -Cd (r = 0.398) and Cu -Cd (r = 0.227) pairs (Correlation matrix is shown in Table 3), indicating the existence of a common source/origin of these metals in the sludge sample. Further confirmation of this hypothesis of 'different heavy metals may have common origin' was secured through multivariate methods of statistical analysis (Hair et al. 1988). In this study, two multivariate techniques were applied: Principal Component Analysis (PCA) and Cluster Analysis (CA). The PCA has emerged as a useful tool for better understanding the relationships among the variables (e.g., metal concentrations in this study) and for revealing groups (or clusters) that are mutually correlated within a data body (Qishlaqi & Moore, 2007). This procedure reduces overall dimensionality of the linearly correlated data by using a smaller number of new independent variables, called principal components (PC), each of which is a linear combination of originally correlated variables. On the other hand, Cluster Analysis (CA) exclusively classifies a set of observations into two or more unknown groups based on combination of internal variables. Therefore, the purpose of CA is to discover a system of organized observations where a number of groups/variables share properties in common, and it is cognitively easier to predict mutual properties based on an overall group membership (Everitt 1993;Jolliffe 1986). This helps define source profiles of variables, such as metal concentrations, and their interpretation in terms of possible sources (Jobson 1991). Principal Component Analysis (PCA) using Varimax normalized rotation was conducted for common source identification. The variables are correlated with two principal components in which 70.3% of the total variance in the data was found. The rotated Principal Component Loadings are given in Table 4. Principle component plot in a rotated space is shown in Fig. 2. The first component with 40.92% of variance comprises Pb, Cu (bold figures in Table 4) with high loadings. This association strongly suggests that these variables have a strong interrelationship. The second component (PC2) contributes Mn and Fe at 29.40 % variance, which also infers the strong correlation between this metal pair. The corresponding cluster analysis Dendrogram is shown in Fig. 3. From the cluster analysis result it can be said that there is a strong correlation between Fe-Mn metal pair, which is a good agreement with PC2, but cluster analysis results did not show a good agreement between Cu-Pb pair. This result suggests that the strong relationship between Cu-Pb pair does not confirm. 1![Fig. 1 : Photographs of some canals in and around the Dhaka city depicting the current pollution scenario](image-2.png "Fig. 1 :") 23![Fig. 2 : Principal component plot in a rotated space](image-3.png "Fig. 2 :Fig. 3 :") 1SerialCanal's NameLocationLengthNo.(Km)1Kalyanpur main khalWestern part of the city32Kalyanpur branch khal -KaWestern part of the city1.53Kalyanpur branch khal -KhaWestern part of the city2.44Kalyanpur branch khal -GhaWestern part of the city1.565Kalyanpur branch Khal -UmoWestern part of the city1.786Kalyanpur branch Khal -ChaWestern part of the city0.987Baunia KhalNorth western part8.88Digun KhalNorth eastern part4.59Mohakhali KhalCentral city2.310Hazaribagh KhalSouth-western part.0.711Shegunbagicha khalCentral-eastern part1.012Manda khalCentral-eastern part1.013Shangbadik ColonyNorth western part1.014Section 2 to Digun Canal through SectionNorth western part3.56 and Rupnagar15Mirpur Housing CanalNorth western part1.016Kashaibari -Boalia to Balu riverNorth eastern part3.017Gerani khalCentral-eastern part5.0 2Dhaka city area 3CdCrCuMnFePbCd1.000Cr-0.2921.000Cu0.227-0.2001.000Mn-0.5470.064-0.1091.000Fe-0.207-0.206-0.0480.8601.000Pb0.398-0.2820.768-0.2730.0311.000 4ComponentPC1PC2Cd0.511-0.525Cr-0.585-0.113Cu0.812-0.066Mn-0.1400.968Fe0.1840.927Pb0.890-0.140 ## Conclusions The present study showed that Cr, Cu and Pb were present as major pollutants in the some canals in the DCC area with high concentration levels, while Cd, Mn and Fe emerged as minor pollutants. Strong positive linear correlations were found between Mn, Cu and Pb from linear regression analysis. Principal component analysis summarizes (reduces) the data set into two major components representing the different interrelationship among the elements. Strong interrelation between Cu -Pb pair and Mn -Fe pair was found from principle component analysis. Corresponding cluster analysis result confirms the relationship between Mn -Fe metal pair but, does not confirm the strong interrelationship between the Cu -Pb metal pair. Comparison with USEPA guideline for sediment showed that Cr, Pb and Cu levels are in far excess of the recommended safe limits for some canals but pollutants concentration in the sludge are below the prescribed hazard limit for land application of sludge. * Chemometrics in monitoring spatial and temporal variations in drinking water quality AAstel MBiziuk APrzyjazny JNamiesnik Water Research 8 2006 * Sources, distribution and variability of air borne trace metals in La Plata city are CBilos JCColombo CNSkorupka PM JRodriguez Argentina. Environ. Pollu 11 2001 * Geostatistical and multivariate statistical method for the assessment of polluted soils; Merits and limitations JWEinax USoldt Chemometrics Intell. Lab 46 1999 * Cluster analysis BSEveritt 1993 Heineman London * Multivariate statistical and GIS-based approach to identify heavy metals sources in soils AFacchinelli ESacchi LMallen Environmental Pollution 114 2001 * A description of how metal pollution occurs in the Tinto-Odiel rias (Huelva-Spain) through the application of cluster analysis JAGrande JBorrego JAMorales MLTorre Marine Pollution Bulletin 46 2003 * Factor and correlation analysis of multivariate environmental data PKHopke Methods of environmental data analysis CNHewitt London, UK Elsevier 1992 * JFHair Jr REAnderson RLTatham WCBack Multivariate data analysis Englewood Cliffs, NJ Prentice-Hall 1988 5th ed. * Principal component analysis ITJolliffe 1986 Springer New York * Applied multivariate data analysis JDJobson 1991 Springer NewYork * Surface water pollution in three urban territories of Nepal, India, and Bangladesh SKKarn HHarada Environ Manage 28 4 2001 * On the contribution of background sources to the heavy metal content of municipal sewage sludge MKoch WRotard Water Sci. Technol 43 2001 * Application chemometrics in river water classification TKowalkowski RZbytniewski JSzpejna BBuszewski Water Research 40 2006 * Multivariate analysis of soil heavy metal pollution and landscape pattern in Changhua county in Taiwan YPLin TPTeng TKChang Landscape and Urban Planning 62 2002 * Metals and micronutrients-food safety issues MJMclaughlin DRParker JMClarke Field Crops Res 60 1999 * A bioavailability-based rationale for controlling metal and metalloid contamination of agricultural land in Australian and New Zealand MJMclaughlin REHammon RMclaren GSpeir TWRogers Aust. J. Soil Res 38 2000 * Pollution source analysis of river water and sewagesludge JNolte Environmental Technology Letters 9 1988 * Source apportionment of trace metals in surface waters of a polluted stream using multivariate statistical analyses HPekey DKarakas LMBakog Marine Pollution Bulletin 49 2004 * Statistical Analysis of Accumulation and Sources of Heavy Metals Occurrence in Agricultural Soils of Khoshk River Banks AQishlaqi FFarid Moore American-Eurasian J. Agric. & Environ. Sci 2 5 2007 * SRSalman YHAbu Rukah 1999 * Multivariate and principal component statistical analysis of contamination in urban and agricultural soils from north Jordan Environmental Geology 38 * Trace inorganics in rural potable water and their correlation to possible sources SSSandhu WJWarren Nelson P Water Research 12 1976 * Discrimination of management effects on soil parameters by using principal component analysis: A multivariate analysis case study MMSena RT SFrighetto PJValarini HTokeshi RJPoppi Soil and Tillage Research 67 2002 * Assessment of the surface water quality in Northern Greece VSimeonov JAStratis CSamara GZachariadis DVoutsa AAnthemidis MSofoniou THKouimtzis Water Research 37 2003 * Multivariate statistical techniques for the evaluation of spatial and temporal variations in water quality of Gomti River (India), a case study KPSingh AMalik DMohan SSinha Water Research 38 2004 * Assessment of surface water quality using multivariate statistical techniques: A case study of the Fuji river basin SShrestha FKazama Japan. Environmental Modelling and Software 22 2007 * Water quality in South Asia BSubramanian Asian J Water Environ Pollut 1 1-2 2004. 2007 Statistical year book (SYB) * Multivariate Analysis of heavy metal contents in soils, sediments and water in the region of Meknes (Central Morocco) MTahri MBounakhla EBilal JJGruffat JMoutte DGarcia Environmental Monitoring and Assessment 102 2005 * Metal pollution in the Golden Horn, Turkey-Contribution of natural and anthropogenic components since 1913 GTTuncer SGTuncel TIBalkas Water Science and Technology 28 1993 * Assessment of seasonal and polluting effects on the quality of river water by exploratory data analysis MVega RPardo EBarraado LDeban Water Research 32 1998 * Pattern recognition techniques for the evaluation of spatial and temporal variations in water quality. A case study: Suquia river basin DAWunderlin MPDiaz MVAme SFPesce ACHued MABistoni Water Research 35 2001 * Multivariate analysis of the Hong Kong Victoria harbour water quality data IM HYeung Environmental Monitoring and Assessment 59 1999