Numerical solutions of the equations of fluid flow and heat transport are used to quantify the effects of groundwater flow on the subsurface thermal regime. Simulations are carried out for a vertical section through a basin with a distance of 40 km separating the regional topographic high and low. Emphasis is placed on understanding the conditions under which advective effects significantly perturb the thermal field. The transition from conduction-dominated to advection-dominated thermal regimes is sharp and depends primarily on the topographic configuration of the water table, the magnitude and spatial distribution of permeability, hydraulic anisotropy, and the depth of active flow. Deviations of surface heat flow from the background heat flux are a measurable effect of groundwater flow and depend on the same factors. Our results show that from 0% to almost 100% of the section may have surface heat flow significantly different from background heat flow, depending upon the nature of the hydrogeologic environment. A limited spatial variability in a distributed set of heat flow measurements and/or linear temperature-depth profiles does not ensure that surface heat flow measurements are not disturbed. The results of our simulations suggest that knowledge of the complete environment of a site, including the water table configuration and subsurface flow system, combined with more closely spaced heat flow measurements may be necessary to unravel the true background heat flux in active flow regions.