It's actually a really good question nand. :) It turns out that this is just small misunderstanding of what is meant by "harmonics" in this context, and there is a very simple explanation.
You are thinking in terms of the harmonics of a time domain signal (and by this I really mean the harmonics you'd get by transforming a time domain signal like the supply voltage or current into the frequency domain). The harmonics that they are referring to in this context however are those of a "spatial domain" signal, specifically the spatial distribution of flux around the machine.
Imagine for example that you had a machine with one pole pair and a very high degree of saliency, such that the flux under the poles was a uniform maximum, and that between the poles was zero. If you looked at the flux as a function of angular position it would basically be a square wave, and so have a Fourier decomposition in the same way that a square wave voltage does.
Notice that in this case however, the third harmonic (for example) doesn't correspond to a something that is changing three times as fast in the time domain, it instead corresponds sinusoidal to a flux component that undergoes three completely cycles as we move one complete rotation around the motor! Do you see the difference.
Hopefully it's now becoming clear, this 3rd harmonic flux component is precisely that which we would get in a machine with 3 times as many pole pairs! Hence 1/3 the speed for the same applied voltage/frequency.