Proper $r$-harmonic functions on the Thurston geometries

For any positive natural number $r\in\N^+$ we construct new explicit proper $r$-harmonic functions on the celebrated $3$-dimensional Thurston geometries $\Sol$, $\Nil$, $\SL2$, \H^2\times\rn and $\s^2\times\rn$.