I am reading Do Carmo's Riemannian Geometry. I got stuck on Problem 14 of Chapter 3. Here is the original problem with the hints:
(Liouville's Theorem) Prove that if $G$ is the geodesic field on $TM$ then $\mathrm{div}G=0$. Conclude from this that the geodesic flow preserves the volume of $TM$.
Hint: Let $p\in M$ and consider a system $(u_1, \ldots, u_n)$ of normal coordinates at $p$. Such coordinates are defined in a normal neighborhood $U$ of $p$ by considering an orthonormal basis $\{e_i\}$ of $T_pM$ and taking $(u_1, \ldots, u_n), q = \mathrm{exp}_p(\sum_iu_ie_i), i=1, \ldots, n$ as coordinates of $q$. In such a coordinate system, $\Gamma_{ij}^k(p)=0$, since the geodesics that pass through $p$ are given by linear equations. Therefore if $X=\sum x_i\frac{\partial}{\partial u_i}$, then $\mathrm{div}X(p)=\sum\frac{\partial x_i}{\partial u_i}$.
Now let $(u_i)$ be normal coordinates in a neighborhood $U\subset M$ around $p\in M$ and let $(u_i, v_j), v=\sum_j v_j\frac{\partial}{\partial u_j}, i, j = 1, \ldots, n$ be coordinates on $TM$ at $(q, v)$. Calculate the volume element of the natural metric of $TM$ at $(q, v), q\in U, v\in T_qM$, and show that it is the volume element of the product metric on $U\times U$ at the point $(q, q)$. Since the divergence of $G$ only depends on the volume element, and $G$ is horizontal, we can calculate $\mathrm{div}G$ in the product metric. Observe that in the coordinates $(u_i, v_j)$ we have $$G(u_i)=v_i, G(v_j)=-\sum_{ik}\Gamma^j_{ik}v_iv_k, k = 1,\ldots, n$$. Since the Christoffel sysbols of the product metric on $U\times U$ vanish at $(p, p)$, we can obtain finally, at $p$, $$\mathrm{div}G=\sum_i\frac{\partial v_i}{\partial u_i}-\sum_j\frac{\partial}{\partial v_j}\Big(\sum_{ik}\Gamma^j_{ik}v_iv_k\Big)=0.$$
Questions:
Where I got stuck is the bold part: Calculate the volume element of the natural metric of $TM$ at $(q, v), q\in U, v\in T_qM$, and show that it is the volume element of the product metric on $U\times U$ at the point $(q, q)$. I do not understand what "natural metric" is. If the natural metric is defined by $$\langle V, W\rangle_{(q, v)}=\langle d\pi(V), d\pi(W)\rangle_p + \langle \frac{Dv}{dt}(0), \frac{Dw}{ds}(0)\rangle_p$$ as is in Problem 2, then it seems too messy for me to calculate and simplify the volume element. Can any one help me on this? Thanks.
Assume that $TM$ has Sasaki metric $G$. Hence note that $M$ is $n$-dimensional totally geodesic submanifold in $TM$. In further $T_xM$ for all $x\in M$ is $n$-dimensional totally geodesic flat submanifold.
And $$ \pi : (TM,G)\rightarrow (M,g)$$ is a Riemannian submersion.
(1) Fix $(x,v)\in TM$. Consider tangent space of $TM$ at $(x,v)$.
Here $n$-dimensional subspace is horizontal lift of $T_xM$ by $\pi$, which is isometric to $(T_xM,g(x))$. And $n$-dimensional subspace orthogonal to that wrt $G$ is a fiber $(\pi^{-1}(x)=T_xM,g(x))$. Hence volume form $V(x,v)$ at $(x,v)$ wrt $G$ is a product of volume form $V(x)$ at $x$ wrt $g$.
Assume that $x$ is coordinate for $M$ Then we have a coordinate vector field $E_i$ Hence $$(x,v)\mapsto (x,v_iE_i)$$ is coordinate for $TM$.
If $g_{ij}:=g(E_i,E_j)$, then $$ V(x)=\sqrt{{\rm det}\ g_{ij} }\ dx^1\cdots dx^n $$
$$ V(x,v)={\rm det}\ g_{ij}\ dx^1\cdots dx^n dv^1\cdots dv^n $$
(2) If $c(t)=x(t)$ is a geodesic in $M$ with $$c'=x_i'(t)E_i(x(t)) $$ then $$ x_i''(t) +x_i' x_j' \Gamma_{jk}^i=0$$
Hence trajectory $(c(t),x_i'(t) E_i )$ has a tangent $$(x_i'E_i,x_j''E_j)=(x_i'E_i,-x_j'x_k'\Gamma_{jk}^iE_i ) $$
So if $x$ is normal coordinate, \begin{align*}&\ \ di_{(x_i'E_i,x_j''E_j)} V(x,v)\\&= d \bigg( {\rm det}\ g_{ab}\ x_i' \bigg)(-1)^i\ dx^1\cdots \widehat{dx^i}\cdots dx^n dv^1\cdots dv^n \\&- d \bigg( {\rm det}\ g_{ab}\ x_j'x_k'\Gamma_{jk}^l\bigg)(-1)^{n+l}\ dx^1\cdots dx^n dv^1\cdots \widehat{dv^l}\cdots dv^n\\&=0\end{align*}
Hence geodesic field has divergence $0$.