Abstract: This paper takes north of the Tianshan Mountains as the study area.Monthly maximum snow cover of the study area was obtained by using MOD10A2.To improve the uniformity of spatial distribution of sample points,more spatial sample observations were gotten by creating virtual meteorological stations in snow-free area.Then,spatial interpolations for maximum snow depth from December to February were made by Kriging and co-Kriging method,and interpolation accuracy was evaluated by cross-testing method.Judged by average standard error and root mean square prediction error,simple co-kriging considering the effects of altitude is the best snow depth interpolation method.With the virtual meteorological stations,the interpolation result clearly shows the spatial distribution of snow depth in study area.