This question emerged when I was thinking about Liouville's theorem of classical mechanics. As far as I understand, the change of any function along the integral curves of Hamiltonian vector field is
$$\frac{d \rho}{d t} = \{\rho, H \} = \mathcal L_X \rho $$
where $\{\,,\}$ is Poisson bracket, $X$ is a vector field associated with $dH$ and $\mathcal L$ is Lie derivative. I have the following understanding of $\frac{d f}{d t}$ --- as we move on the manifold $M$ with a speed $X$, the landscape ($f$ visioned as mountains on $M$) changes. The mountains are still, it us who is moving.
On the other hand, in certain formulations of Liouville's theorem $\frac{\partial}{\partial t}$ is used:
$$\frac{\partial \rho}{\partial t} = - \{\rho, H \} $$
(The question is about the definition of $\frac{\partial}{\partial t}$, not about the contents of Liouville's theorem). I am familiar with the notion of material derivative and partial time derivative from usual formulations of fluid dynamics, but I am not sure of their counterparts in differential geometry.
What is partial time derivative $\frac{\partial}{\partial t}$ and how is it defined invariantly for functions, vector fields and differential forms? How is it related to $X$ and $\mathcal L_X$?
The thing I have missed is that phase space distribution is by definition time-dependent function. Indeed, it is defined as a propagation of the initial distribution with the Hamiltonian flow:
$$\rho_t = \varphi^*_{-t} \, \rho_0$$
thus, partial time derivative makes sense at any point regardless of any flows:
$$\frac{\partial \rho_t}{\partial t} = \lim_{s \to 0} \frac{\rho_{t+s} - \rho_t}{s}$$
On the other hand, the change of a time-dependent function along some flow, that is $\frac{d}{dt}$, is given by:
$$\frac{d}{dt} \left(\varphi^*_t f_t \right) = \varphi^*_t \left[\frac{\partial f_t}{\partial t} + \mathcal L_X \, f_t \right]$$
see, for example, Online Supplement for Jeffery Lee's "Manifolds and Differential Geometry". This expression, by the way, is valid for tensors of all ranks.
Returning to the Liouville's theorem, we get:
$$0 = \frac{d \rho_0}{d t} = \frac{d}{dt} \left(\varphi^*_t \rho_t \right) = \varphi^*_t \left[\frac{\partial \rho_t}{\partial t} + \mathcal L_X \, \rho_t \right],$$
or, in concise notation
$$\frac{d \rho_t}{d t} = \frac{\partial \rho_t}{\partial t} + \{ \rho_t, H \} = 0.$$