The proof is complete.
Remark. A direct computation shows that if
$\alpha $
is a solution of 8 , then Gaussian curvature of the metric 7 takes the form 1 . Since it is supposed that
$K$
is signpreserving, the relation 1 allows to choose
$\alpha $
as a new parameter on
$u$
curves. With respect to the parameter
$\alpha $
we have
$$du=\frac{d\alpha}{K}=\frac{cos\alpha}{a+1cos\alpha}d\alpha $$
and the line element 7 takes the form
$$\begin{array}{c}d{s}^{2}={\left(\frac{cos\alpha}{a+1cos\alpha}\right)}^{2}d{\alpha}^{2}+{sin}^{2}\alpha d{v}^{2}.\end{array}$$ 
(9)

Remark. If
$\xi $
is a unit vector field on the Riemannian manifold
${M}^{n}$
, then the induced metric on
$\xi \left({M}^{n}\right)$
is
$d{\stackrel{~}{s}}^{2}={g}_{ik}d{u}^{i}d{u}^{k}+<{\nabla}_{i}\xi ,{\nabla}_{k}\xi >d{u}^{i}d{u}^{k}.$
If
$\xi $
is a totally geodesic vector field on
${M}^{2}$
, then the metric of
${M}^{2}$
has the standard form 7 and
${\nabla}_{{\partial}_{u}}\xi ={\nabla}_{{e}_{0}}\xi =0,{\nabla}_{{\partial}_{v}}\xi =sin\alpha {\nabla}_{{e}_{1}}\xi =sin\alpha \lambda \eta =2{sin}^{2}(\alpha /2)\eta .$
Thus, we have
$$d{\stackrel{~}{s}}^{2}=d{u}^{2}+{sin}^{2}\alpha d{v}^{2}+4{sin}^{4}(\alpha /2)d{v}^{2}=d{u}^{2}+4{sin}^{2}(\alpha /2)d{v}^{2}.$$
Taking into account 1 we can easily find the Gaussian curvature of the totally geodesic submanifold
$\xi \left({M}^{2}\right)$
, namely
$$\stackrel{~}{K}=\frac{1}{4}K(K2cot(\alpha /2\left){K}_{\alpha}^{\prime}\right),$$
where
$K\left(\alpha \right)$
is the Gaussian curvature of
${M}^{2}$
given by relations 1 and 8 .
The equations 8 and 1 completely determine the class of Riemannian 2dimensional manifolds admitting a totally geodesic unit vector field.
Proposition 1.1
Let
${M}^{2}$
be a Riemannian manifold with a line element of the form
$$d{s}^{2}=d{u}^{2}+{sin}^{2}\alpha \left(u\right)d{v}^{2}.$$
Denote by
$K$
the Gaussian curvature of
${M}^{2}$
. Then
$K=\frac{d\alpha}{du}$
if and only if the function
$\alpha \left(u\right)$
satisfies
$$\frac{d\alpha}{du}=1+\frac{m}{cos\alpha}(m=const).$$
Proof. The sufficient part is already proved. Suppose now that
$$\frac{d\alpha}{du}=K(\ne 0).$$
Then we have
$${\alpha}^{\prime}=K=\frac{{\partial}_{uu}(sin\alpha )}{sin\alpha}=({\alpha}^{\prime}{)}^{2}cot\alpha {\alpha}^{\prime \prime}.$$
Therefore,
${\alpha}^{\prime \prime}={\alpha}^{\prime}(1{\alpha}^{\prime})tan\alpha ,$
or
$\frac{{\alpha}^{\prime \prime}}{{\alpha}^{\prime}1}={\alpha}^{\prime}tan\alpha ,$
or
$$(ln{\alpha}^{\prime}1{)}^{\prime}=(lncos\alpha {)}^{\prime}.$$
Evidently, now
${\alpha}^{\prime}1=\frac{\leftm\right}{cos\alpha }$
where
$m=const$
is a constant of integration. Finally,
$\frac{d\alpha}{du}=1+\frac{m}{cos\alpha}.$
Corollary 1.1
Let
${M}^{2}$
be a Riemannian manifold of constant curvature
$c\ne 0$
. Then
${M}^{2}$
admits a totally geodesic unit vector field if and only if
$c=1$
. This vector field is parallel along meridians and moves along parallels with unit angle speed.
Proof. If
$K=c=const$
, then 1 can be satisfied if and only if
$c=1,a=1$
.
The equation 1 implies an elementary nonexistence result.
Corollary 1.2
Let
${M}^{2}$
be a Riemannian manifold with Gaussian curvature
$K$
. Then
${M}^{2}$
does not admit a totally geodesic unit vector field
$\xi $
with angle speed
$a$
if
$K1<a+1$
.
Proof. Indeed, one can easily see that
$cos\alpha =\frac{a+1}{1K}.$
If
$a+1>K1$
, then we come to a contradiction.
2 Integral tra jectories of the totally geodesic vector field
The integral trajectories of the totally geodesic vector field
$\xi $
can be found easily as follows. Let
$\gamma =\left\{u\right(s),v(s\left)\right\}$
be an integral trajectory. Since
$$\xi =cos\omega {e}_{0}+sin\omega {e}_{1}=cos\omega {\partial}_{u}+\frac{sin\omega}{sin\alpha}{\partial}_{v},$$
we can set
$$\frac{du}{ds}=cos\omega ,\frac{dv}{ds}=\frac{sin\omega}{sin\alpha}$$
and then
$$\frac{du}{dv}=cot\omega sin\alpha .$$
Since
$\alpha =\alpha \left(u\right)$
and
$\omega =av+{\omega}_{0}$
, we come to the equation with separable variables
$$\frac{du}{sin\alpha}=cot\omega dv.$$
Using 8 , we can find
$$\frac{du}{d\alpha}=\frac{cos\alpha}{a1+cos\alpha}$$
and make a parameter change in the left hand side of the equation above.
Then we come to the equation
$$\frac{cos\alpha d\alpha}{sin\alpha (a1+cos\alpha )}=cot\omega dv.$$
Taking primitives, we have
$$\begin{array}{cc}tan(\alpha /2)sin(av+{\omega}_{0})=c(a+(a+2\left){tan}^{2}\right(\alpha /2){)}^{\frac{a+1}{a+2}}& \text{for}a\ne 0,2,\end{array}$$  
$$\begin{array}{cc}tan(\alpha /2)sin(2v+{\omega}_{0})=c{e}^{\frac{1}{2}{tan}^{2}(\alpha /2)}& \text{for}a=2,\end{array}$$  
$$\begin{array}{cc}\frac{1}{2}tan{\omega}_{0}\left(\frac{1}{1cos\alpha}+lntan(\alpha /2\left)\right\right)=vc& \text{for}a=0.\end{array}$$  
Taking into account 3 , we remark that
$tan(\alpha /2)sin\omega =k$
and
${tan}^{2}(\alpha /2)={k}^{2}+{\mathbb{\varkappa}}^{2}$
. Therefore, we have an intrinsic equation on the integral curves of the totally geodesic vector field
$$\begin{array}{cc}k=c{\left[a+(a+2)({k}^{2}+{\mathbb{\varkappa}}^{2})\right]}^{\frac{a+1}{a+2}}& \text{for}a\ne 0,2,\end{array}$$  
$$\begin{array}{cc}k=c{e}^{\frac{1}{2}({k}^{2}+{\mathbb{\varkappa}}^{2})}& \text{for}a=2,\end{array}$$  
$$\begin{array}{cc}k=sin{\omega}_{0}exp\left[2cot{\omega}_{0}(vc)\frac{1}{2}\frac{1+{k}^{2}+{\mathbb{\varkappa}}^{2}}{{k}^{2}+{\mathbb{\varkappa}}^{2}}\right]& \text{for}a=0,\end{array}$$  
where
$c$
is a constant of integration.
Moreover, in any case
$$\xi \left(k\right)=cos\omega {\partial}_{u}[tan(\alpha /2)sin\omega ]+\frac{sin\omega}{sin\alpha}{\partial}_{v}[tan(\alpha /2)sin\omega ]=\frac{cos\omega sin\omega {\alpha}_{u}^{\prime}}{2{cos}^{2}(\alpha /2)}+\frac{asin\omega cos\omega tan(\alpha /2)}{sin\alpha}=\frac{cos\omega sin\omega}{2{cos}^{2}(\alpha /2)}({\alpha}_{u}^{\prime}+a).$$
The equation 8 yields
$$\xi \left(k\right)=\frac{(a+1)cos\omega sin\omega}{2{cos}^{2}(\alpha /2)}\left(1\frac{1}{cos\alpha}\right).$$
Thus, if
$a=1$
, then the integral trajectories of the field
$\xi $
form a family of circles. The metric of
${M}^{2}$
is
$$d{s}^{2}=d{u}^{2}+{sin}^{2}ud{v}^{2}$$
and we are dealing with the unit sphere parameterized by
$$r=\left\{sinucosv,sinusinv,cosu\right\}.$$
These circles satisfy
$$\begin{array}{c}tan(u/2)sinv=c.\end{array}$$ 
(10)

Let
$(\rho ,\phi )$
be polar coordinates in a Cartesian plane which passes through the center of the sphere such that
$(0,0,1)$
is the north pole on the sphere.
Then the parameters
$(\rho ,\phi )$
and
$(u,v)$
are connected via stereographic projection from the south pole as
$$\begin{array}{c}\rho =tan(u/2),\end{array}$$  
$$\begin{array}{c}\phi =v.\end{array}$$  
Therefore, the equation 10 defines a family of parallel straight lines on the Cartesian plane. The family of integral curves of a totally geodesic vector field on the unit sphere can be obtained as inverse images under stereogrphic projection of this family.
An explicit equation of this family is
$$r\left(v\right)=\left\{\frac{2csinvcosv}{{c}^{2}+{sin}^{2}v},\frac{2c{sin}^{2}v}{{c}^{2}+{sin}^{2}v},\frac{{c}^{2}{sin}^{2}v}{{c}^{2}+{sin}^{2}v}\right\}$$
where
$c$
is the geodesic curvature of the corresponding circle. All of these circles pass through the south pole
$(0,0,1)$
when
$v=0,\pi $
. We can find this by using the expression
$tan(u/2)=c/sinv$
and trigonometric expressions for
$sinu$
and
$cosu$
via
$tan(u/2)$
.
The unit sphere is not the unique surface that realizes the metric
9 . Let
$(x,y,z)$
be standard Cartesian coordinates in
${E}^{3}$
. We can find an isometric immersion of the metric 9 into
${E}^{3}$
in a class of a surfaces of revolution. To do this, set
$$\begin{array}{c}x\left(\alpha \right)=sin\alpha ,\end{array}$$  
$$\begin{array}{c}({x}_{\alpha}^{\prime}{)}^{2}+({z}_{\alpha}^{\prime}{)}^{2}={\left(\frac{cos\alpha}{a+1cos\alpha}\right)}^{2}\end{array}$$  
and we easily find
$$\begin{array}{c}x\left(\alpha \right)=sin\alpha ,\end{array}$$  
$$\begin{array}{c}z\left(\alpha \right)={\int}_{{\alpha}_{0}}^{\alpha}\frac{cost}{a+1cost}\sqrt{1(a+1cost{)}^{2}}dt,\end{array}$$  
where the interval of integration is limited by the restrictions
$$\{\begin{array}{c}1+a<cos\alpha <2+a,\\ 2<a<1,\end{array}\text{or}\{\begin{array}{c}a<cos\alpha <1+a,\\ 1<a<0.\end{array}$$
The restrictions mean that if
$a+1\ge 1$
, then the metric 9 does not admit an isometric immersion into
${E}^{3}$
in a class of surfaces of revolution. References

Aminov Yu. The geometry of vector fields. Gordon&Breach Publ., 2000.

Boeckx E., Vanhecke L. Harmonic and minimal radial vector fields. Acta Math. Hungar. 90 (2001), 317331.

Boeckx E., Vanhecke L. Harmonic and minimal vector fields on tangent and unit tangent bundles. Differential Geom. Appl. 13 (2000), 7793.

Boeckx E., Vanhecke L. Isoparametric functions and harmonic and minimal unit vector fields, Contemp. Math. 288 (2001), 20–31.

Borisenko A., Yampolsky A. Riemannian geometry of bundles. Uspehi Mat. Nauk, 26/6 (1991), 5195; Engl. transl.: Russian Math. Surveys, 46/6 (1991), 55–106.

Brito F., Chacon P., Naveira A. On the volume of vector fields on spaces of constant sectional curvature, Preprint, 2001.

GilMedrano O.,LlinaresFuster E. Minimal unit vector fields, Tôhoku Math. J. 54 (2002), 71 – 84.

GilMedrano O.,LlinaresFuster E. Second variation of volume and energy of vector fields. Stability of Hopf vector fields, Math. Ann. 320 (2001), 531545.

GilMedrano O. Relationship between volume and energy of unit vector fields, Diff. Geom. Appl. 15 (2001), 137152.

GilMedrano O., GonzálezDávila J.C., Vanhecke L. Harmonic and minimal invariant unit vector fields on homogeneous Riemannian manifolds, Houston J. Math. 27 (2001), 377409.

Gluck H., Ziller W. On the volume of a unit vector field on the threesphere. Comm. Math. Helv. 61 (1986), 177192.

GonzálezDávila J.C., Vanhecke L. Examples of minimal unit vector fields. Ann Global Anal. Geom. 18 (2000), 385404.

GonzálezDávila J.C., Vanhecke L. Minimal and harmonic characteristic vector fields on threedimensional contact mentic manifolds. J. Geom. 72 (2001), 6576.

Pedersen S. L. Volumes of vector fields on spheres, Trans. Amer. Math. Soc. 336 (1993), 6978.

Sasaki S. On the differential geometry of tangent bundles of Riemannian manifolds. Tôhoku Math. J. 10 (1958), 338354.

Tsukada K.,Vanhecke L. Minimality and harmonicity for Hopf vector fields, Illinois J. Geom. 45 (2001), 441 – 451.

Tsukada K.,Vanhecke L. Invariant minimal unit vector fields on Lie groups, Period. Math. Hungar. 40 (2000), 123133.

Tsukada K.,Vanhecke L. Minimal and harmonic vector fields on
${G}_{2}\left({C}^{m+2}\right)$
and its dual space. Monatsh. Math. 130 (2000), 143154.

Yampolsky A. On the mean curvature of a unit vector field, Math. Publ. Debrecen, 60/12 (2002), 131 – 155.

Yampolsky A. On the intrinsic geometry of a unit vector field, Comment. Math. Univ. Carol. 43/2 (2002), 299317.

Yampolsky A. A totally geodesic property of Hopf vector fields, Acta Math. Hungar. 101/12 (2003), 93112.

Yampolsky A. On extrinsic geometry of unit normal vector field of Riemannian hyperfoliation, Math. Publ. Debrecen 63/4 (2003), 555 – 567.
Department of Geometry, Faculty of Mechanics and Mathematics, Kharkiv National University, Svobody Sq. 4, 61077, Kharkiv, Ukraine. email: yamp@univer.kharkov.ua