We estimate the oscillation of harmonic reproducing kernels. As an application of this estimation we obtain a double integral characterization of harmonic Besov spaces.