We study harmonic Besov spaces b p α on the unit ball of R n , where 0 < p < 1 and α ∈ R. We provide characterizations in terms of partial and radial derivatives and certain radial differential operators that are more compatible with reproducing kernels of harmonic Bergman-Besov spaces. We show that the dual of harmonic Besov space b p α is weighted Bloch space b ∞ β under certain volume integral pairing for 0 < p < 1 and α, β ∈ R. Our other results are about growth at the boundary and atomic decomposition.