Jacobian of a smooth flow on a smooth manifold This is a question about the derivative of Jacobian of the tangent map of a smooth flow. I think it might have an intrinsic formula, but I don't have a clue what it might be.
Let $M$ be an $n$-dimensional Riemannian manifold (say $C^\infty$, orientable, etc if it can save us from some trouble), $X:M\to\mathfrak{X}(M)$ be a smooth vector field (say nowhere vanishing), $\phi_t:M\to M$ be the flow induced by $X$. Now let's consider the tangent map 
$D_x\phi_t:T_xM\to T_{\phi_tx}M$.
Since each tangent space has an inner product structure induced by the metric, we can define the Jacobian of $\phi_t$ at $x$ as $J(t,x)=\det(D_x\phi_t)$.
Now we want to take the derivative with respect to $t$: $J_t(t,x)$. Is there any formula for this derivative? In fact a formula for $J_t(0,x)$ will be good enough.
We can choose a local orthonormal frame around $x$ and then write $D_x\phi_t$ as a time-dependent matrix, say $A(t,x)=(a_{jk}(t,x))=({\bf a}_1(t,x),\cdots,{\bf a}_n(t,x))$. Then $J(t,x)=\det A(t,x)$ and $\displaystyle J_t(t,x)=\sum_{1\le j\le n}\det({\bf a}_1(t,x),\cdots,{\bf \dot{a}}_i(t,x),\cdots,{\bf a}_n(t,x))$.
Then letting $t\to0$, we get that ${\bf a}_j(0,x)={\bf e}_j$ (since $\phi_0=\text{Id}$) and $\displaystyle J_t(0,x)=\sum_{1\le j\le n}\det({\bf e}_1,\cdots,{\bf \dot{a}}_j(0,x),\cdots,{\bf e}_n)=\sum_{1\le j\le n}\dot{a}_{jj}(0,x)$.
I feel like it should be a fundamental fact in Ordinary Differential Equations. But I don't know where to start. Thank you! 

A warmup: consider $M=\mathbb{R}$ and the equation $\dot{x}=f(x)$. Let $\phi(t,x)=\phi_t(x)$ be the solution with initial condition $\phi(0,x)=x$.
Then $J(t,x)=\partial_x \phi(t,x)$ and $J_t(t,x)=\partial_t\partial_x \phi(t,x)=\partial_x\partial_t\phi(t,x)=\partial_x f(x(t))$. In particular $J_t(0,x)=f'(x)$.
Then let $M=\mathbb{R}^2$ and the equation $\dot{x}_i=f_i({\bf x})$. Let $\vec{\phi}(t,{\bf x})=\vec{\phi}_t({\bf x})$ be the solution with initial condition $\vec{\phi}(0,{\bf x})={\bf x}$.
Then $J(t,x)=\partial_1\phi_1\cdot\partial_2\phi_2-\partial_2\phi_1\cdot\partial_1\phi_2$. Taking derivative w.r.t. $t$, $\dot{\phi}_i(t,x)=f_i({\bf x}(t))$ and hence $J_t(0,x)=\sum_j \partial_jf_j(x)=\mathrm{div}(f_1,f_2)$. 
So in general $J_t(0,x)=(\mathrm{div}X)(x)$.
 A: Following the first and third commends above, this question turns out to be a simple exercise in ODE. Let $F:\mathbb{R}^n\to\mathbb{R}^n$ be a smooth vector field and consider the ODE $\dot{x}=F(x)$  ---- $(\ast)$
For simplicity we assume $F$ is bounded and hence there exists a smooth flow $\phi:\mathbb{R}\times\mathbb{R}^n\to\mathbb{R}^n$, $(t,x)\mapsto \phi_t(x)$. For each fixed $t$, $\phi_t$ is a diffeomorphism and let $J(t,x)$ be the Jacobian of $\phi_t$ at $x$, that is, $J(t,x)=\det(D_x\phi_t)$. 
To find an expression for $J(t,\cdot)$, we first rewrite $(\ast)$ as 


*

*$\dot{\phi}(t,x)=F(\phi(t,x))$.


Then taking differential of $x$, the left side gives $D_x(\dot{\phi}(t,x))=\frac{d}{dt}D_x\phi(t,x)$, while the right side gives $D_x(F\circ \phi(t,x))=(D_xF)(\phi(t,x))\cdot D_x\phi(t,x)$. So we get a differential matrix equation:


*

*$\frac{d}{dt}D_x\phi(t,x)=(D_xF)(\phi(t,x))\cdot D_x\phi(t,x)$.


Using the same proof of Liouville's Theorem (here), we calculate 


*

*$\dot{J}(t,x)=\frac{d}{dt}\det D_x\phi(t,x)=\mathrm{Tr}(D_xF)(\phi(t,x))\cdot J(t,x)=\mathrm{div}F(\phi(t,x))\cdot J(t,x)$.


Letting $t=0$, we see $J(0,x)=1$ since $\phi_0=Id$ and hence $\dot{J}(0,x)=\mathrm{div}F(x)$.
More generally by taking integration, we have $J(t,x)=e^{\int_0^t\mathrm{div}F(\phi(s,x))ds}$.
