Solution of 18.07.16
Using OP idea.
Let
$$u=\dfrac{x}2,\quad v=\dfrac{y}2,$$
then translate the condition in form
$$4u^2+4v^2+z^2+4uvz=4,$$
$$(z+2uv)^2 = 4(1-u^2)(1-v^2),$$
$$z=2\sqrt{(1-u^2)(1-v^2)}-2uv \leq (1-u^2)+(1-v^2)-2uv,$$
$$z\leq2-(u+v)^2.$$
So we can prove inequality
$$\sqrt[2]{2(u+v)}+\sqrt[3]{2v+2-(u+v)^2}+\sqrt[4]{2u+2-(u+v)^2} < 4.$$
Using Cauchy-Schwarz inequality, we get
$$\left(\sqrt[2]{2(u+v)}+\sqrt[3]{2v+2-(u+v)^2}+\sqrt[4]{2u+2-(u+v)^2}\right)^2$$
$$\leq (\sqrt[2]4+\sqrt[3]4+\sqrt[4]4)\left(u+v+\left(v+1-\dfrac12(u+v)^2\right)^{2/3}+\left(u+1-\dfrac12(u+v)^2\right)^{1/2}\right) < 16.$$
Maximal value of
$$f(u,v) = u+v+\left(v+1-\dfrac12(u+v)^2\right)^{2/3}+\left(u+1-\dfrac12(u+v)^2\right)^{1/2}$$
can be achieved in local extremum within the area $u,v\in[0,1],$ or in the edges of the area.
The nesessary conditions of local extremum are
$$f'_u=0,\quad f'_v=0,$$
or
$$\begin{cases}
1+\dfrac23\dfrac{-u-v}{\sqrt[3]{v+1-\dfrac12(u+v)^2}}+\dfrac12\dfrac{1-u-v}{\sqrt{u+1-\dfrac12(u+v)^2}} = 0\\
1+\dfrac23\dfrac{1-u-v}{\sqrt[3]{v+1-\dfrac12(u+v)^2}}+\dfrac12\dfrac{-u-v}{\sqrt{u+1-\dfrac12(u+v)^2}} = 0,\\
\end{cases}$$
$$\begin{pmatrix}
-u-v & 1-u-v\\
1-u-v & -u-v
\end{pmatrix}
\begin{pmatrix}
\dfrac2{3\sqrt[3]{v+1-\dfrac12(u+v)^2}}\\
\dfrac1{2\sqrt{u+1-\dfrac12(u+v)^2}}
\end{pmatrix}
=\begin{pmatrix}
-1\\-1
\end{pmatrix}
$$
Solving this linear system by Cramer's rule, we obtain:
$$
\Delta =
\begin{vmatrix}
-u-v & 1-u-v\\
1-u-v & -u-v
\end{vmatrix}
= 2u+2v-1,
$$
$$
\Delta_1 =
\begin{vmatrix}
-1 & 1-u-v\\
-1 & -u-v
\end{vmatrix}
= 1,
$$
$$
\Delta_2 =
\begin{vmatrix}
-u-v & -1\\
1-u-v & -1
\end{vmatrix}
= 1,
$$
$$\begin{cases}
3\sqrt[3]{v+1-\dfrac12(u+v)^2}= 2(2u+2v-1)\\
2\sqrt{u+1-\dfrac12(u+v)^2}= 2u+2v-1
\end{cases}$$
If $$2(u+v)\geq1,$$ then
$$\begin{cases}
27\left(2v+2-(u+v)^2\right)=16(2u+2v-1)^3\\
2\left(2u+2-(u+v)^2\right)=(2u+2v-1)^2,
\end{cases}$$
$$\begin{cases}
108\left(u+v+2-(u+v)^2\right)=32(2u+2v-1)^3+27(2u+2v-1)^2\qquad[2\times(1)+27\times(2)]\\
2\left(2u+2-(u+v)^2\right)=(2u+2v-1)^2.
\end{cases}$$
$$\begin{cases}
256(u+v)^3-168(u+v)^2-24(u+v)-221 = 0\\
u = \dfrac14(6(u+v)^2-4(u+v)-3)
\end{cases}$$
Cubic equation
$$256s^3-168s^2-24s-221 = 0$$
has one real root
$$s = \dfrac1{32}\left(7+3^{5/3}\sqrt[3]{61-8\sqrt{58}}+3^{5/3}\sqrt[3]{61+8\sqrt{58}}\right),$$
so $f(u,v)$ has the unique extremum
$$(u,v,f) \approx(0.392428, 0.875106, 3.08237).$$
The edges of the field achieves when $u\in\{0,1\}$ or $v\in\{0,1\}.$
$$f(0,v) = v + (v+1-\dfrac12v^2)^(2/3) + (1-\dfrac12v^2)^(1/2)$$
has the global maximum $\approx3.03678,$ so
$$f(0,v) < 3.08237.$$
$$f(1,v) = v+1 + (v+1-\dfrac12(1+v)^2)^(2/3) + (2-\dfrac12(1+v)^2)^(1/2)$$
has the global maximum $\approx 2.9714,$ so
$$f(1,v) < 3.08237.$$
$$f(u,0) = u + (1-\dfrac12u^2)^(2/3) + (1+u-\dfrac12u^2)^(1/2)$$
has the global maximum $\approx 2.86142,$ so
$$f(u,0) < 3.08237.$$
$$f(u,1) = u+1 + (2-\dfrac12(u+1)^2)^(2/3) + (u+1-\dfrac12(u+1)^2)^(1/2)$$
has the global maximum $\approx 3.07604,$ so
$$f(u,1) < 3.08237.$$
Note that
$$f(0,0)=2,\quad f(0,1)\approx 3.017477,\quad f(1,0)\approx 2.854705, \quad f(1,1)=2.$$
Given the conditions and the accuracy of the calculations this means that
$$f(u,v) < 3.08238,$$
$$\left(\sqrt[2]{2(u+v)}+\sqrt[3]{2v+2-(u+v)^2}+\sqrt[4]{2u+2-(u+v)^2}\right)^2 < 3.08238\left(\sqrt[2]4+\sqrt[3]4+\sqrt[4]4\right) < 15.4169,$$
$$\sqrt[2]{2(u+v)}+\sqrt[3]{2v+2-(u+v)^2}+\sqrt[4]{2u+2-(u+v)^2} < 3.926,$$
and for $\quad x,y,z \geq 0,\quad x^2+y^2+z^2+xyz=4$
$$\boxed{\sqrt[2]{x+y}+\sqrt[3]{y+z}+\sqrt[4]{z+x}<4.}$$