In this study, selected statistical methods (Pearson correlation, factor analysis and geostatistics) were used to determine correlation between elements and pollution sources in their spatial variability in different parts of Anzali wetland. During of study Major elements and heavy metals (Al, Fe, K, Mg, Mn, Ti, S, P, Na, Ca, Sr, Ba, As, Bi, Cd, Co, Cr, Cu, Mo, V, Ni, Pb, Zn) concentration in 18 points of wetland and rivers flowing into it were measured. Factor analysis and Pearson correlation results demonstrated that Ni has natural and anthropogenic source whereas As, Bi, Cd, Co, Cr, Cu, Mo, Pb, Zn and V have anthropogenic source. The ordinary kriging method applied for Spatial variability and interpolation of data in geostatistical study. Firstly normality of data was determined by using Kolmogorov - Smirnov test and after semi-variogram calculation, Circular, Spherical, Exponential and Gaussian models were fitted to each of heavy metals. RMS parameter of Cross-Validation was used for model validity. For As, Cu and Pb is exponential and for Cd, Zn, Ni circular are the best models. Results for As, Cu, Ni and Pb reveal the high estimation accuracy. Pollution zonation map that is produced by kriging and GIS showed Hendekhale and Shijan are the most polluting districts and Siakishom is the most none pollution part of Anzali wetland.