Abstract: Objective In this article,we derive a spherical strain Kriging formula based on the basic theory of Kriging,and applied it to simulated and real GPS data.We analyzed its difference with the least-square collocationmethod.Crosscheck results indicate that Kriging interpolation is feasible and valid in GPS velocity smoothingand gridding.The Kriging strain results reveal low robustness and obvious edge effects,but the smoothed andgridded results for GPS velocity data during 1999-2004from Kriging interpolation methods are in agreementwith the results calculated by the Least-square collocation method.The strain rate results from the two meth-ods are similar in the whole distribution characteristic,however,the kriging results show low self-consistency.In a word,the Kriging strain method is not as good as the least-square collocation method for robustness andedge effect.