Gradient estimates for Leibenson's equation on Riemannian manifolds

We consider on Riemannian manifolds solutions of the Leibenson equation \begin{equation*} \partial {t}u=Δ{p}u^{q}. \end{equation*} This equation is also known as doubly nonlinear evolution equation. We prove gradient estimates for positive solutions $u$ under the condition that the Ricci curvature on $M$ is bounded from below by a non-positive constant. We distinguish between the case $q(p-1)>1$ (slow diffusion case) and the case $q(p-1)<1$ (fast diffusion case).