In this study, the spatial variations of soil water and heat under bare land (BL), natural snow (NS), compacted snow (CS) and thick snow (TS) treatments were analyzed. The relationship curve between soil temperature and water content conforms to the exponential filtering model, by means of the functional form of the model, it was defined as soil water and heat relation function model. On this basis, soil water and heat function models of 10, 20, 40, 60, 100, and 140 cm were established. Finally, a spatial variation law of the relationship effect was described based on analysising of the differences between the predicted and measured results. During freezing period, the effects of external factors on soil were hindered by snow cover. As the snow increased, the accuracy of the function model gradually improved. During melting period, infiltration by snowmelt affected the relationship between the soil temperature and moisture. With the increasing of snow, the accuracy of the function models gradually decreased. The relationship effects of soil water and heat increased with increasing depth within the frozen zone. In contrast, below the frozen layer, the relationship of soil water and heat was weaker, and the function models were less accurate.