Besov spaces of harmonic functions on the unit ball of R '' are defined by requiring Sufficiently high-order derivatives of functions lie in harmonic Bergman spaces. We compute the reproducing kernels of those Besov spaces that are Hilbert spaces. The kernels turn out to be weighted infinite sums of zonal harmonics and natural radial fractional derivatives of the Poisson kernel. To cite this article: S. Gergun et al., C. R. Acad. Sci. Paris, Ser. I 347 (2009). (C) 2009 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.