We attempt to compute the Surface Average Heat Flow （SAHF） from long-term temperature observations of one hundred seventy-seven observational points at the depths of 0.8, 1.6, and 3.2 m, which were relatively evenly distributed in China's Mainland. We first employ Fourier transformation to remove the influence of atmospheric temperature variations from the observation series, which are classified into the type of the steady-state temperature monotonously increasing with depth （type I） and other three types. Then we compare our results obtained from the data of type I, of which the values are thought to equal to those of the mean borehole heat flow, with those obtained from traditional heat flow observations mainly distributed in North China Craton. In computations of the SAHF at the observation stations, we deduce the thermal diffusivity and volumetric specific heat of the soil by employing harmonic solutions of the heat conduction equation for the same moisture group as the first step, and then we determine the SAHF using Fourier＇s law. Our results indicate that the SAHF derived from shallow earth geothermal data can reflect the heat flow field to a large extent.
