Abstract: The present research was conducted to analyze spatial changes in soil salinity distribution as an aspect of soil degradation and to compare the efficacy of different Geostatistical methods in its estimation and the preparation of maps of the spatial distribution of soil salinity. To estimate soil salinity of non-sampled areas, the methods of Kriging, Co-Kriging and Weighted Moving Average were applied in Geographical Information System (GIS) medium. To evaluate the efficacy of the methods, the cross-evaluation approach with two statistical parameters of mean bias error and mean absolute error was taken in practice. Results indicated the high precise of Kriging method with regression coefficient of 0.98 for the estimation of salinity rates in the areas, for where no data were available before. Estimation error for this method was 1.31 and biass was -0.34 dS m-1 which indicates high accuracy of Kriging method to estimate topsoil salinity and its precise.