Let Ω be a smooth compact oriented 3-dimensional Riemannian manifold with boundary. A quaternion field is a pair q = {α, u} of a function α and a vector field u on Ω. A field q is harmonic if α, u are continuous in Ω and ∇α = rot u, div u = 0 holds into Ω. The space 𝒞(Ω) of harmonic fields is a subspace of the Banach algebra 𝒬 (Ω) of continuous quaternion fields with the point-wise multiplication qq′ = {αα′ − u · u ′ , αu′ + α ′u + u ∧ u ′ }. We prove a Stone-Weierstrass type theorem: the subalgebra ∨𝒞(Ω) generated by harmonic fields is dense in 𝒬 (Ω). Some results on 2-jets of harmonic functions and the uniqueness sets of harmonic fields are provided. Comprehensive study of harmonic fields is motivated by possible applications to inverse problems of mathematical physics.