  We derive an adaptation of Li & Yau estimates for positive solutions of
semilinear heat equations on Riemannian manifolds with nonnegative Ricci
tensor. We then apply these estimates to obtain a Harnack inequality and to
discuss monotonicity, convexity, decay estimates and triviality of ancient and
eternal solutions.
