Surface soil heat flux(G0) is an indispensable component of the surface energy balance and plays an important role in the estimation of surface evapotranspiration(ET). This study calculated G0 in the Heihe River Basin based on the thermal diffusion equation, using the observed soil temperature and moisture profiles, with the aim to analyze the spatial-temporal variations of G0 over the heterogeneous area(with alpine grassland, farmland, and forest). The soil ice content was estimated by the difference in liquid soil water content before and after the melting of the frozen soil and its impact on the calculation of G0 was further analyzed. The results show that:(1) the diurnal variation of G0 is obvious under different underlying surfaces in the Heihe River Basin, and the time when the daily maximum value of G0 occurs is a few minutes to several hours earlier than that of the net radiation flux, which is related to the soil texture, soil moisture, soil thermal properties, and the vegetation coverage;(2) the net radiation flux varies with season and reaches the maximum in summer and the minimum in winter, whereas G0 reaches the maximum in spring rather than in summer, because more vegetation in summer hinders energy transfer into the soil;(3) the proportions of G0 to the net radiation flux are different with seasons and surface types, and the mean values in January are 25.6% at the Arou site, 22.9% at the Yingke site and 4.3% at the Guantan site, whereas the values in July are 2.3%, 1.6% and 0.3%, respectively; and(4) G0 increases when the soil ice content is included in thermal diffusion equation, which improves the surface energy balance closure by 4.3%. Surface soil heat flux (G0) is an indispensable component of the surface energy balance and plays an important role in the estimation of surface evapotranspiration (ET). This study calculated G0 in the Heihe River Basin based on the thermal diffusion equation, using the observed soil temperature and moisture profiles, with the aim to analyze the spatial-temporal variations of G0 over the heterogeneous area (with alpine grassland, farmland, and forest). The soil ice content was estimated by the difference in liquid soil water content before and after the melting of the frozen soil and its impact on the calculation of G0 was further analyzed. The results show that: (1) the diurnal variation of G0 is obvious under different underlying surfaces in the Heihe River Basin, and the time when the daily maximum value of G0 occurs is a few minutes to several hours earlier than that of the net radiation flux, which is related to the soil texture, soil moisture, soil thermal properties, and the vegetat ion coverage; (2) the net radiation flux varies with season and reaches the maximum in summer and the minimum in winter, while G0 reaches the maximum in spring rather than in summer, because more vegetation in summer hinders energy transfer into the soil; ( 3) the proportions of G0 to the net radiation flux are different with seasons and surface types, and the mean values ​​in January are 25.6% at the Arou site, 22.9% at the Yingke site and 4.3% at the Guantan site, while the values in July are 2.3%, 1.6% and 0.3% respectively; and (4) G0 increases when the soil ice content is included in the thermal diffusion equation, which improves the surface energy balance closure by 4.3%.