The Stokes-Einstein rotational diffusion relation tells us that we can write down a rotational diffusion coefficient for a sphere as:
$$D_r \approx \frac{k_B T}{\zeta_f} \approx \frac{k_B T}{(8 \pi \eta)(r)^3}$$
Where $k_B$ is Boltzmann's constant, $T$ is the temperature in Kelvin, $\zeta_f \approx (8 \pi \eta)(r)^3$ is the friction, $\eta$ is the viscosity of the medium (e.g. $\approx 1 \space cP$ in pure water), and $r$ is the radius of the sphere.
For an example calculation of $D_r \approx 0.2 \space rad^2/s$ (or $Hz$) of a $\approx 1 \space \mu m$ radius sphere in pure water at room temperature, please see this site (this shortened URL is going to WolframAlpha).
However, as we can directly see in the calculation, the units for $D_r$ are clearly expressible as rad/s or Hz. Why are the units for $D_r$ always reported as rad$^2$/s in the literature?