In this Wikipedia page it is said that the square roots of -1 in the quaternion ring are the elements of the imaginary sphere. I don't understand why this is so. I don't understand the system that's written there. If I do
$$(a+bi+cj+dk)^2=a^2-b^2-c^2-d^2+2(ab+cd)i+2(ac-bd)j+2(ad+bc)k=-1$$
the system we get is
$$\begin{cases} a^2-b^2-c^2-d^2=-1 \\ ab+cd=0 \\ ac-bd=0 \\ ad+bc=0 \end{cases}$$
I suppose this could be solved to get the same solution, but it is quite tiresome and I'd like to understand how to get to the prettier system the Wikipedia gives.