$\mathbf {ijk}=-1$ is part of the famously stone-carved formula by Sir William Rowan Hamilton, allowing the multiplication of triplets. There is already an intuition question on the topic of quaternions, and a beautifully illustrated post on 3D stereographic projections of quaternion rotations in 4D in the 3Blue1Brown youtube channel.
The question is very specific about how to picture in a geometric or intuitive way the equality $\mathbf {ijk}=-1$ part in the formula
$$\mathbf {i}^2 = \mathbf {j}^2 = \mathbf {k}^2 =\mathbf {ijk}=-1 $$
Attempt: $\mathbf i$, $\mathbf j$ and $\mathbf k$ act on the left of a quaternion as a pure rotation; hence,
$$\mathbf {ijk}\overset{?}=\mathbf {(ij)k}$$
and $\mathbf {i}$ acts on $\mathbf {j},$ rotating $\mathbf {j}$ into $\mathbf {k}$ in keeping with 3Brown1Blue's Grant Sanderson's "right-hand rule", illustrated with the 3D diagram below, corresponding to the stereographic projection of a 4D hypersphere of quaternions with norm $1$ and $0$ real component, where the $-\infty,$ $-\mathbf i,$ $1,$ $\mathbf i,$ $+\infty$ yellow line is the projection of a circle in 4D running through $-1,$ and centered at $0,$ passing through $-\mathbf i$ and $+\mathbf i;$ while the red-blue $\mathbf j,$ $\mathbf k,$ $-\mathbf j,$ $-\mathbf k$ circle corresponds to the 3D projection of a sphere in the 4D hypersphere passing through $-1,$ and centered at $0,$ and reaching both $\pm \mathbf k$ and $\pm \mathbf j.$ Rotating in the direction of $\mathbf i$ (i.e. multiplying on the left by $\mathbf i$) would amount to sliding the yellow line in the direction of the thumb, while rotating the red-blue circle following the clenching of the rest of the fingers:
and resulting in
$$\begin{align}\mathbf {ijk}&\overset{?}=\mathbf {(ij)k}\\ &=\mathbf {kk}\\ &=\mathbf k^2 =-1 \end{align}$$






You are right this can be viewed as a composition of (4D) rotations and, in fact, this particular example provides an example of a very interesting phenomenon in 4D Euclidean geometry. Namely, if we start with an arbitrary point $q$ in $\mathbb{H}$,
$$q = a + bi + cj + dk$$
then
$$kq = ak + bki + ckj + dkk$$
(the scalars commute so this is okay) which becomes through the usual formulas,
$$\begin{align} kq &= ak + bj - ci - d\\ &= -d - ci + bj + ak\end{align}$$
This is a rotation in 4D space as can be seen by noting that its matrix form is
$$k \equiv \mathbf{K} := \begin{bmatrix} 0 && 0 && 0 && -1 \\ 0 && 0 && -1 && 0 \\ 0 && 1 && 0 && 0 \\ 1 && 0 && 0 && 0 \end{bmatrix}$$
which is easily seen to be orthogonal, i.e. $\mathbf{K}^T \mathbf{K} = \mathbf{I}$. Moreover, the expansion by cofactors and minors (which is very easy) shows $\mathrm{det}(\mathbf{K}) = 1$, so it is indeed a proper rotation (i.e. is not a rotoreflection).
Knowing that it's a rotation, the natural question is: what is it a rotation about? In 4D, we don't rotate around axes that are lines, but rather we rotate around planes, because a line gives "too much freedom" for a single angle of rotation to capture in the same way that we don't "rotate around a point" in 3D with a single angle of rotation. (This may be hard to visualized, but that's 4D for 'ya.)
So which plane is it? We may be tempted to note that while this permutes all 4 coordinates, you might see that it fixes the plane where $a = -d$ and $b = -c$, i.e. the quaternions of the form
$$a + bi - bj - ak$$
are unchanged under the transformation $\mathbf{K}$. Now you might think that a bit strange, given that is not a simple coordinate plane - and it is, indeed, an omen. In fact, if you see something like this happening in maths, chances are something very strange is going on - it's no guarantee, but your spidey sense should tingle at it and you should feel disturbed. In particular, we can decompose this a bit further. Taking another look we can see that it is actually the composition of two rotations: one of which does
$$(a + bi + cj + dk) \mapsto (-d + bi + cj + ak)$$
which fixes the $ij$-plane, and is a 90-degree, followed by
$$(a + bi + cj + dk) \mapsto (a - ci + bj + dk)$$
which fixes the $1k$-plane. This is also a 90-degree rotation. These two planes are orthogonal, and when you have a rotation in 4D space that is formed by the composition of two rotations of the same angle through two orthogonal planes, you get what is called an isoclinic or Clifford rotation that in fact fixes infinitely many planes, one other of which we just found. This is a very strange phenomenon with no simple analogue in 3D: in 3D a rotation will fix only one line, its axis, unless it's a degenerate rotation by zero degrees.
Likewise, the rotations by action of $j$ and $i$ are similar Clifford rotations, and they all together act to bring $1$ to $-1$ via $-1 = i(j(k(1))) = ijk$. If you want, you can think of this approximately as being like a hyperdimensional version of a Rubik's cube where you twist the sides through several different positions to move one colored patch from one side to another (not exactly, though, since here the whole "cube" is rotating, but an approximate analogy - the gist is the chained transformations.).