  Let $(M^{n},g)$ be a complete Riemannian manifold. In this paper, we
establish a space-time gradient estimates for positive solutions of nonlinear
parabolic equations $$\partial_{t}u(x,t)=\Delta u(x,t)+a u(x,t)(\log u(x,t))^b
+ q(x,t)A(u(x,t)),$$ on geodesic balls $B(O,r)$ in $M$ with $0<r\leq r$ for
$p>\frac{n}{2}$ when integral Ricci curvature $k(p,1)$ is small enough. By
integrating the gradient estimates, we find the corresponding Harnack
inequalities.
