It is shown that, for the heat equation on \({\mathbb R^d\times \mathbb R}\), d ≥ 1, any convex combination of harmonic (= caloric) measures \({\mu_x^{U_1},\dots,\mu_x^{U_k}}\) , where U1, . . . , Uk are relatively compact open neighborhoods of a given point x, can be approximated by a sequence \({(\mu_x^{W_n})_{n\in\mathbb{N}}}\) of harmonic measures such that each Wn is an open neighborhood of x in \({U_1\cup\dots\cup U_k}\) . Moreover, it is proven that, for every open set U in \({\mathbb R^{d+1}}\) containing x, the extremal representing measures for x with respect to the convex cone of potentials on U (these measures are obtained by balayage, with respect to U, of the Dirac measure at x on Borel subsets of U) are dense in the compact convex set of all representing measures. Since essential ingredients for a proof of corresponding results in the classical case (or more general elliptic situations; see Hansen and Netuka in Adv. Math. 218(4):1181–1223, 2008) are not available for the heat equation, an approach heavily relying on the transit character of the hyperplanes \({\mathbb R^d \times \{c\}}\) , \({c\in\mathbb R}\) , is developed. In fact, the new method is suitable to obtain convexity results for limits of harmonic measures and the density of extremal representing measures on \({X:=X'\times \mathbb R}\) for practically every space–time structure which is given by a sub-Markov semigroup (Pt)t>0 on a space X′ such that there are strictly positive continuous densities \({(t,x,y) \mapsto p_t(x,y)}\) with respect to a (non-atomic) measure on X′. In particular, this includes many diffusions and corresponding symmetric processes given by heat kernels on manifolds and fractals. Moreover, the results may be applied to restrictions of the space–time structure on arbitrary open subsets.