Several properties of balayage of measures in harmonic spaces are studied. In particular, characterisations of thinness of subsets are given. For the heat equation the following result is obtained: suppose that $E={\mathbf{R}}^{m+1}$ is given the presheaf of solutions of$$\sum _{i=1}^m\frac{\partial u}{\partial x_i}=\frac{\partial u}{\partial x_{m+1}}$$and $B$ is a subset of ${\mathbf{R}}^m\times [-\infty ,0]$ satisfying$$\left\{\right(\alpha x,{\alpha }^{2}t):(x,t)\in B,x\in {\mathbf{R}}^m,t\in \mathbf{R}\}\subset B$$for $\alpha \>0$ arbitrarily small. Then $B$ is thin at 0 if and only if $B$ is polar. Similar result for the Laplace equation. At last the reduced of measures is defined and several approximation theorems on reducing and balayage...