A generic approach: the first thing to check is whether your three vectors are linearly independent. They are, so they span a 3D space. Thus, a well-chosen fourth one could complement them into a four-dimensional basis. There would be a infinity of choices: any vector that is not in the 3D space will do the job.
Indeed, the first three vectors are orthogonal, with unit norm as $\sum_1^4 \left(\pm \frac{1}{2}\right)^2 = 1$. Hence, they are pairwise orthonormal. As above, they define a subspace of dimension $3$. Its orthogonal supplement is a 1D vector space, uniquely defined by one non-null vector, which can be scaled by any non-zero scalar.
So if you just want orthogonality, you have an infinity of choices. But in some cases, people uses orthogonal as a proxy for orthonormal. Indeed, your three vectors $v_0$, $v_1$, $v_2$ are of unit norm too.
So in this case, there are only two vectors $v_3$ with unit norm answering your question:
$$ \left[\begin{array}{arr}
\frac{1}{2}\\
-\frac{1}{2}\\
-\frac{1}{2}\\
\frac{1}{2}\end
{array} \right]$$
and its opposite:
$$ \left[\begin{array}{arr}
-\frac{1}{2}\\
\frac{1}{2}\\
\frac{1}{2}\\
-\frac{1}{2}\end
{array} \right]\,.$$
As you can see, $v_0$ exhibits no sign changes, $v_1$ has one, and three for $v_2$. And $v_3$ has two sign changes. You just have rediscovered, up to a factor, the $4$-dimensional Hadamard (or Walsh, or Paley) orthogonal basis:
$$H_4=\left[
\begin{array}{rrrr}
1 & 1 & 1 & 1\\
1 & -1 & -1 & 1\\
1 & 1 & -1 & -1\\
1 & -1 & -1 & 1\\
\end{array}\right]$$
You can find more information about them in: