Concerning the rates and changes in the properties of elements and physico-chemical characteristics of groundwater and describing prepared maps, they have played a major role in the quality management of groundwater resources, therefore based on this issue, appropriate measurement is needed to present for preventing excessive pollution of the aquifer area. There are various methods existing for monitoring and predicting the density of the physico-chemical characteristics of groundwater the selection of appropriate methods depends on the objective circumstances of the region and its data. In this research which was carried out in Plour plain located on part of Haraz–Polour river basin, sampling from qualitative chemical parameters was done and then nutrient concentrations such as nitrate, phosphate and sulfate in 37 different wells in the area were measured. Then using spatial interpolation methods such as ordinary Kriging (OK) geo-statistical methods and other certain techniques such as Inverse Distance Weight (IDW), Radial Basis Function (RBF), Local Polynomial Interpolation (LPI) and Global Polynomial Interpolation (GPI), we assessed the   concentration changes in spatial zones for elements of phosphate, nitrate and sulfate.