Stabiliser codes map to linear subspaces of $\mathbb F_2^{2n}$ which are isotropic w.r.t. the standard symplectic form on this vector space. This is very useful for analysing these codes directly and explicitly computing plenty of properties.
However, it does not align very well with classical coding theory which deals with Euclidean spaces. Computer scientists had already become interested in codes over other finite fields, in particular $\mathbb F_4 = GF(4)$. Thus it seems like a good idea to map said isotropic subspaces over $\mathbb F_2$ to additive self-orthogonal codes over $\mathbb F_4$. If you consider the background of the authors, you might understand why this felt like a natural mapping to them.
And indeed, considering the papers of Calderbank, Rains, Shor, Sloane, Nebe, ... around and after 2000, they have thrown a lot of classical coding theory on the problem to establish some fundamental theorems on stabiliser codes (e.g. http://arxiv.org/abs/math/0001038). The papers, however, are quiet mathematical so a lot of implications might not be immediately clear.
Nevertheless, stabiliser codes are more than just classical additive self-orthogonal codes over $\mathbb F_4$ since this only reflects the structure of the stabiliser group up to phases. Often that's enough but not always.